FAST-PT: a novel algorithm to calculate convolution integrals in cosmological perturbation theory

Journal Article · Thu Sep 01 04:00:00 EDT 2016 · Journal of Cosmology and Astroparticle Physics

DOI:https://doi.org/10.1088/1475-7516/2016/09/015· OSTI ID:22679469

McEwen, Joseph E.; Fang, Xiao; Hirata, Christopher M.

We present a novel algorithm, FAST-PT, for performing convolution or mode-coupling integrals that appear in nonlinear cosmological perturbation theory. The algorithm uses several properties of gravitational structure formation—the locality of the dark matter equations and the scale invariance of the problem—as well as Fast Fourier Transforms to describe the input power spectrum as a superposition of power laws. This yields extremely fast performance, enabling mode-coupling integral computations fast enough to embed in Monte Carlo Markov Chain parameter estimation. We describe the algorithm and demonstrate its application to calculating nonlinear corrections to the matter power spectrum, including one-loop standard perturbation theory and the renormalization group approach. We also describe our public code (in Python) to implement this algorithm. The code, along with a user manual and example implementations, is available at https://github.com/JoeMcEwen/FAST-PT.

• OSTI ID: 22679469
• Journal Information: Journal of Cosmology and Astroparticle Physics, Journal Name: Journal of Cosmology and Astroparticle Physics Journal Issue: 09 Vol. 2016; ISSN 1475-7516
• Country of Publication: United States
• Language: English

Related Subjects

79 ASTRONOMY AND ASTROPHYSICS
ALGORITHMS
COMPUTERIZED SIMULATION
CORRECTIONS
COSMOLOGICAL MODELS
COUPLING
DISTURBANCES
FOURIER TRANSFORMATION
MARKOV PROCESS
MONTE CARLO METHOD
NONLINEAR PROBLEMS
NONLUMINOUS MATTER
PERTURBATION THEORY
RENORMALIZATION
SCALE INVARIANCE
SPECTRA