Fast expansion into harmonics on the ball

Abstract

We devise fast and provably accurate algorithms to transform between an $N\times N\times N$ Cartesian voxel representation of a three-dimensional function and its expansion into the {ball harmonics}, that is, the eigenbasis of the Dirichlet Laplacian on the unit ball in ${\mathbb{R}}^3$. Given $\epsilon >0$, our algorithms achieve relative ${\ell }^1$ - ${\ell }^{\infty }$ accuracy $\epsilon$ in time $O(N^3(\log N{)}^2+N^3|\log \epsilon {|}^2)$, while the naïve direct application of the expansion operators has time complexity $O(N^6)$. We illustrate our methods on numerical examples.

arXiv comment31 pages, 4 figures, 2 tables

References in corpus (9)

• Exact Wavelets on the Ball
• Method of moments for 3-D single particle ab initio modeling with non-uniform distribution of viewing angles
• 3DEX: a code for fast spherical Fourier-Bessel decomposition of 3D surveys
• 3D Spherical Analysis of Baryon Acoustic Oscillations
• A spectral scheme for Kohn-Sham density functional theory of clusters
• Spherical 3D Isotropic Wavelets
• CMB in a box: causal structure and the Fourier-Bessel expansion
• Autocorrelation analysis for cryo-EM with sparsity constraints: Improved sample complexity and projection-based algorithms
• Computing with functions in the ball