From 9b73a937e7dbfca0ff0c49b5003c0416e572bb82 Mon Sep 17 00:00:00 2001
From: Felix Erben <felix.erben@ed.ac.uk>
Date: Mon, 18 Jan 2021 18:57:05 +0000
Subject: [PATCH] bugfix

---
 Grid/qcd/utils/BaryonUtils.h | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/Grid/qcd/utils/BaryonUtils.h b/Grid/qcd/utils/BaryonUtils.h
index 69bf8959..edc5c8d5 100644
--- a/Grid/qcd/utils/BaryonUtils.h
+++ b/Grid/qcd/utils/BaryonUtils.h
@@ -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
   }