Re: [ublas] binding to LAPACK routine gees()

Subject: Re: [ublas] binding to LAPACK routine gees()
From: Thomas Klimpel (Thomas.Klimpel_at_[hidden])
Date: 2008-11-10 11:28:47

Hi,

I don't really understand the error message from the compiler, but I do
see three issues in your code.

>    double workspace[10000];

You should include "boost/numeric/bindings/traits/c_array.hpp" if you
intent to use plain c-arrays with the bindings.

>    ublas::vector<double> v1 (3), v2 (3);

You have to use a complex vector for the eigenvalues if you intent to
use lapack::gees.

>    lapack::gees(m,v1,v2,workspace);

The call sequence with 4 arguments expects a matrix as 3rd argument.

Regards,
Thomas