We implement our method for galaxy clustering and galaxy-galaxy lensing power spectra. It is easily extended to handle integrals involving derivatives of Bessel functions, making it equally applicable to numerically more challenging cases such as contributions from redshift-space distortions and Doppler effects. The new method significantly simplifies the calculation and improves the numerical speed and stability. Here, we present a new method based on our generalized FFTLog algorithm for the efficient computation of angular power spectra beyond the Limber approximation. To reach the high accuracy required by the statistical precision of future surveys, some of these statistics may need to be computed without the commonly employed Limber approximation the exact computation however requires integration over Bessel functions, and a brute-force evaluation is slow to converge. Angular two-point statistics of large-scale structure observables are important cosmological probes.
0 Comments
Leave a Reply. |