portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 05:30:46 +01:00
Code Releases Activity

MSource::Gauss fix missing factor

Browse Source
This commit is contained in:
Nils Asmussen 2019-05-24 20:16:09 +01:00
parent 6b27369ade
commit 4e52e46a2c
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
Hadrons/Modules/MSource/Gauss.hpp
View File

@ -10,9 +10,10 @@ BEGIN_HADRONS_NAMESPACE
/****************************************************************************** /******************************************************************************
* Gauss * * Gauss *
* result[n] = 1/(sqrt(2*pi)*width)^dim * * result[n] = 1/(sqrt(2*pi)*width)^dim *
* * exp(-|n-position|^2/(2*width^2) + i n.mom) * * * exp(-|n-position|^2/(2*width^2) + i p.n) *
* where: * * where: *
* n=(n[0],n[1],...,n[dim-1]) (lattice coordinate) * * n=(n[0],n[1],...,n[dim-1]) (lattice coordinate) *
* p[i]=2*pi/L[i]*mom[i] *
* dim=Nd-1 * * dim=Nd-1 *
******************************************************************************/ ******************************************************************************/
BEGIN_MODULE_NAMESPACE(MSource) BEGIN_MODULE_NAMESPACE(MSource)
@ -103,13 +104,13 @@ void TGauss<FImpl>::execute(void)
envGetTmp(LatticeComplex, component); envGetTmp(LatticeComplex, component);
const int dim=env().getNd()-1; const int dim=env().getNd()-1;
const double fact=-0.5/std::pow(par().width,2); const double fact=-0.5/std::pow(par().width,2);
Complex i(0.0, 1.0); const Complex i(0.0, 1.0);
//exp(fact*|n|^2 +i n.mom) //exp(fact*|n|^2 +i n.mom)
rho=zero; rho=zero;
for(int mu=0; mu<dim; mu++) { for(int mu=0; mu<dim; mu++) {
LatticeCoordinate(component, mu); LatticeCoordinate(component, mu);
rho+=(i*(static_cast<double>(mom_[mu])/env().getDim(mu)))*component; rho+=(i*(mom_[mu]*2*M_PI/env().getDim(mu)))*component;
//FIXME: the next three lines are very inefficient... //FIXME: the next three lines are very inefficient...
// should not need any communication (Cshift) here // should not need any communication (Cshift) here
assert(env().getDim(mu)%2==0); assert(env().getDim(mu)%2==0);