Karim BELABAS on Tue, 10 Sep 2002 16:49:58 +0200 (MEST) |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: bnfisnorm() broken again |
On Mon, 9 Sep 2002, Igor Schein wrote: > On Tue, Sep 10, 2002 at 03:38:40AM +0200, Karim BELABAS wrote: > > On Mon, 9 Sep 2002, Igor Schein wrote: > > > ? bnfisnorm(bnfinit(x^3-2),5) > > > *** division by zero in dvmdii > > > > A few oversights in the new implementation of rnfisnorm in the non-Galois > > case [had forgotten to test it over Q except for quadratic fields...] > > > > Fixed. > > How about this: > > ? bnfisnorm(bnfinit(x^3-x-1),13) > *** incorrect type in lift. Right, the special case where the solution is trivial (rnfisnorm(x) --> [1,x]). The code over Q in bnfisnorm() did not handle it properly. Note: in the CVS code, bnfisnorm is quite useless [ it used to be justified by the weird input format for rnfisnorm ]. N = rnfisnorminit(y, x^3-x-1); [ or rnfisnorminit(y, bnf) if you wish ] rnfisnorm(N, 13); works just as well. And is faster since bnfisnorm is a trivial wrapper for the above code anyway [ N is (a little) costly to compute, esp. if you give bnfisnorm a bnf which is Galois over Q, and do not use the right flag to mention that fact ]. 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/