nan values in perturbations.c

[john-houghteling]

Hi all,

I am currently running CLASS-PT on its own and noticed that ppw->theta_cb is NaN at perturbations.c:5268 for small values of tau. This does not throw an error but I believe it would change quite a few values in nonlinear_pt.c and would lead to the value of P_1loop being slightly off.

Best,
John

[john-houghteling]

I've done a bit of looking around and noticed that in perturb_total_stress_energy rho_plus_p_theta_m and rho_plus_p_m are calculated after they're used to find ppw->theta_cb. So I was able to fix the nan issue by moving the cb portion of the if block at perturbations.c:5707 to directly in front of the if block at perturbations.c:5653.