addressing comments from Danielle and Matt
@@ -257,7 +257,8 @@ def get_isitgr_pk_lin(cosmo):
delta_neff = cosmo['Neff'] - 3.046 # used for BBN YHe comps
- # ISiTGR built on CAMB which defines a neutrino degeneracy factor as T_i = g^(1/4)*T_nu
+ # ISiTGR built on CAMB which defines a neutrino degeneracy
+ # factor as T_i = g^(1/4)*T_nu
# where T_nu is the standard neutrino temperature from first order
# computations
# CLASS defines the temperature of each neutrino species to be
@@ -337,10 +338,9 @@ def get_isitgr_pk_lin(cosmo):
extrap_order_lok=1,
extrap_order_hik=2,
cosmo=cosmo)
-
-
return pk_lin
+
def get_class_pk_lin(cosmo):
"""Run CLASS and return the linear power spectrum.
@@ -457,4 +457,5 @@ def get_class_pk_lin(cosmo):
extrap_order_lok=1,
extrap_order_hik=2,
cosmo=cosmo)
+
return pk_lin
