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

bugfix

Browse Source
This commit is contained in:
Felix Erben 2021-01-18 18:57:05 +00:00
parent fa12b9a329
commit 9b73a937e7
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -951,7 +951,7 @@ void BaryonUtils<FImpl>::BaryonGamma3pt(
spinor result=Zero();
BaryonGamma3ptGroup1Site(Dq_ti,Dq_spec_p[0],Dq_spec_p[1],Dq_tf,GammaJ,GammaBi,GammaBf,wick_contraction,result);
//coalescedWrite(vcorr[ss],vcorr[ss]+result); //diff by factor 10???
coalescedWrite(vcorr[ss],vcorr[ss]+result);
coalescedWrite(vcorr[ss],coalescedRead(vcorr[ss])+result);
});//end loop over lattice sites
} else if (group == 2) {
@ -961,7 +961,7 @@ void BaryonUtils<FImpl>::BaryonGamma3pt(
typedef decltype(coalescedRead(vcorr[0])) spinor;
spinor result=Zero();
BaryonGamma3ptGroup2Site(Dq_spec_p[0],Dq_ti,Dq_spec_p[1],Dq_tf,GammaJ,GammaBi,GammaBf,wick_contraction,result);
coalescedWrite(vcorr[ss],vcorr[ss]+result);
coalescedWrite(vcorr[ss],coalescedRead(vcorr[ss])+result);
});//end loop over lattice sites
} else if (group == 3) {
accelerator_for(ss, grid->oSites(), grid->Nsimd(), {
@ -970,7 +970,7 @@ void BaryonUtils<FImpl>::BaryonGamma3pt(
typedef decltype(coalescedRead(vcorr[0])) spinor;
spinor result=Zero();
BaryonGamma3ptGroup3Site(Dq_spec_p[0],Dq_spec_p[1],Dq_ti,Dq_tf,GammaJ,GammaBi,GammaBf,wick_contraction,result);
coalescedWrite(vcorr[ss],vcorr[ss]+result);
coalescedWrite(vcorr[ss],coalescedRead(vcorr[ss])+result);
});//end loop over lattice sites
}