6,354 research outputs found
Fast algorithms for spherical harmonic expansions, III
We accelerate the computation of spherical harmonic transforms, using what is
known as the butterfly scheme. This provides a convenient alternative to the
approach taken in the second paper from this series on "Fast algorithms for
spherical harmonic expansions." The requisite precomputations become manageable
when organized as a "depth-first traversal" of the program's control-flow
graph, rather than as the perhaps more natural "breadth-first traversal" that
processes one-by-one each level of the multilevel procedure. We illustrate the
results via several numerical examples.Comment: 14 pages, 1 figure, 6 table
Fourier Based Fast Multipole Method for the Helmholtz Equation
The fast multipole method (FMM) has had great success in reducing the
computational complexity of solving the boundary integral form of the Helmholtz
equation. We present a formulation of the Helmholtz FMM that uses Fourier basis
functions rather than spherical harmonics. By modifying the transfer function
in the precomputation stage of the FMM, time-critical stages of the algorithm
are accelerated by causing the interpolation operators to become
straightforward applications of fast Fourier transforms, retaining the
diagonality of the transfer function, and providing a simplified error
analysis. Using Fourier analysis, constructive algorithms are derived to a
priori determine an integration quadrature for a given error tolerance. Sharp
error bounds are derived and verified numerically. Various optimizations are
considered to reduce the number of quadrature points and reduce the cost of
computing the transfer function.Comment: 24 pages, 13 figure
A Fast and Accurate Algorithm for Spherical Harmonic Analysis on HEALPix Grids with Applications to the Cosmic Microwave Background Radiation
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 pixels
(points), the computational complexity of the method is , with an initial set-up cost of . This compares
favorably with 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
- …