PyCosmo: An Integrated Cosmological Boltzmann Solver

Understand

As wide-field surveys yield ever more precise measurements, cosmology has entered a phase of high precision requiring highly accurate and fast theoretical predictions.

• At the heart of most cosmological model predictions is a numerical solution of the Einstein-Boltzmann equations governing the evolution of linear perturbations in the Universe.
• We present PyCosmo, a new Python-based framework to solve this set of equations using a special pur- pose solver based on symbolic manipulations, automatic generation of C++ code and sparsity optimisation.
• The code uses a consistency relation of the field equations to adapt the time step and does not rely on physical approximations for speed-up.