portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

Hadrons: EMLepton: Wall source

Browse Source
This commit is contained in:
Vera Guelpers 2019-05-30 10:07:08 +01:00
parent e35e8da111
commit 166feb6483
1 changed files with 4 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -43,7 +43,7 @@ BEGIN_HADRONS_NAMESPACE
* *
* L(x) = \sum_y S(x,y) i*\gamma_mu*A_mu S(y,xl) \delta_{(tl-x0),dt} * L(x) = \sum_y S(x,y) i*\gamma_mu*A_mu S(y,xl) \delta_{(tl-x0),dt}
* *
* with xl = (0,0,0,tl) * with a wall source for the lepton at tl
* *
* In addition outputs the propagator without photon vertex * In addition outputs the propagator without photon vertex
* *
@ -195,12 +195,9 @@ void TEMLepton<FImpl>::execute(void)
unsigned int tl=0; unsigned int tl=0;
//point source at (0,0,0,tl) //wallsource at tl
position.clear(); sourcetmp = 1.;
for(int tt=0;tt<Nd-1;tt++) position.push_back(0); sourcetmp = where((tlat == tl), sourcetmp, 0.*sourcetmp);
position.push_back(tl);
sourcetmp = zero;
pokeSite(id, sourcetmp, position);
//free propagator from pt source //free propagator from pt source
for (unsigned int s = 0; s < Ns; ++s) for (unsigned int s = 0; s < Ns; ++s)