Add tests

src/diff.jl

@@ -35,6 +35,8 @@ function differentiate(ps::AbstractArray{PT, N}, xs) where {N, PT<:ARPL}
     end
     qs
 end
+differentiate(ps::AbstractArray{<:ARPL}, v::AbstractVariable) = differentiate.(ps, v)
+
 differentiate(p::ARPL, xs) = [differentiate(p, x) for x in xs]
 
 # differentiate(p, [x, y]) with TypedPolynomials promote x to a Monomial

src/rational.jl

@@ -46,10 +46,10 @@ function (-)(r::RationalPoly, s::RationalPoly)
 end
 _minus(p, s::RationalPoly) = (p * s.den - s.num) / s.den
 _minus(s::RationalPoly, p) = (s.num - p * s.den) / s.den
-(-)(p::APL, r::RationalPoly) = _minus(r, p)
+(-)(p::APL, r::RationalPoly) = _minus(p, r)
 (-)(r::RationalPoly, p::APL) = _minus(r, p)
 (-)(r::RationalPoly, α) = _minus(r, α)
-(-)(α, r::RationalPoly) = _minus(r, α)
+(-)(α, r::RationalPoly) = _minus(α, r)
 (-)(r::RationalPoly) = (-r.num) / r.den
 
 (*)(r::RationalPoly, s::RationalPoly) = (r.num*s.num) / (r.den*s.den)

test/comp.jl

@@ -70,6 +70,8 @@ MultivariatePolynomials.terms(p::CustomPoly) = terms(p.p)
             @test 3 != 4x / 2x
             @test 4x / 2x == 2
             @test 3 != 4x / 2x
+            @test 1 - 1/x == (x - 1) / x
+            @test 1 - 1/x != (1 - x) / x
             @test x + x/x == 1 + x^2/x
             @test x - x/x == -(1 - x^2/x)
             @test (1+x)/x - 1 == 1/x

test/diff.jl

@@ -1,5 +1,7 @@
 @testset "Differentiation" begin
     Mod.@polyvar x y
+    @test differentiate(3, y) == 0
+    @test differentiate([x, y], y) == [0, 1]
     @test differentiate(true*x+true*x^2, y) == 0
     @inferred differentiate(true*x+true*x^2, y)
     @test differentiate(x*y + 3y^2 , [x, y]) == [y, x+6y]

test/mono.jl

@@ -82,5 +82,9 @@
         @test σ[2] in (1, 3)
         @test σ[3] in (2, 5)
         @test mergemonovec([[x*y, x, x*y], [x^2*y, x]]) == [x^2*y, x*y, x]
+        @test_throws ArgumentError monovec([1, 2], [x^2])
+        σ, X = sortmonovec((y, x))
+        @test σ == [2, 1]
+        @test X == [x, y]
     end
 end