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.  
  43. -INC PPARAM
  44. -INC CCOPTIO
  45. -INC SMMATRI
  46. -INC SMRIGID
  47.  
  48. REAL*8 V
  49. integer insym
  50. insym = 0
  51.  
  52. * -- DETERMINATION DU NOMBRE D'INCONNUES DU PROBLEME TRAITE --
  53.  
  54. MRIGID = IPRIGI
  55. SEGACT,MRIGID
  56. NRG = IRIGEL(/1)
  57. NBR = IRIGEL(/2)
  58. IPMATR = ICHOLE
  59. IF(NORINC.GT.0 .AND. NORIND.GT.0) THEN
  60. INSYM = 1
  61. ENDIF
  62. IF (NRG.GE.7) THEN
  63. DO 9 IN = 1,NBR
  64. IANTI=IRIGEL(7,IN)
  65. IF(IANTI.GT.0) THEN
  66. INSYM = 1
  67. ENDIF
  68. 9 CONTINUE
  69. ENDIF
  70. SEGDES,MRIGID
  71. *
  72. IF (IPMATR .EQ. 0) THEN
  73. IF (INSYM .EQ. 0) THEN
  74. CALL TRIANG (IPRIGI,1.D-18,0)
  75. ELSE
  76. CALL LDMT1(IPRIGI,1.D-18)
  77. ENDIF
  78. IF (IERR .NE. 0) RETURN
  79. MRIGID = IPRIGI
  80. SEGACT,MRIGID
  81. IPMATR = ICHOLE
  82. SEGDES,MRIGID
  83. END IF
  84. *
  85. MMATRI = IPMATR
  86. SEGACT,MMATRI
  87. MILIGN=IILIGN
  88. SEGDES,MMATRI
  89. SEGACT,MILIGN
  90. N=IPNO(/1)
  91. SEGDES,MILIGN
  92.  
  93.  
  94.  
  95. END
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  

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