portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

First attempt at minimising smearing

Browse Source
This commit is contained in:
Michael Marshall 2019-04-26 15:54:05 +01:00
parent 4333d97958
commit 2e220456d3
1 changed files with 9 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -142,11 +142,15 @@ public:
iQ2 = iQ * iQ;
iQ3 = iQ * iQ2;
//We should check sgn(c0) here already and then apply eq (34) from 0311018
set_uw(u, w, iQ2, iQ3);
set_fj(f0, f1, f2, u, w);
e_iQ = f0 * unity + timesMinusI(f1) * iQ - f2 * iQ2;
if(abs(real(trace(iQ2))) < 0.000000001) // Felix, please check this
e_iQ = unity;
else {
//We should check sgn(c0) here already and then apply eq (34) from 0311018
set_uw(u, w, iQ2, iQ3);
set_fj(f0, f1, f2, u, w);
e_iQ = f0 * unity + timesMinusI(f1) * iQ - f2 * iQ2;
}
};
void set_uw(LatticeComplex& u, LatticeComplex& w, GaugeLinkField& iQ2,