class TestPower(object):
def __init__(self):
# Make a pert class
self.pert = Perturbations(omegab=0.05, omegac=0.25, omegav=0.7, sigma_8=0.8, n=1, H0=70.0,
k_bounds=[np.exp(-21), np.exp(21)], transfer__kmax=10, transfer__k_per_logint=50)
#Get the camb transfer
self.camb_transfer = np.genfromtxt("data/camb_transfer")
#Get the camb power spec
self.camb_power = np.genfromtxt('data/camb_power')
#Get the genmf power (ie. camb power normalised by genmf)
self.genmf_power = np.genfromtxt('data/genmf_power')
def test_transfer(self):
""" Testing whether the transfer function calculated by pycamb is the same as that of camb (via CLI)"""
my_T_vec = self.pert._transfer_function_callable(np.log(self.camb_transfer[:, 0]))
assert rms_diff(my_T_vec, np.log(self.camb_transfer[:, 1]), 1E-3)
