portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-11-04 05:54:32 +00:00
Code Releases Activity

Make MG checking function work level-wise

Browse Source
This commit is contained in:
Daniel Richtmann
2018-01-29 17:18:20 +01:00
parent bfc1411c1f
commit 194e4b94bb
1 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
tests/solver/Test_wilson_ddalphaamg.cc
View File

@@ -539,7 +539,7 @@ public:
std::cout << GridLogMessage << "SAP resid(post) " << std::sqrt(r / Ni) << " " << r << " " << Ni << std::endl;
}
void runChecks(CoarseGrids<nbasis> &cGrids) {
void runChecks(CoarseGrids<nbasis> &cGrids, int whichCoarseGrid) {
/////////////////////////////////////////////
// Some stuff we need for the checks below //
@@ -578,7 +578,7 @@ public:
std::cout << GridLogMessage << "MG correctness check: 0 == (1 - R P) v_c" << std::endl;
std::cout << GridLogMessage << "**************************************************" << std::endl;
random(cGrids.PRNGs[0], cTmps[0]);
random(cGrids.PRNGs[whichCoarseGrid], cTmps[0]);
_Aggregates.PromoteFromSubspace(cTmps[0], fTmps[0]); // P v_c
_Aggregates.ProjectToSubspace(cTmps[1], fTmps[0]); // R P v_c
@@ -599,7 +599,7 @@ public:
std::cout << GridLogMessage << "MG correctness check: 0 == (R D P - D_c) v_c" << std::endl;
std::cout << GridLogMessage << "**************************************************" << std::endl;
random(cGrids.PRNGs[0], cTmps[0]);
random(cGrids.PRNGs[whichCoarseGrid], cTmps[0]);
_Aggregates.PromoteFromSubspace(cTmps[0], fTmps[0]); // P v_c
_FineOperator.Op(fTmps[0], fTmps[1]); // D P v_c
@@ -623,7 +623,7 @@ public:
std::cout << GridLogMessage << "MG correctness check: 0 == |(Im(v_c^dag D_c^dag D_c v_c)|" << std::endl;
std::cout << GridLogMessage << "**************************************************" << std::endl;
random(cGrids.PRNGs[0], cTmps[0]);
random(cGrids.PRNGs[whichCoarseGrid], cTmps[0]);
MdagMOp.Op(cTmps[0], cTmps[1]); // D_c v_c
MdagMOp.AdjOp(cTmps[1], cTmps[2]); // D_c^dag D_c v_c
@@ -923,7 +923,7 @@ int main(int argc, char **argv) {
// FineOperator &Smooth,Matrix &SmootherMatrix)
TrivialPrecon<LatticeFermion> Simple;
Precon.runChecks(cGrids);
Precon.runChecks(cGrids, 0);
std::cout << GridLogMessage << "**************************************************" << std::endl;
std::cout << GridLogMessage << "Building two level VPGCR and FGMRES solvers" << std::endl;