1
0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00

Test wilson flow, added maxTau for adaptive flow

This commit is contained in:
Guido Cossu 2017-06-02 17:02:29 +01:00
parent 9b44189d5a
commit 4a8c4ccfba

View File

@ -33,7 +33,8 @@ namespace Grid{
GRID_SERIALIZABLE_CLASS_MEMBERS(WFParameters, GRID_SERIALIZABLE_CLASS_MEMBERS(WFParameters,
int, steps, int, steps,
double, step_size, double, step_size,
int, meas_interval); int, meas_interval,
double, maxTau); // for the adaptive algorithm
template <class ReaderClass > template <class ReaderClass >
@ -97,11 +98,11 @@ int main(int argc, char **argv) {
WilsonFlow<PeriodicGimplR> WF(WFPar.steps, WFPar.step_size, WFPar.meas_interval); WilsonFlow<PeriodicGimplR> WF(WFPar.steps, WFPar.step_size, WFPar.meas_interval);
WF.smear(Uflow, Umu); WF.smear_adaptive(Uflow, Umu, WFPar.maxTau);
RealD WFlow_plaq = WilsonLoops<PeriodicGimplR>::avgPlaquette(Uflow); RealD WFlow_plaq = WilsonLoops<PeriodicGimplR>::avgPlaquette(Uflow);
RealD WFlow_TC = WilsonLoops<PeriodicGimplR>::TopologicalCharge(Uflow); RealD WFlow_TC = WilsonLoops<PeriodicGimplR>::TopologicalCharge(Uflow);
RealD WFlow_T0 = WF.energyDensityPlaquette(WFPar.steps, Uflow); RealD WFlow_T0 = WF.energyDensityPlaquette(Uflow);
std::cout << GridLogMessage << "Plaquette "<< conf << " " << WFlow_plaq << std::endl; std::cout << GridLogMessage << "Plaquette "<< conf << " " << WFlow_plaq << std::endl;
std::cout << GridLogMessage << "T0 "<< conf << " " << WFlow_T0 << std::endl; std::cout << GridLogMessage << "T0 "<< conf << " " << WFlow_T0 << std::endl;
std::cout << GridLogMessage << "TopologicalCharge "<< conf << " " << WFlow_TC << std::endl; std::cout << GridLogMessage << "TopologicalCharge "<< conf << " " << WFlow_TC << std::endl;
@ -135,7 +136,8 @@ JSON
"WilsonFlow":{ "WilsonFlow":{
"steps": 200, "steps": 200,
"step_size": 0.01, "step_size": 0.01,
"meas_interval": 50 "meas_interval": 50,
"maxTau": 2.0
}, },
"Configurations":{ "Configurations":{
"conf_prefix": "ckpoint_lat", "conf_prefix": "ckpoint_lat",