portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-06-15 22:37:05 +01:00
Code Releases Activity

Just commenting

Browse Source
This commit is contained in:
Guido Cossu
2016-12-12 03:49:01 +00:00
parent 85c055fa30
commit 70068cff51
1 changed files with 1 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
lib/algorithms/iterative/ConjugateGradient.h
View File

@ -161,7 +161,7 @@ class ConjugateGradient : public OperatorFunction<Field> {
axpy(r, -a, mmp, r);// new residual r = r_old - a * Ap 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) { if (ReproTest && !CGState.do_repro) {
CGState.residuals.push_back(cp); // save residuals state CGState.residuals.push_back(cp); // save residuals state
std::cout << GridLogIterative << "ReproTest: Saving state" << std::endl; std::cout << GridLogIterative << "ReproTest: Saving state" << std::endl;