X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Futils.cpp;h=49edabee177ae8a6a05aae870f5e20ce267a4062;hp=e113c04bbdef31bbc72bb09063cdc074b9f13038;hb=ba0fa799375b03d3e8d512b0a6b2b6f1ef14e6e5;hpb=955cb185a85535ab328ffedbfccdc508ce80fa91 diff --git a/ginac/utils.cpp b/ginac/utils.cpp index e113c04b..49edabee 100644 --- a/ginac/utils.cpp +++ b/ginac/utils.cpp @@ -4,7 +4,7 @@ * 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 @@ -25,438 +25,734 @@ #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 (ab) { - return 1; - } - return 0; + if (ab) { + 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(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(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(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(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(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(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(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(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(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(e.bp); - return *n; + const static ex e = ex(numeric(-12)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex_12(void) +const ex & _ex_12(void) { - static ex * e = new ex(_num_12()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(-11)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex_11(void) +const ex & _ex_11(void) { - static ex * e = new ex(_num_11()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(-10)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex_10(void) +const ex & _ex_10(void) { - static ex * e = new ex(_num_10()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(-9)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex_9(void) +const ex & _ex_9(void) { - static ex * e = new ex(_num_9()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(-8)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex_8(void) +const ex & _ex_8(void) { - static ex * e = new ex(_num_8()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(-7)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex_7(void) +const ex & _ex_7(void) { - static ex * e = new ex(_num_7()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(-6)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex_6(void) +const ex & _ex_6(void) { - static ex * e = new ex(_num_6()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(-5)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex_5(void) +const ex & _ex_5(void) { - static ex * e = new ex(_num_5()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(-4)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex_4(void) +const ex & _ex_4(void) { - static ex * e = new ex(_num_4()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(-3)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex_3(void) +const ex & _ex_3(void) { - static ex * e = new ex(_num_3()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(-2)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex_2(void) +const ex & _ex_2(void) { - static ex * e = new ex(_num_2()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(-1)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex_1(void) +const ex & _ex_1(void) { - static ex * e = new ex(_num_1()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(-1,2)); + const static numeric * n = static_cast(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; + 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(e.bp); + return *n; +} + +const ex & _ex_1_3(void) { - const static ex e = ex(numeric(-1,3)); - const static numeric * n = static_cast(e.bp); - return *n; + 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(e.bp); + return *n; } -ex const & _ex_1_3(void) +const ex & _ex_1_4(void) { - static ex * e = new ex(_num_1_3()); - return *e; + 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(e.bp); + return *n; +} + +const ex & _ex0(void) { - const static ex e = ex(numeric(0)); - const static numeric * n = static_cast(e.bp); - return *n; + static ex * e = new ex(_num0()); + return *e; } -ex const & _ex0(void) +// numeric 1/4 +const numeric & _num1_4(void) { - static ex * e = new ex(_num0()); - return *e; + const static ex e = ex(numeric(1,4)); + const static numeric * n = static_cast(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(e.bp); - return *n; + const static ex e = ex(numeric(1,3)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex1_3(void) +const ex & _ex1_3(void) { - static ex * e = new ex(_num1_3()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(1,2)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex1_2(void) +const ex & _ex1_2(void) { - static ex * e = new ex(_num1_2()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(1)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex1(void) +const ex & _ex1(void) { - static ex * e = new ex(_num1()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(2)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex2(void) +const ex & _ex2(void) { - static ex * e = new ex(_num2()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(3)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex3(void) +const ex & _ex3(void) { - static ex * e = new ex(_num3()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(4)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex4(void) +const ex & _ex4(void) { - static ex * e = new ex(_num4()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(5)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex5(void) +const ex & _ex5(void) { - static ex * e = new ex(_num5()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(6)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex6(void) +const ex & _ex6(void) { - static ex * e = new ex(_num6()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(7)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex7(void) +const ex & _ex7(void) { - static ex * e = new ex(_num7()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(8)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex8(void) +const ex & _ex8(void) { - static ex * e = new ex(_num8()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(9)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex9(void) +const ex & _ex9(void) { - static ex * e = new ex(_num9()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(10)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex10(void) +const ex & _ex10(void) { - static ex * e = new ex(_num10()); - return *e; + 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(e.bp); - return *n; + const static ex e = ex(numeric(11)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex11(void) +const ex & _ex11(void) { - static ex * e = new ex(_num11()); - return *e; + 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(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(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(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(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(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(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(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(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(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(12)); - const static numeric * n = static_cast(e.bp); - return *n; + const static ex e = ex(numeric(120)); + const static numeric * n = static_cast(e.bp); + return *n; } -ex const & _ex12(void) +const ex & _ex120(void) { - static ex * e = new ex(_num12()); - return *e; + static ex * e = new ex(_num120()); + return *e; } // comment skeleton for header files @@ -464,24 +760,24 @@ ex const & _ex12(void) // 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 @@ -534,6 +830,6 @@ ex const & _ex12(void) // none -#ifndef NO_GINAC_NAMESPACE +#ifndef NO_NAMESPACE_GINAC } // namespace GiNaC -#endif // ndef NO_GINAC_NAMESPACE +#endif // ndef NO_NAMESPACE_GINAC