Computing the Barnes $G$-function and the gamma function in the entire complex plane

TL;DR

This paper introduces an efficient algorithm for approximating the Barnes G-function and gamma function across the entire complex plane using elementary functions and Padé approximations, achieving high accuracy.

Contribution

The authors develop a novel two-point Padé approximation-based algorithm for computing Barnes G-function and gamma function in the entire complex plane.

Findings

Achieves approximation accuracy of 3×10^{-16} and 3×10^{-31} in the half-plane Re(z)≥3/2.

Provides a reflection formula to extend computations to the entire complex plane.

Produces accurate gamma function approximations as a by-product.

Abstract

We present an algorithm for generating approximations for the logarithm of Barnes $G$-function in the half-plane $Re(z)\ge 3/2$. These approximations involve only elementary functions and are easy to implement. The algorithm is based on a two-point Pad\'e approximation and we use it to provide two approximations to $\ln(G(z))$, accurate to $3 \times 10^{-16}$ and $3 \times 10^{-31}$ in the half-plane $Re(z)\ge 3/2$; a reflection formula is then used to compute Barnes $G$-function in the entire complex plane. A by-product of our algorithm is that it also produces accurate approximations to the gamma function.