Odlyzko-Sch Onhage Algorithm in Conductor Aspect: Bstract
Odlyzko-Sch Onhage Algorithm in Conductor Aspect: Bstract
Odlyzko-Sch Onhage Algorithm in Conductor Aspect: Bstract
IN CONDUCTOR ASPECT
arXiv:math/0306101v2 [math.NT] 26 Jun 2003
JEFFREY STOPPLE
1. I NTRODUCTION .
In [6], Odlyzko and Schönhage developed an algorithm to com-
pute the Riemann zeta function ζ(s) efficiently for values of s very
high up in the critical strip. Their method depends on precompu-
tation of Taylor series expansions of ζ(s) at regularly spaced points,
which in turn can be done efficiently by a clever application of the
Fast Fourier Transform. Rumely later implemented a version of this
for Dirichlet L-functions in [10].
In analytic number theory it is often the case that there is a sym-
metry between what one can prove for large values of t = Im(s) in
the critical strip, and what one can prove for L-functions with large
conductor q. With this in mind, it seems reasonable to try to find an
analog of the Odlyzko-Schönhage algorithm to efficiently compute
values of L-functions with very large conductor.
Roughly speaking, we want to view the L-function as a Mellin
transform of an automorphic form f , and split the integral at the
symmetry point. This gives the extended L-function as an infinite
sum
X
Λ(s, f ) = a(n){G(s, 2πn/q 1/2 ) + G(1 − s, 2πn/q 1/2 )}
n
of values of an incomplete Gamma function
Z ∞
dy
−s
G(s, x) = x Γ(s, x) = exp(−yx)y s .
1 y
1991 Mathematics Subject Classification. 11Y16; 11Y35.
1
2 JEFFREY STOPPLE
def.
X X X
Θ(z, χ) = χ(Q) rQ (n) exp(2πinz) = rχ (n) exp(2πinz).
[Q] n>0 n>0
ODLYZKO ALGORITHM IN CONDUCTOR ASPECT 3
Then
def.
Λ(s, χ) = (q 1/2 /2π)s Γ(s)L(s, χ) =
X
rχ (n) G(s, 2πn/q 1/2 ) + G(1 − s, 2πn/q 1/2 ) .
n>0
We will truncate the infinite series after N terms, and arrange the n <
N into intervals Ij , j = 1, 2, . . . T centered at points xj and with width
∆j . For n in the interval Ij we compute a Taylor series expansion,
truncated to B terms
XB
1/2 1/2
G(s, 2πn/q ) + G(1 − s, 2πn/q ) ≈ Gj,k (s)(2πn/q 1/2 − xj )k ,
k=0
where
Gj,k (s) = G(k) (s, xj ) + G(k) (1 − s, xj ) /k!.
(1)
This gives Z(t, χ) as
Here we have inserted canceling terms ∆kj , ∆−k j to control the small
size of Gj,k (s) and the large size of (2πn/q − xj )k .
1/2
Since the rχ (n) are the coefficients of a weight one cusp form, rχ (n) ≪
n1/2 . Using (5), the left side of (7) is
X
≪ (q/n)1/2 exp(−2πn/q 1/2 )
n>N
Z ∞
1/2
≈q y −1/2 exp(−2πy/q 1/2 )dy
N
Z ∞
1/2
<(q/N) exp(−2πy/q 1/2 )dy
N
1/2
≪q/N exp(−2πN/q 1/2 ).
For this to be less than 10−D we need
log(q) + D log(10) < 2πN/q 1/2 + log(N)/2;
it suffices that
(8) N = q 1/2 log(q · 10D )/2π.
and
j−1
e+1 x1
(11) ∆j = .
e−1 2e
The simplest choice for x1 is 2π/q 1/2 so that the first interval I1
contains the n = 1 term of the series, and this was implemented in
practice. However, because of the exponential growth, for small j
some of the intervals Ij will contain only a single term 2πn/q 1/2 , and
some will contain none at all. Thus for the very small values of j
where the intervals contain fewer than B terms one can do better
computing each term instead of using the Taylor expansion. This
would improve the efficiency about 40% for q near 106 , 30% near 109,
and 20% near 1020 .
We can now determine how many intervals T are needed, since
we need to compute terms in the series out to
With
1 + log(e/(π(e + 1)))/ log((e + 1)/(e − 1)) < 0,
1/ log((e + 1)/(e − 1)) < 1.3
we can simply choose
We will assume the latter from now on. We want ǫ · q 1/4 < 10−D , or
(13) ǫ < q −1/4 10−D ,
which will determine how many terms B we need to take in each
Taylor expansion.
It follows from the choice made in §3 that
∆j 1
= and xj − Rj = x1 /2.
Rj e
Thus the tail (9) of the Taylor series reduces to
2 exp(−x1 /2) −B
(14) ǫ= e (1 − 1/e)−1 .
x1
Again, with x1 = 2π/q 1/2 we combine (13) and (14) to find that we
require
q 1/2 exp(−π/q 1/2 )e−B (e − 1)/(πe) < q −1/4 10−D .
Since
exp(−π/q 1/2 )(e − 1)/(πe) < 1
we ignore it, and so we need q 3/4 10D < eB , and let
(15) B = log(q 3/4 10D ).
Even under the assumption of maximal error we would still only
need B = log(q10D ).
ODLYZKO ALGORITHM IN CONDUCTOR ASPECT 9
takes N ·B operations, which is O(q 1/2 log(q)2 ) by (8) and (15). Subse-
quently, evaluations of Z(t, χ) using (2) require only T ·B operations,
which is O(log(q)2 ) by (12) and (15).
This algorithm was implemented in Mathematica and run on a 400
MHz. Apple Powerbook G4 under OS X 10.2. Mathematica is a good
choice if one wants to re-invent the wheel as few times as possible.
The incomplete Gamma function Γ(s, x) is available; it is computed
via hypergeometric functions and continued fractions according to
10 JEFFREY STOPPLE
0.8
0.6
0.4
0.2
1.2
0.8
0.6
0.4
0.2
0.5 1 1.5
1.2
0.8
0.6
0.4
0.2
√
F IGURE 3. Mean of first zero vs. h(−q)/ q.
1 − sin(2πx)/(2πx),
ODLYZKO ALGORITHM IN CONDUCTOR ASPECT 13
where
[α, β) = [0, .05), [.05, .1), . . . , [1.75, 1.80)
(Here we denote the L-function L(s, φ), so that χ[α,β) can denote a
characteristic function, as usual.) The zeros γ̃ have been re-renormalized
so that the mean above 0. of the first zero is .78.
Finally, for each of the 29 discriminants considered, the mean above
0. of the first zero for the L-functions with that discriminant was
√
compared to the relative size h(−q)/ q of the class number. Fig-
ure 3 shows the plot. Surprisingly, the discriminants with smaller
class number tend to have a larger mean first zero. The correlation
between the two quantities is −.89.
R EFERENCES
[1] G. Dahlquist, Numerical Methods, Prentice Hall, 1974.
[2] N. Katz, and P. Sarnak, Zeros of zeta functions and symmetry, Bull. Amer. Math.
Soc., 36, no. 1, 1999, pp. 1-26.
[3] J. Lagarias, and A. Odlyzko, On computing Artin L-functions in the critical strip,
Math. Comp. 33, 147, 1979, pp.1081-1095.
[4] S. Wolfram, The Mathematica Book, 4th ed., 1999.
[5] H. Montgomery, and P. Weinberger, Notes on small class numbers, Acta Arith.
XXIV 1974, pp. 529-542.
[6] A. Odlyzko, and A. Schönhage, Fast algorithms for multiple evaluations of the
Riemann zeta function, Trans. Amer. Math. Soc. 309, 2, 1988 pp. 797-809.
[7] PARI-GP. By C. Batut, D. Bernardi, H. Cohen , and M. Olivier, currently main-
tained by K. Belabas. http://www.parigp-home.de
[8] PARI-GP mailing list, available at
http://www.parigp-home.de/lists/200212aav
[9] M. Rubinstein, Evidence for a spectral interpretation of the zeros of L-functions,
Ph.D. thesis, Princeton 1998.
[10] R. Rumely, Numerical computations concerning the ERH, Math. Comp., 61, 203,
1993, pp. 415-440.
[11] P. Weinberger, On small zeros of Dirichlet L-functions, Math. Comp., 29 no. 129,
1975, pp. 319-328.
E-mail address: [email protected]