Karim BELABAS on Mon, 10 Mar 2003 19:29:22 +0100 (MET) |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: idealprimedec() bug |
On Mon, 10 Mar 2003, Igor Schein wrote: > 1 more: > > ? idealprimedec(nfinit(x^16+17*x^14+119*x^12+442*x^10+935*x^8+1122*x^6+714*x^4+204*x^2+17),4) > *** bug in GP (Segmentation Fault), please report This one occurs in a completely different place. FpX_normalize couldn't handle the 0 polynomial: this is actually a bug ! Now, the computation enters an infinite loop, ended by stack overflow. I don't think I want to go further here: it would require installing counters in many routines that expect prime numbers, and trigger a pseudo-primality test whenever the counters grow suspiciously larger. Overkill. 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]