def Pk_bbks(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 BBKS
Ref: Bardeen et. al., Astrophys. J., 304, 15 (1986)
"""
return self.A0bbks() * k**self.ns * tf.Tbbks(
k, self.Omega_m_z(0), self.Ob0, self.h)**2. * self.D_plus_z(z)**2.
def integrand_bbks(self, k):
return k**(self.ns + 2.) * tf.Tbbks(k, self.Omega_m_z(0), self.Ob0,
self.h)**2. * (self.Wf(k))**2.