Esope source files

   1 : C LUIW      SOURCE    CHAT      05/01/13    01:26:37     5004           
   2 : CCC
   3 : C Producto v=(A)inv*w con A almacenada como LU
   4 : CCC
   5 :       SUBROUTINE LUIW (A,W,V,N)
   6 :       IMPLICIT INTEGER(I-N)
   7 :       integer n,i,j
   8 :       real*8    v(n),w(n),A(n,n)
   9 : c forward subst L v = w
  10 :       v(1)=w(1)/A(1,1)
  11 :       do i=2,n
  12 :          v(i)=w(i)
  13 :          do j=1,i-1
  14 :             v(i)=v(i)-A(i,j)*v(j)
  15 :          enddo
  16 :          v(i)=v(i)/A(i,i)
  17 :       enddo
  18 : c backward subs U v = v
  19 :       do i=n-1,1,-1
  20 :          do j=i+1,n
  21 :             v(i)=v(i)-A(i,j)*v(j)
  22 :          enddo
  23 :       enddo
  24 :       return
  25 :       end
  26 :  
  27 :