portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-09-20 09:15:38 +01:00
Code Releases Activity

Merge pull request #64 from jch1g10/feature/feynman-rules

Browse Source
This commit is contained in:
Antonin Portelli 2016-11-01 19:53:08 +00:00 committed by GitHub
commit c2d78493c8
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -78,14 +78,14 @@ namespace QCD{
const unsigned int nd = grid->_ndimension; const unsigned int nd = grid->_ndimension;
std::vector<int> &l = grid->_fdimensions; std::vector<int> &l = grid->_fdimensions;
std::vector<int> zm(nd,0); std::vector<int> zm(nd,0);
TComplexD Tone = ComplexD(1.0,0.0); TComplex Tone = Complex(1.0,0.0);
TComplexD Tzero= ComplexD(0.0,0.0); TComplex Tzero= Complex(0.0,0.0);
one = ComplexD(1.0,0.0); one = Complex(1.0,0.0);
out = zero; out = zero;
for(int mu = 0; mu < nd; mu++) for(int mu = 0; mu < nd; mu++)
{ {
RealD twoPiL = M_PI*2./l[mu]; Real twoPiL = M_PI*2./l[mu];
LatticeCoordinate(kmu,mu); LatticeCoordinate(kmu,mu);
kmu = 2.*sin(.5*twoPiL*kmu); kmu = 2.*sin(.5*twoPiL*kmu);
@ -93,7 +93,7 @@ namespace QCD{
} }
pokeSite(Tone, out, zm); pokeSite(Tone, out, zm);
out = one/out; out = one/out;
pokeSite(Tzero, out,zm); pokeSite(Tzero, out, zm);
} }
template<class Gimpl> template<class Gimpl>
@ -107,7 +107,7 @@ namespace QCD{
case ZmScheme::QedTL: case ZmScheme::QedTL:
{ {
std::vector<int> zm(nd,0); std::vector<int> zm(nd,0);
TComplexD Tzero = ComplexD(0.0,0.0); TComplex Tzero = Complex(0.0,0.0);
pokeSite(Tzero, out, zm); pokeSite(Tzero, out, zm);