Télécharger sodiri.eso

Retour à la liste

Numérotation des lignes :

  1. C SODIRI SOURCE CHAT 05/01/13 03:21:20 5004
  2. SUBROUTINE SODIRI(D,X,B,MM,LL,KERRE)
  3. *
  4. *
  5. * resolution d'un du systeme A.X=B
  6. * par la methode du pivot de Gauss
  7. *
  8. IMPLICIT INTEGER(I-N)
  9. IMPLICIT REAL*8(A-H,O-Z)
  10. -INC CCOPTIO
  11. *
  12. REAL*8 X(*), B(*)
  13. DIMENSION D(1)
  14. * DO 12 I=1,LL
  15. * DO 13 J=I,LL
  16. * D(I+J)=A(I,J)
  17. *13 CONTINUE
  18. *12 CONTINUE
  19.  
  20. KERRE=0
  21. PREC=1.e-12
  22. DO 200 N=1,MM
  23. NN=LL*(N-1)+N
  24. DIAG=D(NN)
  25. IF(ABS(DIAG).LE.PREC) THEN
  26. KERRE=49
  27. RETURN
  28. ENDIF
  29. DIAG=1./DIAG
  30. DO 100 J=1,MM
  31. NJ=LL*(N-1)+J
  32. 100 D(NJ)=-D(NJ)*DIAG
  33. DO 150 I=1,MM
  34. IN=LL*(I-1)+N
  35. IF(I.EQ.N) GO TO 150
  36. DO 140 J=1,MM
  37. IF(J.EQ.N) GO TO 140
  38. NJ=LL*(N-1)+J
  39. IJ=LL*(I-1)+J
  40. D(IJ)=D(IJ)+D(IN)*D(NJ)
  41. 140 CONTINUE
  42. 150 D(IN)=D(IN)*DIAG
  43. 200 D(NN)=DIAG
  44. DO 20 I=1,MM
  45. X(I)=0
  46. DO 21 J=1,MM
  47. IND=LL*(I-1)+J
  48. X(I)=X(I)+D(IND)*B(J)
  49. 21 CONTINUE
  50.  
  51. 20 CONTINUE
  52. RETURN
  53. END
  54.  
  55.  
  56.  

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