Bill Hart
unread,Nov 14, 2009, 5:44:20 PM11/14/09Sign in to reply to author
Sign in to forward
You do not have permission to delete messages in this group
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to flint-devel
I've been busy coding up division routines for F_mpz_poly. At this
stage I've only been working on division in Z[x], not pseudo division.
I had to write up some mul_trunc_left routines and I've added a
basecase and divide and conquer routine for division in Z[x].
I also added a shiny new F_mpz_poly_div_exact, which does exact
division by working from one end with Hensel division and from the
other end using ordinary division. For sparse polynomials it is about
twice as fast as what we had before and for dense polynomials about
33% faster.
I also applied some optimisations to all the division code for the
case where you are doing an unbalanced division, i.e. dividing A by B
where len(A) < 2*len(B) - 1. It should be faster than the fmpz_poly
module in those cases now, though I didn't do any timings to confirm
that.
I'm not sure what is next. I think probably heuristic GCD. I have to
do pseudo division some time too I guess.
Bill.