The Gamma Function via Interpolation
Numerical Analysis
2021-09-28 v2 Numerical Analysis
Abstract
A new computational framework for evaluation of the gamma function over the complex plane is developed. The algorithm is based on interpolation by rational functions, and generalizes the classical methods of Lanczos \cite{Lanczos} and Spouge \cite{Spouge} (which we show are also interpolatory). This framework utilizes the exact poles of the gamma function. By relaxing this condition and allowing the poles to vary, a near-optimal rational approximation is possible, which is demonstrated using the adaptive Antoulous Anderson (AAA) algorithm, developed in \cite{AAA,AAA_2020}. The resulting approximations are competitive with Stirling's formula in terms of overall efficiency.
Cite
@article{arxiv.2104.00697,
title = {The Gamma Function via Interpolation},
author = {Matthew F Causley},
journal= {arXiv preprint arXiv:2104.00697},
year = {2021}
}
Comments
19 pages