ccl_test_growth Namespace Reference

Functions
def	read_growth_lowz_benchmark_file ()
def	read_growth_allz_benchmark_file ()
def	compare_growth (z, gfac_bench, Omega_v, w0, wa)
def	test_growth_lowz_model_0 ()
def	test_growth_lowz_model_1 ()
def	test_growth_lowz_model_2 ()
def	test_growth_lowz_model_3 ()
def	test_growth_lowz_model_4 ()
def	test_growth_allz_model_0 ()
def	test_growth_allz_model_1 ()
def	test_growth_allz_model_2 ()
def	test_growth_allz_model_3 ()
def	test_growth_allz_model_4 ()
def	test_mgrowth ()

Variables
int	GROWTH_TOLERANCE = 1
float	Omega_c = 0.25
float	Omega_b = 0.05
int	Neff = 0
int	m_nu = 0
float	h = 0.7
float	A_s = 2.1e-9
float	n_s = 0.96
	Omega_v_vals = np.array([0.7, 0.7, 0.7, 0.65, 0.75])
	w0_vals = np.array([-1.0, -0.9, -0.9, -0.9, -0.9])
	wa_vals = np.array([0.0, 0.0, 0.1, 0.1, 0.1])

Function Documentation

def ccl_test_growth.compare_growth	(		z,
			gfac_bench,
			Omega_v,
			w0,
			wa
	)

Compare growth factor calculated by pyccl with the values in the benchmark
file. This test only works if radiation is explicitly set to 0.

Definition at line 53 of file ccl_test_growth.py.

Referenced by test_growth_allz_model_0(), test_growth_allz_model_1(), test_growth_allz_model_2(), test_growth_allz_model_3(), test_growth_allz_model_4(), test_growth_lowz_model_0(), test_growth_lowz_model_1(), test_growth_lowz_model_2(), test_growth_lowz_model_3(), and test_growth_lowz_model_4().

def compare_growth(z, gfac_bench, Omega_v, w0, wa):
    """
    Compare growth factor calculated by pyccl with the values in the benchmark
    file. This test only works if radiation is explicitly set to 0.
    """

    # Set Omega_K in a consistent way
    Omega_k = 1.0 - Omega_c - Omega_b - Omega_v

    cosmo = ccl.Cosmology(Omega_c=Omega_c, Omega_b=Omega_b, Neff=Neff, m_nu=m_nu,
                       h=h, A_s=A_s, n_s=n_s, Omega_k=Omega_k,
                       w0=w0, wa=wa, Omega_g=0)

    # Calculate distance using pyccl
    a = 1. / (1. + z)
    gfac = ccl.growth_factor_unnorm(cosmo, a)

    # Compare to benchmark data
    assert_allclose(gfac, gfac_bench, atol=1e-12, rtol=GROWTH_TOLERANCE)

def ccl_test_growth.read_growth_allz_benchmark_file

(

)

Read the file containing growth factor benchmarks for the whole redshift range.

Definition at line 36 of file ccl_test_growth.py.

References read_growth_lowz_benchmark_file().

def read_growth_allz_benchmark_file():
    """
    Read the file containing growth factor benchmarks for the whole redshift range.
    """
    # Load data from file
    dat = np.genfromtxt(join(dirname(__file__),"benchmark/growth_cosmomad_allz.txt")).T
    assert(dat.shape == (6,10))

    # Split into redshift column and growth(z) columns
    z = dat[0]
    gfac = dat[1:]
    return z, gfac

# Set-up test data
z_lowz, gfac_lowz = read_growth_lowz_benchmark_file()
z_allz, gfac_allz = read_growth_allz_benchmark_file()

def ccl_test_growth.read_growth_lowz_benchmark_file

(

)

Read the file containing growth factor benchmarks for the low redshifts.

Definition at line 23 of file ccl_test_growth.py.

Referenced by read_growth_allz_benchmark_file().

def read_growth_lowz_benchmark_file():
    """
    Read the file containing growth factor benchmarks for the low redshifts.
    """
    # Load data from file
    dat = np.genfromtxt(join(dirname(__file__),"benchmark/growth_model1-5.txt")).T
    assert(dat.shape == (6,6))

    # Split into redshift column and growth(z) columns
    z = dat[0]
    gfac = dat[1:]
    return z, gfac

def ccl_test_growth.test_growth_allz_model_0

(

)

Definition at line 95 of file ccl_test_growth.py.

References compare_growth().

def test_growth_allz_model_0():
    i = 0
    compare_growth(z_allz, gfac_allz[i], Omega_v_vals[i], w0_vals[i], wa_vals[i])

def ccl_test_growth.test_growth_allz_model_1

(

)

Definition at line 99 of file ccl_test_growth.py.

References compare_growth().

def test_growth_allz_model_1():
    i = 1
    compare_growth(z_allz, gfac_allz[i], Omega_v_vals[i], w0_vals[i], wa_vals[i])

def ccl_test_growth.test_growth_allz_model_2

(

)

Definition at line 103 of file ccl_test_growth.py.

References compare_growth().

def test_growth_allz_model_2():
    i = 2
    compare_growth(z_allz, gfac_allz[i], Omega_v_vals[i], w0_vals[i], wa_vals[i])

def ccl_test_growth.test_growth_allz_model_3

(

)

Definition at line 107 of file ccl_test_growth.py.

References compare_growth().

def test_growth_allz_model_3():
    i = 3
    compare_growth(z_allz, gfac_allz[i], Omega_v_vals[i], w0_vals[i], wa_vals[i])

