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

scalar free propagator fix

Browse Source
This commit is contained in:
Antonin Portelli 2017-01-05 14:58:07 +00:00
parent 673994b281
commit 82b3f54697
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
lib/qcd/action/scalar/Scalar.h
View File

@ -148,10 +148,11 @@ namespace QCD{
void Scalar<SImpl>::MomentumSpacePropagator(ScalarField &out, RealD m) void Scalar<SImpl>::MomentumSpacePropagator(ScalarField &out, RealD m)
{ {
GridBase *grid = out._grid; GridBase *grid = out._grid;
ScalarField kmu(grid); ScalarField kmu(grid), one(grid);
const unsigned int nd = grid->_ndimension; const unsigned int nd = grid->_ndimension;
std::vector<int> &l = grid->_fdimensions; std::vector<int> &l = grid->_fdimensions;
one = Complex(1.0,0.0);
out = m*m; out = m*m;
for(int mu = 0; mu < nd; mu++) for(int mu = 0; mu < nd; mu++)
{ {
@ -161,6 +162,7 @@ namespace QCD{
kmu = 2.*sin(.5*twoPiL*kmu); kmu = 2.*sin(.5*twoPiL*kmu);
out = out + kmu*kmu; out = out + kmu*kmu;
} }
out = one/out;
} }
template <typename SImpl> template <typename SImpl>