Rubinstein’s \(L\)-function calculator¶

This interface provides complete access to Rubinstein’s lcalc calculator with extra PARI functionality compiled in and is a standard part of Sage.

Note

Each call to lcalc runs a complete lcalc process. On a typical Linux system, this entails about 0.3 seconds overhead.

AUTHORS:

Michael Rubinstein (2005): released under GPL the C++ program lcalc
William Stein (2006-03-05): wrote Sage interface to lcalc

class sage.lfunctions.lcalc.LCalc[source]¶

Bases: SageObject

Rubinstein’s \(L\)-functions Calculator.

Type lcalc.[tab] for a list of useful commands that are implemented using the command line interface, but return objects that make sense in Sage. For each command the possible inputs for the \(L\)-function are:

" – (default) the Riemann zeta function
'tau' – the L function of the Ramanujan delta function
E – an elliptic curve over \(\QQ\); defines \(L(E,s)\)

You can also use the complete command-line interface of Rubinstein’s \(L\)-functions calculations program via this class. Type lcalc.help() for a list of commands and how to call them.

analytic_rank(L='')[source]¶

Return the analytic rank of the \(L\)-function at the central critical point.

INPUT:

L – defines \(L\)-function (default: Riemann zeta function)

OUTPUT: integer

Note

Of course this is not provably correct in general, since it is an open problem to compute analytic ranks provably correctly in general.

EXAMPLES:

Sage

sage: E = EllipticCurve('37a')
sage: lcalc.analytic_rank(E)
1

Python

>>> from sage.all import *
>>> E = EllipticCurve('37a')
>>> lcalc.analytic_rank(E)
1

Sage Live

E = EllipticCurve('37a')
lcalc.analytic_rank(E)

help()[source]¶

twist_values(s, dmin, dmax, L='')[source]¶

Return values of \(L(s, \chi_k)\) for each quadratic character \(\chi_k\) whose discriminant \(d\) satisfies \(d_{\min} \leq d \leq d_{\max}\).

INPUT:

s – complex numbers
dmin – integer
dmax – integer
L – defines \(L\)-function (default: Riemann zeta function)

OUTPUT: list of pairs \((d, L(s,\chi_d))\)

EXAMPLES:

Sage

sage: values = lcalc.twist_values(0.5, -10, 10)
sage: values[0][0]
-8
sage: values[0][1] # abs tol 1e-8
1.10042141 + 0.0*I
sage: values[1][0]
-7
sage: values[1][1] # abs tol 1e-8
1.14658567 + 0.0*I
sage: values[2][0]
-4
sage: values[2][1] # abs tol 1e-8
0.667691457 + 0.0*I
sage: values[3][0]
-3
sage: values[3][1] # abs tol 1e-8
0.480867558 + 0.0*I
sage: values[4][0]
5
sage: values[4][1] # abs tol 1e-8
0.231750947 + 0.0*I
sage: values[5][0]
8
sage: values[5][1] # abs tol 1e-8
0.373691713 + 0.0*I

Python

>>> from sage.all import *
>>> values = lcalc.twist_values(RealNumber('0.5'), -Integer(10), Integer(10))
>>> values[Integer(0)][Integer(0)]
-8
>>> values[Integer(0)][Integer(1)] # abs tol 1e-8
1.10042141 + 0.0*I
>>> values[Integer(1)][Integer(0)]
-7
>>> values[Integer(1)][Integer(1)] # abs tol 1e-8
1.14658567 + 0.0*I
>>> values[Integer(2)][Integer(0)]
-4
>>> values[Integer(2)][Integer(1)] # abs tol 1e-8
0.667691457 + 0.0*I
>>> values[Integer(3)][Integer(0)]
-3
>>> values[Integer(3)][Integer(1)] # abs tol 1e-8
0.480867558 + 0.0*I
>>> values[Integer(4)][Integer(0)]
5
>>> values[Integer(4)][Integer(1)] # abs tol 1e-8
0.231750947 + 0.0*I
>>> values[Integer(5)][Integer(0)]
8
>>> values[Integer(5)][Integer(1)] # abs tol 1e-8
0.373691713 + 0.0*I

Sage Live

values = lcalc.twist_values(0.5, -10, 10)
values[0][0]
values[0][1] # abs tol 1e-8
values[1][0]
values[1][1] # abs tol 1e-8
values[2][0]
values[2][1] # abs tol 1e-8
values[3][0]
values[3][1] # abs tol 1e-8
values[4][0]
values[4][1] # abs tol 1e-8
values[5][0]
values[5][1] # abs tol 1e-8

twist_zeros(n, dmin, dmax, L='')[source]¶

Return first \(n\) real parts of nontrivial zeros for each quadratic character \(\chi_k\) whose discriminant \(d\) satisfies \(d_{\min} \leq d \leq d_{\max}\).

INPUT:

n – integer
dmin – integer
dmax – integer
L – defines \(L\)-function (default: Riemann zeta function)

OUTPUT: dictionary; keys are the discriminants \(d\), and values are list of corresponding zeros

EXAMPLES:

Sage

