Adding note because of mvlmm performance test failure

1 files changed, 6 insertions, 0 deletions

diff --git a/src/lapack.cpp b/src/lapack.cpp
index 31b6add..e1a63e1 100644
--- a/src/lapack.cpp
+++ b/src/lapack.cpp
@@ -304,6 +304,12 @@ double CholeskySolve(gsl_matrix *Omega, gsl_vector *Xty, gsl_vector *OiXty) {
 }
 
 // LU decomposition.
+//
+// The functions return GSL_SUCCESS for non-singular matrices. If the
+// matrix is singular, the factorization is still completed: U is
+// singular and the decomposition should not be used to solve linear
+// systems
+
 void LUDecomp(gsl_matrix *LU, gsl_permutation *p, int *signum) {
   // debug_msg("entering");
   enforce_gsl(gsl_linalg_LU_decomp(LU, p, signum));