Bill Allombert on Fri, 27 Jan 2012 10:02:39 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: primepi bug


On Thu, Jan 26, 2012 at 07:21:23PM -0800, John W. Nicholson wrote:
> ? primepi(nextprime(4000000453))
> %25 = 189961830
> ? primepi(nextprime(4000000454))
>   ***   at top-level: primepi(nextprime(40
>   ***                 ^--------------------
>   *** primepi: bug in PARI/GP (Segmentation Fault), please report
>   ***   Break loop: type 'break' to go back to GP
> break> break

This is a bug which is already fixed in the development branch.
(LGr12- off-by-1 error in primepi() when argument is = primelimit)
However you will just get a different error:
  *** primepi: not enough precomputed primes, need primelimit ~ 4000000483.

The fast sieve-based primepi algorithm is not implemented.

Cheers,
Bill.