Re: [ublas] [bindings][lapack] cmputational/stegr problems

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