def test_single_vs_array(): #First sigma2 a1 = peaks.sigma2_at_M(Ma, klin, plin, Omega_m) a2 = np.array([peaks.sigma2_at_M(Mi, klin, plin, Omega_m) for Mi in Ma]) npt.assert_array_equal(a1, a2) a1 = peaks.sigma2_at_R(Ra, klin, plin) a2 = np.array([peaks.sigma2_at_R(Ri, klin, plin) for Ri in Ra]) npt.assert_array_equal(a1, a2)
def test_special(): sigma2 = peaks.sigma2_at_M(M, k, p, Omega_m) Mt = M * (1 - 1e-6 * 0.5) Mb = M * (1 + 1e-6 * 0.5) sigma2t = peaks.sigma2_at_M(Mt, k, p, Omega_m) sigma2b = peaks.sigma2_at_M(Mb, k, p, Omega_m) d, e, f, g = 2.16087369917, 1.18309392312, 0.133881834517, -0.0263615354323 n = mf._dndM_sigma2_precomputed(M, sigma2, sigma2t, sigma2b, Omega_m, d, e, f, g) npt.assert_array_less(n[1:], n[:-1])
def test_s2_and_nu_functions(): #Test the mass calls s2 = peaks.sigma2_at_M(Mass, klin, plin, Omega_m) nu = peaks.nu_at_M(Mass, klin, plin, Omega_m) npt.assert_equal(1.686 / np.sqrt(s2), nu) s2 = peaks.sigma2_at_M(Ma, klin, plin, Omega_m) nu = peaks.nu_at_M(Ma, klin, plin, Omega_m) npt.assert_array_equal(1.686 / np.sqrt(s2), nu) #Now test the R calls R = 1.0 #Mpc/h; arbitrary s2 = peaks.sigma2_at_R(R, klin, plin) nu = peaks.nu_at_R(R, klin, plin) npt.assert_equal(1.686 / np.sqrt(s2), nu)
def test_G(): sigma = np.sqrt(peaks.sigma2_at_M(M, k, p, Omega_m)) Gm = mf.G_at_M(M, k, p, Omega_m) Gs = mf.G_at_sigma(sigma) R = Gm / Gs ones = np.ones_like(R) npt.assert_array_equal(ones, R)
def test_derivatives(): M = np.logspace(13, 14, 2000) ds2dM = peaks.dsigma2dM_at_M(M, klin, plin, Omega_m) s2 = peaks.sigma2_at_M(M, klin, plin, Omega_m) dM = M[1:] - M[:-1] ds2 = s2[1:] - s2[:-1] deriv = ds2 / dM pd = ds2dM[:-1] / deriv npt.assert_array_almost_equal(pd, np.ones_like(pd), 1e-3) return
def test_special(): sigma2 = peaks.sigma2_at_M(M, k, p, Omega_m) dsigma2dM = peaks.dsigma2dM_at_M(M, k, p, Omega_m) d, e, f, g = 2.16087369917, 1.18309392312, 0.133881834517, -0.0263615354323 n = mf._dndM_sigma2_precomputed(M, sigma2, dsigma2dM, Omega_m, d, e, f, g) npt.assert_array_less(n[1:], n[:-1])