Publication Date

8-2016

Date of Final Oral Examination (Defense)

5-18-2016

Type of Culminating Activity

Thesis

Degree Title

Master of Science in Mathematics

Department

Mathematics

Major Advisor

Grady Wright, Ph.D.

Major Advisor

Alex Townsend, Ph.D.

Advisor

Jodi L. Mead, Ph.D.

Abstract

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 di erentiation, 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 [21].

Share

COinS