Add operate_to! with polynomial (#45)

* Add operate_to! with APL

* Add tests

* Fix format

Author: blegat

src/arithmetic.jl

@@ -80,6 +80,24 @@ function MA.operate!(op::Union{typeof(+),typeof(-),typeof(*)}, p::_APL, q::_AE)
     return MA.operate!(op, p, MP.polynomial(q))
 end
 
+function MA.operate_to!(
+    res::MP.AbstractPolynomial,
+    op::typeof(*),
+    p::_AE,
+    q::_APL,
+)
+    return MA.operate_to!(res, op, MP.polynomial(p), q)
+end
+
+function MA.operate_to!(
+    res::MP.AbstractPolynomial,
+    op::typeof(*),
+    p::_APL,
+    q::_AE,
+)
+    return MA.operate_to!(res, op, p, MP.polynomial(q))
+end
+
 # These are not implemented yet for arbitrary bases so we
 # fall back to polynomials
 

test/runtests.jl

@@ -116,6 +116,12 @@ function api_test(B::Type{<:MB.AbstractMonomialIndexed}, degree)
     @test const_alg_el(x[1] => x[2], x[2] => x[1]) == const_mono
     @test differentiate(const_alg_el, x) == differentiate(const_mono, x)
     @test differentiate(const_alg_el, x, 2) == differentiate(const_mono, x, 2)
+    @test MA.operate_to!(polynomial(const_mono), *, const_alg_el, const_mono) ==
+          const_mono
+    @test MA.operate_to!(polynomial(const_mono), *, const_mono, const_alg_el) ==
+          const_mono
+    @test dot([const_alg_el], [const_mono]) == const_mono
+    @test dot([const_mono], [const_alg_el]) == const_mono
 end
 
 function univ_orthogonal_test(