utils.cpp

Go to the documentation of this file.

 
/*
 *  GiNaC Copyright (C) 1999-2026 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
 *  the Free Software Foundation; either version 2 of the License, or
 *  (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see <https://www.gnu.org/licenses/>.
 */
 
#include "ex.h"
#include "numeric.h"
#include "utils.h"
#include "version.h"
 
namespace GiNaC {
 
/* 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;
 
 
pole_error::pole_error(const std::string& what_arg, int degree)
    : domain_error(what_arg), deg(degree) { }
 
int pole_error::degree() const
{
    return deg;
}
 
unsigned log2(unsigned n)
{
    unsigned k;
    for (k = 0; n > 1; n >>= 1)
        ++k;
    return k;
}
 
const numeric
multinomial_coefficient(const std::vector<unsigned> & p)
{
    numeric n = 0, d = 1;
    for (auto & it : p) {
        n = n.add(numeric(it));
        d = d.mul(factorial(numeric(it)));
    }
    return factorial(n).div(d);
}
 
 
int library_init::count = 0;
 
library_init::library_init()
{
    if (count++==0) {
        _num_120_p = (const numeric *)&dynallocate<numeric>(-120);
        _num_60_p = (const numeric *)&dynallocate<numeric>(-60);
        _num_48_p = (const numeric *)&dynallocate<numeric>(-48);
        _num_30_p = (const numeric *)&dynallocate<numeric>(-30);
        _num_25_p = (const numeric *)&dynallocate<numeric>(-25);
        _num_24_p = (const numeric *)&dynallocate<numeric>(-24);
        _num_20_p = (const numeric *)&dynallocate<numeric>(-20);
        _num_18_p = (const numeric *)&dynallocate<numeric>(-18);
        _num_15_p = (const numeric *)&dynallocate<numeric>(-15);
        _num_12_p = (const numeric *)&dynallocate<numeric>(-12);
        _num_11_p = (const numeric *)&dynallocate<numeric>(-11);
        _num_10_p = (const numeric *)&dynallocate<numeric>(-10);
        _num_9_p = (const numeric *)&dynallocate<numeric>(-9);
        _num_8_p = (const numeric *)&dynallocate<numeric>(-8);
        _num_7_p = (const numeric *)&dynallocate<numeric>(-7);
        _num_6_p = (const numeric *)&dynallocate<numeric>(-6);
        _num_5_p = (const numeric *)&dynallocate<numeric>(-5);
        _num_4_p = (const numeric *)&dynallocate<numeric>(-4);
        _num_3_p = (const numeric *)&dynallocate<numeric>(-3);
        _num_2_p = (const numeric *)&dynallocate<numeric>(-2);
        _num_1_p = (const numeric *)&dynallocate<numeric>(-1);
        _num_1_2_p = (const numeric *)&dynallocate<numeric>(-1,2);
        _num_1_3_p = (const numeric *)&dynallocate<numeric>(-1,3);
        _num_1_4_p = (const numeric *)&dynallocate<numeric>(-1,4);
        _num0_p = (const numeric *)&dynallocate<numeric>(0);
        _num0_bp  = _num0_p;  // Cf. class ex default ctor.
        _num1_4_p = (const numeric *)&dynallocate<numeric>(1,4);
        _num1_3_p = (const numeric *)&dynallocate<numeric>(1,3);
        _num1_2_p = (const numeric *)&dynallocate<numeric>(1,2);
        _num1_p = (const numeric *)&dynallocate<numeric>(1);
        _num2_p = (const numeric *)&dynallocate<numeric>(2);
        _num3_p = (const numeric *)&dynallocate<numeric>(3);
        _num4_p = (const numeric *)&dynallocate<numeric>(4);
        _num5_p = (const numeric *)&dynallocate<numeric>(5);
        _num6_p = (const numeric *)&dynallocate<numeric>(6);
        _num7_p = (const numeric *)&dynallocate<numeric>(7);
        _num8_p = (const numeric *)&dynallocate<numeric>(8);
        _num9_p = (const numeric *)&dynallocate<numeric>(9);
        _num10_p = (const numeric *)&dynallocate<numeric>(10);
        _num11_p = (const numeric *)&dynallocate<numeric>(11);
        _num12_p = (const numeric *)&dynallocate<numeric>(12);
        _num15_p = (const numeric *)&dynallocate<numeric>(15);
        _num18_p = (const numeric *)&dynallocate<numeric>(18);
        _num20_p = (const numeric *)&dynallocate<numeric>(20);
        _num24_p = (const numeric *)&dynallocate<numeric>(24);
        _num25_p = (const numeric *)&dynallocate<numeric>(25);
        _num30_p = (const numeric *)&dynallocate<numeric>(30);
        _num48_p = (const numeric *)&dynallocate<numeric>(48);
        _num60_p = (const numeric *)&dynallocate<numeric>(60);
        _num120_p = (const numeric *)&dynallocate<numeric>(120);
 
        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();
    }
}
 
 
library_init::~library_init()
{
    if (--count==0) {
        // It's really necessary to clean up, since the program
        // lifetime might not be the same as libginac.{so,dll} one
        // (e.g. consider // dlopen/dlsym/dlclose sequence).
        // Let the ex dtors care for deleting the numerics!
        _ex120.~ex();
        _ex_120.~ex();
        _ex60.~ex();
        _ex_60.~ex();
        _ex48.~ex();
        _ex_48.~ex();
        _ex30.~ex();
        _ex_30.~ex();
        _ex25.~ex();
        _ex_25.~ex();
        _ex24.~ex();
        _ex_24.~ex();
        _ex20.~ex();
        _ex_20.~ex();
        _ex18.~ex();
        _ex_18.~ex();
        _ex15.~ex();
        _ex_15.~ex();
        _ex12.~ex();
        _ex_12.~ex();
        _ex11.~ex();
        _ex_11.~ex();
        _ex10.~ex();
        _ex_10.~ex();
        _ex9.~ex();
        _ex_9.~ex();
        _ex8.~ex();
        _ex_8.~ex();
        _ex7.~ex();
        _ex_7.~ex();
        _ex6.~ex();
        _ex_6.~ex();
        _ex5.~ex();
        _ex_5.~ex();
        _ex4.~ex();
        _ex_4.~ex();
        _ex3.~ex();
        _ex_3.~ex();
        _ex2.~ex();
        _ex_2.~ex();
        _ex1.~ex();
        _ex_1.~ex();
        _ex1_2.~ex();
        _ex_1_2.~ex();
        _ex1_3.~ex();
        _ex_1_3.~ex();
        _ex1_4.~ex();
        _ex_1_4.~ex();
        _ex0.~ex();
    }
}
 
void library_init::init_unarchivers() { }
 
 
// Flyweight chest of numbers is re-initialized here. Note that this works
// because the numeric* have been dynallocated by the library_init ctor before
// (with the first module that has a static library_init object), so the
// assignments here only increment their refcounts.
 
// static numeric -120
const numeric *_num_120_p;
const ex _ex_120 = ex(*_num_120_p);
 
// static numeric -60
const numeric *_num_60_p;
const ex _ex_60 = ex(*_num_60_p);
 
// static numeric -48
const numeric *_num_48_p;
const ex _ex_48 = ex(*_num_48_p);
 
// static numeric -30
const numeric *_num_30_p;
const ex _ex_30 = ex(*_num_30_p);
 
// static numeric -25
const numeric *_num_25_p;
const ex _ex_25 = ex(*_num_25_p);
 
// static numeric -24
const numeric *_num_24_p;
const ex _ex_24 = ex(*_num_24_p);
 
// static numeric -20
const numeric *_num_20_p;
const ex _ex_20 = ex(*_num_20_p);
 
// static numeric -18
const numeric *_num_18_p;
const ex _ex_18 = ex(*_num_18_p);
 
// static numeric -15
const numeric *_num_15_p;
const ex _ex_15 = ex(*_num_15_p);
 
// static numeric -12
const numeric *_num_12_p;
const ex _ex_12 = ex(*_num_12_p);
 
// static numeric -11
const numeric *_num_11_p;
const ex _ex_11 = ex(*_num_11_p);
 
// static numeric -10
const numeric *_num_10_p;
const ex _ex_10 = ex(*_num_10_p);
 
// static numeric -9
const numeric *_num_9_p;
const ex _ex_9 = ex(*_num_9_p);
 
// static numeric -8
const numeric *_num_8_p;
const ex _ex_8 = ex(*_num_8_p);
 
// static numeric -7
const numeric *_num_7_p;
const ex _ex_7 = ex(*_num_7_p);
 
// static numeric -6
const numeric *_num_6_p;
const ex _ex_6 = ex(*_num_6_p);
 
// static numeric -5
const numeric *_num_5_p;
const ex _ex_5 = ex(*_num_5_p);
 
// static numeric -4
const numeric *_num_4_p;
const ex _ex_4 = ex(*_num_4_p);
 
// static numeric -3
const numeric *_num_3_p;
const ex _ex_3 = ex(*_num_3_p);
 
// static numeric -2
const numeric *_num_2_p;
const ex _ex_2 = ex(*_num_2_p);
 
// static numeric -1
const numeric *_num_1_p;
const ex _ex_1 = ex(*_num_1_p);
 
// static numeric -1/2
const numeric *_num_1_2_p;
const ex _ex_1_2= ex(*_num_1_2_p);
 
// static numeric -1/3
const numeric *_num_1_3_p;
const ex _ex_1_3= ex(*_num_1_3_p);
 
// static numeric -1/4
const numeric *_num_1_4_p;
const ex _ex_1_4= ex(*_num_1_4_p);
 
// static numeric 0
const numeric *_num0_p;
const basic *_num0_bp;
const ex _ex0 = ex(*_num0_p);
 
// static numeric 1/4
const numeric *_num1_4_p;
const ex _ex1_4 = ex(*_num1_4_p);
 
// static numeric 1/3
const numeric *_num1_3_p;
const ex _ex1_3 = ex(*_num1_3_p);
 
// static numeric 1/2
const numeric *_num1_2_p;
const ex _ex1_2 = ex(*_num1_2_p);
 
// static numeric 1
const numeric *_num1_p;
const ex _ex1 = ex(*_num1_p);
 
// static numeric 2
const numeric *_num2_p;
const ex _ex2 = ex(*_num2_p);
 
// static numeric 3
const numeric *_num3_p;
const ex _ex3 = ex(*_num3_p);
 
// static numeric 4
const numeric *_num4_p;
const ex _ex4 = ex(*_num4_p);
 
// static numeric 5
const numeric *_num5_p;
const ex _ex5 = ex(*_num5_p);
 
// static numeric 6
const numeric *_num6_p;
const ex _ex6 = ex(*_num6_p);
 
// static numeric 7
const numeric *_num7_p;
const ex _ex7 = ex(*_num7_p);
 
// static numeric 8
const numeric *_num8_p;
const ex _ex8 = ex(*_num8_p);
 
// static numeric 9
const numeric *_num9_p;
const ex _ex9 = ex(*_num9_p);
 
// static numeric 10
const numeric *_num10_p;
const ex _ex10 = ex(*_num10_p);
 
// static numeric 11
const numeric *_num11_p;
const ex _ex11 = ex(*_num11_p);
 
// static numeric 12
const numeric *_num12_p;
const ex _ex12 = ex(*_num12_p);
 
// static numeric 15
const numeric *_num15_p;
const ex _ex15 = ex(*_num15_p);
 
// static numeric 18
const numeric *_num18_p;
const ex _ex18 = ex(*_num18_p);
 
// static numeric 20
const numeric *_num20_p;
const ex _ex20 = ex(*_num20_p);
 
// static numeric 24
const numeric *_num24_p;
const ex _ex24 = ex(*_num24_p);
 
// static numeric 25
const numeric *_num25_p;
const ex _ex25 = ex(*_num25_p);
 
// static numeric 30
const numeric *_num30_p;
const ex _ex30 = ex(*_num30_p);
 
// static numeric 48
const numeric *_num48_p;
const ex _ex48 = ex(*_num48_p);
 
// static numeric 60
const numeric *_num60_p;
const ex _ex60 = ex(*_num60_p);
 
// static numeric 120
const numeric *_num120_p;
const ex _ex120 = ex(*_num120_p);
 
// comment skeleton for header files
 
 
// member functions
 
    // default constructor, destructor, copy constructor and assignment operator
    // none
 
    // other constructors
    // 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
 
// member variables
// none
    
 
 
// comment skeleton for implementation files
 
 
// default constructor, destructor, copy constructor and assignment operator
 
// public
// protected
 
// other constructors
 
// public
// none
 
// functions overriding virtual functions from base classes
 
// public
// protected
// none
 
// new virtual functions which can be overridden by derived classes
 
// public
// protected
// none
 
// non-virtual functions in this class
 
// public
// protected
// none
 
// static member variables
 
// protected
// private
// none
 
 
} // namespace GiNaC