|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.