Fast computation of temperature and polarization coupling matrices

Abstract

We present a fast and exact method for computing CMB mode-coupling matrices based on an optimised evaluation of Wigner-3j symbols. The method exploits analytic structure in the relevant Wigner-3j symbol configurations appearing in temperature and polarization coupling matrices, expressing all required quantities in terms of a small set of recurrence-generated values which are precomputed and stored in lookup tables. This approach reduces the computational cost of constructing the full coupling matrices whilst maintaining numerical accuracy. We demonstrate the performance of the threej_cosmo implementation using realistic survey masks from current CMB experiments. Relative to standard recursion-based approaches used in existing pseudo-C_l pipelines, the method achieves speedups of 6-25x in practical coupling-matrix constructions, with the largest gains occurring at high multipoles. The algorithm admits efficient parallelisation on both CPUs and GPUs, the latter providing additional acceleration, up to a further order of 50 on modern hardware, without altering the underlying formalism. Beyond full matrix construction, the approach is naturally suited to applications in which only a restricted set of l3 modes is required for each (l1,l2) pair, such as in the computation of band-limited coupling matrices and analytic covariance terms. These features make threej_cosmo a practical backend for pseudo-C_l estimation and related calculations in next-generation CMB analysis pipelines.

9 pages, 2 figures

References in corpus (15)

• Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter Results
• The Simons Observatory: Science goals and forecasts
• MASTER of the CMB Anisotropy Power Spectrum: A Fast Method for Statistical Analysis of Large and Complex CMB Data Sets
• CMB Lensing Reconstruction on the Full Sky
• A unified pseudo-$C_{\ell }$ framework
• Planck 2013 results. XVII. Gravitational lensing by large-scale structure
• Planck 2013 results. XXVII. Doppler boosting of the CMB: Eppur si muove
• Myths and Truths Concerning Estimation of Power Spectra
• A Detailed Description of the CamSpec Likelihood Pipeline and a Reanalysis of the Planck High Frequency Maps
• CMB temperature and polarisation pseudo-Cl estimators and covariances
• Pseudo-$C_{\ell }$ estimators which do not mix E and B modes
• Polarized CMB power spectrum estimation using the pure pseudo-cross-spectrum approach
• Fast computation of angular power spectra and covariances of high-resolution cosmic microwave background maps using the Toeplitz approximation
• Fast Projected Bispectra: the filter-square approach
• Pseudo-$C_{\ell }$s for spin-$s$ fields with component-wise weighting