add neutrino cooling routines #628
Conversation
this will allow us to compute ignition curves
zingale
changed the title
|
Here's a plot made by this new functionality:
|
|
Check out this pull request on See visual diffs & provide feedback on Jupyter Notebooks. Powered by ReviewNB |
zingale
changed the title
|
okay there are unit tests and documentation now. I think this is ready for review. Note: I reorganized the docs a bit -- making more sections. I anticipate the NSE and neutrino sections growing as we add more discussions there. |
|
It's probably a good idea to pull the different contributions out into helper functions/closures. That should limit the scope of the temporary variables, and help pylint figure out what's unused. |
|
I think I removed all the terms Zhi flagged I agree that this needs to be refactored, but I wanted to get something in and tested before doing that. I'd like to do the same with the C++ version. |
|
updated |
| def zfermim12(x): | ||
| """apply a rational function expansion to get the fermi-dirac | ||
| integral of order -1/2 evaluated at x. maximum error is 1.23e-12_rt. | ||
| reference: antia apjs 84,101 1993""" |
There was a problem hiding this comment.
I think this entire function can be removed
|
fixed the remaining ones |
this will allow us to compute ignition curves