portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-05 03:35:55 +01:00
Code Releases Activity

Two heads are better than one. Combined effort and hopefully spatial smearing now fixed!

Browse Source
This commit is contained in:
Michael Marshall 2019-04-26 12:18:11 +01:00
parent b1768ba820
commit 4f3d1ea6e8
2 changed files with 7 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
Grid/qcd/smearing/StoutSmearing.h
View File

@ -53,7 +53,7 @@ class Smear_Stout : public Smear<Gimpl> {
std::vector<double> rho3d(Nd*Nd);
for (int mu=0; mu<Nd; mu++)
for (int nu=0; nu<Nd; nu++)
rho3d[mu + Nd * nu] = (mu == nu || mu == orthogdim || nu == orthogdim) ? 0.0 : rho;
rho3d[mu + Nd * nu] = (mu == nu || mu == orthogdim || nu == orthogdim) ? 1.0 : rho;
return rho3d;
};
@ -62,12 +62,14 @@ public:
/*! Stout smearing with base explicitly specified */
Smear_Stout(Smear<Gimpl>* base) : SmearBase{base} {
std::cout << GridLogDebug << "Stout smearing constructor : Smear_Stout(Smear<Gimpl>* base)" << std::endl
assert(Nc == 3 && "Stout smearing currently implemented only for Nc==3");
}
/*! Construct stout smearing object from explicitly specified rho matrix */
Smear_Stout(const std::vector<double>& rho_)
: OwnedBase{new Smear_APE<Gimpl>(rho_)}, SmearBase{OwnedBase.get()} {
std::cout << GridLogDebug << "Stout smearing constructor : Smear_Stout(const std::vector<double>& rho_)" << std::endl
assert(Nc == 3 && "Stout smearing currently implemented only for Nc==3");
}
@ -75,6 +77,7 @@ public:
Smear_Stout(double rho, int orthogdim = -1)
//: OwnedBase{(orthogdim<0 || orthogdim>=Nd) ? new Smear_APE<Gimpl>(rho) : new Smear_APE<Gimpl>(rho3D(rho,orthogdim))},
: SmearRho{ rho3D(rho,orthogdim) }, OwnedBase{ new Smear_APE<Gimpl>(SmearRho) }, SmearBase{OwnedBase.get()} {
std::cout << GridLogDebug << "Stout smearing constructor : Smear_StoutSmear_Stout(double rho, int orthogdim = -1)\nrho3d=" << rho3D << std::endl;
assert(Nc == 3 && "Stout smearing currently implemented only for Nc==3");
}
@ -89,7 +92,7 @@ public:
GaugeField C(U._grid);
GaugeLinkField tmp(U._grid), iq_mu(U._grid), Umu(U._grid);
std::cout << GridLogDebug << "Stout smearing started\n";
std::cout << GridLogDebug << "Stout smearing started" << std::endl;
// Smear the configurations
SmearBase->smear(C, U);
@ -103,7 +106,7 @@ public:
exponentiate_iQ(tmp, iq_mu);
pokeLorentz(u_smr, tmp * Umu, mu); // u_smr = exp(iQ_mu)*U_mu
}
std::cout << GridLogDebug << "Stout smearing completed\n";
std::cout << GridLogDebug << "Stout smearing completed" << std::endl;
};
void derivative(GaugeField& SigmaTerm, const GaugeField& iLambda,

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

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