From 86a002ae27171a3922d4bd9e7b46ff0df95c51ed Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Fri, 31 Aug 2018 12:00:36 +0000
Subject: Continue debugging calcpab

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

(limited to 'src/lapack.cpp')

diff --git a/src/lapack.cpp b/src/lapack.cpp
index 6e43bd9..165a82d 100644
--- a/src/lapack.cpp
+++ b/src/lapack.cpp
@@ -239,6 +239,8 @@ void lapack_eigen_symmv(gsl_matrix *A, gsl_vector *eval, gsl_matrix *evec,
 double EigenDecomp(gsl_matrix *G, gsl_matrix *U, gsl_vector *eval,
                    const size_t flag_largematrix) {
   lapack_eigen_symmv(G, eval, U, flag_largematrix);
+  assert(!has_nan(eval));
+  // write(eval,"eval");
 
   // Calculate track_G=mean(diag(G)).
   double d = 0.0;
-- 
cgit v1.2.3