On Mon, Sep 7, 2009 at 8:42 AM, Yashodhan Nevatia <yashodhan.nevatia@spaceapplications.com> wrote:

Hi all,

The code fragment is:

double min_delta=0.00001;
matrix<double> U(3,3);

Can't comment on any of the details, but you need to make sure you matrices are column major to support interfacing with fortran. ie.:
matrix<double, ublas::column_major> U(3,3);

I think that in the newest versions of the bindings this will be verified at compilation time, but until then you need to be very careful.