Télécharger nbinc.eso

Retour à la liste

Numérotation des lignes :

  1. C NBINC SOURCE PV 16/11/17 22:00:52 9180
  2. SUBROUTINE NBINC (IPRIGI,N)
  3.  
  4. ***********************************************************************
  5. *
  6. * N B I N C
  7. *
  8. * FONCTION:
  9. * ---------
  10. *
  11. * DETERMINATION DU NOMBRE D'INCONNUES A PARTIR D'UN MRIGID
  12. *
  13. * ENCAPSULATION PARTIELLE DE ALEA1.ESO
  14. *
  15. * PARAMETRES: (E)=ENTREE (S)=SORTIE
  16. * -----------
  17. *
  18. *
  19. * IPRIGI ENTIER (E) POINTEUR D'UNE RIGIDITE
  20. *
  21. *
  22. * N ENTIER (S) DIMENSION DU PROBLEME
  23. *
  24. *
  25. *
  26. *
  27. * AUTEUR, DATE DE CREATION:
  28. * -------------------------
  29. *
  30. * PASCAL BOUDA 2 JUIN 2015
  31. *
  32. * LANGAGE:
  33. * --------
  34. *
  35. * FORTRAN 77 & 90
  36. *
  37. ************************************************************************
  38.  
  39. IMPLICIT INTEGER(I-N)
  40. IMPLICIT REAL*8 (A-H,O-Z)
  41.  
  42. -INC CCOPTIO
  43. -INC SMMATRI
  44. -INC SMRIGID
  45.  
  46. REAL*8 V
  47. integer insym
  48. insym = 0
  49.  
  50. * -- DETERMINATION DU NOMBRE D'INCONNUES DU PROBLEME TRAITE --
  51.  
  52. MRIGID = IPRIGI
  53. SEGACT,MRIGID
  54. NRG = IRIGEL(/1)
  55. NBR = IRIGEL(/2)
  56. IPMATR = ICHOLE
  57. IF(NORINC.GT.0 .AND. NORIND.GT.0) THEN
  58. INSYM = 1
  59. ENDIF
  60. IF (NRG.GE.7) THEN
  61. DO 9 IN = 1,NBR
  62. IANTI=IRIGEL(7,IN)
  63. IF(IANTI.GT.0) THEN
  64. INSYM = 1
  65. ENDIF
  66. 9 CONTINUE
  67. ENDIF
  68. SEGDES,MRIGID
  69. *
  70. IF (IPMATR .EQ. 0) THEN
  71. IF (INSYM .EQ. 0) THEN
  72. CALL TRIANG (IPRIGI,1.D-18,0)
  73. ELSE
  74. CALL LDMT1(IPRIGI,1.D-18)
  75. ENDIF
  76. IF (IERR .NE. 0) RETURN
  77. MRIGID = IPRIGI
  78. SEGACT,MRIGID
  79. IPMATR = ICHOLE
  80. SEGDES,MRIGID
  81. END IF
  82. *
  83. MMATRI = IPMATR
  84. SEGACT,MMATRI
  85. MILIGN=IILIGN
  86. SEGDES,MMATRI
  87. SEGACT,MILIGN
  88. N=IPNO(/1)
  89. SEGDES,MILIGN
  90.  
  91.  
  92.  
  93. END
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales