FuGa3D: Fast full-sky analysis of Galaxy catalogs in 3D

We present FuGa3D, a code for fast computation of correlation functions and power spectra for galaxy survey observables, including galaxy clustering and cosmic shear. We define the redshift-space correlation function (RCF) as the correlation function defined in the parameter space of two redshifts and an angular separation angle. Assuming that there is no preferred direction in the sky, these parameters fully define the relative position of two galaxies, independently of the assumed cosmological model. Once the RCF is constructed, it is easy to compute derived correlation metrics, such as the real-space clustering correlation function and its multipoles. We further define the redshift-space power spectrum as the harmonic counterpart of the RCF, and show that it can be computed efficiently using the discrete galaxy coordinates. We validate the code with simulated mock catalogs. Computing the RCF and the two-point correlation function at 1.5 Mpc (3 Mpc) resolution for a MICE simulation with 46 million galaxies, took 47 node-min for clustering only, and 7.3 node-hours with shear analysis included.