sage: lcalc.twist_zeros(3, -3, 6)
{-3: [8.03973716, 11.2492062, 15.7046192], 5: [6.64845335, 9.83144443, 11.9588456]}

Python

>>> from sage.all import *
>>> lcalc.twist_zeros(Integer(3), -Integer(3), Integer(6))
{-3: [8.03973716, 11.2492062, 15.7046192], 5: [6.64845335, 9.83144443, 11.9588456]}

Sage Live

lcalc.twist_zeros(3, -3, 6)

value(s, L='')[source]¶

Return \(L(s)\) for \(s\) a complex number.

INPUT:

s – complex number
L – defines \(L\)-function (default: Riemann zeta function)

EXAMPLES:

Sage

sage: I = CC.0
sage: lcalc.value(0.5 + 100*I)
2.69261989 - 0.0203860296*I

Python

>>> from sage.all import *
>>> I = CC.gen(0)
>>> lcalc.value(RealNumber('0.5') + Integer(100)*I)
2.69261989 - 0.0203860296*I

Sage Live

I = CC.0
lcalc.value(0.5 + 100*I)

Note, Sage can also compute zeta at complex numbers (using the PARI C library):

Sage

sage: (0.5 + 100*I).zeta()
2.69261988568132 - 0.0203860296025982*I

Python

>>> from sage.all import *
>>> (RealNumber('0.5') + Integer(100)*I).zeta()
2.69261988568132 - 0.0203860296025982*I

Sage Live

(0.5 + 100*I).zeta()

values_along_line(s0, s1, number_samples, L='')[source]¶

Return values of \(L(s)\) at number_samples equally-spaced sample points along the line from \(s_0\) to \(s_1\) in the complex plane.

INPUT:

s0, s1 – complex numbers
number_samples – integer
L – defines \(L\)-function (default: Riemann zeta function)

OUTPUT: list of pairs \((s, L(s))\), where the \(s\) are equally spaced sampled points on the line from \(s_0\) to \(s_1\)

EXAMPLES:

Sage

sage: I = CC.0
sage: values = lcalc.values_along_line(0.5, 0.5+20*I, 5)
sage: values[0][0] # abs tol 1e-8
0.5
sage: values[0][1] # abs tol 1e-8
-1.46035451 + 0.0*I
sage: values[1][0] # abs tol 1e-8
0.5 + 4.0*I
sage: values[1][1] # abs tol 1e-8
0.606783764 + 0.0911121400*I
sage: values[2][0] # abs tol 1e-8
0.5 + 8.0*I
sage: values[2][1] # abs tol 1e-8
1.24161511 + 0.360047588*I
sage: values[3][0] # abs tol 1e-8
0.5 + 12.0*I
sage: values[3][1] # abs tol 1e-8
1.01593665 - 0.745112472*I
sage: values[4][0] # abs tol 1e-8
0.5 + 16.0*I
sage: values[4][1] # abs tol 1e-8
0.938545408 + 1.21658782*I

Python

>>> from sage.all import *
>>> I = CC.gen(0)
>>> values = lcalc.values_along_line(RealNumber('0.5'), RealNumber('0.5')+Integer(20)*I, Integer(5))
>>> values[Integer(0)][Integer(0)] # abs tol 1e-8
0.5
>>> values[Integer(0)][Integer(1)] # abs tol 1e-8
-1.46035451 + 0.0*I
>>> values[Integer(1)][Integer(0)] # abs tol 1e-8
0.5 + 4.0*I
>>> values[Integer(1)][Integer(1)] # abs tol 1e-8
0.606783764 + 0.0911121400*I
>>> values[Integer(2)][Integer(0)] # abs tol 1e-8
0.5 + 8.0*I
>>> values[Integer(2)][Integer(1)] # abs tol 1e-8
1.24161511 + 0.360047588*I
>>> values[Integer(3)][Integer(0)] # abs tol 1e-8
0.5 + 12.0*I
>>> values[Integer(3)][Integer(1)] # abs tol 1e-8
1.01593665 - 0.745112472*I
>>> values[Integer(4)][Integer(0)] # abs tol 1e-8
0.5 + 16.0*I
>>> values[Integer(4)][Integer(1)] # abs tol 1e-8
0.938545408 + 1.21658782*I

Sage Live

I = CC.0
values = lcalc.values_along_line(0.5, 0.5+20*I, 5)
values[0][0] # abs tol 1e-8
values[0][1] # abs tol 1e-8
values[1][0] # abs tol 1e-8
values[1][1] # abs tol 1e-8
values[2][0] # abs tol 1e-8
values[2][1] # abs tol 1e-8
values[3][0] # abs tol 1e-8
values[3][1] # abs tol 1e-8
values[4][0] # abs tol 1e-8
values[4][1] # abs tol 1e-8

Sometimes warnings are printed (by lcalc) when this command is run:

Sage

