固有QR分解結果は2つの方法と異なります

Question

私はQigen分解をEigenを使用して試していますが、結果は次のようなトウの方法から得られますが、エラーを見つけるのを手伝ってください！おかげさまで

Answer 1

// Initialize the sparse matrix 
A.setFromTriplets(triplets.begin(), triplets.end()); 
A.makeCompressed(); 

//Dense matrix method 
MatrixXd MatrixA = A; 
HouseholderQR<MatrixXd> qr(MatrixA); 
MatrixXd Rr = qr.matrixQR().triangularView<Upper>(); 

//Sparse matrix method 
SparseQR < SparseMatrix <double>, COLAMDOrdering<int> > qr; 
qr.compute(A); 
SparseMatrix<double, RowMajor> Rr = qr.matrixR(); 

SparseQR行う列の並べ替えは、フィルインを低減し、ColPivHouseholderQRと同様略ランク明らか分解を達成するために両方のためです。より正確には、HouseholderQRは、A = Q*Rを計算するが、SparseQRは、A*P = Q*Rを計算する。したがって、2つの三角要素が異なることが予想されます。