From 82b3f546970fde47b1e1220679f3c6c772d5eff0 Mon Sep 17 00:00:00 2001
From: Antonin Portelli <antonin.portelli@me.com>
Date: Thu, 5 Jan 2017 14:58:07 +0000
Subject: [PATCH] scalar free propagator fix

---
 lib/qcd/action/scalar/Scalar.h | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/lib/qcd/action/scalar/Scalar.h b/lib/qcd/action/scalar/Scalar.h
index 194f6767..c053e15e 100644
--- a/lib/qcd/action/scalar/Scalar.h
+++ b/lib/qcd/action/scalar/Scalar.h
@@ -148,10 +148,11 @@ namespace QCD{
   void Scalar<SImpl>::MomentumSpacePropagator(ScalarField &out, RealD m)
   {
     GridBase           *grid = out._grid;
-    ScalarField        kmu(grid);
+    ScalarField        kmu(grid), one(grid);
     const unsigned int nd    = grid->_ndimension;
     std::vector<int>   &l    = grid->_fdimensions;
     
+    one = Complex(1.0,0.0);
     out = m*m;
     for(int mu = 0; mu < nd; mu++)
     {
@@ -161,6 +162,7 @@ namespace QCD{
       kmu = 2.*sin(.5*twoPiL*kmu);
       out = out + kmu*kmu;
     }
+    out = one/out;
   }
   
   template <typename SImpl>