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])
