Boost logo

Ublas :

Subject: Re: [ublas] [bindings][lapack] cmputational/stegr problems
From: Rutger ter Borg (rutger_at_[hidden])
Date: 2009-04-13 05:37:03


I've improved argument selection for type selection, which causes stegr (and
a dozen other routines) to choose a corrected argument. I've also added
minimal workspace functions for stegr. It's committed to subversion, please
let me now if it works for you.

Kind regards,

Rutger ter Borg

Affected routines:
boost/numeric/bindings/lapack/computational/bdsqr.hpp
boost/numeric/bindings/lapack/computational/gebak.hpp
boost/numeric/bindings/lapack/computational/ggbak.hpp
boost/numeric/bindings/lapack/computational/hsein.hpp
boost/numeric/bindings/lapack/computational/pteqr.hpp
boost/numeric/bindings/lapack/computational/ptrfs.hpp
boost/numeric/bindings/lapack/computational/pttrs.hpp
boost/numeric/bindings/lapack/computational/stedc.hpp
boost/numeric/bindings/lapack/computational/stegr.hpp
boost/numeric/bindings/lapack/computational/stein.hpp
boost/numeric/bindings/lapack/computational/steqr.hpp
boost/numeric/bindings/lapack/computational/tgevc.hpp
boost/numeric/bindings/lapack/computational/tgsen.hpp
boost/numeric/bindings/lapack/computational/tgsna.hpp
boost/numeric/bindings/lapack/computational/trevc.hpp
boost/numeric/bindings/lapack/computational/trsna.hpp
boost/numeric/bindings/lapack/driver/lalsd.hpp
boost/numeric/bindings/lapack/driver/ptsv.hpp
boost/numeric/bindings/lapack/driver/ptsvx.hpp

Hidekazu Ikeno wrote:

> Hi,
>
> I'm trying to use a 'computational/stegr' binding in my program, but I
> always get
> STATIC_ASSERTION_FAILUE when the complex values are used.
> I have read the lapack/computational/stegr.hpp and have found that , in
> the template
> function to call stegr, value_type is declared as
>
> typedef typename traits::vector_traits< VectorD >::value_type
> value_type;
>
> Because VectorD is a real vector even in the complex version of stegr, the
> template
> function of stegr always trying to instantiate stegr_impl with real types.
> A proper declaration of value_type should be
>
> typedef typename traits::matrix_traits< MatrixZ >::value_type
> value_type;
>
> I think some modifications for lapack_generator.py are required to fix
> this problem.
> Is it easy to fix?
>
>
> Hidekazu Ikeno