From 92ec3404f8a404e7d6420ebfa0f113af5eb6ec6d Mon Sep 17 00:00:00 2001
From: James Harrison <jch1g10@soton.ac.uk>
Date: Mon, 14 Nov 2016 17:59:02 +0000
Subject: [PATCH] Set imaginary part of stochastic QED field to zero after FFT
 into position space

---
 lib/qcd/action/gauge/Photon.h | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/lib/qcd/action/gauge/Photon.h b/lib/qcd/action/gauge/Photon.h
index 852ecb3e..ca0a8d40 100644
--- a/lib/qcd/action/gauge/Photon.h
+++ b/lib/qcd/action/gauge/Photon.h
@@ -172,6 +172,8 @@ namespace QCD{
       pokeLorentz(aTilde, r, mu);
     }
     fft.FFT_all_dim(out, aTilde, FFT::backward);
+
+    out = 0.5*(out + conjugate(out));
   }
 //  template<class Gimpl>
 //  void Photon<Gimpl>::FeynmanGaugeMomentumSpacePropagator_L(GaugeField &out,