Esope source files

   1 : C DESCLU    SOURCE    CHAT      05/01/12    22:43:21     5004           
   2 : CCC
   3 : C Descomposicion LU
   4 : CCC
   5 :       SUBROUTINE DESCLU (A,N)
   6 :       IMPLICIT INTEGER(I-N)
   7 :       integer n,i,j,k
   8 :       real*8    A(n,n),aux
   9 :       do k=1,n-1
  10 :           a(1,k+1)=a(1,k+1)/a(1,1)
  11 :           do i=2,k
  12 :               aux=0.D0
  13 :               do j=1,i-1
  14 :                   aux=aux+a(i,j)*a(j,k+1)
  15 :               enddo
  16 :               a(i,k+1)=(a(i,k+1)-aux)/a(i,i)
  17 :           enddo
  18 :           do i=2,k+1
  19 :               aux=0.D0
  20 :               do j=1,i-1
  21 :                   aux=aux+a(k+1,j)*a(j,i)
  22 :               enddo
  23 :               a(k+1,i)=a(k+1,i)-aux
  24 :           enddo
  25 :       enddo
  26 :       return
  27 :       end
  28 :  
  29 :