That looks great!
As a performance test I tried this:
(79) -> h2323:=((h2*h3*h2*h3)^2);
Type:
XDistributedPolynomial(OrderedVariableList([x,y,z,w,x1,x2,x3,x4,x5]),Fraction(Integer))
Time: 0.00 (IN) + 2.71 (EV) + 0.00 (OT) = 2.71 sec
(80) -> dc_fact h2323
(80)
[
Bill Page wrote:
>
> On Sat, Nov 10, 2018 at 12:02 PM Waldek Hebisch
> wrote:
> ...
> >
> > I have now implemented the lift part of Davenport-Caruso method.
> > You fetch code at:
> >
> > http://www.math.uni.wroc.pl/~hebisch/fricas/dcfact2.input
> > http://www.math.uni.wroc.pl/~hebisch/fricas/nc_