| 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]