Sources Esope | Cast3M

trili
C TRILI     SOURCE    PV        11/03/07    21:18:37     6885
      SUBROUTINE TRILI(A,B,M,LU)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION A(M,M),B(M)
C
      REAL*8 TEMP,PIV,DET
      DET=1.D0
C
      DO 10 K=1,LU
      CALL PIVLI(A,B,M,K,LU)
      PIV=A(K,K)
      DET=DET*PIV
C
      DO 20 I=K+1,LU
      TEMP=A(I,K)/A(K,K)
      B(I)=B(I)-TEMP*B(K)
C
      DO 30 J=K+1,LU
      A(I,J)=A(I,J)-TEMP*A(K,J)
   30 CONTINUE
C
   20 CONTINUE
C
      DET=DET*A(LU,LU)
   10 CONTINUE
C
      RETURN
      END