X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Futils.cpp;h=07c0a1f6ff86f532634491c3eccb12687d47e75e;hp=e113c04bbdef31bbc72bb09063cdc074b9f13038;hb=22abfbe8c78e339188096a5bf749a7c2d4f0a368;hpb=955cb185a85535ab328ffedbfccdc508ce80fa91 diff --git a/ginac/utils.cpp b/ginac/utils.cpp index e113c04b..07c0a1f6 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-2005 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 @@ -18,445 +18,380 @@ * * You should have received a copy of the GNU General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ #include "ex.h" #include "numeric.h" #include "utils.h" +#include "version.h" -#ifndef NO_GINAC_NAMESPACE namespace GiNaC { -#endif // ndef NO_GINAC_NAMESPACE -/** Integer binary logarithm */ -unsigned log2(unsigned n) -{ - unsigned k; - for (k = 0; n > 1; n >>= 1) ++k; - return k; -} - -/** 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) -{ - if (ab) { - return 1; - } - return 0; -} - -////////// -// `construct on first use' chest of numbers -////////// - -// numeric -12 -numeric const & _num_12(void) -{ - const static ex e = ex(numeric(-12)); - const static numeric * n = static_cast(e.bp); - return *n; -} - -ex const & _ex_12(void) -{ - static ex * e = new ex(_num_12()); - return *e; -} - -// numeric -11 -numeric const & _num_11(void) -{ - const static ex e = ex(numeric(-11)); - const static numeric * n = static_cast(e.bp); - return *n; -} - -ex const & _ex_11(void) -{ - static ex * e = new ex(_num_11()); - return *e; -} - -// numeric -10 -numeric const & _num_10(void) -{ - const static ex e = ex(numeric(-10)); - const static numeric * n = static_cast(e.bp); - return *n; -} - -ex const & _ex_10(void) -{ - static ex * e = new ex(_num_10()); - return *e; -} - -// numeric -9 -numeric const & _num_9(void) -{ - const static ex e = ex(numeric(-9)); - const static numeric * n = static_cast(e.bp); - return *n; -} - -ex const & _ex_9(void) -{ - static ex * e = new ex(_num_9()); - return *e; -} - -// numeric -8 -numeric const & _num_8(void) -{ - const static ex e = ex(numeric(-8)); - const static numeric * n = static_cast(e.bp); - return *n; -} - -ex const & _ex_8(void) -{ - static ex * e = new ex(_num_8()); - return *e; -} - -// numeric -7 -numeric const & _num_7(void) -{ - const static ex e = ex(numeric(-7)); - const static numeric * n = static_cast(e.bp); - return *n; -} - -ex const & _ex_7(void) -{ - static ex * e = new ex(_num_7()); - return *e; -} - -// numeric -6 -numeric const & _num_6(void) -{ - const static ex e = ex(numeric(-6)); - const static numeric * n = static_cast(e.bp); - return *n; -} - -ex const & _ex_6(void) -{ - static ex * e = new ex(_num_6()); - return *e; -} - -// numeric -5 -numeric const & _num_5(void) -{ - const static ex e = ex(numeric(-5)); - const static numeric * n = static_cast(e.bp); - return *n; -} - -ex const & _ex_5(void) -{ - static ex * e = new ex(_num_5()); - return *e; -} - -// numeric -4 -numeric const & _num_4(void) -{ - const static ex e = ex(numeric(-4)); - const static numeric * n = static_cast(e.bp); - return *n; -} - -ex const & _ex_4(void) -{ - static ex * e = new ex(_num_4()); - return *e; -} - -// numeric -3 -numeric const & _num_3(void) -{ - const static ex e = ex(numeric(-3)); - const static numeric * n = static_cast(e.bp); - return *n; -} - -ex const & _ex_3(void) -{ - static ex * e = new ex(_num_3()); - return *e; -} - -// numeric -2 -numeric const & _num_2(void) -{ - const static ex e = ex(numeric(-2)); - const static numeric * n = static_cast(e.bp); - return *n; -} +/* Version information buried into the library */ +const int version_major = GINACLIB_MAJOR_VERSION; +const int version_minor = GINACLIB_MINOR_VERSION; +const int version_micro = GINACLIB_MICRO_VERSION; -ex const & _ex_2(void) -{ - static ex * e = new ex(_num_2()); - return *e; -} -// numeric -1 -numeric const & _num_1(void) -{ - const static ex e = ex(numeric(-1)); - const static numeric * n = static_cast(e.bp); - return *n; -} +/** ctor for pole_error exception class. */ +pole_error::pole_error(const std::string& what_arg, int degree) + : domain_error(what_arg), deg(degree) { } -ex const & _ex_1(void) +/** Return the degree of the pole_error exception class. */ +int pole_error::degree() const { - static ex * e = new ex(_num_1()); - return *e; + return deg; } -// numeric -1/2 -numeric const & _num_1_2(void) +/** Integer binary logarithm */ +unsigned log2(unsigned n) { - const static ex e = ex(numeric(-1,2)); - const static numeric * n = static_cast(e.bp); - return *n; + unsigned k; + for (k = 0; n > 1; n >>= 1) + ++k; + return k; } -ex const & _ex_1_2(void) -{ - static ex * e = new ex(_num_1_2()); - return *e; -} -// numeric -1/3 -numeric const & _num_1_3(void) -{ - const static ex e = ex(numeric(-1,3)); - const static numeric * n = static_cast(e.bp); - return *n; -} +////////// +// flyweight chest of numbers is initialized here: +////////// -ex const & _ex_1_3(void) -{ - static ex * e = new ex(_num_1_3()); - return *e; -} +/** How many static objects were created? Only the first one must create + * the static flyweights on the heap. */ +int library_init::count = 0; -// numeric 0 -numeric const & _num0(void) -{ - const static ex e = ex(numeric(0)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -120 +const numeric *_num_120_p; +const ex _ex_120 = _ex_120; -ex const & _ex0(void) -{ - static ex * e = new ex(_num0()); - return *e; -} +// static numeric -60 +const numeric *_num_60_p; +const ex _ex_60 = _ex_60; -// numeric 1/3 -numeric const & _num1_3(void) -{ - const static ex e = ex(numeric(1,3)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -48 +const numeric *_num_48_p; +const ex _ex_48 = _ex_48; -ex const & _ex1_3(void) -{ - static ex * e = new ex(_num1_3()); - return *e; -} +// static numeric -30 +const numeric *_num_30_p; +const ex _ex_30 = _ex_30; -// numeric 1/2 -numeric const & _num1_2(void) -{ - const static ex e = ex(numeric(1,2)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -25 +const numeric *_num_25_p; +const ex _ex_25 = _ex_25; -ex const & _ex1_2(void) -{ - static ex * e = new ex(_num1_2()); - return *e; -} +// static numeric -24 +const numeric *_num_24_p; +const ex _ex_24 = _ex_24; -// numeric 1 -numeric const & _num1(void) -{ - const static ex e = ex(numeric(1)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -20 +const numeric *_num_20_p; +const ex _ex_20 = _ex_20; -ex const & _ex1(void) -{ - static ex * e = new ex(_num1()); - return *e; -} +// static numeric -18 +const numeric *_num_18_p; +const ex _ex_18 = _ex_18; -// numeric 2 -numeric const & _num2(void) -{ - const static ex e = ex(numeric(2)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -15 +const numeric *_num_15_p; +const ex _ex_15 = _ex_15; -ex const & _ex2(void) -{ - static ex * e = new ex(_num2()); - return *e; -} +// static numeric -12 +const numeric *_num_12_p; +const ex _ex_12 = _ex_12; -// numeric 3 -numeric const & _num3(void) -{ - const static ex e = ex(numeric(3)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -11 +const numeric *_num_11_p; +const ex _ex_11 = _ex_11; -ex const & _ex3(void) -{ - static ex * e = new ex(_num3()); - return *e; -} +// static numeric -10 +const numeric *_num_10_p; +const ex _ex_10 = _ex_10; -// numeric 4 -numeric const & _num4(void) -{ - const static ex e = ex(numeric(4)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -9 +const numeric *_num_9_p; +const ex _ex_9 = _ex_9; -ex const & _ex4(void) -{ - static ex * e = new ex(_num4()); - return *e; -} - -// numeric 5 -numeric const & _num5(void) -{ - const static ex e = ex(numeric(5)); - const static numeric * n = static_cast(e.bp); - return *n; -} - -ex const & _ex5(void) -{ - static ex * e = new ex(_num5()); - return *e; -} - -// numeric 6 -numeric const & _num6(void) -{ - const static ex e = ex(numeric(6)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -8 +const numeric *_num_8_p; +const ex _ex_8 = _ex_8; -ex const & _ex6(void) -{ - static ex * e = new ex(_num6()); - return *e; -} +// static numeric -7 +const numeric *_num_7_p; +const ex _ex_7 = _ex_7; -// numeric 7 -numeric const & _num7(void) -{ - const static ex e = ex(numeric(7)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -6 +const numeric *_num_6_p; +const ex _ex_6 = _ex_6; -ex const & _ex7(void) -{ - static ex * e = new ex(_num7()); - return *e; -} - -// numeric 8 -numeric const & _num8(void) -{ - const static ex e = ex(numeric(8)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -5 +const numeric *_num_5_p; +const ex _ex_5 = _ex_5; -ex const & _ex8(void) -{ - static ex * e = new ex(_num8()); - return *e; -} +// static numeric -4 +const numeric *_num_4_p; +const ex _ex_4 = _ex_4; -// numeric 9 -numeric const & _num9(void) -{ - const static ex e = ex(numeric(9)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -3 +const numeric *_num_3_p; +const ex _ex_3 = _ex_3; -ex const & _ex9(void) -{ - static ex * e = new ex(_num9()); - return *e; -} +// static numeric -2 +const numeric *_num_2_p; +const ex _ex_2 = _ex_2; -// numeric 10 -numeric const & _num10(void) -{ - const static ex e = ex(numeric(10)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -1 +const numeric *_num_1_p; +const ex _ex_1 = _ex_1; -ex const & _ex10(void) -{ - static ex * e = new ex(_num10()); - return *e; -} +// static numeric -1/2 +const numeric *_num_1_2_p; +const ex _ex_1_2= _ex_1_2; -// numeric 11 -numeric const & _num11(void) -{ - const static ex e = ex(numeric(11)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric -1/3 +const numeric *_num_1_3_p; +const ex _ex_1_3= _ex_1_3; -ex const & _ex11(void) -{ - static ex * e = new ex(_num11()); - return *e; -} +// static numeric -1/4 +const numeric *_num_1_4_p; +const ex _ex_1_4= _ex_1_4; -// numeric 12 -numeric const & _num12(void) -{ - const static ex e = ex(numeric(12)); - const static numeric * n = static_cast(e.bp); - return *n; -} +// static numeric 0 +const numeric *_num0_p; +const basic *_num0_bp; +const ex _ex0 = _ex0; -ex const & _ex12(void) -{ - static ex * e = new ex(_num12()); - return *e; +// static numeric 1/4 +const numeric *_num1_4_p; +const ex _ex1_4 = _ex1_4; + +// static numeric 1/3 +const numeric *_num1_3_p; +const ex _ex1_3 = _ex1_3; + +// static numeric 1/2 +const numeric *_num1_2_p; +const ex _ex1_2 = _ex1_2; + +// static numeric 1 +const numeric *_num1_p; +const ex _ex1 = _ex1; + +// static numeric 2 +const numeric *_num2_p; +const ex _ex2 = _ex2; + +// static numeric 3 +const numeric *_num3_p; +const ex _ex3 = _ex3; + +// static numeric 4 +const numeric *_num4_p; +const ex _ex4 = _ex4; + +// static numeric 5 +const numeric *_num5_p; +const ex _ex5 = _ex5; + +// static numeric 6 +const numeric *_num6_p; +const ex _ex6 = _ex6; + +// static numeric 7 +const numeric *_num7_p; +const ex _ex7 = _ex7; + +// static numeric 8 +const numeric *_num8_p; +const ex _ex8 = _ex8; + +// static numeric 9 +const numeric *_num9_p; +const ex _ex9 = _ex9; + +// static numeric 10 +const numeric *_num10_p; +const ex _ex10 = _ex10; + +// static numeric 11 +const numeric *_num11_p; +const ex _ex11 = _ex11; + +// static numeric 12 +const numeric *_num12_p; +const ex _ex12 = _ex12; + +// static numeric 15 +const numeric *_num15_p; +const ex _ex15 = _ex15; + +// static numeric 18 +const numeric *_num18_p; +const ex _ex18 = _ex18; + +// static numeric 20 +const numeric *_num20_p; +const ex _ex20 = _ex20; + +// static numeric 24 +const numeric *_num24_p; +const ex _ex24 = _ex24; + +// static numeric 25 +const numeric *_num25_p; +const ex _ex25 = _ex25; + +// static numeric 30 +const numeric *_num30_p; +const ex _ex30 = _ex30; + +// static numeric 48 +const numeric *_num48_p; +const ex _ex48 = _ex48; + +// static numeric 60 +const numeric *_num60_p; +const ex _ex60 = _ex60; + +// static numeric 120 +const numeric *_num120_p; +const ex _ex120 = _ex120; + +/** Ctor of static initialization helpers. The fist call to this is going + * to initialize the library, the others do nothing. */ +library_init::library_init() +{ + if (count++==0) { + (_num_120_p = new numeric(-120))->setflag(status_flags::dynallocated); + (_num_60_p = new numeric(-60))->setflag(status_flags::dynallocated); + (_num_48_p = new numeric(-48))->setflag(status_flags::dynallocated); + (_num_30_p = new numeric(-30))->setflag(status_flags::dynallocated); + (_num_25_p = new numeric(-25))->setflag(status_flags::dynallocated); + (_num_24_p = new numeric(-24))->setflag(status_flags::dynallocated); + (_num_20_p = new numeric(-20))->setflag(status_flags::dynallocated); + (_num_18_p = new numeric(-18))->setflag(status_flags::dynallocated); + (_num_15_p = new numeric(-15))->setflag(status_flags::dynallocated); + (_num_12_p = new numeric(-12))->setflag(status_flags::dynallocated); + (_num_11_p = new numeric(-11))->setflag(status_flags::dynallocated); + (_num_10_p = new numeric(-10))->setflag(status_flags::dynallocated); + (_num_9_p = new numeric(-9))->setflag(status_flags::dynallocated); + (_num_8_p = new numeric(-8))->setflag(status_flags::dynallocated); + (_num_7_p = new numeric(-7))->setflag(status_flags::dynallocated); + (_num_6_p = new numeric(-6))->setflag(status_flags::dynallocated); + (_num_5_p = new numeric(-5))->setflag(status_flags::dynallocated); + (_num_4_p = new numeric(-4))->setflag(status_flags::dynallocated); + (_num_3_p = new numeric(-3))->setflag(status_flags::dynallocated); + (_num_2_p = new numeric(-2))->setflag(status_flags::dynallocated); + (_num_1_p = new numeric(-1))->setflag(status_flags::dynallocated); + (_num_1_2_p = new numeric(-1,2))->setflag(status_flags::dynallocated); + (_num_1_3_p = new numeric(-1,3))->setflag(status_flags::dynallocated); + (_num_1_4_p = new numeric(-1,4))->setflag(status_flags::dynallocated); + (_num0_p = new numeric(0))->setflag(status_flags::dynallocated); + _num0_bp = _num0_p; // Cf. class ex default ctor. + (_num1_4_p = new numeric(1,4))->setflag(status_flags::dynallocated); + (_num1_3_p = new numeric(1,3))->setflag(status_flags::dynallocated); + (_num1_2_p = new numeric(1,2))->setflag(status_flags::dynallocated); + (_num1_p = new numeric(1))->setflag(status_flags::dynallocated); + (_num2_p = new numeric(2))->setflag(status_flags::dynallocated); + (_num3_p = new numeric(3))->setflag(status_flags::dynallocated); + (_num4_p = new numeric(4))->setflag(status_flags::dynallocated); + (_num5_p = new numeric(5))->setflag(status_flags::dynallocated); + (_num6_p = new numeric(6))->setflag(status_flags::dynallocated); + (_num7_p = new numeric(7))->setflag(status_flags::dynallocated); + (_num8_p = new numeric(8))->setflag(status_flags::dynallocated); + (_num9_p = new numeric(9))->setflag(status_flags::dynallocated); + (_num10_p = new numeric(10))->setflag(status_flags::dynallocated); + (_num11_p = new numeric(11))->setflag(status_flags::dynallocated); + (_num12_p = new numeric(12))->setflag(status_flags::dynallocated); + (_num15_p = new numeric(15))->setflag(status_flags::dynallocated); + (_num18_p = new numeric(18))->setflag(status_flags::dynallocated); + (_num20_p = new numeric(20))->setflag(status_flags::dynallocated); + (_num24_p = new numeric(24))->setflag(status_flags::dynallocated); + (_num25_p = new numeric(25))->setflag(status_flags::dynallocated); + (_num30_p = new numeric(30))->setflag(status_flags::dynallocated); + (_num48_p = new numeric(48))->setflag(status_flags::dynallocated); + (_num60_p = new numeric(60))->setflag(status_flags::dynallocated); + (_num120_p = new numeric(120))->setflag(status_flags::dynallocated); + + new((void*)&_ex_120) ex(*_num_120_p); + new((void*)&_ex_60) ex(*_num_60_p); + new((void*)&_ex_48) ex(*_num_48_p); + new((void*)&_ex_30) ex(*_num_30_p); + new((void*)&_ex_25) ex(*_num_25_p); + new((void*)&_ex_24) ex(*_num_24_p); + new((void*)&_ex_20) ex(*_num_20_p); + new((void*)&_ex_18) ex(*_num_18_p); + new((void*)&_ex_15) ex(*_num_15_p); + new((void*)&_ex_12) ex(*_num_12_p); + new((void*)&_ex_11) ex(*_num_11_p); + new((void*)&_ex_10) ex(*_num_10_p); + new((void*)&_ex_9) ex(*_num_9_p); + new((void*)&_ex_8) ex(*_num_8_p); + new((void*)&_ex_7) ex(*_num_7_p); + new((void*)&_ex_6) ex(*_num_6_p); + new((void*)&_ex_5) ex(*_num_5_p); + new((void*)&_ex_4) ex(*_num_4_p); + new((void*)&_ex_3) ex(*_num_3_p); + new((void*)&_ex_2) ex(*_num_2_p); + new((void*)&_ex_1) ex(*_num_1_p); + new((void*)&_ex_1_2) ex(*_num_1_2_p); + new((void*)&_ex_1_3) ex(*_num_1_3_p); + new((void*)&_ex_1_4) ex(*_num_1_4_p); + new((void*)&_ex0) ex(*_num0_p); + new((void*)&_ex1_4) ex(*_num1_4_p); + new((void*)&_ex1_3) ex(*_num1_3_p); + new((void*)&_ex1_2) ex(*_num1_2_p); + new((void*)&_ex1) ex(*_num1_p); + new((void*)&_ex2) ex(*_num2_p); + new((void*)&_ex3) ex(*_num3_p); + new((void*)&_ex4) ex(*_num4_p); + new((void*)&_ex5) ex(*_num5_p); + new((void*)&_ex6) ex(*_num6_p); + new((void*)&_ex7) ex(*_num7_p); + new((void*)&_ex8) ex(*_num8_p); + new((void*)&_ex9) ex(*_num9_p); + new((void*)&_ex10) ex(*_num10_p); + new((void*)&_ex11) ex(*_num11_p); + new((void*)&_ex12) ex(*_num12_p); + new((void*)&_ex15) ex(*_num15_p); + new((void*)&_ex18) ex(*_num18_p); + new((void*)&_ex20) ex(*_num20_p); + new((void*)&_ex24) ex(*_num24_p); + new((void*)&_ex25) ex(*_num25_p); + new((void*)&_ex30) ex(*_num30_p); + new((void*)&_ex48) ex(*_num48_p); + new((void*)&_ex60) ex(*_num60_p); + new((void*)&_ex120) ex(*_num120_p); + + // Initialize print context class info (this is not strictly necessary + // but we do it anyway to make print_context_class_info::dump_hierarchy() + // output the whole hierarchy whether or not the classes are actually + // used) + print_context::get_class_info_static(); + print_dflt::get_class_info_static(); + print_latex::get_class_info_static(); + print_python::get_class_info_static(); + print_python_repr::get_class_info_static(); + print_tree::get_class_info_static(); + print_csrc::get_class_info_static(); + print_csrc_float::get_class_info_static(); + print_csrc_double::get_class_info_static(); + print_csrc_cl_N::get_class_info_static(); + } +} + + +/** Dtor of static initialization helpers. The last call to this is going + * to shut down the library, the others do nothing. */ +library_init::~library_init() +{ + if (--count==0) { + // In theory, we would have to clean up here. But since we were + // only initializing memory in the ctor and that memory is reclaimed + // anyways by the OS when the program exits, we skip this. + } } // comment skeleton for header files @@ -464,31 +399,31 @@ ex const & _ex12(void) // member functions - // default constructor, destructor, copy constructor assignment operator and helpers - // none + // default constructor, destructor, copy constructor and assignment operator + // 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 base 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 ////////// -// default constructor, destructor, copy constructor assignment operator and helpers +// default constructor, destructor, copy constructor and assignment operator ////////// // public @@ -502,7 +437,7 @@ ex const & _ex12(void) // none ////////// -// functions overriding virtual functions from bases classes +// functions overriding virtual functions from base classes ////////// // public @@ -534,6 +469,4 @@ ex const & _ex12(void) // none -#ifndef NO_GINAC_NAMESPACE } // namespace GiNaC -#endif // ndef NO_GINAC_NAMESPACE