portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-11-04 05:54:32 +00:00
Code Releases Activity

To HMC/Mobius2p1fIDSDRGparityEOFA_40ID, added input param to change trajectory length and increased integrator steps for DSDR

Browse Source
This commit is contained in:
Christopher Kelly
2021-12-10 09:06:06 -08:00
parent ddf7540510
commit 551b93ba8e
1 changed files with 5 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
HMC/Mobius2p1fIDSDRGparityEOFA_40ID.cc
View File

@@ -95,6 +95,7 @@ struct EvolParameters: Serializable {
Integer, Trajectories, Integer, Trajectories,
Integer, SaveInterval, Integer, SaveInterval,
Integer, Steps, Integer, Steps,
RealD, TrajectoryLength,
bool, MetropolisTest, bool, MetropolisTest,
std::string, StartingType, std::string, StartingType,
std::vector<Integer>, GparityDirs, std::vector<Integer>, GparityDirs,
@@ -111,6 +112,7 @@ struct EvolParameters: Serializable {
StartingType = "ColdStart"; StartingType = "ColdStart";
GparityDirs.resize(3, 1); //1 for G-parity, 0 for periodic GparityDirs.resize(3, 1); //1 for G-parity, 0 for periodic
Steps = 5; Steps = 5;
TrajectoryLength = 1.0;
} }
}; };
@@ -560,7 +562,7 @@ int main(int argc, char **argv) {
// MD.name = std::string("ForceGradient"); // MD.name = std::string("ForceGradient");
MD.MDsteps = user_params.Steps; MD.MDsteps = user_params.Steps;
MD.trajL = 1.0; MD.trajL = user_params.TrajectoryLength;
typedef HMCWrapper::ImplPolicy GaugeImplPolicy; typedef HMCWrapper::ImplPolicy GaugeImplPolicy;
@@ -651,8 +653,8 @@ int main(int argc, char **argv) {
// Collect actions // Collect actions
//////////////////////////////////// ////////////////////////////////////
ActionLevel<HMCWrapper::Field> Level1(1); //light quark + strange quark ActionLevel<HMCWrapper::Field> Level1(1); //light quark + strange quark
ActionLevel<HMCWrapper::Field> Level2(2); //DSDR ActionLevel<HMCWrapper::Field> Level2(4); //DSDR
ActionLevel<HMCWrapper::Field> Level3(4); //gauge ActionLevel<HMCWrapper::Field> Level3(2); //gauge
///////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////