Fast calculation of the nonlinear redshift-space galaxy power spectrum including selection bias

We present a fast implementation of the next-to-leading order (one-loop) redshift-space galaxy power spectrum by using FFTLog-based methods. Desjacques et al. [J. Cosmol. Astropart. Phys. 12 (2018) 035JCAPBP1475-751610.1088/1475-7516/2018/12/035] have shown that the one-loop galaxy power spectrum in redshift space can be computed with 28 independent loop integrals with 22 bias parameters. Analytical calculation of the angular part of the loop integrals leaves the radial part in the form of a spherical Bessel transformation that is ready to be integrated numerically by using the FFTLog transformation. We find that the original 28 loop integrals can be solved with a total of 85 unique FFTLog transformations, yet leading to a few orders of magnitude speed up over traditional multidimensional integration.