From f1b1179bfdfd30290a1b51397ca96c8b32421d53 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Wed, 17 Jul 2024 12:55:09 -0400 Subject: [PATCH] fix test_jac numerical Jacobian we should not zero out the energy term before calling numerical_jacobian --- unit_test/test_jac/jac_zones.H | 10 ++-------- 1 file changed, 2 insertions(+), 8 deletions(-) diff --git a/unit_test/test_jac/jac_zones.H b/unit_test/test_jac/jac_zones.H index c00cc9348a..afb4616e23 100644 --- a/unit_test/test_jac/jac_zones.H +++ b/unit_test/test_jac/jac_zones.H @@ -26,17 +26,11 @@ bool do_jac (int i, int j, int k, amrex::Array4 const& state, const normalize_abundances_burn(burn_state); - eos_extra_t eos_state; - - burn_to_eos(burn_state, eos_state); - - eos(eos_input_rt, eos_state); - - eos_to_burn(eos_state, burn_state); + eos(eos_input_rt, burn_state); // the integrator doesn't actually care about the initial internal // energy. - burn_state.e = 0.0_rt; + //burn_state.e = 0.0_rt; JacNetArray2D jac_analytic; JacNetArray2D jac_numerical;