In the Grobner package, there is a function poly_pseudo_divide
for multivariable polynomial division. It seems that poly_pseudo_divide
requires polynomials with *integer* (not rational) coefficients:
(%i87) poly_pseudo_divide(x*y + x + 1,[x+y, x+5] , [x,y]);
(%o87) [[y+1,0],-y^2-y+1,1,2] <--- OK
(%i88) poly_pseudo_divide(x*y/ 2 + x + 1,[x+y, x+5] , [x,y]);
(Maxima seems to hang)
Questions:
(1) Is there an option variable that allows poly_pseudo_divide
to work for polynomials with rational coefficients?
(2) If not, does somebody have Maxima code for multivariable
polynomial division that works for polynomials with rational
coefficients? Preferably (for speed), this code work entirely
with CRE expressions.
Notes:
(a) I copied this note to the author of the Maxima Grobner package.
(b) I wrote a pre-processor for poly_pseudo_divide that multiplies
all the polynomials by a constant to transform all coefficients to
integers. I think it works OK, but it's not ideal.
Barton