Karim BELABAS on Wed, 5 Feb 2003 14:41:35 +0100 (MET)


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

Re: galoisinit() bug


On Wed, 5 Feb 2003, Igor Schein wrote:
> On Wed, Feb 05, 2003 at 12:37:04AM -0500, Igor Schein wrote:
>> ? galoisinit(Pol(subst([1,0,4*x,0,2*x^2+24,0,16*x,0,16],x,10^2589)))
>>   ***   division by zero in R/R
>
> Sorry, it's actually polroots() bug, not galoisinit().

A 'double' underflow in polroots.c:split2() [ exp(-745.) = 0. ]

Fixed. [ of course, at defaut precision, most of the roots are 0. ]

    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/