Télécharger la plaquette Cast3M


fr en

Télécharger nbinc.eso

Retour à la liste

Numérotation des lignes : 

nbinc

  1. C NBINC     SOURCE    MB234859  26/08/26    21:15:16     12626          
  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.         CALL TRIANG(IPRIGI,xspetl,0,0,INSYM)
  79.         IF (IERR .NE. 0) RETURN
  80.         MRIGID = IPRIGI
  81.         SEGACT,MRIGID
  82.         IPMATR = ICHOLE
  83.         SEGDES,MRIGID
  84.       ENDIF
  85. *
  86.       MMATRI = IPMATR
  87.       SEGACT,MMATRI
  88.       MILIGN=IILIGN
  89.       SEGDES,MMATRI
  90.       SEGACT,MILIGN
  91.       N=IPNO(/1)
  92.       SEGDES,MILIGN
  93.  
  94.       END
  95.  
  96.  

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