[ublas] novice question. how to use ublas bindings

Subject: [ublas] novice question. how to use ublas bindings
From: Luca Sbardella (luca_at_[hidden])
Date: 2009-08-19 17:22:06

I got the binding from the boost sandbox

http://svn.boost.org/svn/boost/sandbox/numeric_bindings/boost/numeric/bindings/

I'm trying to see if I can use them the to solve an eigenvalue problem
within uBlas.
How do I call the bindings?

I tried to compile the file in
http://svn.boost.org/svn/boost/sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_gees.cpp

but it fails. It seems that is not compatible with the binding for two
reasons

1) the import

<boost/numeric/bindings/lapack/gees.hpp>

is wrong. Should that be

<boost/numeric/bindings/lapack/driver/gees.hpp>   ?

2) How do I call the function gees? In the example above it fails.

lapack::gees( a, e1, z, workspace ) ;

It looks like the functions arguments are wrong. I assume the test is based
on an old version of the bindings.
In fact looking at the header file

http://svn.boost.org/svn/boost/sandbox/numeric_bindings/boost/numeric/bindings/lapack/driver/gees.hpp

the functions requires two extra arguments (char const jobvs, char const
sort).

Any help would be appreciated.

---
Dr Luca Sbardella
Director
Dynamic Quant Ltd