From 194e4b94bb03a91060e0c6049092d6a97c07ef50 Mon Sep 17 00:00:00 2001 From: Daniel Richtmann Date: Mon, 29 Jan 2018 17:18:20 +0100 Subject: [PATCH] Make MG checking function work level-wise --- tests/solver/Test_wilson_ddalphaamg.cc | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/tests/solver/Test_wilson_ddalphaamg.cc b/tests/solver/Test_wilson_ddalphaamg.cc index 95297deb..3ba933e8 100644 --- a/tests/solver/Test_wilson_ddalphaamg.cc +++ b/tests/solver/Test_wilson_ddalphaamg.cc @@ -539,7 +539,7 @@ public: std::cout << GridLogMessage << "SAP resid(post) " << std::sqrt(r / Ni) << " " << r << " " << Ni << std::endl; } - void runChecks(CoarseGrids &cGrids) { + void runChecks(CoarseGrids &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 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;