Télécharger nbinc.eso

Retour à la liste

Numérotation des lignes :

  1. C NBINC SOURCE PB245956 20/12/21 21:15:10 10747
  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. N=0
  52.  
  53. * -- DETERMINATION DU NOMBRE D'INCONNUES DU PROBLEME TRAITE --
  54.  
  55. MRIGID = IPRIGI
  56. SEGACT,MRIGID
  57. NRG = IRIGEL(/1)
  58. NBR = IRIGEL(/2)
  59. * cas d'une matrice vide
  60. IF (NBR.EQ.0) RETURN
  61. IPMATR = ICHOLE
  62. IF(NORINC.GT.0 .AND. NORIND.GT.0) THEN
  63. INSYM = 1
  64. ENDIF
  65. IF (NRG.GE.7) THEN
  66. DO 9 IN = 1,NBR
  67. IANTI=IRIGEL(7,IN)
  68. IF(IANTI.GT.0) THEN
  69. INSYM = 1
  70. ENDIF
  71. 9 CONTINUE
  72. ENDIF
  73. SEGDES,MRIGID
  74. *
  75. IF (IPMATR .EQ. 0) THEN
  76. IF (INSYM .EQ. 0) THEN
  77. CALL TRIANG (IPRIGI,1.D-18,0)
  78. ELSE
  79. CALL LDMT1(IPRIGI,1.D-18)
  80. ENDIF
  81. IF (IERR .NE. 0) RETURN
  82. MRIGID = IPRIGI
  83. SEGACT,MRIGID
  84. IPMATR = ICHOLE
  85. SEGDES,MRIGID
  86. END IF
  87. *
  88. MMATRI = IPMATR
  89. SEGACT,MMATRI
  90. MILIGN=IILIGN
  91. SEGDES,MMATRI
  92. SEGACT,MILIGN
  93. N=IPNO(/1)
  94. SEGDES,MILIGN
  95.  
  96.  
  97.  
  98. END
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  

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