Implement mutable removeleadingterm

Author: blegat

src/monovec.jl

@@ -45,6 +45,11 @@ end
 Base.getindex(x::MonomialVector, i::Integer) = Monomial(x.vars, x.Z[i])
 Base.getindex(x::MonomialVector, i::CartesianIndex{1}) = x[i.I[1]]
 
+function Base.deleteat!(x::MonomialVector, i)
+    deleteat!(x.Z, i)
+    return x
+end
+
 Base.firstindex(x::MonomialVector) = firstindex(x.Z)
 Base.lastindex(x::MonomialVector) = lastindex(x.Z)
 Base.size(x::MonomialVector) = (length(x),)

src/poly.jl

@@ -121,6 +121,11 @@ MP.leadingterm(p::Polynomial) = iszero(p) ? zeroterm(p) : first(terms(p))
 function MP.removeleadingterm(p::Polynomial)
     Polynomial(p.a[2:end], p.x[2:end])
 end
+function MA.mutable_operate!(::typeof(MP.removeleadingterm), p::Polynomial)
+    deleteat!(p.a, 1)
+    deleteat!(p.x, 1)
+    return p
+end
 function MP.removemonomials(p::Polynomial, x::MonomialVector)
     # use the fact that monomials are sorted to do this O(n) instead of O(n^2)
     j = 1