X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Futils.cpp;h=49edabee177ae8a6a05aae870f5e20ce267a4062;hp=7592e4033e8ef6ae7d0eadcb1fffe10102ab7db2;hb=ba0fa799375b03d3e8d512b0a6b2b6f1ef14e6e5;hpb=afdd7fa8c6c0a587f7c80789198551383e8beb7b

diff --git a/ginac/utils.cpp b/ginac/utils.cpp
index 7592e403..49edabee 100644
--- a/ginac/utils.cpp
+++ b/ginac/utils.cpp
@@ -1,9 +1,10 @@
 /** @file utils.cpp
  *
- *  Implementation of several small and furry utilities. */
+ *  Implementation of several small and furry utilities needed within GiNaC
+ *  but not of any interest to the user of the library. */
 
 /*
- *  GiNaC Copyright (C) 1999 Johannes Gutenberg University Mainz, Germany
+ *  GiNaC Copyright (C) 1999-2000 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
@@ -20,56 +21,763 @@
  *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  */
 
+#include "ex.h"
+#include "numeric.h"
 #include "utils.h"
 
-#ifndef NO_GINAC_NAMESPACE
+#ifndef NO_NAMESPACE_GINAC
 namespace GiNaC {
-#endif // ndef NO_GINAC_NAMESPACE
+#endif // ndef NO_NAMESPACE_GINAC
 
+/** ctor for pole_error exception class. */
+pole_error::pole_error(const std::string& what_arg, int degree)
+  : domain_error(what_arg)
+{
+	deg = degree;
+}
+
+/** Return the degree of the pole_error exception class. */
+int pole_error::degree(void) const
+{
+	return deg;
+}
+
+// some compilers (e.g. cygwin) define a macro log2, causing confusion
+#ifndef log2
 /** Integer binary logarithm */
 unsigned log2(unsigned n)
 {
-    unsigned k;
-    for (k = 0; n > 1; n >>= 1) ++k;
-    return k;
+	unsigned k;
+	for (k = 0; n > 1; n >>= 1) ++k;
+	return k;
 }
+#endif
 
 /** Compare two pointers (just to establish some sort of canonical order).
  *  @return -1, 0, or 1 */
-int compare_pointers(void const * a, void const * b)
+int compare_pointers(const void * a, const void * b)
+{
+	if (a<b) {
+		return -1;
+	} else if (a>b) {
+		return 1;
+	}
+	return 0;
+}
+
+//////////
+// `construct on first use' chest of numbers
+//////////
+
+// numeric -120
+const numeric & _num_120(void)
+{
+	const static ex e = ex(numeric(-120));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_120(void)
+{
+	static ex * e = new ex(_num_120());
+	return *e;
+}
+
+// numeric -60
+const numeric & _num_60(void)
+{
+	const static ex e = ex(numeric(-60));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_60(void)
+{
+	static ex * e = new ex(_num_60());
+	return *e;
+}
+
+// numeric -48
+const numeric & _num_48(void)
+{
+	const static ex e = ex(numeric(-48));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_48(void)
+{
+	static ex * e = new ex(_num_48());
+	return *e;
+}
+
+// numeric -30
+const numeric & _num_30(void)
+{
+	const static ex e = ex(numeric(-30));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_30(void)
+{
+	static ex * e = new ex(_num_30());
+	return *e;
+}
+
+// numeric -25
+const numeric & _num_25(void)
+{
+	const static ex e = ex(numeric(-25));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_25(void)
+{
+	static ex * e = new ex(_num_25());
+	return *e;
+}
+
+// numeric -24
+const numeric & _num_24(void)
+{
+	const static ex e = ex(numeric(-24));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_24(void)
+{
+	static ex * e = new ex(_num_24());
+	return *e;
+}
+
+// numeric -20
+const numeric & _num_20(void)
+{
+	const static ex e = ex(numeric(-20));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_20(void)
+{
+	static ex * e = new ex(_num_20());
+	return *e;
+}
+
+// numeric -18
+const numeric & _num_18(void)
+{
+	const static ex e = ex(numeric(-18));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_18(void)
+{
+	static ex * e = new ex(_num_18());
+	return *e;
+}
+
+// numeric -15
+const numeric & _num_15(void)
+{
+	const static ex e = ex(numeric(-15));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_15(void)
+{
+	static ex * e = new ex(_num_15());
+	return *e;
+}
+
+// numeric -12
+const numeric & _num_12(void)
+{
+	const static ex e = ex(numeric(-12));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_12(void)
+{
+	static ex * e = new ex(_num_12());
+	return *e;
+}
+
+// numeric -11
+const numeric & _num_11(void)
+{
+	const static ex e = ex(numeric(-11));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_11(void)
 {
-    if (a<b) {
-        return -1;
-    } else if (a>b) {
-        return 1;
-    }
-    return 0;
+	static ex * e = new ex(_num_11());
+	return *e;
 }
 
+// numeric -10
+const numeric & _num_10(void)
+{
+	const static ex e = ex(numeric(-10));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_10(void)
+{
+	static ex * e = new ex(_num_10());
+	return *e;
+}
+
+// numeric -9
+const numeric & _num_9(void)
+{
+	const static ex e = ex(numeric(-9));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_9(void)
+{
+	static ex * e = new ex(_num_9());
+	return *e;
+}
+
+// numeric -8
+const numeric & _num_8(void)
+{
+	const static ex e = ex(numeric(-8));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_8(void)
+{
+	static ex * e = new ex(_num_8());
+	return *e;
+}
+
+// numeric -7
+const numeric & _num_7(void)
+{
+	const static ex e = ex(numeric(-7));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_7(void)
+{
+	static ex * e = new ex(_num_7());
+	return *e;
+}
+
+// numeric -6
+const numeric & _num_6(void)
+{
+	const static ex e = ex(numeric(-6));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_6(void)
+{
+	static ex * e = new ex(_num_6());
+	return *e;
+}
+
+// numeric -5
+const numeric & _num_5(void)
+{
+	const static ex e = ex(numeric(-5));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_5(void)
+{
+	static ex * e = new ex(_num_5());
+	return *e;
+}
+
+// numeric -4
+const numeric & _num_4(void)
+{
+	const static ex e = ex(numeric(-4));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_4(void)
+{
+	static ex * e = new ex(_num_4());
+	return *e;
+}
+
+// numeric -3
+const numeric & _num_3(void)
+{
+	const static ex e = ex(numeric(-3));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_3(void)
+{
+	static ex * e = new ex(_num_3());
+	return *e;
+}
+
+// numeric -2
+const numeric & _num_2(void)
+{
+	const static ex e = ex(numeric(-2));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_2(void)
+{
+	static ex * e = new ex(_num_2());
+	return *e;
+}
+
+// numeric -1
+const numeric & _num_1(void)
+{
+	const static ex e = ex(numeric(-1));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_1(void)
+{
+	static ex * e = new ex(_num_1());
+	return *e;
+}
+
+// numeric -1/2
+const numeric & _num_1_2(void)
+{
+	const static ex e = ex(numeric(-1,2));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_1_2(void)
+{
+	static ex * e = new ex(_num_1_2());
+	return *e;
+}    
+
+// numeric -1/3
+const numeric & _num_1_3(void)
+{
+	const static ex e = ex(numeric(-1,3));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_1_3(void)
+{
+	static ex * e = new ex(_num_1_3());
+	return *e;
+}    
+
+// numeric -1/4
+const numeric & _num_1_4(void)
+{
+	const static ex e = ex(numeric(-1,4));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex_1_4(void)
+{
+	static ex * e = new ex(_num_1_4());
+	return *e;
+}    
+
+// numeric  0
+const numeric & _num0(void)
+{
+	const static ex e = ex(numeric(0));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex0(void)
+{
+	static ex * e = new ex(_num0());
+	return *e;
+}
+
+// numeric  1/4
+const numeric & _num1_4(void)
+{
+	const static ex e = ex(numeric(1,4));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex1_4(void)
+{
+	static ex * e = new ex(_num1_4());
+	return *e;
+}    
+
+// numeric  1/3
+const numeric & _num1_3(void)
+{
+	const static ex e = ex(numeric(1,3));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex1_3(void)
+{
+	static ex * e = new ex(_num1_3());
+	return *e;
+}    
+
+// numeric  1/2
+const numeric & _num1_2(void)
+{
+	const static ex e = ex(numeric(1,2));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex1_2(void)
+{
+	static ex * e = new ex(_num1_2());
+	return *e;
+}    
+
+// numeric  1
+const numeric & _num1(void)
+{
+	const static ex e = ex(numeric(1));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex1(void)
+{
+	static ex * e = new ex(_num1());
+	return *e;
+}
+
+// numeric  2
+const numeric & _num2(void)
+{
+	const static ex e = ex(numeric(2));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex2(void)
+{
+	static ex * e = new ex(_num2());
+	return *e;
+}
+
+// numeric  3
+const numeric & _num3(void)
+{
+	const static ex e = ex(numeric(3));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex3(void)
+{
+	static ex * e = new ex(_num3());
+	return *e;
+}
+
+// numeric  4
+const numeric & _num4(void)
+{
+	const static ex e = ex(numeric(4));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex4(void)
+{
+	static ex * e = new ex(_num4());
+	return *e;
+}
+
+// numeric  5
+const numeric & _num5(void)
+{
+	const static ex e = ex(numeric(5));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex5(void)
+{
+	static ex * e = new ex(_num5());
+	return *e;
+}
+
+// numeric  6
+const numeric & _num6(void)
+{
+	const static ex e = ex(numeric(6));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex6(void)
+{
+	static ex * e = new ex(_num6());
+	return *e;
+}
+
+// numeric  7
+const numeric & _num7(void)
+{
+	const static ex e = ex(numeric(7));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex7(void)
+{
+	static ex * e = new ex(_num7());
+	return *e;
+}
+
+// numeric  8
+const numeric & _num8(void)
+{
+	const static ex e = ex(numeric(8));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex8(void)
+{
+	static ex * e = new ex(_num8());
+	return *e;
+}
+
+// numeric  9
+const numeric & _num9(void)
+{
+	const static ex e = ex(numeric(9));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex9(void)
+{
+	static ex * e = new ex(_num9());
+	return *e;
+}
+
+// numeric  10
+const numeric & _num10(void)
+{
+	const static ex e = ex(numeric(10));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex10(void)
+{
+	static ex * e = new ex(_num10());
+	return *e;
+}
+
+// numeric  11
+const numeric & _num11(void)
+{
+	const static ex e = ex(numeric(11));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex11(void)
+{
+	static ex * e = new ex(_num11());
+	return *e;
+}
+
+// numeric  12
+const numeric & _num12(void)
+{
+	const static ex e = ex(numeric(12));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex12(void)
+{
+	static ex * e = new ex(_num12());
+	return *e;
+}
+
+// numeric  15
+const numeric & _num15(void)
+{
+	const static ex e = ex(numeric(15));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex15(void)
+{
+	static ex * e = new ex(_num15());
+	return *e;
+}
+
+// numeric  18
+const numeric & _num18(void)
+{
+	const static ex e = ex(numeric(18));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex18(void)
+{
+	static ex * e = new ex(_num18());
+	return *e;
+}
+
+// numeric  20
+const numeric & _num20(void)
+{
+	const static ex e = ex(numeric(20));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex20(void)
+{
+	static ex * e = new ex(_num20());
+	return *e;
+}
+
+// numeric  24
+const numeric & _num24(void)
+{
+	const static ex e = ex(numeric(24));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex24(void)
+{
+	static ex * e = new ex(_num24());
+	return *e;
+}
+
+// numeric  25
+const numeric & _num25(void)
+{
+	const static ex e = ex(numeric(25));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex25(void)
+{
+	static ex * e = new ex(_num25());
+	return *e;
+}
+
+// numeric  30
+const numeric & _num30(void)
+{
+	const static ex e = ex(numeric(30));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex30(void)
+{
+	static ex * e = new ex(_num30());
+	return *e;
+}
+
+// numeric  48
+const numeric & _num48(void)
+{
+	const static ex e = ex(numeric(48));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex48(void)
+{
+	static ex * e = new ex(_num48());
+	return *e;
+}
+
+// numeric  60
+const numeric & _num60(void)
+{
+	const static ex e = ex(numeric(60));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex60(void)
+{
+	static ex * e = new ex(_num60());
+	return *e;
+}
+
+// numeric  120
+const numeric & _num120(void)
+{
+	const static ex e = ex(numeric(120));
+	const static numeric * n = static_cast<const numeric *>(e.bp);
+	return *n;
+}
+
+const ex & _ex120(void)
+{
+	static ex * e = new ex(_num120());
+	return *e;
+}
 
 // comment skeleton for header files
 
 
 // member functions
 
-    // default constructor, destructor, copy constructor assignment operator and helpers
-    // none
+	// default constructor, destructor, copy constructor assignment operator and helpers
+	// none
 
-    // other constructors
-    // none
+	// other constructors
+	// none
 
-    // functions overriding virtual functions from bases classes
-    // none
-    
-    // new virtual functions which can be overridden by derived classes
-    // none
+	// functions overriding virtual functions from bases classes
+	// none
+	
+	// new virtual functions which can be overridden by derived classes
+	// none
 
-    // non-virtual functions in this class
-    // none
+	// non-virtual functions in this class
+	// none
 
 // member variables
 // none
-    
+	
 
 
 // comment skeleton for implementation files
@@ -121,6 +829,7 @@ int compare_pointers(void const * a, void const * b)
 // private
 // none
 
-#ifndef NO_GINAC_NAMESPACE
+
+#ifndef NO_NAMESPACE_GINAC
 } // namespace GiNaC
-#endif // ndef NO_GINAC_NAMESPACE
+#endif // ndef NO_NAMESPACE_GINAC