Microphysics

Microphysics copied to clipboard

template out T derivatives

1

If we are doing a numerical Jacobian, then we don't need to compute the temperature derivatives in sneut, screening, and the rates. We should skip their computation and also not...

zingale

create a test_rhs / test_jac unit test for SDC

2

Currently these call `actual_jac` for the Jacobian, but we want to go through our Jacobian infrastructure here.

zingale

create a temp burn_state for NSE bailout

1

something in the update of the burn_state to check for NSE causes the integration to fail at times when running with SDC+NSE. This works around that. This can be tested...

zingale

see https://github.com/AMReX-Codes/amrex/pull/2601

zingale

Add a templating equivalent to EXTRA_THERMO

2

The functionality is demonstrated in Helmholtz. In order to make test_eos work with USE_EXTRA_THERMO=FALSE, the code that actually does the work needs to be in a templated function (since if...

maxpkatz

Implement analytic Jacobian for powerlaw network

1

In addition to needing an analytic Jacobian implementation for the powerlaw network, we should also add a test of this to the regression test suite. This will allow us to...

maxpkatz

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

zingale

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

zingale

maxpkatz

Compute some Jacobian terms at compilation time

5

The current species Jacobian implementation inefficiently loops through every rate at runtime, computes whether it is used in a given Jacobian element, and is then computed. But we know at...

maxpkatz