From cb285681989768016f86e91492230d75f487ba90 Mon Sep 17 00:00:00 2001 From: Peter Boyle Date: Thu, 6 May 2021 23:12:57 +0200 Subject: [PATCH] Tuning integrator --- HMC/Mobius2p1fEOFA.cc | 20 ++++++++++---------- 1 file changed, 10 insertions(+), 10 deletions(-) diff --git a/HMC/Mobius2p1fEOFA.cc b/HMC/Mobius2p1fEOFA.cc index b1294da5..68a3a5b5 100644 --- a/HMC/Mobius2p1fEOFA.cc +++ b/HMC/Mobius2p1fEOFA.cc @@ -188,15 +188,15 @@ int main(int argc, char **argv) { IntegratorParameters MD; // typedef GenericHMCRunner HMCWrapper; // MD.name = std::string("Leap Frog"); - typedef GenericHMCRunner HMCWrapper; - MD.name = std::string("Force Gradient"); - // typedef GenericHMCRunner HMCWrapper; - // MD.name = std::string("MinimumNorm2"); - MD.MDsteps = 6; + // typedef GenericHMCRunner HMCWrapper; + // MD.name = std::string("Force Gradient"); + typedef GenericHMCRunner HMCWrapper; + MD.name = std::string("MinimumNorm2"); + MD.MDsteps = 12; MD.trajL = 1.0; HMCparameters HMCparams; - HMCparams.StartTrajectory = 590; + HMCparams.StartTrajectory = 211; HMCparams.Trajectories = 1000; HMCparams.NoMetropolisUntil= 0; // "[HotStart, ColdStart, TepidStart, CheckpointStart]\n"; @@ -209,9 +209,9 @@ int main(int argc, char **argv) { TheHMC.Resources.AddFourDimGrid("gauge"); // use default simd lanes decomposition CheckpointerParameters CPparams; - CPparams.config_prefix = "ckpoint_EODWF_lat"; - CPparams.rng_prefix = "ckpoint_EODWF_rng"; - CPparams.saveInterval = 10; + CPparams.config_prefix = "ckpoint_EOFA_lat"; + CPparams.rng_prefix = "ckpoint_EOFA_rng"; + CPparams.saveInterval = 1; CPparams.format = "IEEE64BIG"; TheHMC.Resources.LoadNerscCheckpointer(CPparams); @@ -263,7 +263,7 @@ int main(int argc, char **argv) { FermionActionF::ImplParams ParamsF(boundary); double ActionStoppingCondition = 1e-10; - double DerivativeStoppingCondition = 1e-6; + double DerivativeStoppingCondition = 1e-8; double MaxCGIterations = 30000; ////////////////////////////////////