An efficient algorithm for accelerating the convergence of oscillatory series, useful for computing the polylogarithm and Hurwitz zeta functions

Abstract

This paper sketches a technique for improving the rate of convergence of a general oscillatory sequence, and then applies this series acceleration algorithm to the polylogarithm and the Hurwitz zeta function. As such, it may be taken as an extension of the techniques given by Borwein's "An efficient algorithm for computing the Riemann zeta function", to more general series. The algorithm provides a rapid means of evaluating Li_s(z) for general values of complex s and the region of complex z values given by |z^2/(z-1)|<4. Alternatively, the Hurwitz zeta can be very rapidly evaluated by means of an Euler-Maclaurin series. The polylogarithm and the Hurwitz zeta are related, in that two evaluations of the one can be used to obtain a value of the other; thus, either algorithm can be used to evaluate either function. The Euler-Maclaurin series is a clear performance winner for the Hurwitz zeta, while the Borwein algorithm is superior for evaluating the polylogarithm in the kidney-shaped region. Both algorithms are superior to the simple Taylor's series or direct summation. The primary, concrete result of this paper is an algorithm allows the exploration of the Hurwitz zeta in the critical strip, where fast algorithms are otherwise unavailable. A discussion of the monodromy group of the polylogarithm is included.

37 pages, 6 graphs, 14 full-color phase plots. v3: Added discussion of a fast Hurwitz algorithm; expanded development of the monodromy v4:Correction and clarifiction of monodromy

References in corpus (2)

• Convergence Acceleration via Combined Nonlinear-Condensation Transformations
• Resurgence of the fractional polylogarithms

Cited by in corpus (12)

• Borel and Stokes Nonperturbative Phenomena in Topological String Theory and c=1 Matrix Models
• Systematic time expansion for the Kardar-Parisi-Zhang equation, linear statistics of the GUE at the edge and trapped fermions
• The distinct flavors of Zipf's law in the rank-size and in the size-distribution representations, and its maximum-likelihood fitting
• Kuramoto model with uniformly spaced frequencies:Finite-N asymptotics of the locking threshold
• A practical recipe to fit discrete power-law distributions
• The Lerch zeta function III. Polylogarithms and special values
• Quantisation of the effective string with TBA
• Table of Dirichlet L-Series and Prime Zeta Modulo Functions for Small Moduli
• The Polylogarithm Function in Julia
• The Riemann Hypothesis: A Qualitative Characterization of the Nontrivial Zeros of the Riemann Zeta Function Using Polylogarithms
• A new non-negative distribution with both finite and infinite support
• Fourier series for eclipses on exoplanet binaries