Re: [boost] boost::ublas bindings ? (eigenvalues)

From: Stjepan Rajko (stipe_at_[hidden])
Date: 2007-05-03 14:34:31

Philippe,

The following code snippets might be helpful - I used singular value
decomposition using lapack bindings to get eigenstuff.  I don't know
if it's the same thing as eig() in your case.

ublas::matrix<double > A(rows, cols);

ublas::matrix<double > ATA;
ATA = prod(trans(A),A);

ublas::matrix<double, boost::numeric::ublas::column_major > U(cols,
cols), Vt(cols,cols);
ublas::vector<double> S(cols);
boost::numeric::bindings::lapack::gesvd(ATA, S, U, Vt);

cout << U << endl << endl;
cout << Vt << endl << endl;
cout << S << endl << endl;

HTH,

Stjepan

On 5/3/07, Philippe Vaucher <philippe.vaucher_at_[hidden]> wrote:
> Nevermind I found them in the Sandbox.
> If someone still wants to help me on my eig() problem he's welcome ! :)
>
> Philippe
> _______________________________________________
> Unsubscribe & other changes: http://lists.boost.org/mailman/listinfo.cgi/boost
>