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

Merge pull request #356 from felixerben/bugfix/stoutSmearing

Browse Source 
Jamie's fix
This commit is contained in:
Peter Boyle 2021-04-27 14:10:49 -04:00 committed by GitHub
commit c332d9f08b
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 8 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -85,22 +85,18 @@ public:
std::cout << GridLogDebug << "Stout smearing started\n"; std::cout << GridLogDebug << "Stout smearing started\n";
// Smear the configurations // C contains the staples multiplied by some rho
u_smr = U ; // set the smeared field to the current gauge field
SmearBase->smear(C, U); SmearBase->smear(C, U);
for (int mu = 0; mu < Nd; mu++) { for (int mu = 0; mu < Nd; mu++) {
if( mu == OrthogDim ) continue ;
// u_smr = exp(iQ_mu)*U_mu apart from Orthogdim
Umu = peekLorentz(U, mu); Umu = peekLorentz(U, mu);
if( mu == OrthogDim ){ tmp = peekLorentz(C, mu);
tmp = 1.0; // Don't smear in the orthogonal direction iq_mu = Ta( tmp * adj(Umu));
} exponentiate_iQ(tmp, iq_mu);
else { pokeLorentz(u_smr, tmp * Umu, mu);
tmp = peekLorentz(C, mu);
iq_mu = Ta(
tmp *
adj(Umu)); // iq_mu = Ta(Omega_mu) to match the signs with the paper
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\n";
}; };