From fa12b9a3295f1fdf392b0367de4136e576b5401e Mon Sep 17 00:00:00 2001 From: Felix Erben Date: Wed, 13 Jan 2021 10:04:17 +0000 Subject: [PATCH] bugfix --- Grid/qcd/utils/BaryonUtils.h | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/Grid/qcd/utils/BaryonUtils.h b/Grid/qcd/utils/BaryonUtils.h index 828b4085..69bf8959 100644 --- a/Grid/qcd/utils/BaryonUtils.h +++ b/Grid/qcd/utils/BaryonUtils.h @@ -961,7 +961,7 @@ void BaryonUtils::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],result); + coalescedWrite(vcorr[ss],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::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],result); + coalescedWrite(vcorr[ss],vcorr[ss]+result); });//end loop over lattice sites }