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

bugfix in Baryonutils

Browse Source
This commit is contained in:
ferben 2019-07-18 14:44:55 +01:00
parent cded7670d0
commit 11a8668d19
2 changed files with 7 additions and 21 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Grid/qcd/utils/BaryonUtils.h
View File

@ -66,7 +66,7 @@ void BaryonUtils<FImpl>::ContractBaryons(const PropagatorField &q1,
auto gD1a = GammaA * GammaA * D1; auto gD1a = GammaA * GammaA * D1;
auto gD1b = GammaA * g4 * GammaA * D1; auto gD1b = GammaA * g4 * GammaA * D1;
auto pD1 = 0.5* (gD1a + (double)parity * dD1b); auto pD1 = 0.5* (gD1a + (double)parity * gD1b);
auto gD3 = GammaB * D3; auto gD3 = GammaB * D3;
vobj result=zero; vobj result=zero;

26
Hadrons/Modules/MContraction/Baryon.hpp
View File

@ -33,6 +33,7 @@ See the full license in the file "LICENSE" in the top level distribution directo
#include <Hadrons/Global.hpp> #include <Hadrons/Global.hpp>
#include <Hadrons/Module.hpp> #include <Hadrons/Module.hpp>
#include <Hadrons/ModuleFactory.hpp> #include <Hadrons/ModuleFactory.hpp>
#include "/home/ferben/code/grid_devel/Grid/Grid/qcd/utils/BaryonUtils.h"
BEGIN_HADRONS_NAMESPACE BEGIN_HADRONS_NAMESPACE
@ -134,7 +135,7 @@ void TBaryon<FImpl1, FImpl2, FImpl3>::execute(void)
std::vector<TComplex> buf; std::vector<TComplex> buf;
// C = i gamma_2 gamma_4 => C gamma_5 = - i gamma_1 gamma_3 // C = i gamma_2 gamma_4 => C gamma_5 = - i gamma_1 gamma_3
Gamma GammaA(Gamma::Algebra::Identity); //Still hardcoded 1 /* Gamma GammaA(Gamma::Algebra::Identity); //Still hardcoded 1
Gamma GammaB(Gamma::Algebra::SigmaXZ); //Still hardcoded Cg5 Gamma GammaB(Gamma::Algebra::SigmaXZ); //Still hardcoded Cg5
Gamma g4(Gamma::Algebra::GammaT); //needed for parity P_\pm = 0.5*(1 \pm \gamma_4) Gamma g4(Gamma::Algebra::GammaT); //needed for parity P_\pm = 0.5*(1 \pm \gamma_4)
@ -170,25 +171,6 @@ void TBaryon<FImpl1, FImpl2, FImpl3>::execute(void)
auto Dbc = peekColour(q3,b_snk,c_src); //D_{beta' gamma} auto Dbc = peekColour(q3,b_snk,c_src); //D_{beta' gamma}
auto Dca = peekColour(q1,c_snk,a_src); //D_{gamma' alpha} auto Dca = peekColour(q1,c_snk,a_src); //D_{gamma' alpha}
auto Dcb = peekColour(q1,c_snk,b_src); //D_{gamma' beta} auto Dcb = peekColour(q1,c_snk,b_src); //D_{gamma' beta}
// This needs lees peekColours for some baryons, but does not compile - worth the effort?
/*if (wick_contraction[0] || wick_contraction[4])
auto Daa = peekColour(q2,a_snk,a_src); //D_{alpha' alpha}
if (wick_contraction[0] || wick_contraction[5])
auto Dbb = peekColour(q3,b_snk,b_src); //D_{beta' beta}
if (wick_contraction[0] || wick_contraction[3])
auto Dcc = peekColour(q1,c_snk,c_src); //D_{gamma' gamma}
if (wick_contraction[1] || wick_contraction[3])
auto Dab = peekColour(q2,a_snk,b_src); //D_{alpha' beta}
if (wick_contraction[2] || wick_contraction[5])
auto Dac = peekColour(q2,a_snk,c_src); //D_{alpha' gamma}
if (wick_contraction[2] || wick_contraction[3])
auto Dba = peekColour(q3,b_snk,a_src); //D_{beta' alpha}
if (wick_contraction[1] || wick_contraction[4])
auto Dbc = peekColour(q3,b_snk,c_src); //D_{beta' gamma}
if (wick_contraction[1] || wick_contraction[5])
auto Dca = peekColour(q1,c_snk,a_src); //D_{gamma' alpha}
if (wick_contraction[2] || wick_contraction[4])
auto Dcb = peekColour(q1,c_snk,b_src); //D_{gamma' beta}*/
// This is the \delta_{123}^{123} part // This is the \delta_{123}^{123} part
if (wick_contraction[0]){ if (wick_contraction[0]){
diquark = trace(GammaB * Daa * GammaB * Dbb); //1st GammaB and Daa transposed???? diquark = trace(GammaB * Daa * GammaB * Dbb); //1st GammaB and Daa transposed????
@ -229,8 +211,12 @@ void TBaryon<FImpl1, FImpl2, FImpl3>::execute(void)
} }
} }
} }
*/
ContractBaryons(q1,q2,q3,c);
sliceSum(c,buf,Tp); sliceSum(c,buf,Tp);
for (unsigned int t = 0; t < buf.size(); ++t) for (unsigned int t = 0; t < buf.size(); ++t)
{ {
result.corr[t] = TensorRemove(buf[t]); result.corr[t] = TensorRemove(buf[t]);