Abstract:We propose fast, exact and efficient algorithms for the convolution of two arbitrary functions on the sphere which speed up computations by a factor \order{\sqrt{N}} compared to present methods where $N$ is the number of pixels. No simplifying assumptions are made other than bandlimitation. This reduces typical computation times for convolving the full sky with the asymmetric beam pattern of a megapixel Cosmic Microwave Background (CMB) mission from months to minutes. Our methods enable realistic simulation and careful analysis of data from such missions, taking into account the effects of asymmetric ``point spread functions'' and far side lobes of the physical beam. While motivated by CMB studies, our methods are general and hence applicable to the convolution or filtering of any scalar field on the sphere with an arbitrary, asymmetric kernel. We show in an appendix that the same ideas can be applied to the inverse problems of map-making and beam reconstruction by similarly accelerating the {\em transpose convolution} which is needed for the iterative solution of the normal equations.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to significantly reduce the computational cost of convolving two arbitrary functions on the sphere, especially for the processing of cosmic microwave background (CMB) data. Specifically: 1. **Limitations of Existing Methods**: Current methods are extremely computationally intensive when dealing with physical beams with asymmetric point - spread functions and far - side lobes, resulting in a typical all - sky convolution that may take months. 2. **Proposed Solution**: The authors propose a fast, accurate, and efficient algorithm that can accelerate the convolution calculation of any two band - limited functions on the sphere without simplifying assumptions. The new method rewrites the problem into a form that can utilize the Cooley - Tukey fast Fourier transform (FFT) algorithm, thereby reducing the computation time from months to minutes. 3. **Application Background**: Although this method was initially developed for CMB research, it is applicable to any convolution or filtering problem of scalar fields on the sphere with arbitrary asymmetric kernels. In addition, this method can also be applied to inverse problems in map - making and beam reconstruction. ### Specific Problem Description The key problems mentioned in the paper are: - **Simulation and Correction of Systematic Errors**: It is necessary to simulate and correct systematic errors caused by the imperfections of actual microwave telescopes (such as beam asymmetry and far - side lobes). - **Computational Complexity**: Traditional pixel - based methods face huge computational challenges when dealing with angular scales less than one degree. Even with optimized adaptive grid algorithms, it may still require hundreds of hours of computation time. ### Core Idea of the Solution The authors achieve computational acceleration through the following steps: 1. **Decomposition of Rotation Operations**: Decompose the rotation operation into two auxiliary rotations, thereby simplifying the calculation of the integral. 2. **Spherical Harmonic Basis Representation**: Represent functions and rotation operators in the spherical harmonic basis and use the properties of spherical harmonic functions for efficient computation. 3. **Fast Fourier Transform (FFT)**: Utilize the FFT algorithm to accelerate multi - dimensional Fourier transforms, thereby significantly reducing the computational complexity. ### Improvement in Computational Complexity The computational complexity of traditional methods is $O(L^5)$, while the new method can reduce the computational complexity to $O(L^4\sin\theta_E^2\theta / \pi)$ or lower, depending on the application scenario. In summary, this paper aims to enable scientists to process CMB data more quickly and simulate and analyze complex beam effects more accurately by proposing a new efficient algorithm.

Fast Convolution on the Sphere

Fast and High-Resolution Acoustic Beamforming: A Convolution Accelerated Deconvolution Implementation

Fast computation of radar cross-section by fast multipole method in conjunction with lifting wavelet-like transform

Frequency-independent Scattering for the Large Sphere

Improved cosmic microwave background (de-)lensing using general spherical harmonic transforms

Fast Theoretical Predictions for Spherical Fourier Analysis of Large-Scale Structures

The spherical Fast Multipole Method (sFMM) for Gravitational Lensing Simulation

Efficient polar convolution based on the discrete Fourier-Bessel transform for application in computational biophotonics

Optimal and fast E/B separation with a dual messenger field

Multichannel Sparse Blind Deconvolution on the Sphere

Fast Algorithms and Efficient GPU Implementations for the Radon Transform and the Back-Projection Operator Represented as Convolution Operators

Fast Fourier transforms for the evaluation of convolution products: CPU versus GPU implementation

Differentiable and accelerated spherical harmonic and Wigner transforms

Fast multiscale Gaussian beam methods for wave equations in bounded convex domains

Gabor Transforms on the Sphere with Applications to CMB Power Spectrum Estimation

Optimal filters on the sphere

A novel sampling theorem on the sphere

The Fast Fourier Transform Telescope

Barycentric interpolation formulas for the sphere and the disk