rewrite the test with smaller tol

opencollab · Sep 17, 2024 · f267553 · f267553
1 parent 6183acf
commit f267553
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 5 deletions.
diff --git a/EXAMPLES/MATRIX_MARKET/arpackSolver.hpp b/EXAMPLES/MATRIX_MARKET/arpackSolver.hpp
@@ -297,7 +297,7 @@ class arpackSolver {
       return 0;
     };
 
-    int checkEigVec(EM const & A, EM const * B = nullptr, double const maxResNorm = 1.) const {
+    int checkEigVec(EM const & A, EM const * B = nullptr, double const maxResNorm = 1.e-3) const {
       // Check eigen vectors.
 
       string rs = schur ? "Schur" : "Ritz";
@@ -324,7 +324,7 @@ class arpackSolver {
         right *= lambda;
         EigVecZ residual = left - right;
         if (residual.norm() > maxResNorm) {
-          cerr << endl << "Error: bad eigen value " << i << ":" << endl;
+          cerr << endl << "Error: bad eigen value " << i << " (norm " << std::norm(lambda) << "):" << endl;
           cerr << endl << lambda << endl;
           cerr << endl << "Error: bad eigen vector " << i << " (norm " << V.norm() << "):" << endl;
           cerr << endl << V << endl;

diff --git a/EXAMPLES/MATRIX_MARKET/arpackmm.cpp b/EXAMPLES/MATRIX_MARKET/arpackmm.cpp
@@ -31,7 +31,7 @@ class options {
     invert =
         false;  // Eigen value invertion: look for 1./lambda instead of lambda.
     tol = 1.e-06;
-    maxResNorm = 1.;
+    maxResNorm = 1.e-3;
     maxIt = 100;
     schur = false;  // Compute Ritz vectors.
     slv = "BiCG";
@@ -439,7 +439,7 @@ class options {
     cout << "  --tol T:          tolerance T." << endl;
     cout << "                    default: 1.e-06" << endl;
     cout << "  --maxResNorm R:   maximum residual norm R." << endl;
-    cout << "                    default: 1." << endl;
+    cout << "                    default: 1.e-3" << endl;
     cout << "  --maxIt M:        maximum iterations M." << endl;
     cout << "                    default: 100" << endl;
     cout << "  --schur:          compute Schur vectors." << endl;

diff --git a/EXAMPLES/MATRIX_MARKET/arpackmm.sh b/EXAMPLES/MATRIX_MARKET/arpackmm.sh
@@ -26,7 +26,7 @@ do
     # Note: this is expected as "power-like" methods are designed to find largest eigen values.
     for magOpt in "--mag LM" "--mag SM --noCheck"
     do
-      export mrn="--maxResNorm 1." # Relax residual check to get stable tests.
+      export mrn="--maxResNorm 1.e-2" # Relax residual check to get stable tests.
 
       # Shift slightly to avoid the zero-vector starting problem.
       export shiftZV=""