Karim BELABAS on Sat, 18 Jan 2003 17:38:00 +0100 (MET)


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

Re: gp: algdep(1+.1^10,1) division by zero


On Mon, 13 Jan 2003, Michael Somos wrote:
> pari-dev,
>       Perhaps this is beating a dead horse, but I found this :
>
> ? algdep(1+.1^10,1)
>   ***   division by zero in gdiv, gdivgs or ginv
>
> which can be worked around, and 'algdep()' has definite problems
> otherwise, but still, division by zero should not happen.

Fixed. [ test explicitly for x = 0. instead of relying on gexpo(x) < -M ]

    Karim.
-- 
Karim Belabas                    Tel: (+33) (0)1 69 15 57 48
Dép. de Mathématiques, Bât. 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/