sage: E = EllipticCurve('389a')
sage: values = E.lseries().values_along_line(0.5, 3, 5)
sage: values[0][0] # abs tol 1e-8
0.0
sage: values[0][1] # abs tol 1e-8
0.209951303  + 0.0*I
sage: values[1][0] # abs tol 1e-8
0.5
sage: values[1][1] # abs tol 1e-8
0.0  + 0.0*I
sage: values[2][0] # abs tol 1e-8
1.0
sage: values[2][1] # abs tol 1e-8
0.133768433 - 0.0*I
sage: values[3][0] # abs tol 1e-8
1.5
sage: values[3][1] # abs tol 1e-8
0.360092864 - 0.0*I
sage: values[4][0] # abs tol 1e-8
2.0
sage: values[4][1] # abs tol 1e-8
0.552975867 + 0.0*I

Python

>>> from sage.all import *
>>> E = EllipticCurve('389a')
>>> values = E.lseries().values_along_line(RealNumber('0.5'), Integer(3), Integer(5))
>>> values[Integer(0)][Integer(0)] # abs tol 1e-8
0.0
>>> values[Integer(0)][Integer(1)] # abs tol 1e-8
0.209951303  + 0.0*I
>>> values[Integer(1)][Integer(0)] # abs tol 1e-8
0.5
>>> values[Integer(1)][Integer(1)] # abs tol 1e-8
0.0  + 0.0*I
>>> values[Integer(2)][Integer(0)] # abs tol 1e-8
1.0
>>> values[Integer(2)][Integer(1)] # abs tol 1e-8
0.133768433 - 0.0*I
>>> values[Integer(3)][Integer(0)] # abs tol 1e-8
1.5
>>> values[Integer(3)][Integer(1)] # abs tol 1e-8
0.360092864 - 0.0*I
>>> values[Integer(4)][Integer(0)] # abs tol 1e-8
2.0
>>> values[Integer(4)][Integer(1)] # abs tol 1e-8
0.552975867 + 0.0*I

Sage Live

E = EllipticCurve('389a')
values = E.lseries().values_along_line(0.5, 3, 5)
values[0][0] # abs tol 1e-8
values[0][1] # abs tol 1e-8
values[1][0] # abs tol 1e-8
values[1][1] # abs tol 1e-8
values[2][0] # abs tol 1e-8
values[2][1] # abs tol 1e-8
values[3][0] # abs tol 1e-8
values[3][1] # abs tol 1e-8
values[4][0] # abs tol 1e-8
values[4][1] # abs tol 1e-8

zeros(n, L='')[source]¶

Return the imaginary parts of the first \(n\) nontrivial zeros of the \(L\)-function in the upper half plane, as 32-bit reals.

INPUT:

n – integer
L – defines \(L\)-function (default: Riemann zeta function)

This function also checks the Riemann Hypothesis and makes sure no zeros are missed. This means it looks for several dozen zeros to make sure none have been missed before outputting any zeros at all, so takes longer than self.zeros_of_zeta_in_interval(...).

EXAMPLES:

Sage

sage: lcalc.zeros(4)                           # long time
[14.1347251, 21.0220396, 25.0108576, 30.4248761]
sage: lcalc.zeros(5, L='--tau')                # long time
[9.22237940, 13.9075499, 17.4427770, 19.6565131, 22.3361036]
sage: lcalc.zeros(3, EllipticCurve('37a'))     # long time
[0.000000000, 5.00317001, 6.87039122]

Python

>>> from sage.all import *
>>> lcalc.zeros(Integer(4))                           # long time
[14.1347251, 21.0220396, 25.0108576, 30.4248761]
>>> lcalc.zeros(Integer(5), L='--tau')                # long time
[9.22237940, 13.9075499, 17.4427770, 19.6565131, 22.3361036]
>>> lcalc.zeros(Integer(3), EllipticCurve('37a'))     # long time
[0.000000000, 5.00317001, 6.87039122]

Sage Live

lcalc.zeros(4)                           # long time
lcalc.zeros(5, L='--tau')                # long time
lcalc.zeros(3, EllipticCurve('37a'))     # long time

zeros_in_interval(x, y, stepsize, L='')[source]¶

Return the imaginary parts of (most of) the nontrivial zeros of the \(L\)-function on the line \(\Re(s)=1/2\) with positive imaginary part between \(x\) and \(y\), along with a technical quantity for each.

INPUT:

x, y, stepsize – positive floating point numbers
L – defines \(L\)-function (default: Riemann zeta function)

OUTPUT: list of pairs (zero, S(T)).

Rubinstein writes: The first column outputs the imaginary part of the zero, the second column a quantity related to \(S(T)\) (it increases roughly by 2 whenever a sign change, i.e. pair of zeros, is missed). Higher up the critical strip you should use a smaller stepsize so as not to miss zeros.

EXAMPLES:

Sage

sage: lcalc.zeros_in_interval(10, 30, 0.1)
[(14.1347251, 0.184672916), (21.0220396, -0.0677893290), (25.0108576, -0.0555872781)]

Python

>>> from sage.all import *
>>> lcalc.zeros_in_interval(Integer(10), Integer(30), RealNumber('0.1'))
[(14.1347251, 0.184672916), (21.0220396, -0.0677893290), (25.0108576, -0.0555872781)]

Sage Live

lcalc.zeros_in_interval(10, 30, 0.1)