Re: rnflllgram() regression

[Igor Schein <igor@txc.com>]

On Mon, Mar 10, 2003 at 10:41:09PM +0100, Karim BELABAS wrote:
> I have fixed the numerical instability (not the algorithm). It also fixes an
> old TODO item:
> 
>   4  rnfpolred is numerically unstable:
>      ? rnfpolred(nfinit(quadpoly(904,y)),quadray(904,1))
>        ***   division by zero in gdiv, gdivgs or ginv

This got broken though:

? rnfpolred(nfinit(quadpoly(1129,y)),quadray(1129,1))
  ***   the PARI stack overflows !
  current stack size: 64000000 (61.035 Mbytes)
  [hint] you can increase GP stack with allocatemem()

Didn't need more than a default stack before.

Thanks

Igor