From b9ae0d78147d6329cc88f7a6facbfb624e315b79 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Wed, 9 Sep 2026 12:55:12 -0400 Subject: [PATCH 1/2] have integrator diagnostics respect subtract_internal_energy --- integration/integrator_setup_strang.H | 12 ++++++++++-- 1 file changed, 10 insertions(+), 2 deletions(-) diff --git a/integration/integrator_setup_strang.H b/integration/integrator_setup_strang.H index 6793630e18..ee6d9af857 100644 --- a/integration/integrator_setup_strang.H +++ b/integration/integrator_setup_strang.H @@ -179,7 +179,11 @@ void integrator_cleanup (IntegratorT& int_state, BurnT& state, // Print out some integration statistics, if desired. std::cout << "integration summary: " << std::endl; std::cout << "dens: " << state.rho << " temp: " << state.T << std::endl; - std::cout << " energy released: " << state.e << std::endl; + if (integrator_rp::subtract_internal_energy) { + std::cout << " energy released: " << state.e << std::endl; + } else { + std::cout << " energy released: " << state.e - state_save.e_in << std::endl; + } std::cout << "number of steps taken: " << state.n_step << std::endl; std::cout << "number of f evaluations: " << state.n_rhs << std::endl; } @@ -210,7 +214,11 @@ void integrator_cleanup (IntegratorT& int_state, BurnT& state, std::cout << std::scientific << std::setprecision(OUTDIGITS) << X << " "; } std::cout << std::endl; - std::cout << "energy generated = " << state.e << std::endl; + if (integrator_rp::subtract_internal_energy) { + std::cout << "energy generated = " << state.e << std::endl; + } else { + std::cout << "energy generated = " << state_save.e_in << std::endl; + } #elif defined(ALLOW_GPU_PRINTF) AMREX_DEVICE_PRINTF("[ERROR] integration failed in net, istate = %d, time = %g, dt = %g, dens = %g, temp_start = %g\n", istate, int_state.t, dt, state.rho, state_save.T_in); From c17e6ca3f1a3c015eb5c44282e0febfbebca3e3a Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Fri, 11 Sep 2026 12:14:27 -0400 Subject: [PATCH 2/2] fix --- integration/integrator_setup_strang.H | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/integration/integrator_setup_strang.H b/integration/integrator_setup_strang.H index ee6d9af857..efe6425af8 100644 --- a/integration/integrator_setup_strang.H +++ b/integration/integrator_setup_strang.H @@ -217,7 +217,7 @@ void integrator_cleanup (IntegratorT& int_state, BurnT& state, if (integrator_rp::subtract_internal_energy) { std::cout << "energy generated = " << state.e << std::endl; } else { - std::cout << "energy generated = " << state_save.e_in << std::endl; + std::cout << "energy generated = " << state.e - state_save.e_in << std::endl; } #elif defined(ALLOW_GPU_PRINTF) AMREX_DEVICE_PRINTF("[ERROR] integration failed in net, istate = %d, time = %g, dt = %g, dens = %g, temp_start = %g\n",