Télécharger la plaquette Cast3M


fr en

Télécharger pivli.eso

Retour à la liste

Numérotation des lignes : 

pivli

  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