tighter tolerances

lanl · Jan 3, 2025 · ceec265 · ceec265
1 parent 1d1873f
commit ceec265
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/singularity-eos/eos/eos_base.hpp b/singularity-eos/eos/eos_base.hpp
@@ -877,8 +877,9 @@ class EosBase {
         rhoguess = 0.5 * (rhomin + rhomax);
       }
     }
-    auto status = findRoot(PofRT, press, rhoguess, rhomin, rhomax, 10 * robust::EPS(),
-                           10 * robust::EPS(), rho);
+    // 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.
     // If root find failed to converge, density will at least be
     // within bounds.