diff --git a/EXAMPLES/MATRIX_MARKET/arpackSolver.hpp b/EXAMPLES/MATRIX_MARKET/arpackSolver.hpp index d8fc31f..34fcb14 100644 --- a/EXAMPLES/MATRIX_MARKET/arpackSolver.hpp +++ b/EXAMPLES/MATRIX_MARKET/arpackSolver.hpp @@ -229,6 +229,7 @@ class arpackSolver { int solve(EM & A, EM const * B = nullptr) { stdPb = !B ? true : false; + if (nbCV > A.cols()) nbCV = A.cols(); // Cut-off arpack workspace dim. dumpParameters(); if (verbose == 3) { diff --git a/EXAMPLES/MATRIX_MARKET/arpackmm.cpp b/EXAMPLES/MATRIX_MARKET/arpackmm.cpp index c44c4be..ae4b282 100644 --- a/EXAMPLES/MATRIX_MARKET/arpackmm.cpp +++ b/EXAMPLES/MATRIX_MARKET/arpackmm.cpp @@ -622,8 +622,6 @@ int itrSolve(options& opt, output& out, double const& slvItrILUDropTol, cout << endl; cout << "INP: create A " << readATime << " s" << endl; - if (opt.nbCV > A.cols()) opt.nbCV = A.cols(); /* Cut-off */ - EM B; if (!opt.stdPb) { start = chrono::high_resolution_clock::now(); @@ -711,8 +709,6 @@ int drtSolve(options& opt, output& out) { cout << endl; cout << "INP: create A " << readATime << " s" << endl; - if (opt.nbCV > A.cols()) opt.nbCV = A.cols(); /* Cut-off */ - EM B; if (!opt.stdPb) { start = chrono::high_resolution_clock::now();