Télécharger la plaquette Cast3M


fr en

Télécharger graco8.eso

Retour à la liste

Numérotation des lignes : 

graco8

  1. C GRACO8    SOURCE    PV090527  23/03/08    12:22:21     11617          
  2. C
  3.       SUBROUTINE GRACO8(MMATRX,MVECTX,NOEN,ifacre,ifatra)
  4. C
  5. C  **** EXECUTE LA SOLUTION X DE  (Lt D L) X=F
  6. C
  7. CMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMB
  8. CMB
  9. CMB   Plutot la solution de L.D.Lt ou L.D.Mt (cas non symétrique)
  10. CMB   Elle devrait dons s'appeller DESMON et non MONDES.
  11. CMB
  12. CMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMBCMB
  13. C
  14.       IMPLICIT INTEGER(I-N)
  15.       IMPLICIT REAL*8(A-H,O-Z)
  16.  
  17. -INC SMMATRI
  18. -INC SMELEME
  19. -INC SMVECTD
  20.  
  21. -INC PPARAM
  22. -INC CCOPTIO
  23. -INC CCREEL
  24.       segment ilicre
  25. *  stockage matrice factorisee en creux
  26. *  ilideb position debut de ligne dans ligcre
  27.         integer ilideb(nbinc+1)
  28.         integer ligcrp
  29.       endsegment
  30.       segment ligcre
  31. *  lmatr: longueur reelle ligne
  32. *  posm: numero inconnue
  33. *  valm: valeur terme
  34.         integer posm(lmat)
  35.         real*8 valm(lmat)
  36.       endsegment
  37.  
  38.       MMATRI=MMATRX
  39.  
  40.       MILIGN=IILIGN
  41.       segact milign
  42.       INC=IPNO(/1)
  43.  
  44.       MVECTD=MVECTX
  45.  
  46.       MDNOR=IDNORM
  47.       SEGACT MDNOR
  48.       IF(IDNORD.GT.0) THEN
  49.          MDNO1=IDNORD
  50.          SEGACT MDNO1
  51.       ELSE
  52.          MDNO1=MDNOR
  53.       ENDIF
  54.  
  55.       MDIAG=IDIAG
  56.       SEGACT,MDIAG
  57.  
  58. C
  59. C  ****  DESCENTE: ON RESOU L*C=B.  EN FAIT ON STOCKE C DANS B
  60. C
  61.       ilicre=ifacre
  62.       ligcre=ligcrp
  63.       call graco13(ilideb(1),posm(1),valm(1),vectbb(1),inc)
  64.  
  65. C
  66. C  ****  DIVISION PAR LE TERME DIAGONAL  ****
  67. C
  68.          DO 12 I=1,INC
  69.             J=I
  70.             VECTBB(J)=VECTBB(J)*DIAG(I)
  71.    12    CONTINUE
  72.  
  73. C
  74. C  ****  MONTEE  ****
  75. C
  76.  
  77.       ilicre=ifatra
  78.       ligcre=ligcrp
  79.       call graco14(ilideb(1),posm(1),valm(1),vectbb(1),inc)
  80.  
  81.       RETURN
  82.       END
  83. C
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  

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