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

Moved write outside of loop, some physics corrections

Browse Source
This commit is contained in:
Lanny91 2017-02-08 12:29:33 +00:00
parent b9f7ea47c3
commit e5a7ed4362
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
extras/Hadrons/Modules/MContraction/Meson.hpp
View File

@ -181,7 +181,7 @@ void TMeson<FImpl1, FImpl2>::execute(void)
LatticeCoordinate(coor, mu); LatticeCoordinate(coor, mu);
ph = ph + p[mu]*coor*((1./(env().getGrid()->_fdimensions[mu]))); ph = ph + p[mu]*coor*((1./(env().getGrid()->_fdimensions[mu])));
} }
ph = exp(-(Real)(2*M_PI)*i*ph); ph = exp((Real)(2*M_PI)*i*ph);
parseGammaString(gammaList); parseGammaString(gammaList);
@ -190,7 +190,7 @@ void TMeson<FImpl1, FImpl2>::execute(void)
{ {
Gamma gSnk(gammaList[i].first); Gamma gSnk(gammaList[i].first);
Gamma gSrc(gammaList[i].second); Gamma gSrc(gammaList[i].second);
c = trace((g5*gSnk)*q1*(gSrc*g5)*adj(q2))*ph; c = trace((g5*gSnk)*q1*(adj(gSrc)*g5)*adj(q2))*ph;
sliceSum(c, buf, Tp); sliceSum(c, buf, Tp);
result[i].gamma_snk = gammaList[i].first; result[i].gamma_snk = gammaList[i].first;
@ -200,8 +200,8 @@ void TMeson<FImpl1, FImpl2>::execute(void)
{ {
result[i].corr[t] = TensorRemove(buf[t]); result[i].corr[t] = TensorRemove(buf[t]);
} }
write(writer, "meson", result[i]);
} }
write(writer, "meson", result);
} }
END_MODULE_NAMESPACE END_MODULE_NAMESPACE