Abstract
We devise fast and provably accurate algorithms to transform between an N × N × 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 ℝ3. Given ℰ > 0, our algorithms achieve relative ℓ1-ℓ∞ accuracy ℰ in time O(N3(log N)2 + N3|log ℰ|2), while the naive direct application of the expansion operators has time complexity O(N6). We illustrate our methods on numerical examples.
| Original language | English (US) |
|---|---|
| Pages (from-to) | A1117-A1144 |
| Journal | SIAM Journal on Scientific Computing |
| Volume | 47 |
| Issue number | 2 |
| DOIs | |
| State | Published - 2025 |
All Science Journal Classification (ASJC) codes
- Computational Mathematics
- Applied Mathematics
Keywords
- Laplacian eigenfunctions
- fast transforms
- spherical Bessel
- spherical harmonics