diff --git a/RbcUkqcd.hpp b/RbcUkqcd.hpp index 7c74e94..205b28b 100644 --- a/RbcUkqcd.hpp +++ b/RbcUkqcd.hpp @@ -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(prefix + "_dwf", actionPar);