X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=check%2Fexam_powerlaws.cpp;h=c269e9747e93cb139333e17b9afac443abe1ecfc;hp=f685ffb12e26f1782999322eb9073fa2e5ed8048;hb=d5c259ec7de1d0b78b8a6c8e18178520f2a05f3e;hpb=af922d5eb36ed70e4a9e3ffaf4c24492cf89a1a6

diff --git a/check/exam_powerlaws.cpp b/check/exam_powerlaws.cpp
index f685ffb1..c269e974 100644
--- a/check/exam_powerlaws.cpp
+++ b/check/exam_powerlaws.cpp
@@ -4,7 +4,7 @@
  *  this code, it is a sanity check rather deeply rooted in GiNaC's classes. */
 
 /*
- *  GiNaC Copyright (C) 1999-2000 Johannes Gutenberg University Mainz, Germany
+ *  GiNaC Copyright (C) 1999-2001 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
@@ -31,67 +31,67 @@ static unsigned exam_powerlaws1(void)
 	symbol a("a");
 	symbol b("b");
 	
-	ex e1=power(power(x,a),b);
+	ex e1 = power(power(x,a), b);
 	if (!(is_ex_exactly_of_type(e1,power) &&
-		  is_ex_exactly_of_type(e1.op(0),power) &&
-		  is_ex_exactly_of_type(e1.op(0).op(0),symbol) &&
-		  is_ex_exactly_of_type(e1.op(0).op(1),symbol) &&
-		  is_ex_exactly_of_type(e1.op(1),symbol) &&
-		  e1.is_equal(power(power(x,a),b)) )) {
+	      is_ex_exactly_of_type(e1.op(0),power) &&
+	      is_ex_exactly_of_type(e1.op(0).op(0),symbol) &&
+	      is_ex_exactly_of_type(e1.op(0).op(1),symbol) &&
+	      is_ex_exactly_of_type(e1.op(1),symbol) &&
+	      e1.is_equal(power(power(x,a),b)) )) {
 		clog << "(x^a)^b, x,a,b symbolic wrong" << endl;
 		clog << "returned: " << e1 << endl;
 		return 1;
 	}
 	
-	ex e2=e1.subs(a==1);
+	ex e2 = e1.subs(a==1);
 	if (!(is_ex_exactly_of_type(e2,power) &&
-		  is_ex_exactly_of_type(e2.op(0),symbol) &&
-		  is_ex_exactly_of_type(e2.op(1),symbol) &&
-		  e2.is_equal(power(x,b)) )) {
+	      is_ex_exactly_of_type(e2.op(0),symbol) &&
+	      is_ex_exactly_of_type(e2.op(1),symbol) &&
+	      e2.is_equal(power(x,b)) )) {
 		clog << "(x^a)^b, x,b symbolic, a==1 wrong" << endl;
 		clog << "returned: " << e2 << endl;
 		return 1;
 	}
 	
-	ex e3=e1.subs(a==-1);
+	ex e3 = e1.subs(a==-1);
 	if (!(is_ex_exactly_of_type(e3,power) &&
-		  is_ex_exactly_of_type(e3.op(0),power) &&
-		  is_ex_exactly_of_type(e3.op(0).op(0),symbol) &&
-		  is_ex_exactly_of_type(e3.op(0).op(1),numeric) &&
-		  is_ex_exactly_of_type(e3.op(1),symbol) &&
-		  e3.is_equal(power(power(x,-1),b)) )) {
+	      is_ex_exactly_of_type(e3.op(0),power) &&
+	      is_ex_exactly_of_type(e3.op(0).op(0),symbol) &&
+	      is_ex_exactly_of_type(e3.op(0).op(1),numeric) &&
+	      is_ex_exactly_of_type(e3.op(1),symbol) &&
+	      e3.is_equal(power(power(x,-1),b)) )) {
 		clog << "(x^a)^b, x,b symbolic, a==-1 wrong" << endl;
 		clog << "returned: " << e3 << endl;
 		return 1;
 	}
 	
-	ex e4=e1.subs(lst(a==-1,b==2.5));
+	ex e4 = e1.subs(lst(a==-1, b==2.5));
 	if (!(is_ex_exactly_of_type(e4,power) &&
-		  is_ex_exactly_of_type(e4.op(0),power) &&
-		  is_ex_exactly_of_type(e4.op(0).op(0),symbol) &&
-		  is_ex_exactly_of_type(e4.op(0).op(1),numeric) &&
-		  is_ex_exactly_of_type(e4.op(1),numeric) &&
-		  e4.is_equal(power(power(x,-1),2.5)) )) {
+	      is_ex_exactly_of_type(e4.op(0),power) &&
+	      is_ex_exactly_of_type(e4.op(0).op(0),symbol) &&
+	      is_ex_exactly_of_type(e4.op(0).op(1),numeric) &&
+	      is_ex_exactly_of_type(e4.op(1),numeric) &&
+	      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;
 	}
 	
-	ex e5=e1.subs(lst(a==-0.9,b==2.5));
+	ex e5 = e1.subs(lst(a==-0.9, b==2.5));
 	if (!(is_ex_exactly_of_type(e5,power) &&
-		  is_ex_exactly_of_type(e5.op(0),symbol) &&
-		  is_ex_exactly_of_type(e5.op(1),numeric) &&
-		  e5.is_equal(power(x,numeric(-0.9)*numeric(2.5))) )) {
+	      is_ex_exactly_of_type(e5.op(0),symbol) &&
+	      is_ex_exactly_of_type(e5.op(1),numeric) &&
+	      e5.is_equal(power(x,numeric(-0.9)*numeric(2.5))) )) {
 		clog << "(x^a)^b, x symbolic, a==-0.9, b==2.5 wrong" << endl;
 		clog << "returned: " << e5 << endl;
 		return 1;
 	}
 	
-	ex e6=e1.subs(lst(a==numeric(3)+numeric(5.3)*I,b==-5));
+	ex e6 = e1.subs(lst(a==numeric(3)+numeric(5.3)*I, b==-5));
 	if (!(is_ex_exactly_of_type(e6,power) &&
-		  is_ex_exactly_of_type(e6.op(0),symbol) &&
-		  is_ex_exactly_of_type(e6.op(1),numeric) &&
-		  e6.is_equal(power(x,numeric(-15)+numeric(5.3)*numeric(-5)*I)) )) {
+	      is_ex_exactly_of_type(e6.op(0),symbol) &&
+	      is_ex_exactly_of_type(e6.op(1),numeric) &&
+	      e6.is_equal(power(x,numeric(-15)+numeric(5.3)*numeric(-5)*I)) )) {
 		clog << "(x^a)^b, x symbolic, a==3+5.3*I, b==-5 wrong" << endl;
 		clog << "returned: " << e6 << endl;
 		return 1;
@@ -108,88 +108,88 @@ static unsigned exam_powerlaws2(void)
 	symbol a("a");
 	symbol b("b");
 	
-	ex e1=power(a*x,b);
+	ex e1 = power(a*x,b);
 	if (!(is_ex_exactly_of_type(e1,power) &&
-		  is_ex_exactly_of_type(e1.op(0),mul) &&
-		  (e1.op(0).nops()==2) &&
-		  is_ex_exactly_of_type(e1.op(0).op(0),symbol) &&
-		  is_ex_exactly_of_type(e1.op(0).op(1),symbol) &&
-		  is_ex_exactly_of_type(e1.op(1),symbol) &&
-		  e1.is_equal(power(a*x,b)) )) {
+	      is_ex_exactly_of_type(e1.op(0),mul) &&
+	      (e1.op(0).nops()==2) &&
+	      is_ex_exactly_of_type(e1.op(0).op(0),symbol) &&
+	      is_ex_exactly_of_type(e1.op(0).op(1),symbol) &&
+	      is_ex_exactly_of_type(e1.op(1),symbol) &&
+	      e1.is_equal(power(a*x,b)) )) {
 		clog << "(a*x)^b, x,a,b symbolic wrong" << endl;
 		clog << "returned: " << e1 << endl;
 		return 1;
 	}
 	
-	ex e2=e1.subs(a==3);
+	ex e2 = e1.subs(a==3);
 	if (!(is_ex_exactly_of_type(e2,power) &&
-		  is_ex_exactly_of_type(e2.op(0),mul) &&
-		  (e2.op(0).nops()==2) &&
-		  is_ex_exactly_of_type(e2.op(0).op(0),symbol) &&
-		  is_ex_exactly_of_type(e2.op(0).op(1),numeric) &&
-		  is_ex_exactly_of_type(e2.op(1),symbol) &&
-		  e2.is_equal(power(3*x,b)) )) {
+	      is_ex_exactly_of_type(e2.op(0),mul) &&
+	      (e2.op(0).nops()==2) &&
+	      is_ex_exactly_of_type(e2.op(0).op(0),symbol) &&
+	      is_ex_exactly_of_type(e2.op(0).op(1),numeric) &&
+	      is_ex_exactly_of_type(e2.op(1),symbol) &&
+	      e2.is_equal(power(3*x,b)) )) {
 		clog << "(a*x)^b, x,b symbolic, a==3 wrong" << endl;
 		clog << "returned: " << e2 << endl;
 		return 1;
 	}
 	
-	ex e3=e1.subs(b==-3);
+	ex e3 = e1.subs(b==-3);
 	if (!(is_ex_exactly_of_type(e3,mul) &&
-		  (e3.nops()==2) &&
-		  is_ex_exactly_of_type(e3.op(0),power) &&
-		  is_ex_exactly_of_type(e3.op(1),power) &&
-		  e3.is_equal(power(a,-3)*power(x,-3)) )) {
+	      (e3.nops()==2) &&
+	      is_ex_exactly_of_type(e3.op(0),power) &&
+	      is_ex_exactly_of_type(e3.op(1),power) &&
+	      e3.is_equal(power(a,-3)*power(x,-3)) )) {
 		clog << "(a*x)^b, x,a symbolic, b==-3 wrong" << endl;
 		clog << "returned: " << e3 << endl;
 		return 1;
 	}
 	
-	ex e4=e1.subs(b==4.5);
+	ex e4 = e1.subs(b==4.5);
 	if (!(is_ex_exactly_of_type(e4,power) &&
-		  is_ex_exactly_of_type(e4.op(0),mul) &&
-		  (e4.op(0).nops()==2) &&
-		  is_ex_exactly_of_type(e4.op(0).op(0),symbol) &&
-		  is_ex_exactly_of_type(e4.op(0).op(1),symbol) &&
-		  is_ex_exactly_of_type(e4.op(1),numeric) &&
-		  e4.is_equal(power(a*x,4.5)) )) {
+	      is_ex_exactly_of_type(e4.op(0),mul) &&
+	      (e4.op(0).nops()==2) &&
+	      is_ex_exactly_of_type(e4.op(0).op(0),symbol) &&
+	      is_ex_exactly_of_type(e4.op(0).op(1),symbol) &&
+	      is_ex_exactly_of_type(e4.op(1),numeric) &&
+	      e4.is_equal(power(a*x,4.5)) )) {
 		clog << "(a*x)^b, x,a symbolic, b==4.5 wrong" << endl;
 		clog << "returned: " << e4 << endl;
 		return 1;
 	}
 	
-	ex e5=e1.subs(lst(a==3.2,b==3+numeric(5)*I));
+	ex e5 = e1.subs(lst(a==3.2, b==3+numeric(5)*I));
 	if (!(is_ex_exactly_of_type(e5,mul) &&
-		  (e5.nops()==2) &&
-		  is_ex_exactly_of_type(e5.op(0),power) &&
-		  is_ex_exactly_of_type(e5.op(1),numeric) &&
-		  e5.is_equal(power(x,3+numeric(5)*I)*
+	      (e5.nops()==2) &&
+	      is_ex_exactly_of_type(e5.op(0),power) &&
+	      is_ex_exactly_of_type(e5.op(1),numeric) &&
+	      e5.is_equal(power(x,3+numeric(5)*I)*
 					  power(numeric(3.2),3+numeric(5)*I)) )) {
 		clog << "(a*x)^b, x symbolic, a==3.2, b==3+5*I wrong" << endl;
 		clog << "returned: " << e5 << endl;
 		return 1;
 	}
 	
-	ex e6=e1.subs(lst(a==-3.2,b==3+numeric(5)*I));
+	ex e6 = e1.subs(lst(a==-3.2, b==3+numeric(5)*I));
 	if (!(is_ex_exactly_of_type(e6,mul) &&
-		  (e6.nops()==2) &&
-		  is_ex_exactly_of_type(e6.op(0),power) &&
-		  is_ex_exactly_of_type(e6.op(1),numeric) &&
-		  e6.is_equal(power(-x,3+numeric(5)*I)*
+	      (e6.nops()==2) &&
+	      is_ex_exactly_of_type(e6.op(0),power) &&
+	      is_ex_exactly_of_type(e6.op(1),numeric) &&
+	      e6.is_equal(power(-x,3+numeric(5)*I)*
 					  power(numeric(3.2),3+numeric(5)*I)) )) {
 		clog << "(a*x)^b, x symbolic, a==-3.2, b==3+5*I wrong" << endl;
 		clog << "returned: " << e6 << endl;
 		return 1;
 	}
 	
-	ex e7=e1.subs(lst(a==3+numeric(5)*I,b==3.2));
+	ex e7 = e1.subs(lst(a==3+numeric(5)*I, b==3.2));
 	if (!(is_ex_exactly_of_type(e7,power) &&
-		  is_ex_exactly_of_type(e7.op(0),mul) &&
-		  (e7.op(0).nops()==2) &&
-		  is_ex_exactly_of_type(e7.op(0).op(0),symbol) &&
-		  is_ex_exactly_of_type(e7.op(0).op(1),numeric) &&
-		  is_ex_exactly_of_type(e7.op(1),numeric) &&
-		  e7.is_equal(power((3+numeric(5)*I)*x,3.2)) )) {
+	      is_ex_exactly_of_type(e7.op(0),mul) &&
+	      (e7.op(0).nops()==2) &&
+	      is_ex_exactly_of_type(e7.op(0).op(0),symbol) &&
+	      is_ex_exactly_of_type(e7.op(0).op(1),numeric) &&
+	      is_ex_exactly_of_type(e7.op(1),numeric) &&
+	      e7.is_equal(power((3+numeric(5)*I)*x,3.2)) )) {
 		clog << "(a*x)^b, x symbolic, a==3+5*I, b==3.2 wrong" << endl;
 		clog << "returned: " << e7 << endl;
 		return 1;
@@ -216,8 +216,8 @@ static unsigned exam_powerlaws3(void)
 	
 	ex e3 = power(numeric(5),numeric(1,2));
 	if (!(is_ex_exactly_of_type(e3,power) &&
-		  e3.op(0).is_equal(numeric(5)) &&
-		  e3.op(1).is_equal(numeric(1,2)))) {
+	      e3.op(0).is_equal(numeric(5)) &&
+	      e3.op(1).is_equal(numeric(1,2)))) {
 		clog << "5^(1/2) wrongly returned " << e3 << endl;
 		return 1;
 	}
@@ -277,7 +277,7 @@ static unsigned exam_powerlaws5(void)
 	ex e2 = pow(0,a);
 	if (!(is_ex_exactly_of_type(e2,power))) {
 		clog << "0^a was evaluated to " << e2
-			 << " though nothing is known about a." << endl;
+		     << " though nothing is known about a." << endl;
 		return 1;
 	}