Re: rnfkummer-induced bug

[Karim BELABAS <Karim.Belabas@math.u-psud.fr>]

On Fri, 9 Aug 2002, Igor Schein wrote:
> This resulted in more regression:
>
> ? setrand(1);rnfkummer(bnrinit(bnfinit(quadpoly(1020,y)),31,1),matdiagonal([5,1,1]));
>   ***   division by zero in powgi.

I've smoothed this out [ get_arch_real() didn't obey its specification in a
trivial special case, which was never used before ]

I've committed all my remaining rnfkummer patches. In its current state,
rnfkummer should be able to treat extensions of arbitrary prime degree (in
particular larger than 5...), again provided that bnfinit( K(zeta_ell) )
succeeds.

I don't intend to modify it further. You can crash-test it now :-)

    Karim.
-- 
Karim Belabas                    Tel: (+33) (0)1 69 15 57 48
Dép. de Mathematiques, Bat. 425  Fax: (+33) (0)1 69 15 60 19
Université Paris-Sud             Email: Karim.Belabas@math.u-psud.fr
F-91405 Orsay (France)           http://www.math.u-psud.fr/~belabas/
--
PARI/GP Home Page: http://www.parigp-home.de/