Add the case of c1=-1 and c2>0 to the rewriting rule of exponents: (x^c1)^c2 -> x...

[ginac.git] / check / exam_powerlaws.cpp

diff --git a/check/exam_powerlaws.cpp b/check/exam_powerlaws.cpp
index cab16dfb5852314e1da23a918bc4eceb738b973a..beb4dfb28d6f4da024b543114ff922f713066fc7 100644 (file)
--- a/check/exam_powerlaws.cpp
+++ b/check/exam_powerlaws.cpp
@@ -69,14 +69,14 @@ static unsigned exam_powerlaws1()
                return 1;
        }
        
-       ex e4 = e1.subs(lst(a==-1, b==2.5));
+       ex e4 = e1.subs(lst(a==-1, b==-2.5));
        if (!(is_exactly_a<power>(e4) &&
              is_exactly_a<power>(e4.op(0)) &&
              is_exactly_a<symbol>(e4.op(0).op(0)) &&
              is_exactly_a<numeric>(e4.op(0).op(1)) &&
              is_exactly_a<numeric>(e4.op(1)) &&
-             e4.is_equal(power(power(x,-1),2.5)) )) {
-               clog << "(x^a)^b, x symbolic, a==-1, b==2.5 wrong" << endl;
+             e4.is_equal(power(power(x,-1),-2.5)) )) {
+               clog << "(x^a)^b, x symbolic, a==-1, b==-2.5 wrong" << endl;
                clog << "returned: " << e4 << endl;
                return 1;
        }