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

Need to validate range about which exp^iQ is considered unity

Browse Source
This commit is contained in:
Michael Marshall 2019-04-26 16:00:35 +01:00
parent 2e220456d3
commit e223d0b99f
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -142,7 +142,8 @@ public:
iQ2 = iQ * iQ; iQ2 = iQ * iQ;
iQ3 = iQ * iQ2; iQ3 = iQ * iQ2;
if(abs(real(trace(iQ2))) < 0.000000001) // Felix, please check this Real tr = real(trace(iQ2));
if(tr > -0.0000001 && tr < 0.00000001 ) // Felix, please check this
e_iQ = unity; e_iQ = unity;
else { else {
//We should check sgn(c0) here already and then apply eq (34) from 0311018 //We should check sgn(c0) here already and then apply eq (34) from 0311018