Merge pull request #104 from JuliaSymbolics/s/polynormalize

export mpoly functionality separately as polynormalize

Author: shashi

src/simplify.jl

@@ -47,18 +47,18 @@ function simplify(x, ctx=DefaultCtx();
                   rules=default_rules(x, ctx),
                   fixpoint=true,
                   applyall=true,
-                  mpoly=false,
                   kwargs...)
-    if mpoly
-        x = to_term(to_mpoly(x)...)
-    end
     if fixpoint
         SymbolicUtils.fixpoint(rules, x, ctx; applyall=applyall)
     else
         rules(x, ctx; applyall=applyall, kwargs...)
     end
 end
 
+function polynormalize(x)
+    to_term(to_mpoly(x)...)
+end
+
 
 Base.@deprecate simplify(x, rules::RuleSet; kwargs...)  simplify(x, rules=rules; kwargs...)
 

test/nf.jl

@@ -1,8 +1,8 @@
 using SymbolicUtils, Test
+using SymbolicUtils: polynormalize
 @testset "polyform" begin
     @syms a b c d
-    @test simplify(a * (b + -1 * c) + -1 * (b * a + -1 * c * a), mpoly=true) == 0
-    @eqtest simplify(sin((a+b)^2)^2; mpoly=true) == sin(a^2+b^2+2*a*b)^2
-    # fixme: can this be made faster?
-    @test simplify(sin((a+b)^2)^2 + cos((a+b)^2)^2; mpoly=true) == 1
+    @test polynormalize(a * (b + -1 * c) + -1 * (b * a + -1 * c * a)) == 0
+    @eqtest simplify(polynormalize(sin((a+b)^2)^2)) == simplify(sin(a^2+2*(b*a)+b^2)^2)
+    @test simplify(polynormalize(sin((a+b)^2)^2 + cos((a+b)^2)^2)) == 1
 end