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

Added method to compute and return the Wilson flow energy density over some number of steps

Browse Source
This commit is contained in:
Christopher Kelly 2021-06-30 17:24:00 -04:00
parent ac4f2d9798
commit 75a1f85162
1 changed files with 20 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -74,8 +74,15 @@ public:
}
void smear_adaptive(GaugeField&, const GaugeField&, RealD maxTau);
//Compute t^2 <E(t)> for timestep 'step' from the plaquette
RealD energyDensityPlaquette(unsigned int step, const GaugeField& U) const;
//Compute t^2 <E(t)> for time taus, set by smear_adaptive, from the plaquette
RealD energyDensityPlaquette(const GaugeField& U) const;
//Evolve the gauge field by Nstep steps of epsilon and return the energy density computed every interval steps
std::vector<RealD> flowMeasureEnergyDensityPlaquette(const GaugeField& U) const;
};
@ -154,11 +161,20 @@ RealD WilsonFlow<Gimpl>::energyDensityPlaquette(const GaugeField& U) const {
return 2.0 * taus * taus * SG.S(U)/U.Grid()->gSites();
}
template <class Gimpl>
std::vector<RealD> WilsonFlow<Gimpl>::flowMeasureEnergyDensityPlaquette(const GaugeField& U) const{
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
evolve_step(V);
if( step % measure_interval == 0){
out.push_back( energyDensityPlaquette(step,V) );
}
}
return out;
}
//#define WF_TIMING
template <class Gimpl>
void WilsonFlow<Gimpl>::smear(GaugeField& out, const GaugeField& in) const {
out = in;