diff --git a/lib/algorithms/iterative/ConjugateGradient.h b/lib/algorithms/iterative/ConjugateGradient.h index 45236e2c..acc202e7 100644 --- a/lib/algorithms/iterative/ConjugateGradient.h +++ b/lib/algorithms/iterative/ConjugateGradient.h @@ -161,7 +161,7 @@ class ConjugateGradient : public OperatorFunction { axpy(r, -a, mmp, r);// new residual r = r_old - a * Ap - cp = norm2(r, ReprTest); // bookmarking this norm + cp = norm2(r, ReprTest); // bookkeeping this norm if (ReproTest && !CGState.do_repro) { CGState.residuals.push_back(cp); // save residuals state std::cout << GridLogIterative << "ReproTest: Saving state" << std::endl;