portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-09 23:45:36 +00:00
Code Releases Activity

Gparity 48ID tuning and exposure of trajectory length as input variable

Browse Source
This commit is contained in:
Christopher Kelly 2021-11-16 07:31:41 -08:00
parent f7e9621492
commit a1211cdcce
2 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Grid/qcd/smearing/WilsonFlow.h
View File

@ -166,8 +166,10 @@ std::vector<RealD> WilsonFlow<Gimpl>::flowMeasureEnergyDensityPlaquette(const Ga
std::vector<RealD> out;
GaugeField V(U);
for (unsigned int step = 0; step < Nstep; step++) { //bn tau = epsilon*(step+1) so tau after performing step=0 is epsilon
std::cout << GridLogMessage << "[WilsonFlow] Evolving step " << step << std::endl;
evolve_step(V);
if( step % measure_interval == 0){
std::cout << GridLogMessage << "[WilsonFlow] Computing energy density for step " << step << std::endl;
out.push_back( energyDensityPlaquette(step,V) );
}
}

8
HMC/Mobius2p1fIDSDRGparityEOFA_48ID.cc
View File

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