Fast algorithms for Jacobi expansions via nonoscillatory phase functions
Numerical Analysis
2018-03-13 v1
Abstract
We describe a suite of fast algorithms for evaluating Jacobi polynomials, applying the corresponding discrete Sturm-Liouville eigentransforms and calculating Gauss-Jacobi quadrature rules. Our approach is based on the well-known fact that Jacobi's differential equation admits a nonoscillatory phase function which can be loosely approximated via an affine function over much of its domain. Our algorithms perform better than currently available methods in most respects. We illustrate this with several numerical experiments, the source code for which is publicly available.
Cite
@article{arxiv.1803.03889,
title = {Fast algorithms for Jacobi expansions via nonoscillatory phase functions},
author = {James Bremer and Haizhao Yang},
journal= {arXiv preprint arXiv:1803.03889},
year = {2018}
}