diff --git a/tests/forces/Test_double_ratio.cc b/tests/forces/Test_double_ratio.cc index 0a350692..dc6713f9 100644 --- a/tests/forces/Test_double_ratio.cc +++ b/tests/forces/Test_double_ratio.cc @@ -476,6 +476,20 @@ int main (int argc, char ** argv) // ForceTest(BdyNf2eo,U,DDHMCFilter); //////////////////// One flavour boundary det //////////////////// + RationalActionParams OFRp; // Up/down + OFRp.lo = 6.0e-5; + OFRp.hi = 90.0; + OFRp.inv_pow = 2; + OFRp.MaxIter = SP_iters; // get most shifts by 2000, stop sharing space + OFRp.action_tolerance= 1.0e-8; + OFRp.action_degree = 18; + OFRp.md_tolerance= 1.0e-5; + OFRp.md_degree = 14; + // OFRp.degree = 20; converges + // OFRp.degree = 16; + OFRp.precision= 80; + OFRp.BoundsCheckFreq=0; + /* OneFlavourRationalParams OFRp; // Up/down OFRp.lo = 4.0e-5; OFRp.hi = 90.0; @@ -485,6 +499,23 @@ int main (int argc, char ** argv) OFRp.degree = 18; OFRp.precision= 80; OFRp.BoundsCheckFreq=0; + */ + std::vector ActionTolByPole({ + 1.0e-7,1.0e-8,1.0e-8,1.0e-8, + 1.0e-8,1.0e-8,1.0e-8,1.0e-8, + 1.0e-8,1.0e-8,1.0e-8,1.0e-8, + 1.0e-8,1.0e-8,1.0e-8,1.0e-8, + 1.0e-8,1.0e-8 + }); + std::vector MDTolByPole({ + 1.6e-5,5.0e-6,1.0e-6,3.0e-7, // soften convergence more more + // 1.0e-6,3.0e-7,1.0e-7,1.0e-7, + // 3.0e-6,1.0e-6,1.0e-7,1.0e-7, // soften convergence + 1.0e-8,1.0e-8,1.0e-8,1.0e-8, + 1.0e-8,1.0e-8,1.0e-8,1.0e-8, + 1.0e-8,1.0e-8 + }); + /* std::vector ActionTolByPole({ 1.0e-8,1.0e-8,1.0e-8,1.0e-8, 1.0e-8,1.0e-8,1.0e-8,1.0e-8, @@ -499,9 +530,9 @@ int main (int argc, char ** argv) // 1.0e-8,1.0e-8,1.0e-8,1.0e-8, 1.0e-8,1.0e-8,1.0e-8,1.0e-8, 1.0e-8,1.0e-8,1.0e-8,1.0e-8, - 1.0e-8,1.0e-8,1.0e-8,1.0e-8, 1.0e-8,1.0e-8 }); + */ OneFlavourEvenOddRatioRationalPseudoFermionAction BdySqrt(DdwfDirichlet,DdwfPeriodic,OFRp); BdySqrt.SetTolerances(ActionTolByPole,MDTolByPole); ForceTest(BdySqrt,U,DDHMCFilter);