| Bill Allombert on Wed, 26 Jun 2013 15:08:00 +0200 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
| Re: Slow polynomial division? |
On Wed, Jun 26, 2013 at 10:41:55AM +0100, Richard in Reading wrote:
> ? allocatemem(2^33)
> *** Warning: new stack size = 8589934592 (8192.000 Mbytes).
> ? e=ellinit("69696bf1");
> ? f=elldivpol(e,11*15)/elldivpol(e,55);
>
> The divisor has a degree and coefficients with digit-length around 1500 and
> the dividend has degree and coefficient-digit-length around 15000. After
> running for 111 hours the division seemed not to have completed.
>
> Calculating the polynomials themselves seems to take a minute or so.
>
> Is this expected?
Not really. The division only takes 12 minutes on paridev (with both 2.6.1 and 2.5.4).
I suggest you use \ instead of / since you know the result is a polynomial.
f=elldivpol(e,11*15)\elldivpol(e,55);
Cheers,
Bill.