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

adding C1M and restructuring solvers #1

Merged
portelli merged 5 commits from mdicarlo/HadronsPresets:main into main 2023-06-01 16:39:14 +01:00
Conversation 6 Commits 5 Files Changed 1 +8 -14
1 changed files with 8 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files
Showing only changes of commit ab747d2b1d - Show all commits

22
RbcUkqcd.hpp
View File

@ -62,11 +62,6 @@ struct RbcUkqcd
inline static constexpr DeflationParameters c1m32DeflPar{3.0e-06, 5.5, 101, 100, 110, 120};
// Light solver: load deflation from disk
static inline void addLightLCDSolver(Application &app, const RbcUkqcd::EnsembleParameters &par,
const std::string solverName, const std::string gaugeName,
const std::string gaugeTransform,
const std::string eigenpackPath, const double residual);
static inline void addM0LightLCDSolver(Application &app, const std::string solverName,
const std::string gaugeName,
const std::string gaugeTransform,
@ -125,11 +120,10 @@ struct RbcUkqcd
// Implementations /////////////////////////////////////////////////////////////////////////////////
// Load deflation from disk
void RbcUkqcd::addLightLCDSolver(Application &app, const RbcUkqcd::EnsembleParameters &par,
const std::string solverName, const std::string gaugeName,
const std::string gaugeTransform, const std::string eigenpackPath,
const double residual)
// Light M0 (load deflation from disk)
void RbcUkqcd::addM0LightLCDSolver(Application &app, const std::string solverName,
const std::string gaugeName, const std::string gaugeTransform,
const std::string eigenpackPath, const double residual)
{
const std::string prefix = solverName;
@ -143,10 +137,10 @@ void RbcUkqcd::addLightLCDSolver(Application &app, const RbcUkqcd::EnsembleParam
MAction::ScaledDWF::Par actionPar;
actionPar.gauge = gaugeName;
actionPar.Ls = par.Ls;
actionPar.M5 = par.M5;
actionPar.mass = par.ml;
actionPar.scale = par.scale;
actionPar.Ls = RbcUkqcd::m0LCDPar.Ls;
actionPar.M5 = RbcUkqcd::m0LCDPar.M5;
actionPar.mass = RbcUkqcd::m0LCDPar.ml;
actionPar.scale = RbcUkqcd::m0LCDPar.scale;
actionPar.boundary = "1 1 1 1";
actionPar.twist = "0. 0. 0. 0.";
app.createModule<MAction::ScaledDWF>(prefix + "_dwf", actionPar);