The gamma function via interpolation

DOI10.1007/S11075-021-01204-8zbMATH Open1491.65021arXiv2104.00697OpenAlexW3204357867MaRDI QIDQ2138402FDOQ2138402

Authors: Matthew F. Causley 

Publication date: 11 May 2022

Published in: Numerical Algorithms (Search for Journal in Brave)

Abstract: A new computational framework for evaluation of the gamma function $Gamma(z)$ 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.

Full work available at URL: https://arxiv.org/abs/2104.00697

Recommendations

• On the evaluation of the gamma function
• Computation of the Gamma, Digamma, and Trigamma Functions
• An algorithm for the evaluation of the gamma function and ramifications. I
• An algorithm for the evaluation of the gamma function and ramifications. I
• Computing the Gamma Function Using Contour Integrals and Rational Approximations

zbMATH Keywords

gamma functioninterpolationLanczos approximationStirling seriesSpouge approximation

Mathematics Subject Classification ID

Computation of special functions and constants, construction of tables (65D20) Gamma, beta and polygamma functions (33B15)

Cites Work

• Algorithm 814
• The RKFIT Algorithm for Nonlinear Rational Approximation
• The AAA Algorithm for Rational Approximation
• Title not available (Why is that?)
• Title not available (Why is that?)
• Numerical recipes. The art of scientific computing.
• The special functions and their approximations. Vol. I, II
• The convergence of Padé approximants to functions with branch points
• Best estimates of the generalized Stirling formula
• Padé approximants for functions with branch points -- strong asymptotics of Nuttall-Stahl polynomials
• Title not available (Why is that?)
• Talbot quadratures and rational approximations
• Title not available (Why is that?)
• Title not available (Why is that?)
• Leonhard Euler's Integral: A Historical Profile of the Gamma Function: In Memoriam: Milton Abramowitz
• Computation of the Gamma, Digamma, and Trigamma Functions
• A Precision Approximation of the Gamma Function
• An approximation of the gamma function
• Functions of a Complex Variable: Theory and Technique
• On the coefficients of the asymptotic expansion of \(n!\)
• A Berry-Esseen type bound for the kernel density estimator based on a weakly dependent and randomly left truncated data
• Computing the Gamma Function Using Contour Integrals and Rational Approximations
• Gamma function asymptotics by an extension of the method of steepest descents
• An Algorithm for Real and Complex Rational Minimax Approximation
• Unified approaches to the approximations of the gamma function
• A bootstrap method for sum-of-poles approximations
• The early history of the factorial function
• Wielandt's Theorem About the Γ-Function
• Computing the Principal Branch of log-Gamma
• Error bounds and exponential improvements for the asymptotic expansions of the gamma function and its reciprocal
• Computation of matrix gamma function
• Gamma and Factorial in the Monthly

Cited In (3)

• Computation of the Gamma, Digamma, and Trigamma Functions
• On the evaluation of the gamma function
• Title not available (Why is that?)

Uses Software

• Algorithm 814
• RKToolbox
• AAA

This page was built for publication: The gamma function via interpolation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2138402)