diff options
| -rw-r--r-- | src/lapack.cpp | 6 |
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)); |