also add sanity check for rhomax vs rhomin

lanl · Jan 3, 2025 · ee5e374 · ee5e374
1 parent 7d9c90a
commit ee5e374
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/singularity-eos/eos/eos_base.hpp b/singularity-eos/eos/eos_base.hpp
@@ -863,9 +863,12 @@ class EosBase {
     constexpr Real DEFAULT_RHO_GUESS = 12;
 
     CRTP copy = *(static_cast<CRTP const *>(this));
+
     // P(rho) not monotone. When relevant, bound rhopmin.
     Real rhomin = std::max(copy.RhoPmin(temp), copy.MinimumDensity());
     Real rhomax = copy.MaximumDensity();
+    PORTABLE_REQUIRE(rhomax > rhomin, "max bound > min bound");
+
     auto PofRT = [&](const Real r) {
       return copy.PressureFromDensityTemperature(r, temp, lambda);
     };