From 785bc7a14f15dedfe26977a1af16b16f0c3d01ad Mon Sep 17 00:00:00 2001
From: Mashy Green <mashy@me.com>
Date: Mon, 10 Mar 2025 12:29:04 +0000
Subject: [PATCH] Adding staple zeroing fix

---
 Grid/qcd/utils/WilsonLoops.h | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

diff --git a/Grid/qcd/utils/WilsonLoops.h b/Grid/qcd/utils/WilsonLoops.h
index 655255d6..6e9318c0 100644
--- a/Grid/qcd/utils/WilsonLoops.h
+++ b/Grid/qcd/utils/WilsonLoops.h
@@ -302,7 +302,11 @@ public:
   }
 
   static void Staple(GaugeMat &staple, const std::vector<GaugeMat> &Umu, int mu) {
-    staple = Zero();
+
+    autoView(staple_v, staple, AcceleratorWrite);
+    accelerator_for(i, staple.Grid()->oSites(), Simd::Nsimd(), {
+        staple_v[i] = Zero();
+    });
 
     for (int nu = 0; nu < Nd; nu++) {