|
|
Ublas : |
Subject: Re: [ublas] novice question. how to use ublas bindings
From: Rutger ter Borg (rutger_at_[hidden])
Date: 2009-08-20 01:31:26
Luca Sbardella wrote:
> 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.
>
Hello Luca,
the bindings at sandbox/numeric_bindings are the "trunk" version, whereas
sandbox/numeric_bindings-v1 are "stable"
You could try the more stable version, or help us improve the trunk version.
For the trunk version, you could try 'N' for jobvs, and 'N' for sort. Select
is probably new, too, please try 0 for that one. If unsure about any
arguments of a lapack driver functions, please check the manual of lapack
itself. E.g., consult your function at
http://www.netlib.org/lapack/double/
http://www.netlib.org/lapack/complex/
Cheers,
Rutger