Michael Zingale

for simplified SDC, if T is outside of range do we still need to include advective terms?

1

In `vode_rhs_simplified_sdc.H`, we have a check: ``` if (state.T = MAX_TEMP) { for (int n = 1; n

it basically just dumps out the elements of the Jacobian

In the new way of automatically computing the species Jacobians we have a temporary array `spec_jac` that we will and then we copy it into the main `jac` array to...

This EOS: https://arxiv.org/abs/2104.00691v1 has a more consistent treatment of Coulomb effects The license appears to be LGPL, so I am not sure if we could directly use it, but we'd...

implement Radau integrator

1

It should be straight forward to adapt the BackwardEuler integrator to do a Radau integration, just by changing the single_step routine and the form of the timestep estimator

allow for multiple networks

4

It would be nice to have two different reaction networks built in and be able to chose which net we use in a zone based on thermodynamic conditions. E.g., for...

We should have an eos mode that takes density and entropy as inputs

add new C12 resonant rate as an option

1

We should add a runtime parameter to the aprox rate module that allows us to select the C12+C12 rate from https://www.nature.com/articles/s41586-018-0149-4