Télécharger pivli.eso

Retour à la liste

Numérotation des lignes :

  1. C PIVLI SOURCE CB215821 16/04/21 21:17:57 8920
  2. SUBROUTINE PIVLI(A,B,M,L,LU)
  3. C
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. DIMENSION A(M,M),B(M)
  7. C
  8. C TEST de consistance numérique
  9. C
  10. DO 125 I=1,LU
  11. DO 126 K=1,LU
  12. C
  13. IF (ABS(A(I,K)).LE.1.E-6) THEN
  14. A(I,K)=1.E-6
  15. ENDIF
  16. CONTINUE
  17. IF (ABS(B(I)).LE.1.E-6) THEN
  18. B(I)=1.E-6
  19. ENDIF
  20. C
  21. 126 CONTINUE
  22. 125 CONTINUE
  23. C
  24. SUP=ABS(A(L,L))
  25. NUM=L
  26. EF=1.E-10
  27. C
  28. DO 30 I=L+1,LU
  29. C
  30. IF (ABS(A(I,L)).GT.SUP) THEN
  31. SUP=ABS(A(I,L))
  32. NUM=I
  33. ENDIF
  34. C
  35. 30 CONTINUE
  36. C
  37. IF (SUP.LT.EF) THEN
  38. C PRINT*
  39. C PRINT*, 'IL N EXISTE AUCUNE VALEUR DE RESOLUTION'
  40. C PRINT*
  41. C PRINT*, 'LA MATRICE N EST PAS INVERSIBLE'
  42. C PRINT*
  43. STOP
  44. ENDIF
  45. C
  46. DO 35 I=L,LU
  47. C
  48. T=A(L,I)
  49. A(L,I)=A(NUM,I)
  50. A(NUM,I)=T
  51. C
  52. 35 CONTINUE
  53. C
  54. T=B(L)
  55. B(L)=B(NUM)
  56. B(NUM)=T
  57. C
  58. RETURN
  59. END
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  

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