Karim BELABAS on Wed, 27 Aug 2003 19:09:39 +0200 (MEST) |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: factor() boob |
On Wed, 13 Aug 2003, Phil Carmody wrote: > The following version of GP: > > GP/PARI CALCULATOR Version 2.2.5 (alpha) > i686 running linux (ix86/GMP-4.1.2 kernel) 32-bit version > (readline v4.1 enabled, extended help available) > > > does naughty things when asked to factor simple difference of squares > quadratics: > > > gp > for(i=1,908,if(length(factor(x^2-i*i)~)==1,print(i"\t("x^2-i*i")\t"factor(x^2-i*i)))) > 58 (x^2 - 3364) Mat([x^2 - 3364, 1]) [...] > Has this been noticed before or even already fixed? Yes and yes [ introduced in 2.2.5, fixed in 2.2.6 ( missing factor 2 in root_bound(), when bounding the root of largest modulus ). ]. When using the unstable version, you should try and upgrade once in a while ( 2.2.5 is 6 months old already!. As usual, CVS is the preferred way to keep up. ) 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 http://www.math.u-psud.fr/~belabas/ F-91405 Orsay (France) http://www.parigp-home.de/ [PARI/GP]