More debugging and a performance check

1 files changed, 1 insertions, 1 deletions

diff --git a/src/lapack.cpp b/src/lapack.cpp
index 125e5a0..79d49fd 100644
--- a/src/lapack.cpp
+++ b/src/lapack.cpp
@@ -319,7 +319,7 @@ void LUDecomp(gsl_matrix *LU, gsl_permutation *p, int *signum) {
 void LUInvert(const gsl_matrix *LU, const gsl_permutation *p, gsl_matrix *ret_inverse) {
   // debug_msg("entering");
   auto det = LULndet(LU);
-  assert(det != 1.0);
+  enforce_msg(det != 1.0,"LU determinant is zero -> LU is not invertable");
 
   enforce_gsl(gsl_linalg_LU_invert(LU, p, ret_inverse));
 }