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

Added momentum scaling to scalar HMC theories in order to follow UKQCD/CPS conventions

Browse Source
This commit is contained in:
Henrique B. R. 2020-04-02 17:38:47 +01:00
parent 05ebc458e2
commit 2c22db841a
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/action/scalar/ScalarImpl.h
View File

@ -1,5 +1,13 @@
#pragma once
#define CPS_MD_TIME
#ifdef CPS_MD_TIME
#define HMC_MOMENTUM_DENOMINATOR (2.0)
#else
#define HMC_MOMENTUM_DENOMINATOR (1.0)
#endif
NAMESPACE_BEGIN(Grid);
template <class S>
@ -20,13 +28,17 @@ public:
typedef Field PropagatorField;
static inline void generate_momenta(Field& P, GridParallelRNG& pRNG){
RealD scale = ::sqrt(HMC_MOMENTUM_DENOMINATOR);
// CPS and UKQCD conventions not yet implemented for U(1) scalars.
gaussian(pRNG, P);
P *= scale;
}
static inline Field projectForce(Field& P){return P;}
static inline void update_field(Field& P, Field& U, double ep) {
U += P*ep;
std::cout << "Field updated. Epsilon = " << std::setprecision(10) << ep << std::endl;
}
static inline RealD FieldSquareNorm(Field& U) {
@ -66,7 +78,7 @@ public:
}
static void FreePropagator(const Field &in, Field &out,
const Field &momKernel)
const Field &momKernel)
{
FFT fft((GridCartesian *)in.Grid());
Field inFT(in.Grid());
@ -139,14 +151,17 @@ public:
static inline void generate_momenta(Field &P, GridParallelRNG &pRNG)
{
RealD scale = ::sqrt(HMC_MOMENTUM_DENOMINATOR); // Being consistent with CPS and UKQCD conventions
#ifndef USE_FFT_ACCELERATION
Group::GaussianFundamentalLieAlgebraMatrix(pRNG, P);
#else
Field Pgaussian(P.Grid()), Pp(P.Grid());
ComplexField p2(P.Grid()); p2 = zero;
RealD M = FFT_MASS;
Group::GaussianFundamentalLieAlgebraMatrix(pRNG, Pgaussian);
FFT theFFT((GridCartesian*)P.Grid());
@ -156,8 +171,8 @@ public:
p2 = sqrt(p2);
Pp *= p2;
theFFT.FFT_all_dim(P, Pp, FFT::backward);
#endif //USE_FFT_ACCELERATION
P *= scale;
}
static inline Field projectForce(Field& P) {return P;}
@ -166,7 +181,8 @@ public:
{
#ifndef USE_FFT_ACCELERATION
double t0=usecond();
U += P*ep;
U += P*ep;
std::cout << "Field updated. Epsilon = " << std::setprecision(10) << ep << std::endl;
double t1=usecond();
double total_time = (t1-t0)/1e6;
std::cout << GridLogIntegrator << "Total time for updating field (s) : " << total_time << std::endl;