Karim BELABAS on Sat, 7 Sep 2002 20:47:44 +0200 (MEST)


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

Re: yet another rnfkummer() posting


On Thu, 5 Sep 2002, Igor Schein wrote:
> ? bezoutres(Pol(sin(x)+Pi),Pol(cos(x)));
>   ***   bug in subresext, please report

This one is easy for a change. In the subresultant algorithm, I have replaced
the final divisibility test by an approximate division, issuing a warning if
the remainder is not an exact 0.

(19:45) gp >  bezoutres(Pol(sin(x)+Pi),Pol(cos(x)));
  ***   Warning: inexact computation in subresext.

This is completely analogous to the older:

(19:45) gp >  bezout(Pol(sin(x)+Pi),Pol(cos(x)));
  ***   Warning: inexact computation in bezoutpol.

    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/