A fast and accurate algorithm for spherical harmonic analysis on HEALPix grids with applications to the cosmic microwave background radiation

DOI10.1016/J.JCP.2020.109544MaRDI QIDQ781995zbMATH OpenOpenAlexFDO

Authors Kathryn P. Drake, Grady B. Wright

Publication date 21 July 2020

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1904.10514

zbMATH Keywords

nonuniform fast Fourier transformcosmic microwave background radiationdouble Fourier spherefast spherical harmonic transformHEALPix

Mathematics Subject Classification ID

Astrophysical cosmology (85A40) Numerical methods for discrete and fast Fourier transforms (65T50) Applications to the sciences (65Z05)

Abstract: The Hierarchical Equal Area isoLatitude Pixelation (HEALPix) scheme is used extensively in astrophysics for data collection and analysis on the sphere. The scheme was originally designed for studying the Cosmic Microwave Background (CMB) radiation, which represents the first light to travel during the early stages of the universe's development and gives the strongest evidence for the Big Bang theory to date. Refined analysis of the CMB angular power spectrum can lead to revolutionary developments in understanding the nature of dark matter and dark energy. In this paper, we present a new method for performing spherical harmonic analysis for HEALPix data, which is a central component to computing and analyzing the angular power spectrum of the massive CMB data sets. The method uses a novel combination of a non-uniform fast Fourier transform, the double Fourier sphere method, and Slevinsky's fast spherical harmonic transform (Slevinsky, 2019). For a HEALPix grid with $N$ pixels (points), the computational complexity of the method is $mathcalO(Nlo{g}^{2}N)$, with an initial set-up cost of $mathcalO(N^{3/2}logN)$. This compares favorably with $mathcalO(N^{3/2})$ runtime complexity of the current methods available in the HEALPix software when multiple maps need to be analyzed at the same time. Using numerical experiments, we demonstrate that the new method also appears to provide better accuracy over the entire angular power spectrum of synthetic data when compared to the current methods, with a convergence rate at least two times higher.

Recommendations

• Fast spin \(\pm \)2 spherical harmonics transforms and application in cosmology
• Non-linear data analysis on the sphere. The quest for anomalies in the cosmic microwave background
• GAUSS–LEGENDRE SKY PIXELIZATION (GLESP) FOR CMB MAPS
• scientific article; zbMATH DE number 1803183
• The Gauss-Legendre sky pixelization for the CMB polarization (GLESP-pol) errors due to pixelization of the CMB sky

Cites work

• scientific article; zbMATH DE number 1688612 (Why is no real title available?)
• scientific article; zbMATH DE number 1064576 (Why is no real title available?)
• scientific article; zbMATH DE number 961607 (Why is no real title available?)
• A comparison of popular point configurations on S^2
• A fast algorithm for particle simulations
• A nonuniform fast Fourier transform based on low rank approximation
• A partition of the unit sphere into regions of equal area and small diameter
• Accelerating the Nonuniform Fast Fourier Transform
• An Algorithm for the Machine Calculation of Complex Fourier Series
• An algorithm for the rapid evaluation of special function transforms
• Computing with functions in spherical and polar geometries. I. The sphere
• Distributing many points on a sphere
• Fast and backward stable transforms between spherical harmonic expansions and bivariate Fourier series
• Fast evaluation of quadrature formulae on the sphere
• LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares
• Meshfree approximation methods with Matlab. With CD-ROM.
• Numerical solution of saddle point problems
• Spherical harmonics and approximations on the unit sphere. An introduction

Cited in

(3)

• Approximation properties of the double Fourier sphere method
• A double Fourier sphere method for \(d\)-dimensional manifolds
• Fast spin \(\pm \)2 spherical harmonics transforms and application in cosmology

Describes a project that uses

Uses Software

• CRAIG
• FastTransforms.jl
• IGLOO
• Healpix
• LSQR

This page was built for publication: A fast and accurate algorithm for spherical harmonic analysis on HEALPix grids with applications to the cosmic microwave background radiation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q781995)