Boost logo

Ublas :

Subject: Re: [ublas] compile error on gesv.hpp on 64Bit system
From: Kraus Philipp (philipp.kraus_at_[hidden])
Date: 2010-07-25 12:15:05


Hi Rutger,

Am 25.07.2010 um 16:15 schrieb Rutger ter Borg:

> Kraus Philipp wrote:
>
> [snip]
>
>> Do I have forgot anything? I hope anyone can help me and send me a
>> tip
>> or an idea.
>> Thanks a lot
>>
>
> Do you have an example of your code? Are you using fortran_int_t?

This is a code snip that I call:

  ublas::matrix<T, ublas::column_major> l_matrix(p_matrix);

  ublas::matrix<T, ublas::column_major> l_eigvec(l_matrix.size1(),
l_matrix.size2());
  ublas::vector<T> l_eigval(l_matrix.size1());

  ublas::vector<T> l_tmp1(l_eigval.size());
  ublas::matrix<T, ublas::column_major>
l_tmp2(l_matrix.size1(),l_matrix.size2());
  ublas::matrix<T, ublas::column_major>
l_tmp3(l_matrix.size1(),l_matrix.size2());

  linalg::geev( 'N', 'V', l_matrix, l_eigval, l_tmp1,l_tmp2,
l_eigvec, linalg::optimal_workspace() );

with this namespaces
namespace ublas = boost::numeric::ublas;
namespace blas = boost::numeric::bindings::blas;
namespace linalg = boost::numeric::bindings::lapack;

Is this enought information otherwise I would create tomorrow a full
program.
The code excerpt ist part of a template class.

Thx

Phil