X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?a=blobdiff_plain;f=ginac%2Fbasic.h;h=242bd4569817f04412698b353d613532a0379083;hb=1f6eb80698f9fb8e464209b4167061a901decfdb;hp=a376855d6f7454747d3b09c60a2cef6b8965a227;hpb=487e5659efe401683eee0381b0d23f967ffffc3c;p=ginac.git diff --git a/ginac/basic.h b/ginac/basic.h index a376855d..242bd456 100644 --- a/ginac/basic.h +++ b/ginac/basic.h @@ -3,7 +3,7 @@ * Interface to GiNaC's ABC. */ /* - * 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 @@ -27,9 +27,16 @@ #include #include +// CINT needs to work properly with +#include + #include #include -#include +#include + +#ifndef NO_GINAC_NAMESPACE +namespace GiNaC { +#endif // ndef NO_GINAC_NAMESPACE class basic; class ex; @@ -37,7 +44,8 @@ class symbol; class lst; class numeric; -typedef vector exvector; +//typedef vector exvector; +typedef vector exvector; // CINT does not like vector<...,default_alloc> #define INLINE_BASIC_CONSTRUCTORS @@ -55,8 +63,6 @@ public: #ifdef INLINE_BASIC_CONSTRUCTORS : tinfo_key(TINFO_basic), flags(0), refcount(0) { - debugmsg("basic default constructor",LOGLEVEL_CONSTRUCT); - // nothing to do } #else ; @@ -65,9 +71,8 @@ public: virtual ~basic() #ifdef INLINE_BASIC_CONSTRUCTORS { - debugmsg("basic destructor",LOGLEVEL_DESTRUCT); destroy(0); - ASSERT((!(flags & status_flags::dynallocated))||(refcount==0)); + GINAC_ASSERT((!(flags & status_flags::dynallocated))||(refcount==0)); } #else ; @@ -76,7 +81,6 @@ public: basic(basic const & other) #ifdef INLINE_BASIC_CONSTRUCTORS { - debugmsg("basic copy constructor",LOGLEVEL_CONSTRUCT); copy(other); } #else @@ -99,8 +103,6 @@ protected: #ifdef INLINE_BASIC_CONSTRUCTORS : tinfo_key(ti), flags(0), refcount(0) { - debugmsg("basic constructor with tinfo_key",LOGLEVEL_CONSTRUCT); - // nothing to do } #else ; @@ -112,14 +114,14 @@ protected: // new virtual functions which can be overridden by derived classes public: // only const functions please (may break reference counting) virtual basic * duplicate() const; + virtual void print(ostream & os,unsigned upper_precedence=0) const; virtual void printraw(ostream & os) const; virtual void printtree(ostream & os, unsigned indent) const; - virtual void print(ostream & os,unsigned upper_precedence=0) const; virtual void printcsrc(ostream & os, unsigned type, unsigned upper_precedence=0) const; virtual void dbgprint(void) const; virtual void dbgprinttree(void) const; virtual bool info(unsigned inf) const; - virtual int nops() const; + virtual unsigned nops() const; virtual ex op(int const i) const; virtual ex & let_op(int const i); virtual ex operator[](ex const & index) const; @@ -136,7 +138,7 @@ public: // only const functions please (may break reference counting) virtual ex subs(lst const & ls, lst const & lr) const; virtual ex normal(lst &sym_lst, lst &repl_lst, int level=0) const; virtual numeric integer_content(void) const; - virtual ex smod(numeric const &xi) const; + virtual ex smod(const numeric &xi) const; virtual numeric max_coefficient(void) const; virtual exvector get_indices(void) const; virtual ex simplify_ncmul(exvector const & v) const; @@ -184,52 +186,39 @@ extern int max_recursion_level; // convenience macros +#ifndef NO_GINAC_NAMESPACE + #define is_of_type(OBJ,TYPE) \ - (dynamic_cast(const_cast(&OBJ))!=0) + (dynamic_cast(const_cast(&OBJ))!=0) #define is_exactly_of_type(OBJ,TYPE) \ - ((OBJ).tinfo()==TINFO_##TYPE) + ((OBJ).tinfo()==GiNaC::TINFO_##TYPE) #define is_ex_of_type(OBJ,TYPE) \ - (dynamic_cast(const_cast((OBJ).bp))!=0) + (dynamic_cast(const_cast((OBJ).bp))!=0) #define is_ex_exactly_of_type(OBJ,TYPE) \ - ((*(OBJ).bp).tinfo()==TINFO_##TYPE) + ((*(OBJ).bp).tinfo()==GiNaC::TINFO_##TYPE) -#define are_ex_trivially_equal(EX1,EX2) \ - ((EX1).bp==(EX2).bp) +#else // ndef NO_GINAC_NAMESPACE -// global functions +#define is_of_type(OBJ,TYPE) \ + (dynamic_cast(const_cast(&OBJ))!=0) -inline unsigned rotate_left_31(unsigned n) -{ - // clear highest bit and shift 1 bit to the left - n=(n & 0x7FFFFFFFU) << 1; +#define is_exactly_of_type(OBJ,TYPE) \ + ((OBJ).tinfo()==TINFO_##TYPE) - // overflow? clear highest bit and set lowest bit - if (n & 0x80000000U) { - n=(n & 0x7FFFFFFFU) | 0x00000001U; - } +#define is_ex_of_type(OBJ,TYPE) \ + (dynamic_cast(const_cast((OBJ).bp))!=0) - ASSERT(n<0x80000000U); +#define is_ex_exactly_of_type(OBJ,TYPE) \ + ((*(OBJ).bp).tinfo()==TINFO_##TYPE) - return n; -} +#endif // ndef NO_GINAC_NAMESPACE -inline unsigned golden_ratio_hash(unsigned n) -{ -#if 0 - // This requires ´long long´ (or an equivalent 64 bit type)---which is, - // unfortunately, not ANSI-compliant: - unsigned long long l = n * 0x4f1bbcddLL; - return (l & 0x7fffffffU) ^ (l >> 32); -#else - // This requires ´long double´ to have a mantissa of at least 64 bit--- - // which is not guaranteed by any standard: - const static long double golden_ratio=.618033988749894848204586834370; - long double m=golden_ratio*n; - return unsigned((m-int(m))*0x80000000); -#endif -} +#ifndef NO_GINAC_NAMESPACE +} // namespace GiNaC +#endif // ndef NO_GINAC_NAMESPACE #endif // ndef __GINAC_BASIC_H__ +