From 4bfe678218236256477afc69b5bfb687aa9b3b03 Mon Sep 17 00:00:00 2001
From: Nils Asmussen <n.asmussen@soton.ac.uk>
Date: Sun, 2 Jun 2019 12:36:57 +0100
Subject: [PATCH] MSource::Gauss Integer is unsigned...

---
 Hadrons/Modules/MSource/Gauss.hpp | 17 ++++++-----------
 1 file changed, 6 insertions(+), 11 deletions(-)

diff --git a/Hadrons/Modules/MSource/Gauss.hpp b/Hadrons/Modules/MSource/Gauss.hpp
index 2519569e..347b01f5 100644
--- a/Hadrons/Modules/MSource/Gauss.hpp
+++ b/Hadrons/Modules/MSource/Gauss.hpp
@@ -24,8 +24,8 @@ public:
     GRID_SERIALIZABLE_CLASS_MEMBERS(GaussPar,
                                     std::string, position,
                                     std::string, mom,
-                                    int,         tA,
-                                    int,         tB,
+                                    Integer,     tA,
+                                    Integer,     tB,
                                     double,      width);
 };
 
@@ -115,8 +115,6 @@ void TGauss<FImpl>::execute(void)
     const Real fact=-0.5/std::pow(par().width,2);
     const Complex i(0.0, 1.0);
     const SitePropagator idMat=[](){ SitePropagator s; s=1.; return s; }();
-    const vTInteger vTtA=par().tA;
-    const vTInteger vTtB=par().tB;
 
     ScalarRho=zero;
     for(int mu=0; mu<dim; mu++) {
@@ -126,8 +124,6 @@ void TGauss<FImpl>::execute(void)
         const int Lmu=env().getDim(mu);
         const int LmuHalf=Lmu/2;
         const int posMu=position_[mu];
-        const vTInteger vTLmuHalf=LmuHalf;
-        const vTInteger vTposMu=posMu;
 
         LatticeCoordinate(component, mu);
         LatticeCoordinate(compHelper, mu);
@@ -137,16 +133,15 @@ void TGauss<FImpl>::execute(void)
 
         //Gauss distribution
         component-=Complex(posMu);
-        compHelper-=vTposMu;
-        if(posMu<Lmu/2)
+        if(posMu<LmuHalf)
         {
-            component=where((compHelper>vTLmuHalf),
+            component=where((compHelper>Integer(posMu+LmuHalf)),
                     component-Complex(Lmu),
                     component);
         }
         else
         {
-            component=where((compHelper<=-vTLmuHalf),
+            component=where((compHelper<=Integer(posMu-LmuHalf)),
                     component+Complex(Lmu),
                     component);
         }
@@ -162,7 +157,7 @@ void TGauss<FImpl>::execute(void)
 
     //select time slices
     LatticeCoordinate(compHelper, dim);
-    ScalarRho=where((compHelper>=vTtA && compHelper<=vTtB),
+    ScalarRho=where((compHelper>=par().tA && compHelper<=par().tB),
           ScalarRho,
           0.*ScalarRho);