Télécharger nbinc.eso

Retour à la liste

Numérotation des lignes :

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

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