Karim BELABAS on Wed, 15 Jan 2003 23:43:49 +0100 (MET)


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

Re: recent factorint() regression


On Wed, 15 Jan 2003, Igor Schein wrote:
> factorint(2300240759,3)
>
> gets stuck in an infinite loop.
>
> It was broken 5 days ago.

By me. I have incorporated trial division in the prime testing algorithms,
and catered for the small values. Two unrelated typos could affect the result
for arguments around 2^32. I have fixed both.

pari-2.1.5:
(22:45) gp > t=0;forstep(x = 2^32-10^7+1,2^32,2, t += isprime(x)); t
time = 2mn, 7,260 ms.
%1 = 450525

Current CVS:
time = 32,300 ms.
%1 = 450525

Cheers,

    Karim.
-- 
Karim Belabas                    Tel: (+33) (0)1 69 15 57 48
Dép. de Mathématiques, Bât. 425  Fax: (+33) (0)1 69 15 60 19
Université Paris-Sud             Email: Karim.Belabas@math.u-psud.fr
F-91405 Orsay (France)           http://www.math.u-psud.fr/~belabas/
--
PARI/GP Home Page: http://www.parigp-home.de/