From cfdb56f3148a2730f63fd0a1927d9ec36a90d380 Mon Sep 17 00:00:00 2001 From: Peter Boyle Date: Tue, 13 May 2025 18:35:01 -0400 Subject: [PATCH] Run measurements at t=0 too --- Grid/qcd/smearing/WilsonFlow.h | 10 ++++++++++ 1 file changed, 10 insertions(+) diff --git a/Grid/qcd/smearing/WilsonFlow.h b/Grid/qcd/smearing/WilsonFlow.h index dc135823..188c6973 100644 --- a/Grid/qcd/smearing/WilsonFlow.h +++ b/Grid/qcd/smearing/WilsonFlow.h @@ -252,6 +252,11 @@ void WilsonFlow::smear(GaugeField& out, const GaugeField& in) const{ out = in; RealD taus = 0.; + + // Perform initial t=0 measurements + for(auto const &meas : this->functions) + meas.second(0,taus,out); + for (unsigned int step = 1; step <= Nstep; step++) { //step indicates the number of smearing steps applied at the time of measurement auto start = std::chrono::high_resolution_clock::now(); evolve_step(out, taus); @@ -336,6 +341,11 @@ void WilsonFlowAdaptive::smear(GaugeField& out, const GaugeField& in) con RealD taus = 0.; RealD eps = init_epsilon; unsigned int step = 0; + + // Perform initial t=0 measurements + for(auto const &meas : this->functions) + meas.second(step,taus,out); + do{ int step_success = evolve_step_adaptive(out, taus, eps); step += step_success; //step will not be incremented if the integration step fails