Add rules to distribute conj over sum and prod

JuliaSymbolics · Dec 12, 2024 · dc73ca9 · dc73ca9
1 parent 59e3aa6
commit dc73ca9
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 0 deletions.
diff --git a/src/simplify_rules.jl b/src/simplify_rules.jl
@@ -63,6 +63,8 @@ let
         @rule(~x \ ~y => ~y / (~x))
         @rule(one(~x) => one(symtype(~x)))
         @rule(zero(~x) => zero(symtype(~x)))
+        @rule(conj(+(~~x)) => sum(conj.(~~x)))
+        @rule(conj(*(~~x)) => prod(conj.(~~x)))
         @rule(conj(~x::_isreal) => ~x)
         @rule(real(~x::_isreal) => ~x)
         @rule(imag(~x::_isreal) => zero(symtype(~x)))

diff --git a/test/rulesets.jl b/test/rulesets.jl
@@ -25,6 +25,8 @@ end
     @eqtest simplify(Term{Real}(real, [x])) == x
     @eqtest simplify(Term{Real}(imag, [x])) == 0
     @eqtest simplify(Term{Real}(imag, [y])) == imag(y)
+    @eqtest simplify(conj(x*y^2)) == x*conj(y^2)
+    @eqtest simplify(conj(x+y)) == x+conj(y)
     @eqtest simplify(x - y) == x + -1 * y
     @eqtest simplify(x - sin(y)) == x + -1 * sin(y)
     @eqtest simplify(-sin(x)) == -1 * sin(x)