Michael Zingale

what's the purpose of integrate_energy now?

1

It is unclear what the purpose of setting `integrator.integrate_energy = 0` is now. If we want to decouple the energy release from the temperature, e.g., to do a constant T...

switch Strang integration to work with (rho X, rho e)

3

Right now Strang integrates (X, e) and SDC integrates (rho X, rho e). We should sync them. algebraically, since rho is constant in Strang integration, this should not affect things...

With the integrators now templated, it should be straightforward to create a wrapper that can call multiple different integrators depending on the thermodynamic conditions. For example, we could try RKC...

create a eos_rh_t type

1

For MAESTROeX, it would be useful to have an `eos_rh_t` type

consider changing rhs() interface in integrator to pass y instead of integrator type

2

Right now, we pass the full integrator type (like `vode_t`) to the `rhs()` wrapper and then copy `vode_t integrator_state.y()` into the `burn_t` state to do the burn. But some integrators...

We should do a JOSS paper describing the overall Microphysics repo We'll use this issue to track thoughts. One thing that would be good to demonstrate is how it can...

We should track how many EOS calls it takes for us to converge in Newton-Raphson and see whether we can first invert the EOS in single precision and then do...

We should adopt a code of conduct for our project. I adapted the Astropy code of conduct (which is CC licensed). My intention is for this code of conduct to...

this removes old Fortran callbacks ## PR summary ## PR motivation ## PR checklist - [ ] test suite needs to be run on this PR - [ ] this...