Esope source files

   1 : C TRILI     SOURCE    PV        11/03/07    21:18:37     6885           
   2 :       SUBROUTINE TRILI(A,B,M,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 :       REAL*8 TEMP,PIV,DET
   9 :       DET=1.D0
  10 : C
  11 :       DO 10 K=1,LU
  11 :       CALL PIVLI(A,B,M,K,LU)
  12 :       PIV=A(K,K)
  13 :       DET=DET*PIV
  14 : C
  15 :       DO 20 I=K+1,LU
  16 :       TEMP=A(I,K)/A(K,K)
  17 :       B(I)=B(I)-TEMP*B(K)
  18 : C
  19 :       DO 30 J=K+1,LU
  20 :       A(I,J)=A(I,J)-TEMP*A(K,J)
  21 :    30 CONTINUE
  22 : C
  23 :    20 CONTINUE
  24 : C
  25 :       DET=DET*A(LU,LU)
  26 :    10 CONTINUE
  27 : C
  28 :       RETURN
  29 :       END
  30 : 
  31 :  
  32 :  
  33 :