def ccl_test_growth.test_growth_allz_model_4

(

)

Definition at line 111 of file ccl_test_growth.py.

References compare_growth().

def test_growth_allz_model_4():
    i = 4
    compare_growth(z_allz, gfac_allz[i], Omega_v_vals[i], w0_vals[i], wa_vals[i])

def ccl_test_growth.test_growth_lowz_model_0

(

)

Definition at line 74 of file ccl_test_growth.py.

References compare_growth().

def test_growth_lowz_model_0():
    i = 0
    compare_growth(z_lowz, gfac_lowz[i], Omega_v_vals[i], w0_vals[i], wa_vals[i])

def ccl_test_growth.test_growth_lowz_model_1

(

)

Definition at line 78 of file ccl_test_growth.py.

References compare_growth().

def test_growth_lowz_model_1():
    i = 1
    compare_growth(z_lowz, gfac_lowz[i], Omega_v_vals[i], w0_vals[i], wa_vals[i])

def ccl_test_growth.test_growth_lowz_model_2

(

)

Definition at line 82 of file ccl_test_growth.py.

References compare_growth().

def test_growth_lowz_model_2():
    i = 2
    compare_growth(z_lowz, gfac_lowz[i], Omega_v_vals[i], w0_vals[i], wa_vals[i])

def ccl_test_growth.test_growth_lowz_model_3

(

)

Definition at line 86 of file ccl_test_growth.py.

References compare_growth().

def test_growth_lowz_model_3():
    i = 3
    compare_growth(z_lowz, gfac_lowz[i], Omega_v_vals[i], w0_vals[i], wa_vals[i])

def ccl_test_growth.test_growth_lowz_model_4

(

)

Definition at line 90 of file ccl_test_growth.py.

References compare_growth().

def test_growth_lowz_model_4():
    i = 4
    compare_growth(z_lowz, gfac_lowz[i], Omega_v_vals[i], w0_vals[i], wa_vals[i])

# 0.01 < z < 1000 tests

def ccl_test_growth.test_mgrowth

(

)

Compare the modified growth function computed by CCL against the exact
result for a particular modification of the growth rate.

Definition at line 115 of file ccl_test_growth.py.

References Catch::Generators.range().

def test_mgrowth():
    """
    Compare the modified growth function computed by CCL against the exact
    result for a particular modification of the growth rate.
    """
    # Define differential growth rate arrays
    nz_mg = 128
    z_mg = np.zeros(nz_mg)
    df_mg = np.zeros(nz_mg)
    for i in range(0, nz_mg):
        z_mg[i] = 4. * (i + 0.0) / (nz_mg - 1.)
        df_mg[i] = 0.1 / (1. + z_mg[i])

    # Define two test cosmologies, without and with modified growth respectively
    cosmo1 = ccl.Cosmology(Omega_c=0.25, Omega_b=0.05, Omega_k=0., Neff=0., m_nu=0.,
                        w0=-1., wa=0., h=0.7, A_s=2.1e-9, n_s=0.96)
    cosmo2 = ccl.Cosmology(Omega_c=0.25, Omega_b=0.05, Omega_k=0., Neff=0., m_nu=0.,
                        w0=-1., wa=0., h=0.7, A_s=2.1e-9, n_s=0.96,
                        z_mg=z_mg, df_mg=df_mg)

    # We have included a growth modification \delta f = K*a, with K==0.1
    # (arbitrarily). This case has an analytic solution, given by
    # D(a) = D_0(a)*exp(K*(a-1)). Here we compare the growth computed by CCL
    # with the analytic solution.
    a = 1. / (1. + z_mg)

    d1 = ccl.growth_factor(cosmo1, a)
    d2 = ccl.growth_factor(cosmo2, a)
    f1 = ccl.growth_rate(cosmo1, a)
    f2 = ccl.growth_rate(cosmo2, a)

    f2r = f1 + 0.1*a
    d2r = d1 * np.exp(0.1*(a-1.))

    # Check that ratio of calculated and analytic results is within tolerance
    assert_allclose(d2r, d2, rtol=GROWTH_TOLERANCE)
    assert_allclose(f2r, f2, rtol=GROWTH_TOLERANCE)

Variable Documentation

float ccl_test_growth.A_s = 2.1e-9

Definition at line 15 of file ccl_test_growth.py.

int ccl_test_growth.GROWTH_TOLERANCE = 1

Definition at line 6 of file ccl_test_growth.py.

float ccl_test_growth.h = 0.7

Definition at line 14 of file ccl_test_growth.py.

int ccl_test_growth.m_nu = 0

Definition at line 13 of file ccl_test_growth.py.

float ccl_test_growth.n_s = 0.96

Definition at line 16 of file ccl_test_growth.py.

int ccl_test_growth.Neff = 0

Definition at line 12 of file ccl_test_growth.py.

float ccl_test_growth.Omega_b = 0.05

Definition at line 11 of file ccl_test_growth.py.

float ccl_test_growth.Omega_c = 0.25

Definition at line 10 of file ccl_test_growth.py.

ccl_test_growth.Omega_v_vals = np.array([0.7, 0.7, 0.7, 0.65, 0.75])

Definition at line 19 of file ccl_test_growth.py.

ccl_test_growth.w0_vals = np.array([-1.0, -0.9, -0.9, -0.9, -0.9])

Definition at line 20 of file ccl_test_growth.py.

ccl_test_growth.wa_vals = np.array([0.0, 0.0, 0.1, 0.1, 0.1])

Definition at line 21 of file ccl_test_growth.py.