diff --git a/AMDiS/src/CGSolver.hh b/AMDiS/src/CGSolver.hh index 00c001cea2eda6ed1b20c39846ccc643c85eac32..821204983148a1abf8bb3415768d596537e9f92e 100644 --- a/AMDiS/src/CGSolver.hh +++ b/AMDiS/src/CGSolver.hh @@ -108,15 +108,13 @@ namespace AMDiS { // Check if x is a solution of the linear system. if (SOLVE_INFO(iter, sqrt(alpha1), &old_res)) { - *p = *x; - *p *= -1.0; - - // r = b - Ax - matVec->matVec(NoTranspose, *p, *r); - *r += *b; - - ::std::cout << "Final res: " << norm(r) << ::std::endl; - + *p = *x; + *p *= -1.0; + + // r = b - Ax + matVec->matVec(NoTranspose, *p, *r); + *r += *b; + return(iter); }