English

Computing Stieltjes constants using complex integration

Classical Analysis and ODEs 2018-08-14 v3 Numerical Analysis

Abstract

The generalized Stieltjes constants γ_n(v)\gamma\_n(v) are, up to a simple scaling factor, the Laurent series coefficients of the Hurwitz zeta function ζ(s,v)\zeta(s,v) about its unique pole s=1s = 1. In this work, we devise an efficient algorithm to compute these constants to arbitrary precision with rigorous error bounds, for the first time achieving this with low complexity with respect to the order~nn. Our computations are based on an integral representation with a hyperbolic kernel that decays exponentially fast. The algorithm consists of locating an approximate steepest descent contour and then evaluating the integral numerically in ball arithmetic using the Petras algorithm with a Taylor expansion for bounds near the saddle point. An implementation is provided in the Arb library. We can, for example, compute γ_n(1)\gamma\_n(1) to 1000 digits in a minute for any nn up to n=10100n=10^{100}. We also provide other interesting integral representations for γ_n(v)\gamma\_n(v), ζ(s)\zeta(s), ζ(s,v)\zeta(s,v), some polygamma functions and the Lerch transcendent.

Keywords

Cite

@article{arxiv.1804.01679,
  title  = {Computing Stieltjes constants using complex integration},
  author = {Fredrik Johansson and Iaroslav Blagouchine},
  journal= {arXiv preprint arXiv:1804.01679},
  year   = {2018}
}