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

QedFVol: remove unnecessary copies of free propagator from shifted sources in ScalarVP

Browse Source
This commit is contained in:
James Harrison 2017-04-27 14:17:50 +01:00
parent bd466a55a8
commit 3ac27e5596
2 changed files with 8 additions and 55 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

57
extras/Hadrons/Modules/MScalar/ScalarVP.cc
View File

@ -48,16 +48,12 @@ void TScalarVP::setup(void)
propTadName_ = getName() + "_propTad"; propTadName_ = getName() + "_propTad";
phaseName_.clear(); phaseName_.clear();
muGFSrcName_.clear();
muProp0Name_.clear();
muPropQName_.clear(); muPropQName_.clear();
muPropSunName_.clear(); muPropSunName_.clear();
muPropTadName_.clear(); muPropTadName_.clear();
for (unsigned int mu = 0; mu < env().getNd(); ++mu) for (unsigned int mu = 0; mu < env().getNd(); ++mu)
{ {
phaseName_.push_back("_shiftphase_" + std::to_string(mu)); phaseName_.push_back("_shiftphase_" + std::to_string(mu));
muGFSrcName_.push_back("_" + getName() + "_DinvSrc_" + std::to_string(mu));
muProp0Name_.push_back(getName() + "_prop0_" + std::to_string(mu));
muPropQName_.push_back(getName() + "_propQ_" + std::to_string(mu)); muPropQName_.push_back(getName() + "_propQ_" + std::to_string(mu));
muPropSunName_.push_back(getName() + "_propSun_" + std::to_string(mu)); muPropSunName_.push_back(getName() + "_propSun_" + std::to_string(mu));
muPropTadName_.push_back(getName() + "_propTad_" + std::to_string(mu)); muPropTadName_.push_back(getName() + "_propTad_" + std::to_string(mu));
@ -78,24 +74,10 @@ void TScalarVP::setup(void)
{ {
env().registerLattice<ScalarField>(GFSrcName_); env().registerLattice<ScalarField>(GFSrcName_);
} }
if (!env().hasRegisteredObject(muGFSrcName_[0]))
{
for (unsigned int mu = 0; mu < env().getNd(); ++mu)
{
env().registerLattice<ScalarField>(muGFSrcName_[mu]);
}
}
if (!env().hasRegisteredObject(prop0Name_)) if (!env().hasRegisteredObject(prop0Name_))
{ {
env().registerLattice<ScalarField>(prop0Name_); env().registerLattice<ScalarField>(prop0Name_);
} }
if (!env().hasRegisteredObject(muProp0Name_[0]))
{
for (unsigned int mu = 0; mu < env().getNd(); ++mu)
{
env().registerLattice<ScalarField>(muProp0Name_[mu]);
}
}
env().registerLattice<ScalarField>(propQName_); env().registerLattice<ScalarField>(propQName_);
for (unsigned int mu = 0; mu < env().getNd(); ++mu) for (unsigned int mu = 0; mu < env().getNd(); ++mu)
{ {
@ -168,23 +150,6 @@ void TScalarVP::execute(void)
{ {
GFSrc_ = env().getObject<ScalarField>(GFSrcName_); GFSrc_ = env().getObject<ScalarField>(GFSrcName_);
} }
// cache G*exp(i*k_mu)*F*src
if (!env().hasCreatedObject(muGFSrcName_[0]))
{
for (unsigned int mu = 0; mu < env().getNd(); ++mu)
{
muGFSrc_.push_back(env().createLattice<ScalarField>(muGFSrcName_[mu]));
fft.FFT_all_dim(*(muGFSrc_[mu]), source, FFT::forward);
*(muGFSrc_[mu]) = (*freeMomProp_)*(*phase_[mu])*(*muGFSrc_[mu]);
}
}
else
{
for (unsigned int mu = 0; mu < env().getNd(); ++mu)
{
muGFSrc_.push_back(env().getObject<ScalarField>(muGFSrcName_[mu]));
}
}
// cache position-space free scalar propagators // cache position-space free scalar propagators
if (!env().hasCreatedObject(prop0Name_)) if (!env().hasCreatedObject(prop0Name_))
{ {
@ -195,21 +160,6 @@ void TScalarVP::execute(void)
{ {
prop0_ = env().getObject<ScalarField>(prop0Name_); prop0_ = env().getObject<ScalarField>(prop0Name_);
} }
if (!env().hasCreatedObject(muProp0Name_[0]))
{
for (unsigned int mu = 0; mu < env().getNd(); ++mu)
{
muProp0_.push_back(env().createLattice<ScalarField>(muProp0Name_[mu]));
fft.FFT_all_dim(*(muProp0_[mu]), *(muGFSrc_[mu]), FFT::backward);
}
}
else
{
for (unsigned int mu = 0; mu < env().getNd(); ++mu)
{
muProp0_.push_back(env().getObject<ScalarField>(muProp0Name_[mu]));
}
}
// PROPAGATOR CALCULATION // PROPAGATOR CALCULATION
// Propagator from unshifted source // Propagator from unshifted source
@ -220,15 +170,18 @@ void TScalarVP::execute(void)
// Propagators from shifted sources // Propagators from shifted sources
std::vector<ScalarField *> muPropQ_, muPropSun_, muPropTad_; std::vector<ScalarField *> muPropQ_, muPropSun_, muPropTad_;
ScalarField buf(env().getGrid());
for (unsigned int mu = 0; mu < env().getNd(); ++mu) for (unsigned int mu = 0; mu < env().getNd(); ++mu)
{ {
muPropQ_.push_back(env().createLattice<ScalarField>(muPropQName_[mu])); muPropQ_.push_back(env().createLattice<ScalarField>(muPropQName_[mu]));
muPropSun_.push_back(env().createLattice<ScalarField>(muPropSunName_[mu])); muPropSun_.push_back(env().createLattice<ScalarField>(muPropSunName_[mu]));
muPropTad_.push_back(env().createLattice<ScalarField>(muPropTadName_[mu])); muPropTad_.push_back(env().createLattice<ScalarField>(muPropTadName_[mu]));
buf = adj(*phase_[mu])*(*GFSrc_);
chargedProp(*(muPropQ_[mu]), *(muPropSun_[mu]), *(muPropTad_[mu]), chargedProp(*(muPropQ_[mu]), *(muPropSun_[mu]), *(muPropTad_[mu]),
*(muGFSrc_[mu]), fft); buf, fft);
} }
} }
// Calculate O(q) and O(q^2) terms of momentum-space charged propagator // Calculate O(q) and O(q^2) terms of momentum-space charged propagator

6
extras/Hadrons/Modules/MScalar/ScalarVP.hpp
View File

@ -50,10 +50,10 @@ private:
private: private:
std::string freeMomPropName_, GFSrcName_, prop0Name_, std::string freeMomPropName_, GFSrcName_, prop0Name_,
propQName_, propSunName_, propTadName_; propQName_, propSunName_, propTadName_;
std::vector<std::string> phaseName_, muGFSrcName_, muProp0Name_, std::vector<std::string> phaseName_, muPropQName_, muPropSunName_,
muPropQName_, muPropSunName_, muPropTadName_; muPropTadName_;
ScalarField *freeMomProp_, *GFSrc_, *prop0_; ScalarField *freeMomProp_, *GFSrc_, *prop0_;
std::vector<ScalarField *> phase_, muGFSrc_, muProp0_; std::vector<ScalarField *> phase_;
EmField *A; EmField *A;
}; };