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    MB234859  26/06/10    21:15:37     12569          
  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. -INC PPARAM
  21. -INC CCOPTIO
  22. -INC CCREEL
  23. -INC SILICRE
  24.  
  25.       MMATRI=MMATRX
  26.       MILIGN=IILIGN
  27.       segact milign
  28.       INC=IPNO(/1)
  29.  
  30.       MVECTD=MVECTX
  31.  
  32.       MDNOR=IDNORM
  33.       SEGACT MDNOR
  34.       IF(IDNORD.GT.0) THEN
  35.          MDNO1=IDNORD
  36.          SEGACT MDNO1
  37.       ELSE
  38.          MDNO1=MDNOR
  39.       ENDIF
  40.  
  41.       MDIAG=IDIAG
  42.       SEGACT,MDIAG
  43.  
  44. C
  45. C  ****  DESCENTE: ON RESOU L*C=B.  EN FAIT ON STOCKE C DANS B
  46. C
  47.       ilicre=ifacre
  48.       ligcre=ligcrp
  49.       call graco13(ilideb(1),posm(1),valm(1),vectbb(1),inc)
  50.  
  51. C
  52. C  ****  DIVISION PAR LE TERME DIAGONAL  ****
  53. C
  54.          DO 12 I=1,INC
  55.             J=I
  56.             VECTBB(J)=VECTBB(J)*DIAG(I)
  57.    12    CONTINUE
  58.  
  59. C
  60. C  ****  MONTEE  ****
  61. C
  62.  
  63.       ilicre=ifatra
  64.       ligcre=ligcrp
  65.       call graco14(ilideb(1),posm(1),valm(1),vectbb(1),inc)
  66.  
  67.       RETURN
  68.       END
  69. C
  70.  
  71.  
  72.  
  73.  
  74.  
  75.  
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  

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