X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Fsymbol.cpp;h=2f8afc006700e698779b5bd4bb933d21138b1d10;hp=f84d98cbe996bd9237f46b54012b7e7e7f1baede;hb=979b2636e20e037573c7c178a2e99fba282ed711;hpb=487e5659efe401683eee0381b0d23f967ffffc3c diff --git a/ginac/symbol.cpp b/ginac/symbol.cpp index f84d98cb..2f8afc00 100644 --- a/ginac/symbol.cpp +++ b/ginac/symbol.cpp @@ -3,7 +3,7 @@ * Implementation of GiNaC's symbolic objects. */ /* - * GiNaC Copyright (C) 1999 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2004 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,195 +25,245 @@ #include "symbol.h" #include "lst.h" +#include "archive.h" +#include "tostring.h" #include "utils.h" -#include "idx.h" +#include "inifcns.h" + +namespace GiNaC { + +GINAC_IMPLEMENT_REGISTERED_CLASS_OPT(symbol, basic, + print_func(&symbol::do_print). + print_func(&symbol::do_print_latex). + print_func(&symbol::do_print_tree). + print_func(&symbol::do_print_python_repr)) ////////// -// default constructor, destructor, copy constructor assignment operator and helpers +// default constructor ////////// -symbol::symbol() : basic(TINFO_symbol) +// symbol + +symbol::symbol() + : inherited(TINFO_symbol), asexinfop(new assigned_ex_info), serial(next_serial++), name(autoname_prefix() + ToString(serial)), TeX_name(name), ret_type(return_types::commutative), ret_type_tinfo(TINFO_symbol), domain(domain::complex) { - debugmsg("symbol default constructor",LOGLEVEL_CONSTRUCT); - serial=next_serial++; - name=autoname_prefix()+ToString(serial); - asexinfop=new assigned_ex_info; - setflag(status_flags::evaluated); + setflag(status_flags::evaluated | status_flags::expanded); } -symbol::~symbol() +// realsymbol + +realsymbol::realsymbol() +{ + domain = domain::real; +} + +////////// +// other constructors +////////// + +// public + +// symbol + +symbol::symbol(const std::string & initname, unsigned domain) + : inherited(TINFO_symbol), asexinfop(new assigned_ex_info), serial(next_serial++), name(initname), TeX_name(default_TeX_name()), ret_type(return_types::commutative), ret_type_tinfo(TINFO_symbol), domain(domain) { - debugmsg("symbol destructor",LOGLEVEL_DESTRUCT); - destroy(0); + setflag(status_flags::evaluated | status_flags::expanded); } -symbol::symbol(symbol const & other) +symbol::symbol(const std::string & initname, const std::string & texname, unsigned domain) + : inherited(TINFO_symbol), asexinfop(new assigned_ex_info), serial(next_serial++), name(initname), TeX_name(texname), ret_type(return_types::commutative), ret_type_tinfo(TINFO_symbol), domain(domain) { - debugmsg("symbol copy constructor",LOGLEVEL_CONSTRUCT); - copy(other); + setflag(status_flags::evaluated | status_flags::expanded); } -void symbol::copy(symbol const & other) +symbol::symbol(const std::string & initname, unsigned rt, unsigned rtt, unsigned domain) + : inherited(TINFO_symbol), asexinfop(new assigned_ex_info), serial(next_serial++), name(initname), TeX_name(default_TeX_name()), ret_type(rt), ret_type_tinfo(rtt), domain(domain) { - basic::copy(other); - name=other.name; - serial=other.serial; - asexinfop=other.asexinfop; - ++asexinfop->refcount; + setflag(status_flags::evaluated | status_flags::expanded); } -void symbol::destroy(bool call_parent) +symbol::symbol(const std::string & initname, const std::string & texname, unsigned rt, unsigned rtt, unsigned domain) + : inherited(TINFO_symbol), asexinfop(new assigned_ex_info), serial(next_serial++), name(initname), TeX_name(texname), ret_type(rt), ret_type_tinfo(rtt), domain(domain) { - if (--asexinfop->refcount == 0) { - delete asexinfop; - } - if (call_parent) { - basic::destroy(call_parent); - } + setflag(status_flags::evaluated | status_flags::expanded); } -// how should the following be interpreted? -// symbol x; -// symbol y; -// x=y; -// probably as: x=ex(y); +// realsymbol + +realsymbol::realsymbol(const std::string & initname, unsigned domain) + : symbol(initname, domain) { } -////////// -// other constructors -////////// +realsymbol::realsymbol(const std::string & initname, const std::string & texname, unsigned domain) + : symbol(initname, texname, domain) { } -// public +realsymbol::realsymbol(const std::string & initname, unsigned rt, unsigned rtt, unsigned domain) + : symbol(initname, rt, rtt, domain) { } -symbol::symbol(string const & initname) : basic(TINFO_symbol) -{ - debugmsg("symbol constructor from string",LOGLEVEL_CONSTRUCT); - name=initname; - serial=next_serial++; - asexinfop=new assigned_ex_info; - setflag(status_flags::evaluated); -} +realsymbol::realsymbol(const std::string & initname, const std::string & texname, unsigned rt, unsigned rtt, unsigned domain) + : symbol(initname, texname, rt, rtt, domain) { } ////////// -// functions overriding virtual functions from bases classes +// archiving ////////// -// public +/** Construct object from archive_node. */ +symbol::symbol(const archive_node &n, lst &sym_lst) + : inherited(n, sym_lst), asexinfop(new assigned_ex_info), serial(next_serial++) +{ + if (!n.find_string("name", name)) + name = autoname_prefix() + ToString(serial); + if (!n.find_string("TeXname", TeX_name)) + TeX_name = default_TeX_name(); + if (!n.find_unsigned("domain", domain)) + domain = domain::complex; + if (!n.find_unsigned("return_type", ret_type)) + ret_type = return_types::commutative; + if (!n.find_unsigned("return_type_tinfo", ret_type_tinfo)) + ret_type_tinfo = TINFO_symbol; + setflag(status_flags::evaluated | status_flags::expanded); +} -basic * symbol::duplicate() const +/** Unarchive the object. */ +ex symbol::unarchive(const archive_node &n, lst &sym_lst) { - debugmsg("symbol duplicate",LOGLEVEL_DUPLICATE); - return new symbol(*this); + ex s = (new symbol(n, sym_lst))->setflag(status_flags::dynallocated); + + // If symbol is in sym_lst, return the existing symbol + for (lst::const_iterator it = sym_lst.begin(); it != sym_lst.end(); ++it) { + if (is_a(*it) && (ex_to(*it).name == ex_to(s).name)) + return *it; + } + + // Otherwise add new symbol to list and return it + sym_lst.append(s); + return s; } -bool symbol::info(unsigned inf) const +/** Archive the object. */ +void symbol::archive(archive_node &n) const { - if (inf==info_flags::symbol) return true; - if (inf==info_flags::polynomial || inf==info_flags::integer_polynomial || inf==info_flags::rational_polynomial || inf==info_flags::rational_function) { - return true; - } else { - return basic::info(inf); - } + inherited::archive(n); + n.add_string("name", name); + if (TeX_name != default_TeX_name()) + n.add_string("TeX_name", TeX_name); + if (domain != domain::complex) + n.add_unsigned("domain", domain); + if (ret_type != return_types::commutative) + n.add_unsigned("return_type", ret_type); + if (ret_type_tinfo != TINFO_symbol) + n.add_unsigned("return_type_tinfo", ret_type_tinfo); } -ex symbol::expand(unsigned options) const +////////// +// functions overriding virtual functions from base classes +////////// + +// public + +void symbol::do_print(const print_context & c, unsigned level) const { - return this->hold(); + c.s << name; } -bool symbol::has(ex const & other) const +void symbol::do_print_latex(const print_latex & c, unsigned level) const { - if (is_equal(*other.bp)) return true; - return false; + c.s << TeX_name; } -int symbol::degree(symbol const & s) const +void symbol::do_print_tree(const print_tree & c, unsigned level) const { - return compare_same_type(s)==0 ? 1 : 0; + c.s << std::string(level, ' ') << name << " (" << class_name() << ")" << " @" << this + << ", serial=" << serial + << std::hex << ", hash=0x" << hashvalue << ", flags=0x" << flags << std::dec + << ", domain=" << domain + << std::endl; } -int symbol::ldegree(symbol const & s) const +void symbol::do_print_python_repr(const print_python_repr & c, unsigned level) const { - return compare_same_type(s)==0 ? 1 : 0; + c.s << class_name() << "('" << name; + if (TeX_name != default_TeX_name()) + c.s << "','" << TeX_name; + c.s << "')"; } -ex symbol::coeff(symbol const & s, int const n) const +bool symbol::info(unsigned inf) const { - if (compare_same_type(s)==0) { - return n==1 ? exONE() : exZERO(); - } else { - return n==0 ? *this : exZERO(); - } + if (inf == info_flags::symbol) + return true; + if (inf == info_flags::polynomial || + inf == info_flags::integer_polynomial || + inf == info_flags::cinteger_polynomial || + inf == info_flags::rational_polynomial || + inf == info_flags::crational_polynomial || + inf == info_flags::rational_function) + return true; + if (inf == info_flags::real) + return domain == domain::real; + else + return inherited::info(inf); } ex symbol::eval(int level) const { - if (level == -max_recursion_level) { - throw(std::runtime_error("max recursion level reached")); - } - - if (asexinfop->is_assigned) { - setflag(status_flags::evaluated); - if (level==1) { - return (asexinfop->assigned_expression); - } else { - return (asexinfop->assigned_expression).eval(level); - } - } else { - return this->hold(); - } + if (level == -max_recursion_level) + throw(std::runtime_error("max recursion level reached")); + + if (asexinfop->is_assigned) { + setflag(status_flags::evaluated); + if (level==1) + return (asexinfop->assigned_expression); + else + return (asexinfop->assigned_expression).eval(level); + } else { + return this->hold(); + } } -ex symbol::subs(lst const & ls, lst const & lr) const +ex symbol::conjugate() const { - ASSERT(ls.nops()==lr.nops()); -#ifdef DOASSERT - for (int i=0; idomain == domain::complex) { + return GiNaC::conjugate_function(*this).hold(); + } else { + return *this; + } } // protected -int symbol::compare_same_type(basic const & other) const +/** Implementation of ex::diff() for single differentiation of a symbol. + * It returns 1 or 0. + * + * @see ex::diff */ +ex symbol::derivative(const symbol & s) const { - ASSERT(is_of_type(other,symbol)); - const symbol *o = static_cast(&other); - if (serial==o->serial) return 0; - return serial < o->serial ? -1 : 1; + if (compare_same_type(s)) + return _ex0; + else + return _ex1; } -bool symbol::is_equal_same_type(basic const & other) const +int symbol::compare_same_type(const basic & other) const { - ASSERT(is_of_type(other,symbol)); - const symbol *o = static_cast(&other); - return serial==o->serial; + GINAC_ASSERT(is_a(other)); + const symbol *o = static_cast(&other); + if (serial==o->serial) return 0; + return serial < o->serial ? -1 : 1; } -unsigned symbol::return_type(void) const +bool symbol::is_equal_same_type(const basic & other) const { - return return_types::commutative; -} - -unsigned symbol::return_type_tinfo(void) const -{ - return tinfo_key; + GINAC_ASSERT(is_a(other)); + const symbol *o = static_cast(&other); + return serial==o->serial; } -unsigned symbol::calchash(void) const +unsigned symbol::calchash() const { - // return golden_ratio_hash(tinfo()) ^ serial; - hashvalue=golden_ratio_hash(golden_ratio_hash(0x55555555U ^ serial)); - setflag(status_flags::hash_calculated); - return hashvalue; + hashvalue = golden_ratio_hash(tinfo() ^ serial); + setflag(status_flags::hash_calculated); + return hashvalue; } ////////// @@ -228,28 +278,52 @@ unsigned symbol::calchash(void) const // public -void symbol::assign(ex const & value) +void symbol::assign(const ex & value) { - asexinfop->is_assigned=1; - asexinfop->assigned_expression=value; - clearflag(status_flags::evaluated); + asexinfop->is_assigned = true; + asexinfop->assigned_expression = value; + clearflag(status_flags::evaluated | status_flags::expanded); } -void symbol::unassign(void) +void symbol::unassign() { - if (asexinfop->is_assigned) { - asexinfop->is_assigned=0; - asexinfop->assigned_expression=exZERO(); - } - setflag(status_flags::evaluated); + if (asexinfop->is_assigned) { + asexinfop->is_assigned = false; + asexinfop->assigned_expression = _ex0; + } + setflag(status_flags::evaluated | status_flags::expanded); } // private -string & symbol::autoname_prefix(void) +/** Symbols not constructed with a string get one assigned using this + * prefix and a number. */ +std::string & symbol::autoname_prefix() +{ + static std::string *s = new std::string("symbol"); + return *s; +} + +/** Return default TeX name for symbol. This recognizes some greek letters. */ +std::string symbol::default_TeX_name() const { - static string * s=new string("symbol"); - return *s; + if (name=="alpha" || name=="beta" || name=="gamma" + || name=="delta" || name=="epsilon" || name=="varepsilon" + || name=="zeta" || name=="eta" || name=="theta" + || name=="vartheta" || name=="iota" || name=="kappa" + || name=="lambda" || name=="mu" || name=="nu" + || name=="xi" || name=="omicron" || name=="pi" + || name=="varpi" || name=="rho" || name=="varrho" + || name=="sigma" || name=="varsigma" || name=="tau" + || name=="upsilon" || name=="phi" || name=="varphi" + || name=="chi" || name=="psi" || name=="omega" + || name=="Gamma" || name=="Delta" || name=="Theta" + || name=="Lambda" || name=="Xi" || name=="Pi" + || name=="Sigma" || name=="Upsilon" || name=="Phi" + || name=="Psi" || name=="Omega") + return "\\" + name; + else + return name; } ////////// @@ -258,22 +332,15 @@ string & symbol::autoname_prefix(void) // private -unsigned symbol::next_serial=0; - -// string const symbol::autoname_prefix="symbol"; - -////////// -// global constants -////////// - -const symbol some_symbol; -type_info const & typeid_symbol=typeid(some_symbol); +unsigned symbol::next_serial = 0; ////////// // subclass assigned_ex_info ////////// /** Default ctor. Defaults to unassigned. */ -symbol::assigned_ex_info::assigned_ex_info(void) : is_assigned(0), refcount(1) +symbol::assigned_ex_info::assigned_ex_info() throw() : is_assigned(false) { } + +} // namespace GiNaC