[Coin-ipopt] How to use subroutine EV_HLV(...)

[Su Ba]

Dear Andreas,

1.  What is different between EV_H(...) and EV_HLV(...) ?  Could I save 
memory when I used EV_HLV(..)?

2.  Could you explain to me how to use subroutine EV_HLV(...)?

I have the problem of very large Hessian matrix.  Here is a small example of 
my Hessian:

        H1
            H2         0
                H3
                   .
H=                  .
             0           .
                           H16


where:  Hi (which is one of H1,H2,...H16) are matrices of 6x6.  They have a 
similar form except their values.

Now, if I use subroutine EV_HLV(...) like this, is its definition like this 
correct? (Suppose that algorithm parameters are changed as Ipopt manual)


subroutine EV_HLV(...)
....
     do i=1,16
          calculate Hi =hi(X)                          (to get 6x6 matrix)
          do j=1,6
          calculate indice of k-th variable
          VOUT(k)=VOUT(k) + Hi*VIN(k)
          end do
     end do
end

Are dimension of VOUT and VIN 6x16 =96 in my example ?

Thank you in advance.

Danh

_________________________________________________________________
Express yourself instantly with MSN Messenger! Download today it's FREE! 
http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/