Xavier Roblot on Wed, 30 Jun 1999 13:23:51 -0400


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

Re: bnrL1() bug


Igor Schein wrote:

> ? bnf=bnfinit(x^2+x+258);
> ? bnr1=bnrinit(bnf,1,1);
> ? bnr7=bnrinit(bnf,7,1);
> ? bnrL1(bnr7)
>   ***   the PARI stack overflows !!!

The corresponding ray-class number is 840. This explains why this
computation is not possible without a _huge_ stack.

> ? bnrL1(bnr1)
>   ***   segmentation fault: bug in GP (please report).

Here is a patch. Replace the following line of GetValue1 in stark.c 

L2090  c = mpneg(gdiv(gmul(hk, Rk), wk));
           ^^^^^ 
by 

L2090  c = gneg_i(gdiv(gmul(hk, Rk), wk));
           ^^^^^^

---
Xavier Roblot