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@gmail.com> 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