X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Fcontainer.pl;h=0ef3b896d47908f8c4355b1a86637b75c9120d59;hp=4cb7f8839e0c63faaccabfefc3950a397e6cb9b1;hb=083b0f50275a536be807fa2a34c1e278098e12f5;hpb=a8507b8af1c08d9b27d98d57f95c7ca1a8671e27 diff --git a/ginac/container.pl b/ginac/container.pl index 4cb7f883..0ef3b896 100755 --- a/ginac/container.pl +++ b/ginac/container.pl @@ -52,11 +52,11 @@ if ($reserve) { if ($prepend) { $PREPEND_INTERFACE=<=0); - ASSERT(i=0); + GINAC_ASSERT(i -#include -typedef ${STLHEADER} ${STLT}; +// CINT needs to work properly with and +#include + +#include "basic.h" +#include "ex.h" + +#ifndef NO_NAMESPACE_GINAC +namespace GiNaC { +#endif // ndef NO_NAMESPACE_GINAC + +// typedef ${STLHEADER} ${STLT}; +typedef ${STLHEADER} ${STLT}; // CINT does not like ${STLHEADER}<...,default_alloc> class ${CONTAINER} : public basic { + GINAC_DECLARE_REGISTERED_CLASS(${CONTAINER}, basic) public: ${CONTAINER}(); @@ -140,28 +153,28 @@ protected: public: ${CONTAINER}(${STLT} const & s, bool discardable=0); ${CONTAINER}(${STLT} * vp); // vp will be deleted - explicit ${CONTAINER}(ex const & e1); - explicit ${CONTAINER}(ex const & e1, ex const & e2); - explicit ${CONTAINER}(ex const & e1, ex const & e2, ex const & e3); - explicit ${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4); - explicit ${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4, ex const & e5); - explicit ${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4, ex const & e5, ex const & e6); - explicit ${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4, ex const & e5, ex const & e6, - ex const & e7); - explicit ${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4, ex const & e5, ex const & e6, - ex const & e7, ex const & e8); - explicit ${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4, ex const & e5, ex const & e6, - ex const & e7, ex const & e8, ex const & e9); - explicit ${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4, ex const & e5, ex const & e6, - ex const & e7, ex const & e8, ex const & e9, - ex const &e10); + explicit ${CONTAINER}(const ex & e1); + explicit ${CONTAINER}(const ex & e1, const ex & e2); + explicit ${CONTAINER}(const ex & e1, const ex & e2, const ex & e3); + explicit ${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4); + explicit ${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4, const ex & e5); + explicit ${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4, const ex & e5, const ex & e6); + explicit ${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4, const ex & e5, const ex & e6, + const ex & e7); + explicit ${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4, const ex & e5, const ex & e6, + const ex & e7, const ex & e8); + explicit ${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4, const ex & e5, const ex & e6, + const ex & e7, const ex & e8, const ex & e9); + explicit ${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4, const ex & e5, const ex & e6, + const ex & e7, const ex & e8, const ex & e9, + const ex &e10); public: basic * duplicate() const; @@ -169,23 +182,23 @@ public: void print(ostream & os, unsigned upper_precedence=0) const; void printtree(ostream & os, unsigned indent) const; bool info(unsigned inf) const; - int nops() const; - ex & let_op(int const i); + unsigned nops() const; + ex & let_op(int i); ex expand(unsigned options=0) const; - bool has(ex const & other) const; + bool has(const ex & other) const; ex eval(int level=0) const; ex evalf(int level=0) const; ex normal(lst &sym_lst, lst &repl_lst, int level=0) const; - ex diff(symbol const & s) const; - ex subs(lst const & ls, lst const & lr) const; + ex derivative(const symbol & s) const; + ex subs(const lst & ls, const lst & lr) const; protected: - int compare_same_type(basic const & other) const; - bool is_equal_same_type(basic const & other) const; + int compare_same_type(const basic & other) const; + bool is_equal_same_type(const basic & other) const; unsigned return_type(void) const; // new virtual functions which can be overridden by derived classes public: - virtual ${CONTAINER} & append(ex const & b); + virtual ${CONTAINER} & append(const ex & b); ${PREPEND_INTERFACE} protected: virtual void printseq(ostream & os, char openbracket, char delim, @@ -199,8 +212,8 @@ protected: ${STLT} evalchildren(int level) const; ${STLT} evalfchildren(int level) const; ${STLT} normalchildren(int level) const; - ${STLT} diffchildren(symbol const & s) const; - ${STLT} * subschildren(lst const & ls, lst const & lr) const; + ${STLT} diffchildren(const symbol & s) const; + ${STLT} * subschildren(const lst & ls, const lst & lr) const; protected: ${STLT} seq; @@ -210,11 +223,22 @@ protected: // global constants extern const ${CONTAINER} some_${CONTAINER}; -extern type_info const & typeid_${CONTAINER}; +extern const type_info & typeid_${CONTAINER}; + +// utility functions +inline const ${CONTAINER} &ex_to_${CONTAINER}(const ex &e) +{ + return static_cast(*e.bp); +} -// macros +inline ${CONTAINER} &ex_to_nonconst_${CONTAINER}(const ex &e) +{ + return static_cast<${CONTAINER} &>(*e.bp); +} -#define ex_to_${CONTAINER}(X) (static_cast<${CONTAINER} const &>(*(X).bp)) +#ifndef NO_NAMESPACE_GINAC +} // namespace GiNaC +#endif // ndef NO_NAMESPACE_GINAC #endif // ndef __GINAC_${CONTAINER_UC}_H__ @@ -223,7 +247,9 @@ END_OF_INTERFACE $implementation=< #include -#include "ginac.h" +#include "${CONTAINER}.h" +#include "ex.h" +#include "archive.h" +#include "debugmsg.h" + +#ifndef NO_NAMESPACE_GINAC +namespace GiNaC { +#endif // ndef NO_NAMESPACE_GINAC + +GINAC_IMPLEMENT_REGISTERED_CLASS(${CONTAINER}, basic) ${RESERVE_IMPLEMENTATION} @@ -264,7 +299,7 @@ ${RESERVE_IMPLEMENTATION} // public -${CONTAINER}::${CONTAINER}() : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}() : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} default constructor",LOGLEVEL_CONSTRUCT); } @@ -295,14 +330,14 @@ ${CONTAINER} const & ${CONTAINER}::operator=(${CONTAINER} const & other) void ${CONTAINER}::copy(${CONTAINER} const & other) { - basic::copy(other); + inherited::copy(other); seq=other.seq; } void ${CONTAINER}::destroy(bool call_parent) { seq.clear(); - if (call_parent) basic::destroy(call_parent); + if (call_parent) inherited::destroy(call_parent); } ////////// @@ -311,7 +346,7 @@ void ${CONTAINER}::destroy(bool call_parent) // public -${CONTAINER}::${CONTAINER}(${STLT} const & s, bool discardable) : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}(${STLT} const & s, bool discardable) : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} constructor from ${STLT}", LOGLEVEL_CONSTRUCT); @@ -322,15 +357,15 @@ ${CONTAINER}::${CONTAINER}(${STLT} const & s, bool discardable) : basic(TINFO_$ } } -${CONTAINER}::${CONTAINER}(${STLT} * vp) : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}(${STLT} * vp) : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} constructor from ${STLT} *",LOGLEVEL_CONSTRUCT); - ASSERT(vp!=0); + GINAC_ASSERT(vp!=0); seq.swap(*vp); delete vp; } -${CONTAINER}::${CONTAINER}(ex const & e1) : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}(const ex & e1) : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} constructor from 1 ex", LOGLEVEL_CONSTRUCT); @@ -338,7 +373,7 @@ ${CONTAINER}::${CONTAINER}(ex const & e1) : basic(TINFO_${CONTAINER_UC}) seq.push_back(e1); } -${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2) : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}(const ex & e1, const ex & e2) : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} constructor from 2 ex", LOGLEVEL_CONSTRUCT); @@ -347,8 +382,8 @@ ${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2) : basic(TINFO_${CONTAIN seq.push_back(e2); } -${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3) - : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}(const ex & e1, const ex & e2, const ex & e3) + : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} constructor from 3 ex", LOGLEVEL_CONSTRUCT); @@ -358,8 +393,8 @@ ${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3) seq.push_back(e3); } -${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4) : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4) : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} constructor from 4 ex", LOGLEVEL_CONSTRUCT); @@ -370,8 +405,8 @@ ${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, seq.push_back(e4); } -${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4, ex const & e5) : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4, const ex & e5) : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} constructor from 5 ex", LOGLEVEL_CONSTRUCT); @@ -383,9 +418,9 @@ ${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, seq.push_back(e5); } -${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4, ex const & e5, ex const & e6) - : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4, const ex & e5, const ex & e6) + : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} constructor from 6 ex", LOGLEVEL_CONSTRUCT); @@ -398,9 +433,9 @@ ${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, seq.push_back(e6); } -${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4, ex const & e5, ex const & e6, - ex const & e7) : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4, const ex & e5, const ex & e6, + const ex & e7) : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} constructor from 7 ex", LOGLEVEL_CONSTRUCT); @@ -414,9 +449,9 @@ ${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, seq.push_back(e7); } -${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4, ex const & e5, ex const & e6, - ex const & e7, ex const & e8) : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4, const ex & e5, const ex & e6, + const ex & e7, const ex & e8) : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} constructor from 8 ex", LOGLEVEL_CONSTRUCT); @@ -431,10 +466,10 @@ ${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, seq.push_back(e8); } -${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4, ex const & e5, ex const & e6, - ex const & e7, ex const & e8, ex const & e9) - : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4, const ex & e5, const ex & e6, + const ex & e7, const ex & e8, const ex & e9) + : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} constructor from 9 ex", LOGLEVEL_CONSTRUCT); @@ -450,11 +485,11 @@ ${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, seq.push_back(e9); } -${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, - ex const & e4, ex const & e5, ex const & e6, - ex const & e7, ex const & e8, ex const & e9, - ex const &e10) - : basic(TINFO_${CONTAINER_UC}) +${CONTAINER}::${CONTAINER}(const ex & e1, const ex & e2, const ex & e3, + const ex & e4, const ex & e5, const ex & e6, + const ex & e7, const ex & e8, const ex & e9, + const ex &e10) + : basic(TINFO_${CONTAINER}) { debugmsg("${CONTAINER} constructor from 10 ex", LOGLEVEL_CONSTRUCT); @@ -471,6 +506,40 @@ ${CONTAINER}::${CONTAINER}(ex const & e1, ex const & e2, ex const & e3, seq.push_back(e10); } +////////// +// archiving +////////// + +/** Construct object from archive_node. */ +${CONTAINER}::${CONTAINER}(const archive_node &n, const lst &sym_lst) : inherited(n, sym_lst) +{ + debugmsg("${CONTAINER} constructor from archive_node", LOGLEVEL_CONSTRUCT); + for (unsigned int i=0; true; i++) { + ex e; + if (n.find_ex("seq", e, sym_lst, i)) + seq.push_back(e); + else + break; + } +} + +/** Unarchive the object. */ +ex ${CONTAINER}::unarchive(const archive_node &n, const lst &sym_lst) +{ + return (new ${CONTAINER}(n, sym_lst))->setflag(status_flags::dynallocated); +} + +/** Archive the object. */ +void ${CONTAINER}::archive(archive_node &n) const +{ + inherited::archive(n); + ${STLT}::const_iterator i = seq.begin(), iend = seq.end(); + while (i != iend) { + n.add_ex("seq", *i); + i++; + } +} + ////////// // functions overriding virtual functions from bases classes ////////// @@ -518,7 +587,7 @@ void ${CONTAINER}::printtree(ostream & os, unsigned indent) const // ${CONTAINER}::info() will be implemented by user elsewhere"; -int ${CONTAINER}::nops() const +unsigned ${CONTAINER}::nops() const { return seq.size(); } @@ -538,9 +607,9 @@ ex ${CONTAINER}::expand(unsigned options) const // a ${CONTAINER} 'has' an expression if it is this expression itself or a child 'has' it -bool ${CONTAINER}::has(ex const & other) const +bool ${CONTAINER}::has(const ex & other) const { - ASSERT(other.bp!=0); + GINAC_ASSERT(other.bp!=0); if (is_equal(*other.bp)) return true; for (${STLT}::const_iterator it=seq.begin(); it!=seq.end(); ++it) { if ((*it).has(other)) return true; @@ -570,12 +639,12 @@ ex ${CONTAINER}::normal(lst &sym_lst, lst &repl_lst, int level) const return n.bp->basic::normal(sym_lst,repl_lst,level); } -ex ${CONTAINER}::diff(symbol const & s) const +ex ${CONTAINER}::derivative(const symbol & s) const { return this${CONTAINER}(diffchildren(s)); } -ex ${CONTAINER}::subs(lst const & ls, lst const & lr) const +ex ${CONTAINER}::subs(const lst & ls, const lst & lr) const { ${STLT} * vp=subschildren(ls,lr); if (vp==0) { @@ -586,9 +655,9 @@ ex ${CONTAINER}::subs(lst const & ls, lst const & lr) const // protected -int ${CONTAINER}::compare_same_type(basic const & other) const +int ${CONTAINER}::compare_same_type(const basic & other) const { - ASSERT(is_of_type(other,${CONTAINER})); + GINAC_ASSERT(is_of_type(other,${CONTAINER})); ${CONTAINER} const & o=static_cast<${CONTAINER} const &> (const_cast(other)); int cmpval; @@ -607,9 +676,9 @@ int ${CONTAINER}::compare_same_type(basic const & other) const return 1; } -bool ${CONTAINER}::is_equal_same_type(basic const & other) const +bool ${CONTAINER}::is_equal_same_type(const basic & other) const { - ASSERT(is_of_type(other,${CONTAINER})); + GINAC_ASSERT(is_of_type(other,${CONTAINER})); ${CONTAINER} const & o=static_cast<${CONTAINER} const &> (const_cast(other)); if (seq.size()!=o.seq.size()) return false; @@ -635,7 +704,7 @@ unsigned ${CONTAINER}::return_type(void) const // public -${CONTAINER} & ${CONTAINER}::append(ex const & b) +${CONTAINER} & ${CONTAINER}::append(const ex & b) { ensure_if_modifiable(); seq.push_back(b); @@ -757,7 +826,7 @@ ${STLT} ${CONTAINER}::normalchildren(int level) const return s; } -${STLT} ${CONTAINER}::diffchildren(symbol const & y) const +${STLT} ${CONTAINER}::diffchildren(const symbol & y) const { ${STLT} s; RESERVE(s,seq.size()); @@ -768,7 +837,7 @@ ${STLT} ${CONTAINER}::diffchildren(symbol const & y) const } /* obsolete subschildren -${STLT} ${CONTAINER}::subschildren(lst const & ls, lst const & lr) const +${STLT} ${CONTAINER}::subschildren(const lst & ls, const lst & lr) const { ${STLT} s; RESERVE(s,seq.size()); @@ -779,7 +848,7 @@ ${STLT} ${CONTAINER}::subschildren(lst const & ls, lst const & lr) const } */ -${STLT} * ${CONTAINER}::subschildren(lst const & ls, lst const & lr) const +${STLT} * ${CONTAINER}::subschildren(const lst & ls, const lst & lr) const { // returns a NULL pointer if nothing had to be substituted // returns a pointer to a newly created epvector otherwise @@ -788,7 +857,7 @@ ${STLT} * ${CONTAINER}::subschildren(lst const & ls, lst const & lr) const ${STLT}::const_iterator last=seq.end(); ${STLT}::const_iterator cit=seq.begin(); while (cit!=last) { - ex const & subsed_ex=(*cit).subs(ls,lr); + const ex & subsed_ex=(*cit).subs(ls,lr); if (!are_ex_trivially_equal(*cit,subsed_ex)) { // something changed, copy seq, subs and return it @@ -830,7 +899,11 @@ unsigned ${CONTAINER}::precedence=10; ////////// const ${CONTAINER} some_${CONTAINER}; -type_info const & typeid_${CONTAINER}=typeid(some_${CONTAINER}); +const type_info & typeid_${CONTAINER}=typeid(some_${CONTAINER}); + +#ifndef NO_NAMESPACE_GINAC +} // namespace GiNaC +#endif // ndef NO_NAMESPACE_GINAC END_OF_IMPLEMENTATION