CCL: Neutrino masses error handling
Just running:
cosmo = ccl.Cosmology(Omega_c=0.27, Omega_b=0.045, h=0.67,
sigma8=0.8, n_s=0.96, w0=-1., wa=-0.0,
m_nu=0.05)
gets into an infinite loop that you can only Ctrl+C
out of. This will happen whenever you pass m_nu<0.06.
- Why is this happening? It is a bit weird that m_nu=0 is fine, but the gap between 0 and 0.06 is not. In principle people should be able to put in whatever positive neutrino masses they want, even if they’re unphysical (and I’d like to understand what “unphysical” means exactly - maybe @c-d-leonard knows?).
- In any case, it should not be possible to fall into that infinite loop. Looking at the code, I think it should be very easy to sort that out.
About this issue
- Original URL
- State: closed
- Created 5 years ago
- Comments: 15 (1 by maintainers)
Agreed, these issues keep coming up. I’m going to make a PR with the current thing so the bug fix gets in; seems like having an informative error message may have to wait for those bigger changes.
@aferte - insttructions here: https://confluence.slac.stanford.edu/display/LSSTDESC/Automatic+Maintenance+of+GitHub+Teams