portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-11-04 05:54:32 +00:00
Code Releases Activity

Hadrons: EMLepton: multiple source-sink separations at once

Browse Source
This commit is contained in:
Vera Guelpers
2019-07-03 14:36:34 +01:00
parent b120ef1fe4
commit eac6337466
1 changed files with 38 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -54,8 +54,9 @@ BEGIN_HADRONS_NAMESPACE
* - action: fermion action used for propagator (string) * - action: fermion action used for propagator (string)
* - emField: photon field A_mu (string) * - emField: photon field A_mu (string)
* - mass: input mass for the lepton propagator * - mass: input mass for the lepton propagator
* - boundary: boundary conditions for the lepton propagator, e.g. "1 1 1 -1"
* - twist: twisted boundary for lepton propagator, e.g. "0.0 0.0 0.0 0.5" * - twist: twisted boundary for lepton propagator, e.g. "0.0 0.0 0.0 0.5"
* - deltat: source-sink separation * - deltat: list of source-sink separations
* *
*******************************************************************************/ *******************************************************************************/
@@ -74,7 +75,7 @@ public:
double, mass, double, mass,
std::string , boundary, std::string , boundary,
std::string, twist, std::string, twist,
unsigned int, deltat); std::vector<unsigned int>, deltat);
}; };
template <typename FImpl> template <typename FImpl>
@@ -124,7 +125,12 @@ std::vector<std::string> TEMLepton<FImpl>::getInput(void)
template <typename FImpl> template <typename FImpl>
std::vector<std::string> TEMLepton<FImpl>::getOutput(void) std::vector<std::string> TEMLepton<FImpl>::getOutput(void)
{ {
std::vector<std::string> out = {getName(), getName() + "_free"}; std::vector<std::string> out = {};
for(int i=0; i<par().deltat.size(); i++)
{
out.push_back(std::to_string(par().deltat[i]) + "_" + getName() + "_free");
out.push_back(std::to_string(par().deltat[i]) + "_" + getName());
}
return out; return out;
} }
@@ -134,8 +140,11 @@ template <typename FImpl>
void TEMLepton<FImpl>::setup(void) void TEMLepton<FImpl>::setup(void)
{ {
Ls_ = env().getObjectLs(par().action); Ls_ = env().getObjectLs(par().action);
envCreateLat(PropagatorField, getName()); for(int i=0; i<par().deltat.size(); i++)
envCreateLat(PropagatorField, getName() + "_free"); {
envCreateLat(PropagatorField, std::to_string(par().deltat[i]) + "_" + getName() + "_free");
envCreateLat(PropagatorField, std::to_string(par().deltat[i]) + "_" + getName());
}
envTmpLat(FermionField, "source", Ls_); envTmpLat(FermionField, "source", Ls_);
envTmpLat(FermionField, "sol", Ls_); envTmpLat(FermionField, "sol", Ls_);
envTmpLat(FermionField, "tmp"); envTmpLat(FermionField, "tmp");
@@ -157,9 +166,6 @@ void TEMLepton<FImpl>::execute(void)
RealD mass = par().mass; RealD mass = par().mass;
Complex ci(0.0,1.0); Complex ci(0.0,1.0);
PropagatorField &Aslashlep = envGet(PropagatorField, getName());
PropagatorField &lep = envGet(PropagatorField, getName() + "_free");
envGetTmp(FermionField, source); envGetTmp(FermionField, source);
envGetTmp(FermionField, sol); envGetTmp(FermionField, sol);
envGetTmp(FermionField, tmp); envGetTmp(FermionField, tmp);
@@ -227,6 +233,22 @@ void TEMLepton<FImpl>::execute(void)
} }
} }
for(unsigned int dt=0;dt<par().deltat.size();dt++){
PropagatorField &lep = envGet(PropagatorField, std::to_string(par().deltat[dt]) + "_" + getName() + "_free");
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);
// free propagator for fixed source-sink separation
lep = where(tlat == (tl-par().deltat[dt]+nt)%nt, proptmp, lep);
}
//account for possible anti-periodic boundary in time
lep = where( tlat >= nt-par().deltat[dt], boundary[Tp]*lep, lep);
}
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
@@ -234,9 +256,6 @@ void TEMLepton<FImpl>::execute(void)
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);
// free propagator for fixed source-sink separation
lep = where(tlat == (tl-par().deltat+nt)%nt, proptmp, lep);
// i*A_mu*gamma_mu // i*A_mu*gamma_mu
sourcetmp = zero; sourcetmp = zero;
for(unsigned int mu=0;mu<=3;mu++) for(unsigned int mu=0;mu<=3;mu++)
@@ -276,13 +295,17 @@ 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); for(unsigned int dt=0;dt<par().deltat.size();dt++){
PropagatorField &Aslashlep = envGet(PropagatorField, std::to_string(par().deltat[dt]) + "_" + getName());
Aslashlep = where(tlat == (tl-par().deltat[dt]+nt)%nt, proptmp, Aslashlep);
}
} }
//account for possible anti-periodic boundary in time //account for possible anti-periodic boundary in time
Aslashlep = where( tlat >= nt-par().deltat, boundary[Tp]*Aslashlep, Aslashlep); for(unsigned int dt=0;dt<par().deltat.size();dt++){
lep = where( tlat >= nt-par().deltat, boundary[Tp]*lep, lep); PropagatorField &Aslashlep = envGet(PropagatorField, std::to_string(par().deltat[dt]) + "_" + getName());
Aslashlep = where( tlat >= nt-par().deltat[dt], boundary[Tp]*Aslashlep, Aslashlep);
}
} }
END_MODULE_NAMESPACE END_MODULE_NAMESPACE