Accelerating spherical harmonic transforms for a large number of sky maps

Abstract

The spherical harmonic transform is a powerful tool in the analysis of spherical data sets, such as the cosmic microwave background data. In this work, we present a new scheme for the spherical harmonic transforms that supports both CPU and GPU computations, which is specially efficient on a large number of sky maps. By comparing our implementation with the standard Libsharp-HEALPix program, we demonstrate 2-10 times speedup for the CPU implementation, and up to 30 times speedup when a state-of-the-art GPU is employed. This new scheme's software package is available via an open source GitHub repository.

arXiv commentOpen access. The software package is available at https://github.com/liuhao-cn/fastSHT

References in corpus (16)

• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• The Simons Observatory: Science goals and forecasts
• An All-Sky Analysis of Polarization in the Microwave Background
• Efficient Spherical Harmonic Transforms aimed at pseudo-spectral numerical simulations
• ACTPol: A polarization-sensitive receiver for the Atacama Cosmology Telescope
• Gauss--Legendre Sky Pixelization (GLESP) for CMB maps
• SPTpol: an instrument for CMB polarization measurements with the South Pole Telescope
• Probing Primordial Gravitational Waves: Ali CMB Polarization Telescope
• Libsharp - spherical harmonic transforms revisited
• CMB Polarization Experiments
• Fast and Exact Spin-s Spherical Harmonic Transforms
• Libpsht - algorithms for efficient spherical harmonic transforms
• Methods for pixel domain correction of EB leakage
• Sparse image reconstruction on the sphere: implications of a new sampling theorem
• Blind correction of the EB-leakage in the pixel domain
• General solutions of the leakage in integral transforms and applications to the EB-leakage and detection of the cosmological gravitational wave background