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

Browse Source
This commit is contained in:
ferben 2019-02-08 14:33:06 +00:00
parent cbd2dfe53f
commit c1341b8ed2
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
Hadrons/Modules/MDistil/BContraction.hpp
View File

@ -151,7 +151,7 @@ void TBContraction<FImpl>::execute(void)
Gamma::Algebra::SigmaYT, // C = i gamma_2 gamma_4 Gamma::Algebra::SigmaYT, // C = i gamma_2 gamma_4
Gamma::Algebra::GammaYGamma5, // i gamma_4 C gamma_5 = i gamma_2 gamma_5 Gamma::Algebra::GammaYGamma5, // i gamma_4 C gamma_5 = i gamma_2 gamma_5
}; };
//std::vector<Complex> factor23 = {(0.,-1.),(0.,1.),(0.,1.)}; std::vector<Complex> factor23 = {(0.,-1.),(0.,1.),(0.,1.)};
for (int i1=0 ; i1 < N_1 ; i1++){ for (int i1=0 ; i1 < N_1 ; i1++){
for (int i2=0 ; i2 < N_2 ; i2++){ for (int i2=0 ; i2 < N_2 ; i2++){
@ -160,6 +160,8 @@ void TBContraction<FImpl>::execute(void)
for (int t=0 ; t < Nt ; t++){ for (int t=0 ; t < Nt ; t++){
Bindex = i1 + N_1*(i2 + N_2*(i3 + N_3*(imom+Nmom*t))); Bindex = i1 + N_1*(i2 + N_2*(i3 + N_3*(imom+Nmom*t)));
ExtractSliceLocal(tmp1,one[i1],0,t,3); ExtractSliceLocal(tmp1,one[i1],0,t,3);
ExtractSliceLocal(tmp2,two[i2],0,t,3);
ExtractSliceLocal(tmp3,three[i3],0,t,3);
parallel_for (unsigned int sU = 0; sU < grid3d->oSites(); ++sU) parallel_for (unsigned int sU = 0; sU < grid3d->oSites(); ++sU)
{ {
for (int ie=0 ; ie < 6 ; ie++){ for (int ie=0 ; ie < 6 ; ie++){
@ -186,7 +188,6 @@ void TBContraction<FImpl>::execute(void)
g_diquark()(isl,isr)() = tmp111()(isl)(); g_diquark()(isl,isr)() = tmp111()(isl)();
} }
} }
// Really only the trace? Should check baryons again! laph paper lists c_{alpha,beta,gamma}, gattringer-lang two gamma matrices.
for (int is=0 ; is < 4 ; is++){ for (int is=0 ; is < 4 ; is++){
BField[Bindex]+=(double)epsilon_sgn[ie]*tmp11s()(is)()*g_diquark()(is,is)(); BField[Bindex]+=(double)epsilon_sgn[ie]*tmp11s()(is)()*g_diquark()(is,is)();
} }