Microphysics
Microphysics copied to clipboard
Nuclear Partition Function Test Case
Don and I will go through this together soon.
The main result:
Edot Ratio (PF/NOPF): 2.293683464
The estimated for 1e10k is 2.294
This pull request contains two main things.
- A single rate network.
-
test_pf
This test is basicallytest_rhs
run twice. Once withuse_partition_factors = .true.
and once withuse_partition_factors = .false.
We then construct the rhs and compare the two Edots. They should differ by the partition function factor of mg23 at the specific temperature. For now, I have hardcoded the solution into the code to compare to. The solution was generated usingtest_pf/solution.ipynb
closing this -- we have a test_part_func
unit test now.