CFD/eigen
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Add reference for the default threshold in sparse QR

Browse Source
This commit is contained in:
Desire NUENTSA 2013-02-25 14:26:55 +01:00
parent ced8dfc0d9
commit cc35c44256
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
Eigen/src/SparseQR/SparseQR.h
View File

@ -286,7 +286,10 @@ void SparseQR<MatrixType,OrderingType>::factorize(const MatrixType& mat)
m_pmat.innerNonZeroPtr()[p] = mat.outerIndexPtr()[i+1] - mat.outerIndexPtr()[i]; m_pmat.innerNonZeroPtr()[p] = mat.outerIndexPtr()[i+1] - mat.outerIndexPtr()[i];
} }
// Compute the default threshold. /* Compute the default threshold, see :
* Tim Davis, "Algorithm 915, SuiteSparseQR: Multifrontal Multithreaded Rank-Revealing
* Sparse QR Factorization, ACM Trans. on Math. Soft. 38(1), 2011, Page 8:3
*/
if(m_useDefaultThreshold) if(m_useDefaultThreshold)
{ {
RealScalar max2Norm = 0.0; RealScalar max2Norm = 0.0;