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

Spatial smearing doesn't work yet. Fixed inconsistency in naming of perambulator in PerambLight.hpp

Browse Source
This commit is contained in:
Michael Marshall 2019-04-19 13:54:25 +01:00
parent 143b75956c
commit 7214681e11
2 changed files with 5 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Hadrons/Modules/MDistil/LapEvec.hpp
View File

@ -235,7 +235,7 @@ void TLapEvec<GImpl>::execute(void)
{ {
const StoutParameters &Stout{par().Stout}; const StoutParameters &Stout{par().Stout};
envGetTmp(GaugeField, Umu_stout); envGetTmp(GaugeField, Umu_stout);
Smear_Stout<PeriodicGimplR> LS(Stout.parm, Tdir); Smear_Stout<PeriodicGimplR> LS(Stout.parm );//, Tdir); // should be spatial - doesn't work yet
for (int i = 0; i < Stout.steps; i++) { for (int i = 0; i < Stout.steps; i++) {
LS.smear(Umu_stout, Umu_smear); LS.smear(Umu_stout, Umu_smear);
Umu_smear = Umu_stout; Umu_smear = Umu_stout;

9
Hadrons/Modules/MDistil/PerambLight.hpp
View File

@ -176,11 +176,10 @@ void TPerambLight<FImpl>::execute(void)
const std::string &UniqueIdentifier{par().UniqueIdentifier}; const std::string &UniqueIdentifier{par().UniqueIdentifier};
auto &noise = envGet(std::vector<Complex>, par().noise); auto &noise = envGet(std::vector<Complex>, par().noise);
auto &perambulator = envGet(Perambulator<SpinVector COMMA 6 COMMA sizeof(Real)>, auto &perambulator = envGet(Perambulator<SpinVector COMMA 6 COMMA sizeof(Real)>, getName());
getName() + "_perambulator_light"); auto &epack = envGet(Grid::Hadrons::EigenPack<LatticeColourVector>, par().eigenPack);
auto &epack = envGet(Grid::Hadrons::EigenPack<LatticeColourVector>, par().eigenPack); auto &unsmeared_sink = envGet(std::vector<FermionField>, getName() + "_unsmeared_sink");
auto &unsmeared_sink = envGet(std::vector<FermionField>, getName() + "_unsmeared_sink");
// Load perambulator if it exists on disk instead of creating it // Load perambulator if it exists on disk instead of creating it