#11475: improve prime_pi (speedup + small fixes)
-----------------------------+----------------------------------------------
Reporter: rohana | Owner: was
Type: enhancement | Status: needs_work
Priority: major | Milestone: sage-4.7.1
Component: number theory | Keywords: primes, prime counting, prime_pi
Work_issues: | Upstream: N/A
Reviewer: | Author: R. Andrew Ohana
Merged: | Dependencies:
-----------------------------+----------------------------------------------
Changes (by rohana):
* status: needs_review => needs_work
Comment:
I have verified that `prime_pi` gives correct output for `i*10**12` with
`0 <= i <= 1000` on x86_64 :). I am probably going to remove the bound
flag and instead put a warning about larger values not being well tested
if/when they are called.
Unfortunately, I have also come across an issue with solaris 10 sparc
(t2.math). The root of the cause is that `mpz_export` has the bit ordering
arguments flipped. I am not aware of another way to get unsigned ints out
of gmp integers, nor am I aware of a way to check the platform from
cython. Also, I don't know if other platforms have this same issue
(solaris86 or mac powerpc in particular), so if someone could test that,
please let me know.
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/11475#comment:33>
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.
