The current test of classy power spectra fails, due to the interpolation creating >% level errors with respect to the previous function implementation, depending on the redshift.

We need to decide how to define a 'truth' power spectrum in order to create test data to validate which, if either of these is more correct.

This would probably require:

• creating a baseline P(k,z) where we had consistency at the 10^-4 level between camb and class run natively
• using this P(k,z) as the test data for our implementations which include interpolation