Re: [ublas] To calculate svd ?

From: Karl Meerbergen (Karl.Meerbergen_at_[hidden])
Date: 2005-01-13 08:56:04

Rakesh K Sinha wrote:

>Hi, 
>  I was wondering if there is any way I could get the SVD of the given matrix. 
>I was checking the documentation at the Boost library web page and could not
>find any there. Hence I am posting it here.
>
>  A word of thanks to the uBLAS developers for coming up with a nice C++ modern
>API that makes it a compelling a library to integrate with C++ applications.
>
>--
>  Karthik. 
>
>
>	
>

Hello,

You could use the lapack bindings in boost-sandbox 
(boost/numeric/bindings/lapack/gesvd.hpp).
First you have to include the traits class for the matrix and vector 
bindings. For example, if your matrices are 
ublas::matrix<T,column_major> include 
boost/numeric/bindings/traits/ublas_matrix.hpp.
If the svd values are stored in a ublas vector also include
boost/numeric/bindings/traits/ublas_vector.hpp

Karl