portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-04 19:25:56 +01:00
Code Releases Activity

To DWF+I G-parity evolution code, added ability to specify number of MD steps in params and an optional usage mode that reads the config and checks the plaq/checksum agree then exits

Browse Source
This commit is contained in:
Christopher Kelly 2021-02-16 10:41:52 -05:00
parent daa095c519
commit e0f6a146d8
1 changed files with 16 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
HMC/DWF2p1fIwasakiGparity.cc
View File

@ -72,6 +72,7 @@ struct EvolParameters: Serializable {
Integer, StartTrajectory, Integer, StartTrajectory,
Integer, Trajectories, Integer, Trajectories,
Integer, SaveInterval, Integer, SaveInterval,
Integer, Steps,
bool, MetropolisTest, bool, MetropolisTest,
std::string, StartingType, std::string, StartingType,
std::vector<Integer>, GparityDirs, std::vector<Integer>, GparityDirs,
@ -86,6 +87,7 @@ struct EvolParameters: Serializable {
SaveInterval = 5; SaveInterval = 5;
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;
} }
}; };
@ -237,11 +239,14 @@ int main(int argc, char **argv) {
std::cout << GridLogMessage << "Grid is setup to use " << threads << " threads" << std::endl; std::cout << GridLogMessage << "Grid is setup to use " << threads << " threads" << std::endl;
std::string param_file = "params.xml"; std::string param_file = "params.xml";
bool file_load_check = false;
for(int i=1;i<argc;i++){ for(int i=1;i<argc;i++){
if(std::string(argv[i]) == "--param_file"){ std::string sarg(argv[i]);
if(sarg == "--param_file"){
assert(i!=argc-1); assert(i!=argc-1);
param_file = argv[i+1]; param_file = argv[i+1];
break; }else if(sarg == "--read_check"){ //check the fields load correctly and pass checksum/plaquette repro
file_load_check = true;
} }
} }
@ -291,7 +296,7 @@ int main(int argc, char **argv) {
typedef ConjugateHMCRunnerD<MinimumNorm2> HMCWrapper; //NB: This is the "Omelyan integrator" typedef ConjugateHMCRunnerD<MinimumNorm2> HMCWrapper; //NB: This is the "Omelyan integrator"
typedef HMCWrapper::ImplPolicy GaugeImplPolicy; typedef HMCWrapper::ImplPolicy GaugeImplPolicy;
MD.name = std::string("MinimumNorm2"); MD.name = std::string("MinimumNorm2");
MD.MDsteps = 5; //5 steps of 0.2 for GP* ensembles MD.MDsteps = user_params.Steps;
MD.trajL = 1.0; MD.trajL = 1.0;
HMCparameters HMCparams; HMCparameters HMCparams;
@ -358,6 +363,14 @@ int main(int argc, char **argv) {
GaugeImplPolicy::setDirections(dirs4); //gauge BC GaugeImplPolicy::setDirections(dirs4); //gauge BC
//Run optional gauge field checksum checker and exit
if(file_load_check){
TheHMC.initializeGaugeFieldAndRNGs(Ud);
std::cout << GridLogMessage << " Done" << std::endl;
Grid_finalize();
return 0;
}
//////////////////////////////////// ////////////////////////////////////
// Collect actions // Collect actions