skypyproject / skypy

SkyPy: A package for modelling the Universe.
BSD 3-Clause "New" or "Revised" License
118 stars 38 forks source link

Test power spectra against "truth", both massless and massive nu #570

Open itrharrison opened 1 year ago

itrharrison commented 1 year ago

Description

This PR addresses #567 to improve power spectrum tests. I have made the following changes to the module:

And implemented new / changed the tests such that both CLASS and CAMB are tested against a 'truth' cosmology which is calculated from CAMB, with separate tests for massive and massless neutrinos. We could do a lot of playing with accuracy settings here, both in terms of the rtol in the test and the accuracy settings within each Boltzmann code. Currently these are specified so that the tests pass with the default accuracy settings for each Boltzmann code. Thoughts welcome.

Merging will close #567.

Checklist

itrharrison commented 1 year ago

Here are the two rtol plots for CAMB and CLASS for both the massless and massive (m_nu=[0.0, 0.0, 0.06]) cases.

camb-numodels class-numodels