From 3e712fe64324391e7e91cabd20c43ade9e04ccea Mon Sep 17 00:00:00 2001 From: Peter Boyle Date: Wed, 17 Apr 2019 11:54:17 +0100 Subject: [PATCH] Scale momentum convention to CPS/UKQCD MD time --- Grid/qcd/action/gauge/GaugeImplTypes.h | 12 ++++++++++-- 1 file changed, 10 insertions(+), 2 deletions(-) diff --git a/Grid/qcd/action/gauge/GaugeImplTypes.h b/Grid/qcd/action/gauge/GaugeImplTypes.h index 818e5405..95680900 100644 --- a/Grid/qcd/action/gauge/GaugeImplTypes.h +++ b/Grid/qcd/action/gauge/GaugeImplTypes.h @@ -29,6 +29,14 @@ directory #ifndef GRID_GAUGE_IMPL_TYPES_H #define GRID_GAUGE_IMPL_TYPES_H +#define CPS_MD_TIME + +#ifdef CPS_MD_TIME +#define HMC_MOMENTUM_DENOMINATOR (2.0) +#else +#define HMC_MOMENTUM_DENOMINATOR (1.0) +#endif + namespace Grid { namespace QCD { @@ -105,13 +113,13 @@ public: // // Expect cxmua variance sqrt(2). // - // Must scale the momentum by sqrt(2) up to invoke CPS and UKQCD conventions + // Must scale the momentum by sqrt(2) to invoke CPS and UKQCD conventions // LinkField Pmu(P._grid); Pmu = Zero(); for (int mu = 0; mu < Nd; mu++) { SU::GaussianFundamentalLieAlgebraMatrix(pRNG, Pmu); - RealD scale = ::sqrt(2) ; + RealD scale = ::sqrt(HMC_MOMENTUM_DENOMINATOR) ; Pmu = Pmu*scale; PokeIndex(P, Pmu, mu); }