A Fast Transform for Spherical Harmonics
Martin J. Mohlenkamp The Journal of Fourier Analysis and Applications 5(2/3):159–184, 1999.
A liation. Department of Mathematics, Yale University; and Department of Applied Mathemat- ics, University of Colorado, Campus Box 526, Boulder CO 80309-0526, [email protected].
Math Subject Classi cation. Primary 65T20; Secondary 42C10, 33C55.
Keywords and Phrases. Spherical Harmonics, Fast Transforms, Associated Legendre Functions.
Acknowledgements and Notes. I would like to thank my thesis advisor, R.R. Coifman, for his help and guidance.
Abstract
Spherical Harmonics arise on the sphere S2 in the same way that the (Fourier) exponential functions ik e k∈Z arise on the circle. Spherical Harmonic series have many of the same wonderful properties as{ Fourier} series, but have lacked one important thing: a numerically stable fast transform analogous to the Fast Fourier Transform. Without a fast transform, evaluating (or expanding in) Spherical Harmonic series on the computer is slow—for large computations prohibitively slow. This paper provides a fast transform. For a grid of (N 2) points on the sphere, a direct calculation has computational complexity (N 4), but a simpleO separation of variables and Fast Fourier Transform reduce it to (N 3) time. OHere we present algorithms with times (N 5/2 log N) and (N 2(log N)2). O The problem quickly reduces to the fastO application of matricesO of Associated Legendre Functions of certain orders. The essential insight is that although these matrices are dense and oscillatory, locally they can be represented e ciently in trigonometric series.
1 Introduction
Spherical Harmonics arise on the sphere S2 in the same way that the (Fourier) exponential functions ik e k∈Z arise on the circle. Spherical Harmonic series have many of the same wonderful properties as{ Fourier} series, but have lacked one important thing: a numerically stable fast transform analogous to the Fast Fourier Transform. Without a fast transform, evaluating (or expanding in) Spherical Harmonic series on the computer is slow—for large computations prohibitively slow. This paper provides a fast transform. Using spherical coordinates ( (0, ), [0, 2 )), the Spherical Harmonics for S2 can be m im ∈ ∈ m represented as Pn (cos )e where n m , m and n are integers, and the Pn are Associated Legendre functions.{ When prop} erly normalized, | | this forms a basis for the sphere with measure sin d d .
1 1 INTRODUCTION 2
On a grid of 2N points equispaced in by N gaussian nodes in cos , N 2 basis functions are resolved. If we expand (or evaluate) directly using these basis functions, it will take 2N 4 compu- tations. In the variable, our basis is simply Fourier series and we can use the FFT (Fast Fourier Transform). Applying the FFT for N values of requires N 2N log(2N) computations, and leaves us with a set of 2N problems indexed by m. Each problem is: on an interval of N points, expand m (evaluate) a given function using the basis Pn (cos ) with measure sin d . If we do each of these directly, it will take 2N (N 2) computations.{ Our }overall cost is then (N 3). O m O 2 Here we present a method for expanding (evaluating) using Pn (cos ) in (N log N) time, with constant independent of m. This method allows us to reduce{ the overall} costOto (N 2 log2 N). We also present an (N 5/2 log N) algorithm that, although asymptotically slower, performedO better in our tests. Both methoO ds use precomputed, compressed representations of the Associated Legendre function bases. This paper is based on the thesis [Moh97]. More details may be found there.
1.1 Statement of the Problem We are given a set of Spherical Harmonics:
˜m im Pn (cos )e /√2 (1) (m,n)∈IN n o ˜m de ned on the sphere ( (0, ), [0, 2 )). Pn is the Associated Legendre function of degree n and order m, normalized∈ so that our∈ set of functions is orthonormal on the sphere. The indices allowed are restricted to the set
= (m, n) Z Z : 0 n < N, n m n . (2) IN { ∈ } m For the evaluation problem, we are also given a set of coe cients n (m,n)∈IN and we form the sum { } m ˜m im f( , ) = n Pn (cos )e /√2 . (3) XIN This sum is evaluated on a grid in with 2N equispaced points in and N gaussian nodes in cos : k ( , ) = cos 1(gN ), 2 : j, k Z; 0 j < N, 0 k < 2N . (4) j k j 2N ∈ ½ µ ¶ ¾ (See e.g. [FPTV92] for discussion of the gaussian nodes.) The evaluation problem for Spherical 2 2 m Harmonics is to compute the (N ) output values f( j , k) from the (N ) input values n in a fast, numerically stable way. O O The expansion problem for Spherical Harmonics is dual to the evaluation problem. We are given 2 2 the (N ) sampled values f( j , k) of a function of the form (3) and wish to compute the (N ) O m O output values n in a fast, numerically stable way.
1.2 Summary of Results The problem of the Fast Transform quickly reduces to a set of 2N problems indexed by m ( N < m < N). (See Section 2.3.) The problems are to nd a fast application for the matrices
m ˜m Mn,j = sin jPn (cos j) . (5) n,j ³p ´ These matrices are dense and oscillatory. (See Figure 2 in Section 2.4.) 2 BACKGROUND 3
To understand these matrices, we model them using quasi-classical (WKB) frequency estimates. (Section 2.4.) The model is