mdicarlo/HadronsPresets
1
0
Fork 1
forked from portelli/HadronsPresets
Code Issues Pull Requests Packages Projects Releases Wiki Activity

try double prec CG for M0 charm

Browse Source
This commit is contained in:
Antonin Portelli 2023-06-10 15:44:18 +01:00
parent 3936ed389e
commit 3c30d73912
1 changed files with 6 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
RbcUkqcd.hpp
View File

@ -409,13 +409,6 @@ void RbcUkqcd::addM0CharmSolver(Application &app, const std::string solverName,
smearPar.orthogDim = ""; smearPar.orthogDim = "";
app.createModule<MGauge::StoutSmearing>(prefix + "_gauge_3stout", smearPar); app.createModule<MGauge::StoutSmearing>(prefix + "_gauge_3stout", smearPar);
// Gauge field FP32 cast
MUtilities::GaugeSinglePrecisionCast::Par gaugeCastPar;
gaugeCastPar.field = prefix + "_gauge_3stout";
app.createModule<MUtilities::GaugeSinglePrecisionCast>(prefix + "_gauge_3stout_fp32",
gaugeCastPar);
// Scaled DWF action + FP32 version // Scaled DWF action + FP32 version
MAction::ScaledDWF::Par actionPar; MAction::ScaledDWF::Par actionPar;
@ -427,20 +420,15 @@ void RbcUkqcd::addM0CharmSolver(Application &app, const std::string solverName,
actionPar.boundary = boundary; actionPar.boundary = boundary;
actionPar.twist = "0. 0. 0. 0."; actionPar.twist = "0. 0. 0. 0.";
app.createModule<MAction::ScaledDWF>(prefix + "_dwf", actionPar); app.createModule<MAction::ScaledDWF>(prefix + "_dwf", actionPar);
actionPar.gauge = prefix + "_gauge_3stout_fp32";
app.createModule<MAction::ScaledDWFF>(prefix + "_dwf_fp32", actionPar);
// Mixed-precision red-black preconditionned CG // Red-black preconditionned CG
MSolver::MixedPrecisionRBPrecCG::Par solverPar; MSolver::RBPrecCG::Par solverPar;
solverPar.innerAction = prefix + "_dwf_fp32"; solverPar.action = prefix + "_dwf";
solverPar.outerAction = prefix + "_dwf"; solverPar.maxIteration = 30000;
solverPar.maxInnerIteration = 30000;
solverPar.maxOuterIteration = 100;
solverPar.residual = residual; solverPar.residual = residual;
solverPar.innerGuesser = ""; solverPar.guesser = "";
solverPar.outerGuesser = ""; app.createModule<MSolver::RBPrecCG>(solverName, solverPar);
app.createModule<MSolver::MixedPrecisionRBPrecCG>(solverName, solverPar);
} }
} // namespace hadpresets } // namespace hadpresets