From cc35c4425698f8a9b8293af0c02dfc385442ad6b Mon Sep 17 00:00:00 2001
From: Desire NUENTSA <desire.nuentsa_wakam@inria.fr>
Date: Mon, 25 Feb 2013 14:26:55 +0100
Subject: Add reference for the default threshold in sparse QR

---
 Eigen/src/SparseQR/SparseQR.h | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

(limited to 'Eigen/src/SparseQR')

diff --git a/Eigen/src/SparseQR/SparseQR.h b/Eigen/src/SparseQR/SparseQR.h
index 8fc0a7c3c..b3d5cd208 100644
--- a/Eigen/src/SparseQR/SparseQR.h
+++ b/Eigen/src/SparseQR/SparseQR.h
@@ -286,7 +286,10 @@ void SparseQR<MatrixType,OrderingType>::factorize(const MatrixType& mat)
     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) 
   {
     RealScalar max2Norm = 0.0;
-- 
cgit v1.2.3