portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

Hadrons: Lepton Propagator for kl2, sign swap for antiperiodic boundary

Browse Source
This commit is contained in:
Vera Guelpers 2019-05-10 12:46:18 +01:00
parent 4f0631615f
commit 79e930ba12
1 changed files with 6 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
Hadrons/Modules/MFermion/EMLepton.hpp
View File

@ -233,6 +233,7 @@ void TEMLepton<FImpl>::execute(void)
for(tl=0;tl<nt;tl++){ for(tl=0;tl<nt;tl++){
//shift free propagator to different source positions //shift free propagator to different source positions
//account for possible anti-periodic boundary in time
proptmp = Cshift(freetmp,Tp, -tl); proptmp = Cshift(freetmp,Tp, -tl);
proptmp = where( tlat < tl, boundary[Tp]*proptmp, proptmp); 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 // keep the result for the desired delta t
Aslashlep = where(tlat == (tl-par().deltat+nt)%nt, proptmp, Aslashlep); 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 END_MODULE_NAMESPACE