From 739c2308b5ce9a9464dbbd9057dbe49f6b04cf59 Mon Sep 17 00:00:00 2001
From: James Harrison <jch1g10@soton.ac.uk>
Date: Tue, 15 Nov 2016 13:07:52 +0000
Subject: [PATCH] Set imaginary part of stochastic QED field to zero using
 real() instead of conjugate().

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

diff --git a/lib/qcd/action/gauge/Photon.h b/lib/qcd/action/gauge/Photon.h
index ca0a8d40..b6c1b76f 100644
--- a/lib/qcd/action/gauge/Photon.h
+++ b/lib/qcd/action/gauge/Photon.h
@@ -173,7 +173,7 @@ namespace QCD{
     }
     fft.FFT_all_dim(out, aTilde, FFT::backward);
 
-    out = 0.5*(out + conjugate(out));
+    out = real(out);
   }
 //  template<class Gimpl>
 //  void Photon<Gimpl>::FeynmanGaugeMomentumSpacePropagator_L(GaugeField &out,