Sources Esope | Cast3M

vtluiw
C VTLUIW    SOURCE    CHAT      05/01/13    04:10:14     5004
CCC
C Producto vt*(A)inv*w con A almacenada como LU
CCC
             FUNCTION VTLUIW (V,A,W,N)
      IMPLICIT INTEGER(I-N)
      integer n,i
      real*8    v(n),w(n),A(n,n),auxv(8),aux1,vtLUiw
c L U auxv = w
      call LUiw(A,w,auxv,n)
c vt_LUinv_w = v^t * auxv
      aux1=0.D0
      do i=1,n
         aux1=aux1+v(i)*auxv(i)
      enddo
c      write(*,*)' nA-1m',aux1
      vtLUiw=aux1
      return
      end