Karim BELABAS on Thu, 9 Oct 2003 21:04:42 +0200 (MEST) |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: polredabs failure |
On Thu, 9 Oct 2003, Igor Schein wrote: > ? p1=x^16-4*x^15-20*x^14+80*x^13+138*x^12-712*x^11+736*x^10+1200*x^9-8138*x^8+9288*x^7+1904*x^6-25248*x^5+73988*x^4-19184*x^3+41192*x^2+69456*x+30724; > ? p2=x^16 + 66*x^12 + 3627*x^8 + 48114*x^4 + 531441; > ? p1==polredabs(p1) > 1 > \\ p1 already *reduced* > ? poldegree((v=polcompositum(p1,p2))[#v]) > 16 > \\ isomorphic > ? norml2(polroots(p1)) > 91.81780460041328907655758261 > ? norml2(polroots(p2)) > 87.63560920082657815311516524 > \\ so we should have polredabs(p1)==p2 A complicated typo. The smallest element was found, then deleted due to a typo which assessed to it the norm of another (larger) element, during the "Sorting" phase (when there are too many small vectors, the list overflows, and we need to discard some of them). Fixed in CVS. 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]