- power::derivative(): slight optimization.

diff --git a/ginac/power.cpp b/ginac/power.cpp
index a0c6ea441620719d6b4dfb7d9f857b8a487999a3..14878a9187f40a588b0ea498336793fd5f1de943 100644 (file)
--- a/ginac/power.cpp
+++ b/ginac/power.cpp
@@ -545,9 +545,9 @@ ex power::derivative(const symbol & s) const
                return mul(newseq, exponent);
        } else {
                // D(b^e) = b^e * (D(e)*ln(b) + e*D(b)/b)
                return mul(newseq, exponent);
        } else {
                // D(b^e) = b^e * (D(e)*ln(b) + e*D(b)/b)

-               return mul(power(basis, exponent),
+               return mul(*this,

                           add(mul(exponent.diff(s), log(basis)),
                           add(mul(exponent.diff(s), log(basis)),

-                          mul(mul(exponent, basis.diff(s)), power(basis, -1))));
+                          mul(mul(exponent, basis.diff(s)), power(basis, _ex_1()))));

        }
 }
 
        }
 }