From 4a4203c610143f0f77ab6c65a427b18de8593925 Mon Sep 17 00:00:00 2001
From: Felix Erben <dc-erbe1@tesseract-login1.ib0.sgi.cluster.dirac.ed.ac.uk>
Date: Thu, 18 Apr 2019 19:10:49 +0100
Subject: [PATCH] fixed stout smearing for now

---
 Grid/qcd/smearing/StoutSmearing.h          | 6 +++++-
 Hadrons/Modules/MGauge/StoutSmearing3D.hpp | 3 ++-
 2 files changed, 7 insertions(+), 2 deletions(-)

diff --git a/Grid/qcd/smearing/StoutSmearing.h b/Grid/qcd/smearing/StoutSmearing.h
index e4e11a63..65501c7e 100644
--- a/Grid/qcd/smearing/StoutSmearing.h
+++ b/Grid/qcd/smearing/StoutSmearing.h
@@ -71,9 +71,13 @@ public:
     }
 
   /*! Default constructor. rho is constant in all directions, optionally except for orthogonal dimension */
-  Smear_Stout(double rho = 1.0, int orthogdim = -1)
+  /*Smear_Stout(double rho = 1.0, int orthogdim = -1)
     : OwnedBase{new Smear_APE<Gimpl>(rho3D(rho,orthogdim))}, SmearBase{OwnedBase.get()} {
     assert(Nc == 3 && "Stout smearing currently implemented only for Nc==3");
+  }*/
+
+  Smear_Stout(double rho = 1.0) : SmearBase(new Smear_APE<Gimpl>(rho)) {
+    assert(Nc == 3);//                  "Stout smearing currently implemented only for Nc==3");
   }
 
   ~Smear_Stout() {}  // delete SmearBase...
diff --git a/Hadrons/Modules/MGauge/StoutSmearing3D.hpp b/Hadrons/Modules/MGauge/StoutSmearing3D.hpp
index ef9db9e9..c5724943 100644
--- a/Hadrons/Modules/MGauge/StoutSmearing3D.hpp
+++ b/Hadrons/Modules/MGauge/StoutSmearing3D.hpp
@@ -113,7 +113,8 @@ void TStoutSmearing3D<GImpl>::execute(void)
                  << " of 3D-stout smearing and rho= " << par().rho
 		 << "orthogonal to dimension " << par().orthogdim << std::endl;
 
-    Smear_Stout<GImpl> smearer(par().rho, par().orthogdim);
+    //Smear_Stout<GImpl> smearer(par().rho, par().orthogdim);
+    Smear_Stout<GImpl> smearer(par().rho);
     auto               &U    = envGet(GaugeField, par().gauge);
     auto               &Usmr = envGet(GaugeField, getName());