Merge pull request #115 from JuliaAlgebra/bl/unary_prod

Implement unary product

Author: blegat

src/operators.jl

@@ -45,6 +45,7 @@ Base.:-(m::AbstractMonomialLike) = (-1) * m
 Base.:-(t::AbstractTermLike) = (-coefficient(t)) * monomial(t)
 Base.:-(p::APL) = polynomial((-).(terms(p)))
 Base.:+(p::Union{APL, RationalPoly}) = p
+Base.:*(p::Union{APL, RationalPoly}) = p
 
 # Avoid adding a zero constant that might artificially increase the Newton polytope
 # Need to add polynomial conversion for type stability

test/polynomial.jl

@@ -42,6 +42,8 @@ const MP = MultivariatePolynomials
     Mod.@polyvar y
 
     p = 3x^2*y^4 + 2x
+    @test +(p) === p
+    @test *(p) === p
     @test terms(p)[end] == 2x
     typetests(p)
     typetests([p, x + y])