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

Merge pull request #65 from jch1g10/feature/qed-fvol

Browse Source
This commit is contained in:
Antonin Portelli 2016-11-01 19:53:25 +00:00 committed by GitHub
commit 330a9b3f4c
2 changed files with 8 additions and 8 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;
std::vector<int> &l = grid->_fdimensions;
std::vector<int> zm(nd,0);
TComplexD Tone = ComplexD(1.0,0.0);
TComplexD Tzero= ComplexD(0.0,0.0);
TComplex Tone = Complex(1.0,0.0);
TComplex Tzero= Complex(0.0,0.0);
one = ComplexD(1.0,0.0);
one = Complex(1.0,0.0);
out = zero;
for(int mu = 0; mu < nd; mu++)
{
RealD twoPiL = M_PI*2./l[mu];
Real twoPiL = M_PI*2./l[mu];
LatticeCoordinate(kmu,mu);
kmu = 2.*sin(.5*twoPiL*kmu);
@ -93,7 +93,7 @@ namespace QCD{
}
pokeSite(Tone, out, zm);
out = one/out;
pokeSite(Tzero, out,zm);
pokeSite(Tzero, out, zm);
}
template<class Gimpl>
@ -107,7 +107,7 @@ namespace QCD{
case ZmScheme::QedTL:
{
std::vector<int> zm(nd,0);
TComplexD Tzero = ComplexD(0.0,0.0);
TComplex Tzero = Complex(0.0,0.0);
pokeSite(Tzero, out, zm);

4
programs/qed-fvol/qed-fvol.cc
View File

@ -90,8 +90,8 @@ int main(int argc, char *argv[])
}
}
double vol = grid.gSites();
double faces = (1.0 * Nd * (Nd - 1)) / 2.0;
Real vol = grid.gSites();
Real faces = (1.0 * Nd * (Nd - 1)) / 2.0;
Complex avgPlaqA = sum(trace(plaqA));
avgPlaqA = avgPlaqA / vol / faces;