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/