portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

WilsonMG: Read in MG parameters from xml in test

Browse Source
This commit is contained in:
Daniel Richtmann 2018-04-03 16:03:11 +02:00
parent ff6413a764
commit 57a49ed22f
No known key found for this signature in database
GPG Key ID: B33C490AF0772057
1 changed files with 16 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
tests/solver/Test_wilsonclover_mg.cc
View File

@ -520,7 +520,7 @@ public:
auto coarseSolverMaxIter = _MultiGridParams.coarseSolverMaxOuterIter * _MultiGridParams.coarseSolverMaxInnerIter; auto coarseSolverMaxIter = _MultiGridParams.coarseSolverMaxOuterIter * _MultiGridParams.coarseSolverMaxInnerIter;
// On the coarsest level we only have a fine what I above call the fine level, no coarse one // On the coarsest level we only have what I above call the fine level, no coarse one
TrivialPrecon<FineVector> fineTrivialPreconditioner; TrivialPrecon<FineVector> fineTrivialPreconditioner;
FlexibleGeneralisedMinimalResidual<FineVector> fineFGMRES( FlexibleGeneralisedMinimalResidual<FineVector> fineFGMRES(
_MultiGridParams.coarseSolverTol, coarseSolverMaxIter, fineTrivialPreconditioner, _MultiGridParams.coarseSolverMaxInnerIter, false); _MultiGridParams.coarseSolverTol, coarseSolverMaxIter, fineTrivialPreconditioner, _MultiGridParams.coarseSolverMaxInnerIter, false);
@ -592,7 +592,7 @@ int main(int argc, char **argv) {
WilsonFermionR Dw(Umu, *FGrid, *FrbGrid, mass); WilsonFermionR Dw(Umu, *FGrid, *FrbGrid, mass);
WilsonCloverFermionR Dwc(Umu, *FGrid, *FrbGrid, mass, csw_r, csw_t, wilsonAnisCoeff, wcImplparams); WilsonCloverFermionR Dwc(Umu, *FGrid, *FrbGrid, mass, csw_r, csw_t, wilsonAnisCoeff, wcImplparams);
// Params for two-level MG preconditioner // Default params for two-level MG preconditioner (TODO: use sensible ones)
mgParams.nLevels = 2; mgParams.nLevels = 2;
mgParams.blockSizes = {{2, 2, 2, 2}}; mgParams.blockSizes = {{2, 2, 2, 2}};
mgParams.smootherTol = {1e-14}; mgParams.smootherTol = {1e-14};
@ -606,7 +606,7 @@ int main(int argc, char **argv) {
mgParams.coarseSolverMaxOuterIter = 1; mgParams.coarseSolverMaxOuterIter = 1;
mgParams.coarseSolverMaxInnerIter = 1; mgParams.coarseSolverMaxInnerIter = 1;
// // Params for three-level MG preconditioner // // Default params for three-level MG preconditioner (TODO: use sensible ones)
// mgParams.nLevels = 3; // mgParams.nLevels = 3;
// mgParams.blockSizes = {{2, 2, 2, 2}, {2, 2, 1, 1}}; // mgParams.blockSizes = {{2, 2, 2, 2}, {2, 2, 1, 1}};
// mgParams.smootherTol = {1e-14, 1e-14}; // mgParams.smootherTol = {1e-14, 1e-14};
@ -620,7 +620,7 @@ int main(int argc, char **argv) {
// mgParams.coarseSolverMaxOuterIter = 1; // mgParams.coarseSolverMaxOuterIter = 1;
// mgParams.coarseSolverMaxInnerIter = 1; // mgParams.coarseSolverMaxInnerIter = 1;
// // // Params for four-level MG preconditioner // // Default params for four-level MG preconditioner (TODO: use sensible ones)
// mgParams.nLevels = 4; // mgParams.nLevels = 4;
// mgParams.blockSizes = {{2, 2, 2, 2}, {2, 2, 1, 1}, {1, 1, 2, 1}}; // mgParams.blockSizes = {{2, 2, 2, 2}, {2, 2, 1, 1}, {1, 1, 2, 1}};
// mgParams.smootherTol = {1e-14, 1e-14, 1e-14}; // mgParams.smootherTol = {1e-14, 1e-14, 1e-14};
@ -634,6 +634,18 @@ int main(int argc, char **argv) {
// mgParams.coarseSolverMaxOuterIter = 1; // mgParams.coarseSolverMaxOuterIter = 1;
// mgParams.coarseSolverMaxInnerIter = 1; // mgParams.coarseSolverMaxInnerIter = 1;
{
XmlWriter writer("mg_params_template.xml");
write(writer, "Params", mgParams);
std::cout << GridLogMessage << " Written mg_params_template.xml" << std::endl;
}
{
std::string paramFileName{"./mg_params.xml"};
XmlReader reader(paramFileName);
read(reader, "Params", mgParams);
}
checkParameterValidity(mgParams); checkParameterValidity(mgParams);
std::cout << mgParams << std::endl; std::cout << mgParams << std::endl;