* 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
#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)
{
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;
// `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
-numeric const & _num_12(void)
+const numeric & _num_12(void)
{
const static ex e = ex(numeric(-12));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex_12(void)
+const ex & _ex_12(void)
{
static ex * e = new ex(_num_12());
return *e;
}
// numeric -11
-numeric const & _num_11(void)
+const numeric & _num_11(void)
{
const static ex e = ex(numeric(-11));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex_11(void)
+const ex & _ex_11(void)
{
static ex * e = new ex(_num_11());
return *e;
}
// numeric -10
-numeric const & _num_10(void)
+const numeric & _num_10(void)
{
const static ex e = ex(numeric(-10));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex_10(void)
+const ex & _ex_10(void)
{
static ex * e = new ex(_num_10());
return *e;
}
// numeric -9
-numeric const & _num_9(void)
+const numeric & _num_9(void)
{
const static ex e = ex(numeric(-9));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex_9(void)
+const ex & _ex_9(void)
{
static ex * e = new ex(_num_9());
return *e;
}
// numeric -8
-numeric const & _num_8(void)
+const numeric & _num_8(void)
{
const static ex e = ex(numeric(-8));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex_8(void)
+const ex & _ex_8(void)
{
static ex * e = new ex(_num_8());
return *e;
}
// numeric -7
-numeric const & _num_7(void)
+const numeric & _num_7(void)
{
const static ex e = ex(numeric(-7));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex_7(void)
+const ex & _ex_7(void)
{
static ex * e = new ex(_num_7());
return *e;
}
// numeric -6
-numeric const & _num_6(void)
+const numeric & _num_6(void)
{
const static ex e = ex(numeric(-6));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex_6(void)
+const ex & _ex_6(void)
{
static ex * e = new ex(_num_6());
return *e;
}
// numeric -5
-numeric const & _num_5(void)
+const numeric & _num_5(void)
{
const static ex e = ex(numeric(-5));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex_5(void)
+const ex & _ex_5(void)
{
static ex * e = new ex(_num_5());
return *e;
}
// numeric -4
-numeric const & _num_4(void)
+const numeric & _num_4(void)
{
const static ex e = ex(numeric(-4));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex_4(void)
+const ex & _ex_4(void)
{
static ex * e = new ex(_num_4());
return *e;
}
// numeric -3
-numeric const & _num_3(void)
+const numeric & _num_3(void)
{
const static ex e = ex(numeric(-3));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex_3(void)
+const ex & _ex_3(void)
{
static ex * e = new ex(_num_3());
return *e;
}
// numeric -2
-numeric const & _num_2(void)
+const numeric & _num_2(void)
{
const static ex e = ex(numeric(-2));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex_2(void)
+const ex & _ex_2(void)
{
static ex * e = new ex(_num_2());
return *e;
}
// numeric -1
-numeric const & _num_1(void)
+const numeric & _num_1(void)
{
const static ex e = ex(numeric(-1));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex_1(void)
+const ex & _ex_1(void)
{
static ex * e = new ex(_num_1());
return *e;
}
// numeric -1/2
-numeric const & _num_1_2(void)
+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;
}
-ex const & _ex_1_2(void)
+const ex & _ex_1_2(void)
{
static ex * e = new ex(_num_1_2());
return *e;
}
// numeric -1/3
-numeric const & _num_1_3(void)
+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;
}
-ex const & _ex_1_3(void)
+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
-numeric const & _num0(void)
+const numeric & _num0(void)
{
const static ex e = ex(numeric(0));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex0(void)
+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
-numeric const & _num1_3(void)
+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;
}
-ex const & _ex1_3(void)
+const ex & _ex1_3(void)
{
static ex * e = new ex(_num1_3());
return *e;
}
// numeric 1/2
-numeric const & _num1_2(void)
+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;
}
-ex const & _ex1_2(void)
+const ex & _ex1_2(void)
{
static ex * e = new ex(_num1_2());
return *e;
}
// numeric 1
-numeric const & _num1(void)
+const numeric & _num1(void)
{
const static ex e = ex(numeric(1));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex1(void)
+const ex & _ex1(void)
{
static ex * e = new ex(_num1());
return *e;
}
// numeric 2
-numeric const & _num2(void)
+const numeric & _num2(void)
{
const static ex e = ex(numeric(2));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex2(void)
+const ex & _ex2(void)
{
static ex * e = new ex(_num2());
return *e;
}
// numeric 3
-numeric const & _num3(void)
+const numeric & _num3(void)
{
const static ex e = ex(numeric(3));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex3(void)
+const ex & _ex3(void)
{
static ex * e = new ex(_num3());
return *e;
}
// numeric 4
-numeric const & _num4(void)
+const numeric & _num4(void)
{
const static ex e = ex(numeric(4));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex4(void)
+const ex & _ex4(void)
{
static ex * e = new ex(_num4());
return *e;
}
// numeric 5
-numeric const & _num5(void)
+const numeric & _num5(void)
{
const static ex e = ex(numeric(5));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex5(void)
+const ex & _ex5(void)
{
static ex * e = new ex(_num5());
return *e;
}
// numeric 6
-numeric const & _num6(void)
+const numeric & _num6(void)
{
const static ex e = ex(numeric(6));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex6(void)
+const ex & _ex6(void)
{
static ex * e = new ex(_num6());
return *e;
}
// numeric 7
-numeric const & _num7(void)
+const numeric & _num7(void)
{
const static ex e = ex(numeric(7));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex7(void)
+const ex & _ex7(void)
{
static ex * e = new ex(_num7());
return *e;
}
// numeric 8
-numeric const & _num8(void)
+const numeric & _num8(void)
{
const static ex e = ex(numeric(8));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex8(void)
+const ex & _ex8(void)
{
static ex * e = new ex(_num8());
return *e;
}
// numeric 9
-numeric const & _num9(void)
+const numeric & _num9(void)
{
const static ex e = ex(numeric(9));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex9(void)
+const ex & _ex9(void)
{
static ex * e = new ex(_num9());
return *e;
}
// numeric 10
-numeric const & _num10(void)
+const numeric & _num10(void)
{
const static ex e = ex(numeric(10));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex10(void)
+const ex & _ex10(void)
{
static ex * e = new ex(_num10());
return *e;
}
// numeric 11
-numeric const & _num11(void)
+const numeric & _num11(void)
{
const static ex e = ex(numeric(11));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex11(void)
+const ex & _ex11(void)
{
static ex * e = new ex(_num11());
return *e;
}
// numeric 12
-numeric const & _num12(void)
+const numeric & _num12(void)
{
const static ex e = ex(numeric(12));
const static numeric * n = static_cast<const numeric *>(e.bp);
return *n;
}
-ex const & _ex12(void)
+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
// none
-#ifndef NO_GINAC_NAMESPACE
+#ifndef NO_NAMESPACE_GINAC
} // namespace GiNaC
-#endif // ndef NO_GINAC_NAMESPACE
+#endif // ndef NO_NAMESPACE_GINAC