Sources Esope | Cast3M

luiw
C LUIW      SOURCE    CHAT      05/01/13    01:26:37     5004
CCC
C Producto v=(A)inv*w con A almacenada como LU
CCC
      SUBROUTINE LUIW (A,W,V,N)
      IMPLICIT INTEGER(I-N)
      integer n,i,j
      real*8    v(n),w(n),A(n,n)
c forward subst L v = w
      v(1)=w(1)/A(1,1)
      do i=2,n
         v(i)=w(i)
         do j=1,i-1
            v(i)=v(i)-A(i,j)*v(j)
         enddo
         v(i)=v(i)/A(i,i)
      enddo
c backward subs U v = v
      do i=n-1,1,-1
         do j=i+1,n
            v(i)=v(i)-A(i,j)*v(j)
         enddo
      enddo
      return
      end