Michael Zingale

At the moment, each network has its own binding energy and computes mion it. But we can automate this by having the python script that generates the `network_header.H` use the...

add neutrino losses to aprox13

2

This comes out of discussions with Sam Jones, Aron Michel, and @carlnotsagan We should implement the neutrino losses from the weak reactions. This would mean keeping track of each reaction...

add the Potekhin and Chabrier (2010,2013) EOS

3

The Potekhin and Chabrier (2010,2013) EOS is freely available (in Fortran) here: http://www.ioffe.ru/astro/EIP/ This EOS seems to be preferred by MESA: https://docs.mesastar.org/en/r15140/eos/overview.html There is also support here for a magnetized...

other 3-alpha rates?

1

We are still using an old rate for 3-alpha in most nets. Some updates are found in: Angulo, C., Arnould, M., Rayet, M., Descouvemont, P., Baye, D., et al. (1999):...

this addresses #762

we should be able to query an EOS to see if an output is supported

1

Not all EOSes fill all the thermodynamic data. We should be able to query whether something is supported? Like `is_filled("ne")`.

analytic Jacobian species derivatives are missing terms

1

We want to express our Jacobian in terms of (rho, e, X_k), but the reaction system is more naturally in terms of (ρ, T, X_k). Consider dX/dt = f(ρ, T,...

what if we pass VODE failure message out of burn through burn_t

4

Currently we don't dump out the burn failure when we are using CUDA, which makes debugging hard on GPUs. What if we add a `std::string` to the `burn_t` that can...

It would be good to have a script that ran a single zone burn for a given network using all of the integrators and plotted the mass fraction evolution vs....

use C++ autodiff to do Jacobians

4

Instead of hand-coding the Jacobians, we should explore using a C++ autodiff package, e.g., https://autodiff.github.io/ to construct the analytic Jacobians.