#4670: prime_pi for input ~ 10^10 causes PariError
---------------------------+------------------------------------------------
Reporter: roed | Owner: was
Type: defect | Status: new
Priority: minor | Milestone:
Component: number theory | Keywords:
---------------------------+------------------------------------------------
Calling the primepi function on a large pari integer (10^10) causes an
error. The issue is that in sage/libs/pari/gen.pyx the function
init_primes casts the input to an unsigned long. If we don't want to
allow initialization with input bigger than this, we should give a better
error.
sage: prime_pi(10^10)
---------------------------------------------------------------------------
PariError Traceback (most recent call
last)
/Users/Roed/Math/sage-3.2/<ipython console> in <module>()
/Users/Roed/Math/sage-3.2/local/lib/python2.5/site-
packages/sage/functions/transcendental.pyc in __call__(self, x)
363 from sage.rings.integer import Integer
364 pari.init_primes(pari(x)+Integer(1))
--> 365 return ZZ(pari(x).primepi())
366
367 def plot(self, xmin=0, xmax=100, *args, **kwds):
/Users/Roed/Math/sage-3.2/local/lib/python2.5/site-
packages/sage/libs/pari/gen.so in sage.libs.pari.gen._pari_trap
(sage/libs/pari/gen.c:37972)()
PariError: impossible assignment I-->S (23)
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4670>
Sage <http://sagemath.org/>
Sage - Open Source Mathematical Software: Building the Car Instead of
Reinventing the Wheel
--~--~---------~--~----~------------~-------~--~----~
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
-~----------~----~----~----~------~----~------~--~---
