Karim Belabas on Fri, 21 Sep 2012 21:47:48 +0200


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

Re: Bug in APR-CL


* Bill Allombert [2012-09-21 18:38]:
> On Fri, Sep 21, 2012 at 11:33:32AM -0400, Charles Greathouse wrote:
> > > isprime(5368962301599408606279497323618896374219)
> >   ***   at top-level: isprime(536896230159
> >   ***                 ^--------------------
> >   *** isprime: bug in aprcl: e(t) too small, please report.
> >   ***   Break loop: type 'break' to go back to GP prompt
> 
> Thanks for your report!
> 
> This bug was introduced in rev
> 
> commit 95fd5cf91a7bc08a39c008896f4a1b73257dde96
> Author: Karim Belabas <Karim.Belabas@math.u-bordeaux1.fr>
> Date:   Fri Aug 10 14:10:11 2012 +0200
> 
>     rewrite e / compt -> compute_e / compute_t
> 
>     - recompute table t -> e(t)
>     - add zv_prod_Z
>     - remove calcjac [ no point in precomputing, compute as needed ]
>     - allow larger prime divisors in e(t) [ avoid infinite loops ]
>     - reduce memory needed by flaglp
> 

and is now fixed in

commit b0900773f2ffab3813cb0e26b1709f54b76c85e9
Author: Karim Belabas <Karim.Belabas@math.u-bordeaux1.fr>
Date:   Fri Sep 21 21:40:50 2012 +0200

    fix rounding error in script used to generate compute_t_small()
    
    isprime(5368962301599408606279497323618896374219) => e(t) too small


Thanks for your report !

    K.B.
-- 
Karim Belabas, IMB (UMR 5251)  Tel: (+33) (0)5 40 00 26 17
Universite Bordeaux 1          Fax: (+33) (0)5 40 00 69 50
351, cours de la Liberation    http://www.math.u-bordeaux1.fr/~belabas/
F-33405 Talence (France)       http://pari.math.u-bordeaux1.fr/  [PARI/GP]
`