Télécharger sodiri.eso

Retour à la liste

Numérotation des lignes :

sodiri
  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.  
  11. -INC PPARAM
  12. -INC CCOPTIO
  13. *
  14. REAL*8 X(*), B(*)
  15. DIMENSION D(1)
  16. * DO 12 I=1,LL
  17. * DO 13 J=I,LL
  18. * D(I+J)=A(I,J)
  19. *13 CONTINUE
  20. *12 CONTINUE
  21.  
  22. KERRE=0
  23. PREC=1.e-12
  24. DO 200 N=1,MM
  25. NN=LL*(N-1)+N
  26. DIAG=D(NN)
  27. IF(ABS(DIAG).LE.PREC) THEN
  28. KERRE=49
  29. RETURN
  30. ENDIF
  31. DIAG=1./DIAG
  32. DO 100 J=1,MM
  33. NJ=LL*(N-1)+J
  34. 100 D(NJ)=-D(NJ)*DIAG
  35. DO 150 I=1,MM
  36. IN=LL*(I-1)+N
  37. IF(I.EQ.N) GO TO 150
  38. DO 140 J=1,MM
  39. IF(J.EQ.N) GO TO 140
  40. NJ=LL*(N-1)+J
  41. IJ=LL*(I-1)+J
  42. D(IJ)=D(IJ)+D(IN)*D(NJ)
  43. 140 CONTINUE
  44. 150 D(IN)=D(IN)*DIAG
  45. 200 D(NN)=DIAG
  46. DO 20 I=1,MM
  47. X(I)=0
  48. DO 21 J=1,MM
  49. IND=LL*(I-1)+J
  50. X(I)=X(I)+D(IND)*B(J)
  51. 21 CONTINUE
  52.  
  53. 20 CONTINUE
  54. RETURN
  55. END
  56.  
  57.  
  58.  

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