From c08e633d2f858b35d2939d92c8a1b82d36168944 Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Thu, 27 Sep 2018 09:35:13 +0000
Subject: Inlining GSL calls

---
 src/lapack.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

(limited to 'src/lapack.cpp')

diff --git a/src/lapack.cpp b/src/lapack.cpp
index bf73338..eb5b16b 100644
--- a/src/lapack.cpp
+++ b/src/lapack.cpp
@@ -319,8 +319,8 @@ void LUDecomp(gsl_matrix *LU, gsl_permutation *p, int *signum) {
 // introductory textbook on numerical linear algebra for details).
 void LUInvert(const gsl_matrix *LU, const gsl_permutation *p, gsl_matrix *ret_inverse) {
   // debug_msg("entering");
-  auto det = LULndet(LU);
-  enforce_msg(det != 1.0,"LU determinant is zero -> LU is not invertable");
+  if (is_check_mode())
+    LULndet(LU);
 
   enforce_gsl(gsl_linalg_LU_invert(LU, p, ret_inverse));
 }
-- 
cgit v1.2.3