Nonlinear perturbation theory extension of the Boltzmann code CLASS

We present a new open-source code that calculates one-loop power auto- and cross-power spectra for matter fields and biased tracers in real and redshift space. These spectra incorporate all ingredients required for a direct application to data: non-linear bias, redshift-space distortions, infra-red resummation, counterterms, and the Alcock-Paczynski effect. Our code is based on the Boltzmann solver CLASS and inherits its advantage: user friendliness, ease of modification, high speed, and simple interface with other software. We present detailed descriptions of the theoretical model, the code structure, approximations, and accuracy tests. A typical end-to-end run for one cosmology takes $\sim 0.3$ seconds, which is sufficient for Markov Chain Monte Carlo parameter extraction. As an example, we apply the code to data from the Baryon Oscillation Spectroscopic Survey (BOSS) and infer cosmological parameters from the shape of the galaxy power spectrum.