Karim Belabas on Sat, 11 Aug 2012 08:56:48 +0200


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

Re: poolroots division by zero


* Sören Lennart Berg [2012-08-11 02:07]:
> I'm getting a 'division by zero' from polroots for the following input:
> 
> ? polroots(64.0*z^6-640.0*z^5+2640.0*z^4-5760.0*z^3+7040.0*z^2-4608.0*z+1280.0)
>   ***   at top-level: polroots(64.0*z^6-64
>   ***                 ^--------------------
>   *** polroots: division by zero
>   ***   Break loop: type 'break' to go back to GP
[...]

I cannot reproduce the problem in pari-2.5.2 (nor in the 'master'
branch), with either GMP or the native kernel, and valgrind does not
show anything suspicious. ( All on a 64bit machine. )

- What exact version were you using ?

- If not pari-2.5.2 (or master HEAD), can you upgrade and test again ?

Thanks for your report !

    K.B.
-- 
Karim Belabas, IMB (UMR 5251)  Tel: (+33) (0)5 40 00 26 17
Universite Bordeaux 1          Fax: (+33) (0)5 40 00 69 50
351, cours de la Liberation    http://www.math.u-bordeaux1.fr/~belabas/
F-33405 Talence (France)       http://pari.math.u-bordeaux1.fr/  [PARI/GP]
`