From a24297bb1cb686d3630db3a0da8fa8a2079ad178 Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Sat, 25 Aug 2018 10:18:41 +0000
Subject: Added remarks

---
 src/lapack.cpp | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

(limited to 'src/lapack.cpp')

diff --git a/src/lapack.cpp b/src/lapack.cpp
index b309a1b..0fc2de4 100644
--- a/src/lapack.cpp
+++ b/src/lapack.cpp
@@ -193,7 +193,9 @@ void lapack_eigen_symmv(gsl_matrix *A, gsl_vector *eval, gsl_matrix *evec,
     double WORK_temp[1];
     int IWORK_temp[1];
 
-    if (is_check_mode()) disable_segfpe(); // disable fast NaN checking for now
+    // disable fast NaN checking for now - dsyevr throws NaN errors,
+    // but fixes them (apparently)
+    if (is_check_mode()) disable_segfpe();
 
     // DSYEVR - computes selected eigenvalues and, optionally,
     // eigenvectors of a real symmetric matrix
@@ -208,8 +210,8 @@ void lapack_eigen_symmv(gsl_matrix *A, gsl_vector *eval, gsl_matrix *evec,
 
     if (INFO != 0) cerr << "ERROR: value of INFO is " << INFO;
     enforce_msg(INFO == 0, "lapack_eigen_symmv failed");
-    LWORK = (int)WORK_temp[0];
-    LIWORK = (int)IWORK_temp[0];
+    LWORK = (int)WORK_temp[0];    // The dimension of the array work.
+    LIWORK = (int)IWORK_temp[0];  // The dimension of the array iwork, lwork≥ max(1, 10n).
 
     double *WORK = new double[LWORK];
     int *IWORK = new int[LIWORK];
-- 
cgit v1.2.3