From 009ccd581ede8faf0ba748fa49a1757419106e23 Mon Sep 17 00:00:00 2001
From: Felix Erben <felix.erben@ed.ac.uk>
Date: Mon, 26 Apr 2021 10:36:33 +0100
Subject: [PATCH] bugfix 3D stout smearing

---
 Grid/qcd/smearing/StoutSmearing.h | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/Grid/qcd/smearing/StoutSmearing.h b/Grid/qcd/smearing/StoutSmearing.h
index ed2ccdb6..629f81e2 100644
--- a/Grid/qcd/smearing/StoutSmearing.h
+++ b/Grid/qcd/smearing/StoutSmearing.h
@@ -89,11 +89,12 @@ public:
     SmearBase->smear(C, U);
 
     for (int mu = 0; mu < Nd; mu++) {
-      if( mu == OrthogDim )
+      Umu = peekLorentz(U, mu);
+      if( mu == OrthogDim ){
         tmp = 1.0;  // Don't smear in the orthogonal direction
+      }
       else {
         tmp = peekLorentz(C, mu);
-        Umu = peekLorentz(U, mu);
         iq_mu = Ta(
                    tmp *
                    adj(Umu));  // iq_mu = Ta(Omega_mu) to match the signs with the paper