load_default_power

randomfield.powertools.load_default_power(scaled_by_h=True)

Loads a default power spectrum P(k,z=0).

The default power spectrum was created using:

from randomfield.cosmotools import create_cosmology
cosmology = create_cosmology()
from randomfield.powertools import calculate_power
power = calculate_power(cosmology, 1e-4, 22., scaled_by_h=True)
np.savetxt('default_power.dat', power)

The range of k values used here, 1e-4 h/Mpc <= k <= 22 h/Mpc, is sufficient to cover grids with spacing >= 0.25 Mpc/h and and box dimensions <= 50 Gpc/h. If scaled_by_h is True, the units are h/Mpc for k and (Mpc/h)**3 for P(k). Otherwise, replace Mpc/h with Mpc in the units.