Michael Zingale
Michael Zingale
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...
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
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
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