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

Tests more reductions and stops if failure is found in an iteration

Browse Source
This commit is contained in:
Guido Cossu
2016-12-12 03:33:21 +00:00
parent 90fedbd2af
commit 85c055fa30
1 changed files with 10 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -102,24 +102,24 @@ class ConjugateGradient : public OperatorFunction<Field> {
std::cout << GridLogMessage << "Starting reproducibility test, full check every "
<< ReprTest.interval << " calls" << std::endl;
// Initial residual computation & set up
RealD guess = norm2(psi);
assert(std::isnan(guess) == 0);
Linop.HermOpAndNorm(psi, mmp, d, b);
if(!ReprTest.do_check)
ReprTest.reset();
ReprTest.enable_reprocheck=ReproTest;
// Initial residual computation & set up
RealD guess = norm2(psi, ReprTest);
assert(std::isnan(guess) == 0);
Linop.HermOpAndNorm(psi, mmp, d, b);// eventually split this for the norm check
r = src - mmp;
p = r;
a = norm2(p);
a = norm2(p, ReprTest);
cp = a;
ssq = norm2(src);
ssq = norm2(src, ReprTest);
std::cout << GridLogIterative << "ConjugateGradient: guess " << guess << std::endl;
std::cout << GridLogIterative << "ConjugateGradient: src " << ssq << std::endl;
@ -174,7 +174,7 @@ class ConjugateGradient : public OperatorFunction<Field> {
std::cout << GridLogMessage << " at k=" << k << std::endl;
std::cout << GridLogMessage << "saved residual = " << CGState.residuals[k-1]
<< " cp = " << cp << std::endl;
//exit(-1); do not stop, report all the failures
exit(1); // exit after the first failure
}
}
b = cp / c;