portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-09-20 01:05:38 +01:00
Code Releases Activity

Merge pull request #217 from guelpers/feature/EMlepwall

Browse Source 
Hadrons: EMLepton: Wall source
This commit is contained in:
Antonin Portelli 2019-05-30 11:13:27 +02:00 committed by GitHub
commit b120ef1fe4
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
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)