Karim BELABAS on Fri, 4 Jun 1999 15:31:59 +0200 (MET DST) |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: polredabs(,4) |
[Igor:] > In 2.0.14: > > ? printp(polredabs(x^12-4*x^9+22*x^6-12*x^3+2,4)~)' > > [(x^12 - 8 x^9 + 18 x^6 + 4 x^3 + 2)] > [(x^12 - 4 x^9 + 22 x^6 - 12 x^3 + 2)] > [(x^12 + 16)] > > In 2.0.15: > > ? printp(polredabs(x^12-4*x^9+22*x^6-12*x^3+2,4)~)' > > [[(x^12 - 4 x^9 + 22 x^6 - 12 x^3 + 2), (x)]] > > Looks like the meaning of flag=4 changed, but the documentation doesn't > reflect it. Is the old functionality lost? It did not change. A typo which prevented one flag from acting (and another one computed a too optimistic bound...) Karim. *** src/basemath/bibli1.c.orig Thu May 27 17:08:51 1999 --- src/basemath/bibli1.c Fri Jun 4 15:29:25 1999 *************** *** 2067,2073 **** { GEN p1,y,beta; ! if (nf_RAW) { long i, c = lg(z); y=cgetg(c,t_VEC); --- 2067,2073 ---- { GEN p1,y,beta; ! if (flun & nf_RAW) { long i, c = lg(z); y=cgetg(c,t_VEC); *************** *** 2693,2698 **** --- 2693,2699 ---- if (i == r1) p = gmul2n(p, 1); p = gadd(p, gnorm((GEN)r[i])); } + if (i == r1) p = gmul2n(p, 1); return p; } __ Karim Belabas email: Karim.Belabas@math.u-psud.fr Dep. de Mathematiques, Bat. 425 Universite Paris-Sud Tel: (00 33) 1 69 15 57 48 F-91405 Orsay (France) Fax: (00 33) 1 69 15 60 19 -- PARI/GP Home Page: http://hasse.mathematik.tu-muenchen.de/ntsw/pari/