|
|
Ublas : |
Subject: Re: [ublas] [bindings][lapack] cmputational/stegr problems
From: Rutger ter Borg (rutger_at_[hidden])
Date: 2009-04-14 02:41:43
I knew this would come around at some point -- should be fixed now.
Interesting, I see that stegr has an abstol argument that is not used
anymore by (at least by version 3.2 of) lapack. Perhaps we should do
something about that as well.
Kind regards,
Rutger
Hidekazu Ikeno wrote:
> Thanks. Now traits::is_complex<std::complex<T> > is evaluated as expected.
>
> I found one more bug in stegr template funstions.
> The arguments vl, vu and abstol should be real type, i.e.,
>
> template< typename VectorD, typename VectorE, typename VectorW,
> typename MatrixZ, typename VectorISUPPZ >
> inline integer_t stegr( char const jobz, char const range,
> integer_t const n, VectorD& d, VectorE& e,
> typename traits::matrix_traits< MatrixZ >::value_type const vl,
> //
> <--
> typename traits::matrix_traits< MatrixZ >::value_type const vu, //
> <--
> integer_t const il, integer_t const iu,
> typename traits::matrix_traits< MatrixZ >::value_type const
> abstol,
> // <--
> integer_t& m, VectorW& w, MatrixZ& z, VectorISUPPZ& isuppz );
>
> should be modified to
>
> template< typename VectorD, typename VectorE, typename VectorW,
> typename MatrixZ, typename VectorISUPPZ >
> inline integer_t stegr( char const jobz, char const range,
> integer_t const n, VectorD& d, VectorE& e,
> typename traits::vector_traits< VectorW >::value_type const vl, //
> <--
> typename traits::vector_traits< VectorW >::value_type const vu, //
> <--
> integer_t const il, integer_t const iu,
> typename traits::vector_traits< VectorW >::value_type const
> abstol,
> // <--
> integer_t& m, VectorW& w, MatrixZ& z, VectorISUPPZ& isuppz );
>
> (it's really complicated...).
>
> At the moment, the binding works fines if we call stegr_impl<T>::invoke()
> directly.
>
>
>
> Best regards,
>
> Hidekazu Ikeno
>
>