{"id":"arxiv/2607.15310","url":"https://arxiv.org/abs/2607.15310","tags":["math.GM"],"score":0,"title":"A computational algorithm for the Hardy function $Z(t)$, utilising sub-sequences of generalised cubic Gauss sums, with an overall operational complexity of $O\\bigl((t/\\varepsilon_t)^{[.25,.3]}(\\log t)^{2+o(1)}\\bigr)$, for $t \\in [10^{23},10^{35}]$","domain":"math.GM","source":"arxiv","authors":["David Lewis","Ashley Brereton"],"reviews":0,"subject":"math.GM","summary":"In 2011 G. A. Hiary devised a computational algorithm for the Hardy function $Z(t)$, requiring just $O(t^{1/3} (\\log(t))^κ)$ operations. This compares to $O(\\sqrt t)$ operations necessary for computing $Z(t)$ using the classical Riemann-Siegel formula. The methodology involved the sub-division of the Riemann-Siegel formula into sequences of quadratic Gauss/exponential sums of various lengths $N$. Such sums can be computed rapidly, in order $\\log(N)$ operations, using standard recursive schemes. More recently, the principal author developed a similar algorithm with an…","category":"research","comments":0,"disputed":false,"division":"math","mentions":0,"discussed":0,"published":"2026-07-15","reactions":{}}