Télécharger graco8.eso

Retour à la liste

Numérotation des lignes :

  1. C GRACO8 SOURCE PV 16/11/17 21:59:38 9180
  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 CCOPTIO
  21. -INC CCREEL
  22. segment ilicre
  23. * stockage matrice factorisee en creux
  24. * ilideb position debut de ligne dans ligcre
  25. integer ilideb(nbinc+1)
  26. integer ligcrp
  27. endsegment
  28. segment ligcre
  29. * lmatr: longueur reelle ligne
  30. * posm: numero inconnue
  31. * valm: valeur terme
  32. integer posm(lmat)
  33. real*8 valm(lmat)
  34. endsegment
  35.  
  36. MMATRI=MMATRX
  37.  
  38. MILIGN=IILIGN
  39. segact milign
  40. INC=IPNO(/1)
  41.  
  42. MVECTD=MVECTX
  43.  
  44. MDNOR=IDNORM
  45. SEGACT MDNOR
  46. IF(IDNORD.GT.0) THEN
  47. MDNO1=IDNORD
  48. SEGACT MDNO1
  49. ELSE
  50. MDNO1=MDNOR
  51. ENDIF
  52.  
  53. MDIAG=IDIAG
  54. SEGACT,MDIAG
  55.  
  56. C
  57. C **** DESCENTE: ON RESOU L*C=B. EN FAIT ON STOCKE C DANS B
  58. C
  59. ilicre=ifacre
  60. ligcre=ligcrp
  61. call graco13(ilideb(1),posm(1),valm(1),vectbb(1),inc)
  62.  
  63. C
  64. C **** DIVISION PAR LE TERME DIAGONAL ****
  65. C
  66. DO 12 I=1,INC
  67. J=I
  68. VECTBB(J)=VECTBB(J)/DIAG(I)
  69. 12 CONTINUE
  70.  
  71. C
  72. C **** MONTEE ****
  73. C
  74.  
  75. ilicre=ifatra
  76. ligcre=ligcrp
  77. call graco14(ilideb(1),posm(1),valm(1),vectbb(1),inc)
  78.  
  79. RETURN
  80. END
  81. C
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  

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