Sources Esope | Cast3M

luinv
C LUINV     SOURCE    CHAT      05/01/13    01:26:19     5004
CCC
C Calculo de Ainv=(LUa)^-1
CCC
      SUBROUTINE LUINV (A,AI,N)
      IMPLICIT INTEGER(I-N)
      integer n,i,j
      real*8 A(n,n),Ai(n,n),w(10)
      do i=1,n
        do j=1,n
          Ai(j,i)=0.D0
          w(j)=0.D0
        enddo
        w(i)=1.D0
        call LUiw(A,w,Ai(1,i),n)
      enddo
      return
      end