X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=check%2Fexam_misc.cpp;h=8c28fb5344edcde971182e42616bb75b5f721381;hp=fc9c41ebc48e61d302f544fce1fb04bf662f7220;hb=13905769fe865764a271f0e2aeb22f422e1289d7;hpb=9506030d2059ef868d5da8dcbe714c9c5b50ff50

diff --git a/check/exam_misc.cpp b/check/exam_misc.cpp
index fc9c41eb..8c28fb53 100644
--- a/check/exam_misc.cpp
+++ b/check/exam_misc.cpp
@@ -3,7 +3,7 @@
  */
 
 /*
- *  GiNaC Copyright (C) 1999-2010 Johannes Gutenberg University Mainz, Germany
+ *  GiNaC Copyright (C) 1999-2020 Johannes Gutenberg University Mainz, Germany
  *
  *  This program is free software; you can redistribute it and/or modify
  *  it under the terms of the GNU General Public License as published by
@@ -89,49 +89,6 @@ static unsigned exam_expand_power()
 	return result;
 }
 
-static unsigned exam_sqrfree()
-{
-	unsigned result = 0;
-	symbol x("x"), y("y");
-	ex e1, e2;
-	
-	e1 = (1+x)*pow((2+x),2)*pow((3+x),3)*pow((4+x),4);
-	e2 = sqrfree(expand(e1),lst(x));
-	if (e1 != e2) {
-		clog << "sqrfree(expand(" << e1 << ")) erroneously returned "
-		     << e2 << endl;
-		++result;
-	}
-	
-	e1 = (x+y)*pow((x+2*y),2)*pow((x+3*y),3)*pow((x+4*y),4);
-	e2 = sqrfree(expand(e1));
-	if (e1 != e2) {
-		clog << "sqrfree(expand(" << e1 << ")) erroneously returned "
-		     << e2 << endl;
-		++result;
-	}
-	e2 = sqrfree(expand(e1),lst(x));
-	if (e1 != e2) {
-		clog << "sqrfree(expand(" << e1 << "),[x]) erroneously returned "
-		     << e2 << endl;
-		++result;
-	}
-	e2 = sqrfree(expand(e1),lst(y));
-	if (e1 != e2) {
-		clog << "sqrfree(expand(" << e1 << "),[y]) erroneously returned "
-		     << e2 << endl;
-		++result;
-	}
-	e2 = sqrfree(expand(e1),lst(x,y));
-	if (e1 != e2) {
-		clog << "sqrfree(expand(" << e1 << "),[x,y]) erroneously returned "
-		     << e2 << endl;
-		++result;
-	}
-	
-	return result;
-}
-
 /* Arithmetic Operators should behave just as one expects from built-in types.
  * When somebody screws up the operators this routine will most probably fail
  * to compile.  Unfortunately we can only test the stuff that is allowed, not
@@ -164,7 +121,7 @@ static unsigned exam_operator_semantics()
 		++result;
 	}
 	
-	// Prefix/postfix increment/decrement behaviour:
+	// Prefix/postfix increment/decrement behavior:
 	e1 = 7; e2 = 4;
 	i1 = 7; i2 = 4;
 	e1 = (--e2 = 2)++;
@@ -220,7 +177,7 @@ static unsigned exam_subs()
 	}
 
 	// And this used to fail in GiNaC 1.5.8 because it first substituted
-	// exp(x) -> exp(log(x)) -> x, and then substitued again x -> log(x)
+	// exp(x) -> exp(log(x)) -> x, and then substituted again x -> log(x)
 	e1 = exp(x);
 	e2 = e1.subs(x == log(x));
 	if (!e2.is_equal(x)) {
@@ -298,6 +255,29 @@ static unsigned exam_subs_algebraic()
 	return result;
 }
 
+/* Test suitable cases of the exponent power law: (e^t)^s=e^(ts). */
+static unsigned exam_exponent_power_law()
+{
+	unsigned result = 0;
+	symbol x("x");
+	realsymbol s("s");
+	possymbol t("t");
+
+	exmap pwr_exp =
+		{ {pow(exp(x), 2), exp(2*x)},
+		  {pow(exp(s), t), exp(s*t)},
+		  {exp(x)*pow(exp(x),-1), 1} };
+
+	for (auto e : pwr_exp) {
+		if (! (e.first.is_equal(e.second)) ) {
+			clog << "power of exponent  " << e.first << " produces error.\n";
+			++result;
+		}
+	}
+
+	return result;
+}
+
 unsigned exam_misc()
 {
 	unsigned result = 0;
@@ -307,11 +287,11 @@ unsigned exam_misc()
 	result += exam_expand_subs();  cout << '.' << flush;
 	result += exam_expand_subs2();  cout << '.' << flush;
 	result += exam_expand_power(); cout << '.' << flush;
-	result += exam_sqrfree(); cout << '.' << flush;
 	result += exam_operator_semantics(); cout << '.' << flush;
 	result += exam_subs(); cout << '.' << flush;
 	result += exam_joris(); cout << '.' << flush;
 	result += exam_subs_algebraic(); cout << '.' << flush;
+	result += exam_exponent_power_law(); cout << '.' << flush;
 	
 	return result;
 }