Re: [ublas] lapack geqrf matrix Q

From: Eddie Vedder (Eddie.Vedder_at_[hidden])
Date: 2005-12-22 04:02:25

> You can extract the unitary matrix with DORMQR. There is a binding in 
> the boost_sandbox: ormqr.hpp

You are right. This is much better:

ublas::vector<double> vecTau( k );
lapack::geqrf( matAQR, vecTau );
ublas::identity_matrix<double> matI( k );
ublas::matrix<double, ublas::column_major> matC = matI;
lapack::ormqr( 'R', 'N', matAQR, vecTau, matC, lapack::optimal_workspace()
);

Ed

-- 
GMX DSL-Flatrate 1 Jahr kostenlos* + WLAN-Router ab 0,- Euro*
Bis 31.12.2005 einsteigen! Infos unter: http://www.gmx.net/de/go/dsl