Esope source files

   1 : C LUINV     SOURCE    CHAT      05/01/13    01:26:19     5004           
   2 : CCC
   3 : C Calculo de Ainv=(LUa)^-1
   4 : CCC
   5 :       SUBROUTINE LUINV (A,AI,N)
   6 :       IMPLICIT INTEGER(I-N)
   7 :       integer n,i,j
   8 :       real*8 A(n,n),Ai(n,n),w(10)
   9 :       do i=1,n
  10 :         do j=1,n
  11 :           Ai(j,i)=0.D0
  12 :           w(j)=0.D0
  13 :         enddo
  14 :         w(i)=1.D0
  14 :         call LUiw(A,w,Ai(1,i),n)
  15 :       enddo
  16 :       return
  17 :       end
  18 :  
  19 :