From 9dc885d2974aa9c054e8baf54630389005bc78e7 Mon Sep 17 00:00:00 2001 From: Daniel Richtmann Date: Thu, 18 Jan 2018 17:02:04 +0100 Subject: [PATCH] Fix a bug in Wilson MG The calculation of the lattice size of a second coarse level was incorrect. --- tests/solver/Test_wilson_ddalphaamg.cc | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/tests/solver/Test_wilson_ddalphaamg.cc b/tests/solver/Test_wilson_ddalphaamg.cc index 439df423..95297deb 100644 --- a/tests/solver/Test_wilson_ddalphaamg.cc +++ b/tests/solver/Test_wilson_ddalphaamg.cc @@ -113,7 +113,8 @@ public: std::cout << GridLogMessage << "Constructing " << coarsegrids << " CoarseGrids" << std::endl; for(int cl = 0; cl < coarsegrids; ++cl) { // may be a bit ugly and slow but not perf critical - LattSizes.push_back({GridDefaultLatt()}); + // need to differentiate between first and other coarse levels in size calculation + LattSizes.push_back({cl == 0 ? GridDefaultLatt() : LattSizes[cl - 1]}); Seeds.push_back(std::vector(LattSizes[cl].size())); for(int d = 0; d < LattSizes[cl].size(); ++d) {