From 79e930ba1203852f3405cc1aa32627b82a3dbb2a Mon Sep 17 00:00:00 2001
From: Vera Guelpers <Vera.Guelpers@ed.ac.uk>
Date: Fri, 10 May 2019 12:46:18 +0100
Subject: [PATCH] Hadrons: Lepton Propagator for kl2, sign swap for
 antiperiodic boundary

---
 Hadrons/Modules/MFermion/EMLepton.hpp | 6 ++++++
 1 file changed, 6 insertions(+)

diff --git a/Hadrons/Modules/MFermion/EMLepton.hpp b/Hadrons/Modules/MFermion/EMLepton.hpp
index a67f7495..1a0dfeef 100644
--- a/Hadrons/Modules/MFermion/EMLepton.hpp
+++ b/Hadrons/Modules/MFermion/EMLepton.hpp
@@ -233,6 +233,7 @@ void TEMLepton<FImpl>::execute(void)
     for(tl=0;tl<nt;tl++){
 
 	//shift free propagator to different source positions
+	//account for possible anti-periodic boundary in time
 	proptmp = Cshift(freetmp,Tp, -tl);
 	proptmp = where( tlat < tl, boundary[Tp]*proptmp, proptmp);
 
@@ -280,6 +281,11 @@ void TEMLepton<FImpl>::execute(void)
 	// keep the result for the desired delta t
 	Aslashlep = where(tlat == (tl-par().deltat+nt)%nt, proptmp, Aslashlep);
     }
+
+    //account for possible anti-periodic boundary in time
+    Aslashlep = where( tlat >= nt-par().deltat, boundary[Tp]*Aslashlep, Aslashlep);
+    lep = where( tlat >= nt-par().deltat, boundary[Tp]*lep, lep);
+
 }
 
 END_MODULE_NAMESPACE