Evaluate abs(z)^(2n) -> z^n * conjugate(z)^n.

[ginac.git] / ginac / inifcns.cpp

diff --git a/ginac/inifcns.cpp b/ginac/inifcns.cpp
index 02c909f313180456e7598794500cad82c4a0370f..d4087e22b5003fd837d278f552e799fb765a1a97 100644 (file)
--- a/ginac/inifcns.cpp
+++ b/ginac/inifcns.cpp
@@ -302,10 +302,12 @@ static ex abs_imag_part(const ex& arg)
 
 static ex abs_power(const ex & arg, const ex & exp)
 {
-       if (arg.is_equal(arg.conjugate()) && ((is_a<numeric>(exp) && ex_to<numeric>(exp).is_even())
-                                               || exp.info(info_flags::even)))
-               return power(arg, exp);
-       else
+       if ((is_a<numeric>(exp) && ex_to<numeric>(exp).is_even()) || exp.info(info_flags::even)) {
+               if (arg.info(info_flags::real) || arg.is_equal(arg.conjugate()))
+                       return power(arg, exp);
+               else
+                       return power(arg, exp/2)*power(arg.conjugate(), exp/2);
+       } else
                return power(abs(arg), exp).hold();
 }