No label defined (Q4550206)

The Riemann zeta-function \(\zeta(\sigma+it)\) may be computed to moderate accuracy, requiring \(O(t^{1/2+\varepsilon})\) arithmetic operations, by use of the Riemann-Siegel formula. The Euler-Maclaurin summation formula may be used for computing \(\zeta(\sigma+it)\) to arbitrarily high accuracy, using \(O(t^{1+\varepsilon}+d^{1+\varepsilon})\) operations to compute \(d\) digits. In the present paper it is shown that numerical quadrature may be used to compute \(\zeta(\sigma+it)\) to arbitrarily high accuracy. This method requires \(O(t^{1/2+\varepsilon})\) operations to calculate \(d\) digits of \(\zeta(\sigma+it)\), for fixed \(d\) as \(t\to\infty\). For large \(t\), this method should improve on the Euler-Maclaurin method for sufficiently large \(t\) with efficiency comparable to that of the Riemann-Siegel formula. For the entire collection see [Zbl 0980.00025].