[BUG FIX] arpack may succeed to compute more EV than expected. (#197)

opencollab · Feb 13, 2019 · dfc8abd · dfc8abd
1 parent a98f855
commit dfc8abd
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 4 deletions.
diff --git a/PARPACK/TESTS/MPI/icb_parpack_cpp.cpp b/PARPACK/TESTS/MPI/icb_parpack_cpp.cpp
@@ -75,7 +75,8 @@ void real_symmetric_runner() {
                                    &(workd[ipntr[1] - 1]));
   }
   // check number of ev found by arpack.
-  if (iparam[4] != nev || info != 0) {
+  if (iparam[4] < nev /*arpack may succeed to compute more EV than expected*/ || info != 0) {
+    std::cout << "ERROR: iparam[4] " << iparam[4] << ", nev " << nev << ", info " << info << std::endl;
     throw std::domain_error("Error inside ARPACK routines");
   }
 
@@ -153,7 +154,8 @@ void complex_symmetric_runner() {
   }
 
   // check number of ev found by arpack
-  if (iparam[4] != nev || info != 0) {
+  if (iparam[4] < nev /*arpack may succeed to compute more EV than expected*/ || info != 0) {
+    std::cout << "ERROR: iparam[4] " << iparam[4] << ", nev " << nev << ", info " << info << std::endl;
     throw std::domain_error("Error inside ARPACK routines");
   }
 

diff --git a/TESTS/icb_arpack_cpp.cpp b/TESTS/icb_arpack_cpp.cpp
@@ -73,7 +73,8 @@ void real_symmetric_runner() {
   }
 
   // check number of ev found by arpack.
-  if (iparam[4] != nev || info != 0) {
+  if (iparam[4] < nev /*arpack may succeed to compute more EV than expected*/ || info != 0) {
+    std::cout << "ERROR: iparam[4] " << iparam[4] << ", nev " << nev << ", info " << info << std::endl;
     throw std::domain_error("Error inside ARPACK routines");
   }
 
@@ -151,7 +152,8 @@ void complex_symmetric_runner() {
   }
 
   // check number of ev found by arpack.
-  if (iparam[4] != nev || info != 0) {
+  if (iparam[4] < nev /*arpack may succeed to compute more EV than expected*/ || info != 0) {
+    std::cout << "ERROR: iparam[4] " << iparam[4] << ", nev " << nev << ", info " << info << std::endl;
     throw std::domain_error("Error inside ARPACK routines");
   }