portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-07 04:35:56 +01:00
Code Releases Activity

continued work on baryons

Browse Source
This commit is contained in:
Felix Erben 2019-07-09 17:42:36 +01:00
parent cc3346073e
commit 98cf20cf06
1 changed files with 11 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
Hadrons/Modules/MContraction/Nucleon.hpp
View File

@ -127,11 +127,12 @@ void TNucleon<FImpl1, FImpl2, FImpl3>::execute(void)
auto &q2 = envGet(PropagatorField2, par().q2);
auto &q3 = envGet(PropagatorField3, par().q2);
envGetTmp(LatticeComplex, c);
//envGetTmp(LatticeComplex, quark2);
//envGetTmp(LatticeComplex, quark3);
envGetTmp(LatticeComplex, diquark);
Result result;
int nt = env().getDim(Tp);
result.corr.resize(nt);
std::vector<TComplex> buf;
// C = i gamma_2 gamma_4 => C gamma_5 = - i gamma_1 gamma_3
Gamma Cg5(Gamma::Algebra::SigmaXZ);
Gamma g4(Gamma::Algebra::GammaT); //needed for parity P_\pm = 0.5*(1 \pm \gamma_4)
@ -179,7 +180,13 @@ void TNucleon<FImpl1, FImpl2, FImpl3>::execute(void)
}
}
// saveResult(par().output, "meson", result);
/* sliceSum(c,buf,Tp);
for (unsigned int t = 0; t < buf.size(); ++t)
{
result.corr[t] = TensorRemove(buf[t]);
}*/
// saveResult(par().output, "baryon", result);
}
END_MODULE_NAMESPACE