| Karim BELABAS on Tue, 4 Feb 2003 21:33:08 +0100 (MET) |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
| Re: polredabs() observation |
On Fri, 31 Jan 2003, Igor Schein wrote:
> On Tue, Jan 28, 2003 at 10:09:08AM -0500, Igor Schein wrote:
>> polredabs() on the following polynomial is much slower in latest CVS compared to stable:
>>
>> x^16-224*x^14+19520*x^12-843488*x^10+19172700*x^8-223732800*x^6+1177835136*x^4-1779647552*x^2+171662404
>
> Actually, here's the very recent change that made polredabs()
> considerably slower in many cases:
>
> \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
> revision 1.627
> date: 2003/01/04 19:00:43; author: karim; state: Exp; lines: +2 -0
> 48- simplified fincke_pohst() [don't recompute gram_matrix]
> \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
>
> And here's a smaller yet even more drastic example:
>
> polredabs(x^4+51976*x^3+837697192*x^2+4716826989696*x+5560701377110272);
Typo: the smallvectors() bound was the square of what it should have been...
I was confused by the terse comment [ due to me, from code
reverse-engineering ] at the header of sqred1(), and didn't check which
"Cholesky form" it was about. It turned out it was not the standard one:-(.
Both fixed.
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/