Michael Zingale

just to understand -- there is nothing special about `networks/ase` except that you added neutrons to it, right? This method should work for any pynucastro net that has neutrons? (can...

can you add some documentation to https://amrex-astro.github.io/Microphysics/docs/dev/nse.html#self-consistent-nse explaining how to use this, including what requirements there are on networks.

tests pass: http://groot.astro.sunysb.edu/Microphysics/test-suite/gfortran/2022-11-30-001/index.html

also now `use_number_densities`

going to leave `integrate_energy` for #1215

is this still something we want?

closing this -- we have a `test_part_func` unit test now.

I think we are effectively doing this now for all of the pynucastro nets

interesting question: since screening depends on the plasma state, does this mean that the Jacobian will never be sparse if we differentiate the screening with composition automagically.

yes, I realize that. I just thought it was neat. Never thought about this part before...