API

log of the gamma function

same as above but suitable for use in a piecewise integrator. re-use the moment I(ka) which was I(kb) previously

Convenience function to create an ionisation history from some tables

integrates a quadratic polynomial with moments. input f, f′, f″ are those evaluated at a, over interval (a,b). Uses a preconstructed interpolator for the moments.

peebles_Xₑ(par, Xₑ₀, x_start, x_end)

Solve the Peebles equation over a span of scale factors, and then construct an interpolator mapping scale factor to the resulting ionization fraction.

Arguments:

• par: cosmological parameters
• Xₑ₀: initial ionization fraction
• x_start: scale factor to begin integration
• x_end: scale factor to end integration

Returns:

• generic function: interpolator for Xₑ(x)

saha_peebles_recombination(par::AbstractCosmoParams)

Utility function for generating a decent approximation to Xₑ in ΛCDM recombination, using the Saha equation until z=1587.4 and then the Peebles equation for the rest.