def Pk_wh(self, k, z):
"""
Matter Power Spectra Pk in units if h^{-3} Mpc^{3} as a function of k in units of [h Mpc^{-1}]
and z;
Transfer function is taken to be Eisenstein & Hu
(ref(Eisenstein and Hu, Astrophys. J., 496, 605 (1998)))
"""
return self.A0wh() * k**self.ns * tf.Twh(
k, self.Omega_m_z(0), self.Ob0, self.h)**2. * self.D_plus_z(z)**2.
def integrand_wh(self, k):
return k**(self.ns + 2.) * tf.Twh(k, self.Omega_m_z(0), self.Ob0,
self.h)**2. * (self.Wf(k))**2.