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

started stout smearing for small w

Browse Source
This commit is contained in:
Felix Erben 2019-04-26 14:27:40 +01:00
parent 3ac5a69a57
commit e70e03f560
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -139,6 +139,7 @@ public:
iQ2 = iQ * iQ; iQ2 = iQ * iQ;
iQ3 = iQ * iQ2; iQ3 = iQ * iQ2;
//We should check sgn(c0) here already and then apply eq (34) from 0311018
set_uw(u, w, iQ2, iQ3); set_uw(u, w, iQ2, iQ3);
set_fj(f0, f1, f2, u, w); set_fj(f0, f1, f2, u, w);
@ -198,9 +199,8 @@ public:
} }
LatticeComplex func_xi0(const LatticeComplex& w) const { LatticeComplex func_xi0(const LatticeComplex& w) const {
// Define a function to do the check // Definition from arxiv 0311018
// if( w < 1e-4 ) std::cout << GridLogWarning<< "[Smear_stout] w too small: //if (abs(w) < 0.05) {w2 = w*w; return 1.0 - w2/6.0 * (1.0-w2/20.0 * (1.0-w2/42.0));}
// "<< w <<"\n";
return sin(w) / w; return sin(w) / w;
} }