Computing the truncated theta function via Mordell integral

Abstract

Hiary [3] has presented an algorithm which allows to evaluate the truncated theta function $\sum_{k=0}^n \exp(2πı(zk+τk^2))$ to within $\pm ε$ in $O(\ln(\tfrac{n}ε)^κ)$ arithmetic operations for any real $z$ and $τ$. This remarkable result has many applications in Number Theory, in particular it is the crucial element in Hiary's algorithm for computing $ζ(\tfrac{1}{2}+ıt)$ to within $\pm t^{-λ}$ in $O_λ(t^{\frac{1}{3}}\ln(t)^κ)$ arithmetic operations, see [2]. We present a significant simplification of Hiary's algorithm for evaluating the truncated theta function. Our method avoids the use of the Poisson summation formula, and substitutes it with an explicit identity involving the Mordell integral. This results in an algorithm which is efficient, conceptually simple and easy to implement.