は、QRに大きな行列を適用し、私は固有のライブラリでC++にMATLABコードを変換する必要がありますが、私はQRでいくつかの問題が分解してい

Question

を分解し、MATLABが機能を持っています

[Q,R]=qr(A,0); // A is m-by-n 

それは経済規模を生産m>nの場合は、の最初のnカラムと、Rの最初のnローだけが計算されます。 m<=nの場合は[Q,R]=qr(A)と同じです。

私は固有ライブラリで計算しようとしました。しかし、Aは20000x1000なので、QR分解時にアプリケーションクラッシュが常に発生します。固有モードやその他の方法でエコノミーサイズの分解をどのように生成するかわかりません。

[Q,R]=qr(A,0)をC++/Eigenに変換するにはどうすればよいですか？

Answer 1

MatrixXd A(m,n);  
HouseholderQR<MatrixXd> qr; 
qr.compute(A); 
temp= qr.matrixQR().triangularView<Upper>(); 
Q= qr.householderQ() * Eigen::MatrixXd::Identity(m, n); 
R=temp.topRows(n);