Date of Final Oral Examination (Defense)
Type of Culminating Activity
Master of Science in Mathematics
Grady Wright, Ph.D.
Alex Townsend, Ph.D.
Jodi L. Mead, Ph.D.
A new low rank approximation method for computing with functions in polar and spherical geometries is developed. By synthesizing a classic procedure known as the double Fourier sphere (DFS) method with a structure-preserving variant of Gaussian elimination, approximants to functions on the sphere and disk can be constructed that (1) preserve the bi-periodicity of the sphere, (2) are smooth over the poles of the sphere (and origin of the disk), (3) allow for the use of FFT-based algorithms, and (4) are near-optimal in their underlying discretizations. This method is used to develop a suite of fast, scalable algorithms that exploit the low rank form of approximants to reduce many operations to essentially 1D procedures. This includes algorithms for differentiation, integration, and vector calculus. Combining these ideas with Fourier and ultraspherical spectral methods results in an optimal complexity solver for Poisson's equation, which can be used to solve problems with 108 degrees of freedom in just under a minute on a laptop computer. All of these algorithms have been implemented and are publicly available in the open-source computing system called Chebfun .
Wilber, Heather Denise, "Numerical Computing with Functions on the Sphere and Disk" (2016). Boise State University Theses and Dissertations. 1158.