#11475: improve prime_pi (speedup + small fixes)
-----------------------------------------------------------+----------------
Reporter: rohana | Owner:
was
Type: enhancement | Status:
needs_work
Priority: major | Milestone:
sage-5.1
Component: number theory | Resolution:
Keywords: primes, prime counting, prime_pi sd40.5 | Work issues:
Report Upstream: N/A | Reviewers:
Yann Laigle-Chapuy, Leif Leonhardy, William Stein, Karl-Dieter Crisman
Authors: R. Andrew Ohana | Merged in:
Dependencies: | Stopgaps:
-----------------------------------------------------------+----------------
Changes (by jdemeyer):
* status: positive_review => needs_work
Comment:
This fails on arando (Linux Ubuntu 12.04 i686, 32 bit):
{{{
sage -t --long -force_lib devel/sage/sage/functions/prime_pi.pyx
**********************************************************************
File
"/var/lib/buildbot/build/sage/arando-1/arando_full/build/sage-5.1.beta4/devel
/sage-main/sage/functions/prime_pi.pyx", line 123:
sage: prime_pi(10^10)
Expected:
455052511
Got:
70462980
**********************************************************************
File
"/var/lib/buildbot/build/sage/arando-1/arando_full/build/sage-5.1.beta4/devel
/sage-main/sage/functions/prime_pi.pyx", line 231:
sage: for i in (32..42): prime_pi(2^i) # long time (13s on sage.math,
2011)
Expected:
203280221
393615806
762939111
1480206279
2874398515
5586502348
10866266172
21151907950
41203088796
80316571436
156661034233
Got:
0
0
0
0
0
0
0
0
0
0
0
**********************************************************************
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/11475#comment:64>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
