Igor Schein on Mon, 10 Mar 2003 20:16:30 -0500

[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: rnflllgram() regression

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.