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