# HG changeset patch
# User Burcin Erocal <burcin@erocal.org>
# Date 1273115129 -7200
# Node ID 6b51a9fc888e469965eefe8f75586da72738a3b1
# Parent  1fc650c9f6319f745be590a56c655f24b98ff5b2
Restrict automatic conjugation of power objects. The only automatic evaluation performed is to return the same expression when the basis is positive and exponent is a rational.

diff --git a/ginac/power.cpp b/ginac/power.cpp
--- a/ginac/power.cpp
+++ b/ginac/power.cpp
@@ -567,6 +567,14 @@
 			}
 
 			if (exponent_is_rational) {
+				// ^(-c1, c2) -> ((-1)^c2)*(c1^c2)
+				// for example sqrt(-2) -> I*sqrt(2)
+				if (num_basis->is_negative() &&
+						!num_basis->is_equal(*_num_1_p))
+					return mul(power(_ex_1, *num_exponent),
+						power((*num_basis)*(*_num_1_p),
+							*num_exponent));
+
 				const numeric res = num_basis->power(*num_exponent);
 				if (res.is_crational()) {
 					return res;
@@ -798,12 +806,19 @@
 
 ex power::conjugate() const
 {
+	// Restrict automatic evaluation of this to fix Sage bug #8775
+	if (is_a<numeric>(exponent) && ex_to<numeric>(exponent).is_rational()
+			&& basis.info(info_flags::positive))
+		return *this;
+	return conjugate_function(*this).hold();
+	/*
 	ex newbasis = basis.conjugate();
 	ex newexponent = exponent.conjugate();
 	if (are_ex_trivially_equal(basis, newbasis) && are_ex_trivially_equal(exponent, newexponent)) {
 		return *this;
 	}
 	return (new power(newbasis, newexponent))->setflag(status_flags::dynallocated);
+	*/
 }
 
 ex power::real_part() const