Bill Allombert on Tue, 16 Dec 2014 18:57:56 +0100

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

Re: polynomial partial fractions

On Tue, Dec 16, 2014 at 09:17:33PM +1100, Kevin Ryde wrote:
> Karim Belabas <> writes:
> >
> >   C = vector(#P, i, Mod(B \/ PE[i], PE[i]));
> >   [Q, vector(#P, i, padic(lift(A/C[i]), P[i])), PE];
> Nice.  I don't think I knew to go by poly inverses before.  How does
> that rate against a matsolve()?  About the same but less taking apart
> and putting back together polynomials perhaps?

Inverting POLMOD is done using the extended Euclid algorithm which is likely to
be much faster than matsolve() (O(n^2) instead of O(n^3)).