A (hopefully vaguely useful) python library for cosmology.

Ian Harrison 2012-
Christopher Lovell 2013-




To use the Eisenstein & Hu fitting functions located in the EH/ directory
you must first build the extensions and locate them in your current directory.

Luckily this is all done within setup.py:

      $ python setup.py build_ext --inplace

You can then import the power and tf_fit modules:

e.g.    import EH.power as power




Example usage:

from pycosmo.cosmology import *

wmap7 = Cosmology() # creates a default cosmology

wmap7_st = Cosmology(hmf_type='st') # create a cosmology with a Sheth-Tormen hmf
wmap7_st_fnl.f_nl = 100 # add some primoridal non-Gaussianity

lnm_arr = linspace(log(1.e13), log(1.e18), 100)

# calculate z=0 mass functions for the two cosmologies
hmf1 = wmap7.dndlnm(lnm_arr, 0.e0)
hmf2 = wmap7_st_fnl.dndlnm(lnm_arr, 0.e0)

# calculate EVS of haloes in z=0 hypersurface for the two cosmologies
phi_max1, lnm_arr = evs_hypersurface_pdf(cosm=wmap7)
phi_max2, lnm_arr = evs_hypersurface_pdf(cosm=wmap7_st_fnl)