Re: [ublas] Wrapping a large Fortran array

From: David M Garza (David.M.Garza_at_[hidden])
Date: 2006-01-31 14:38:08

SNOPT.  I'm using the interface where it will expect the matrix in 
compressed column form.  The argument list for the callback function 
only provides the value array--it expects you to remember the index 
arrays separately.  I guess if you know SNOPT you know this already... 
Thanks!
David

eeitan_at_[hidden] wrote:
> 
> Quoting David M Garza <David.M.Garza_at_[hidden]>:
> 
> 
>>I'm working with a numerical optimizer designed for large, sparse
>>problems and written in Fortran.  When the the optimizer calls back for
>>a Jacobian evaluation it passes the array where it expects the Jacobian
>>values to be placed.  Is there a simple way to build a uBLAS sparse
>>matrix object which simply wraps this array?  My first thought is no,
>>but my second thought is "I could make a custom allocator..."  Has
>>someone already tackled this problem?  I expect this array to be big
>>(10,000-100,000 elements) so I'd like to avoid just copying it.
>>
>>Thanks,
>>David
>>
>>_______________________________________________
>>ublas mailing list
>>ublas_at_[hidden]
>>http://lists.boost.org/mailman/listinfo.cgi/ublas
>>
> 
> if you could tell me which optimizer it is, i will try to do so.
> 
> 
> ----------------------------------------------------------------
> This message was sent using IMP, the Internet Messaging Program.
> 
> _______________________________________________
> ublas mailing list
> ublas_at_[hidden]
> http://lists.boost.org/mailman/listinfo.cgi/ublas
> 
>