Skip to content

have integrator diagnostics respect subtract_internal_energy #2089

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
zingale merged 3 commits into from
Sep 14, 2026
Merged

have integrator diagnostics respect subtract_internal_energy #2089

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
b9ae0d7
have integrator diagnostics respect subtract_internal_energy
zingale Sep 9, 2026
07d36fd
Merge branch 'development' into fix_energy_released
zingale Sep 11, 2026
c17e6ca
fix
zingale Sep 11, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
12 changes: 10 additions & 2 deletions integration/integrator_setup_strang.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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;
}
Expand Down Expand Up @@ -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.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",
istate, int_state.t, dt, state.rho, state_save.T_in);
Expand Down
Loading