From e223d0b99f4059db0b5f7370141a9474065f8ca1 Mon Sep 17 00:00:00 2001
From: Michael Marshall <43034299+mmphys@users.noreply.github.com>
Date: Fri, 26 Apr 2019 16:00:35 +0100
Subject: [PATCH] Need to validate range about which exp^iQ is considered unity

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

diff --git a/Grid/qcd/smearing/StoutSmearing.h b/Grid/qcd/smearing/StoutSmearing.h
index fb86483b..590e7821 100644
--- a/Grid/qcd/smearing/StoutSmearing.h
+++ b/Grid/qcd/smearing/StoutSmearing.h
@@ -142,7 +142,8 @@ public:
     iQ2 = iQ * iQ;
     iQ3 = iQ * iQ2;
 
-    if(abs(real(trace(iQ2))) < 0.000000001) // Felix, please check this
+    Real tr = real(trace(iQ2));
+    if(tr > -0.0000001 && tr < 0.00000001 ) // Felix, please check this
       e_iQ = unity;
     else {
       //We should check sgn(c0) here already and then apply eq (34) from 0311018