Karim BELABAS on Thu, 10 Jul 2003 19:11:37 +0200 (MEST)


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

Re: numbpart() bug


On Thu, 10 Jul 2003, Igor Schein wrote:
>           GP/PARI CALCULATOR Version 2.2.7 (development CHANGES-1.775)
>     UltraSparc running solaris (C portable/GMP-4.1.2 kernel) 64-bit version
>                          compiled: Jun 28 2003, gcc-3.3
>               (readline v4.3 enabled, extended help not available)
>
> ? default(parisize,128m)
>    parisize = 128000000
> ? numbpart(10^15+2)
> Segmentation fault \\ Session crashes
>
> However,
>
> ? numbpart(10^15+3)
>   ***   arg to partition function must be <10^15.
>
> So 10^15+2 falls in the no man land.

Fixed in CVS. [ in floating point arithmetic, x < 10^15 is not the same as
log_10(x) < 15 ... ]

    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              http://www.math.u-psud.fr/~belabas/
F-91405 Orsay (France)            http://www.parigp-home.de/  [PARI/GP]