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

QedFVol: set infinite-volume photon propagator to 1 at x=0,

Browse Source 
so that momentum-spage photon propagator is non-negative.
Need to check whether this is sufficient for all volumes.
This commit is contained in:
James Harrison 2018-02-14 16:30:09 +00:00
parent 724cf02d4a
commit 4680a977c3
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
lib/qcd/action/gauge/Photon.h
View File

@ -109,7 +109,7 @@ namespace QCD{
std::vector<int> &l = grid->_fdimensions; std::vector<int> &l = grid->_fdimensions;
std::vector<int> x0(nd,0); std::vector<int> x0(nd,0);
TComplex Tone = Complex(1.0,0.0); TComplex Tone = Complex(1.0,0.0);
TComplex Tzero = Complex(0.0,0.0); // TComplex Tzero = Complex(0.0,0.0);
FFT fft(grid); FFT fft(grid);
one = Complex(1.0,0.0); one = Complex(1.0,0.0);
@ -122,7 +122,7 @@ namespace QCD{
} }
pokeSite(Tone, out, x0); pokeSite(Tone, out, x0);
out = one/out; out = one/out;
pokeSite(Tzero, out, x0); // pokeSite(Tzero, out, x0);
fft.FFT_all_dim(out, out, FFT::forward); fft.FFT_all_dim(out, out, FFT::forward);
} }