diff --git a/singularity-eos/eos/eos_base.hpp b/singularity-eos/eos/eos_base.hpp index 88124677e32..d9fe29d20f1 100644 --- a/singularity-eos/eos/eos_base.hpp +++ b/singularity-eos/eos/eos_base.hpp @@ -877,7 +877,6 @@ class EosBase { rhoguess = 0.5 * (rhomin + rhomax); } } - // JMM: Demand as much tolerance as we can, but don't reset rho below. auto status = findRoot(PofRT, press, rhoguess, rhomin, rhomax, robust::EPS(), robust::EPS(), rho); // JMM: This needs to not fail and instead return something sane. diff --git a/test/eos_unit_test_helpers.hpp b/test/eos_unit_test_helpers.hpp index 39a51873b59..80db84884ba 100644 --- a/test/eos_unit_test_helpers.hpp +++ b/test/eos_unit_test_helpers.hpp @@ -158,7 +158,7 @@ CheckRhoSieFromPT(EOS eos, Real rho, Real T, std::min(std::abs(singularity::robust::ratio(residual, P)), std::abs(residual)); bool results_good = (isClose(rho, rtest, 1e-8) && isClose(sie, etest, 1e-8)) // This is as good as it will get sometimes. - || (std::abs(frac_residual) <= 10 * singularity::robust::EPS()); + || (std::abs(frac_residual) <= 1e-8); if (!results_good) { printf("RhoSie of PT failure!\n" "\trho_true = %.14e\n"