Boost logo

Ublas :

From: Thomas Lemaire (thomas.lemaire_at_[hidden])
Date: 2005-02-07 09:50:17

Dear list,

I already succesfully used lapack::gesv() binding.

When calling lapack::syev(), my program hangs. I get no error from boost
bindings, neither from LAPACK. I both tryed with minimal and optimal
workspace, and 'L' and 'U' part of the matrix.

any help ??

Here is a test code:

#include "boost/numeric/bindings/lapack/syev.hpp"
#include "boost/numeric/bindings/traits/ublas_matrix.hpp"
#include "boost/numeric/bindings/traits/ublas_vector.hpp"

namespace lapack = boost::numeric::bindings::lapack;

void test_syev(const ublas::symmetric_matrix<double>& cov_,
ublas::vector<double>& lambda)
    ublas::matrix<double, ublas::column_major> A(cov_);
    lapack::syev( 'V', 'L', A, lambda, lapack::minimal_workspace() );

PS: I know lapack proposes some better suited routines for symmetric
real matrix, but the bindings do not exist yet in boost, so I first want
to try this more general routine...