make the assertion on rho, e more informative (compressible solvers) (#…

…328)
python-hydro · Jan 24, 2025 · ab26c9c · ab26c9c
1 parent 49f7447
commit ab26c9c
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/pyro/compressible/simulation.py b/pyro/compressible/simulation.py
@@ -65,8 +65,10 @@ def cons_to_prim(U, gamma, ivars, myg):
                   out=np.zeros_like(U[:, :, ivars.iener]),
                   where=(U[:, :, ivars.idens] != 0.0))
 
-    assert e.v().min() > 0.0
-    assert q.v(n=ivars.irho).min() > 0.0
+    e_min = e.v().min()
+    rho_min = q.v(n=ivars.irho).min()
+
+    assert e_min > 0.0 and rho_min > 0.0, f"invalid state, min(rho) = {rho_min}, min(e) = {e_min}"
 
     q[:, :, ivars.ip] = eos.pres(gamma, q[:, :, ivars.irho], e)