Computing Stieltjes constants using complex integration

Abstract : The generalized Stieltjes constants $\gamma_n(v)$ are, up to a simple scaling factor, the Laurent series coefficients of the Hurwitz zeta function $\zeta(s,v)$ about its unique pole $s = 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~$n$. 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 $\gamma_n(1)$ to 1000 digits in a minute for any $n$ up to $n=10^{100}$. We also provide other interesting integral representations for $\gamma_n(v)$, $\zeta(s)$, $\zeta(s,v)$, some polygamma functions and the Lerch transcendent.
Type de document :
Pré-publication, Document de travail
2018
Liste complète des métadonnées

Littérature citée [29 références]  Voir  Masquer  Télécharger

https://hal.inria.fr/hal-01758620
Contributeur : Fredrik Johansson <>
Soumis le : mercredi 30 mai 2018 - 17:14:10
Dernière modification le : vendredi 8 juin 2018 - 13:40:35

Fichiers

stieltjes.pdf
Fichiers produits par l'(les) auteur(s)

Identifiants

  • HAL Id : hal-01758620, version 2
  • ARXIV : 1804.01679

Collections

Citation

Fredrik Johansson, Iaroslav Blagouchine. Computing Stieltjes constants using complex integration. 2018. 〈hal-01758620v2〉

Partager

Métriques

Consultations de la notice

48

Téléchargements de fichiers

28