From e5a7ed43628516e675b1a470702dc625a4edfad4 Mon Sep 17 00:00:00 2001 From: Lanny91 Date: Wed, 8 Feb 2017 12:29:33 +0000 Subject: [PATCH] Moved write outside of loop, some physics corrections --- extras/Hadrons/Modules/MContraction/Meson.hpp | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/extras/Hadrons/Modules/MContraction/Meson.hpp b/extras/Hadrons/Modules/MContraction/Meson.hpp index 1ea46eb6..9581b8d5 100644 --- a/extras/Hadrons/Modules/MContraction/Meson.hpp +++ b/extras/Hadrons/Modules/MContraction/Meson.hpp @@ -181,7 +181,7 @@ void TMeson::execute(void) LatticeCoordinate(coor, 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); @@ -190,7 +190,7 @@ void TMeson::execute(void) { Gamma gSnk(gammaList[i].first); 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); result[i].gamma_snk = gammaList[i].first; @@ -200,8 +200,8 @@ void TMeson::execute(void) { result[i].corr[t] = TensorRemove(buf[t]); } - write(writer, "meson", result[i]); } + write(writer, "meson", result); } END_MODULE_NAMESPACE