Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: PARI/gp wins my Riemann zeta speed contest, so far
Replies: 9   Last Post: Feb 15, 2013 10:43 PM

 Messages: [ Previous | Next ]
 David Bernier Posts: 3,892 Registered: 12/13/04
PARI/gp wins my Riemann zeta speed contest, so far
Posted: Jan 21, 2013 12:23 AM

I had done my own Riemann zeta computations with Bernoulli numbers,
computing in days a few 10,000 (say 55,000) decimals of
the imaginary part of the first non-trivial zero
1/2 + i*14.13... using PARI-gp.

This time, I used PARI-gp's own zeta(.):

? system(date);zeta(t);system(date)
Sun Jan 20 03:36:28 EST 2013
Sun Jan 20 07:10:17 EST 2013

// 20,000 decimals precision
// t is with 1/10^400 of first non-trivial zero.

3 hours and 34 minutes for 20017 significant digits
near 1/2 + i*14.134725141734693790457251983562470270784257

? \p
realprecision = 20017 significant digits (20000 digits displayed)
? a=zeta(t);

? ?
Help topics: for a list of relevant subtopics, type ?n for n in
0: user-defined functions (aliases, installed and user functions)
2: CONVERSIONS and similar elementary functions
3: TRANSCENDENTAL functions
4: NUMBER THEORETICAL functions
5: Functions related to ELLIPTIC CURVES
6: Functions related to general NUMBER FIELDS

12: The PARI community

Also:
=====================

On PARI-gp, I'm trying to get the readline library included.
That means backspace reprints the previous query, etc.

David Bernier

Date Subject Author
1/21/13 David Bernier
1/21/13 James Waldby
1/21/13 David Bernier
1/21/13 David Bernier
1/21/13 James Waldby
2/3/13 David Bernier
1/21/13 James Waldby
2/7/13 David Bernier
2/15/13 David Bernier
2/15/13 David Bernier