X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Fsymbol.h;h=c174ad5f5136b10e83d8dc0dc1cf526b07b142d6;hp=3b187d617bf9d71f90159fe18ab0df53549b4595;hb=1b735625b2759c1627046f9cb1baf834f4d26f5d;hpb=29b420743599f001b8bfc8b63fb6ac81aa3dfba7 diff --git a/ginac/symbol.h b/ginac/symbol.h index 3b187d61..c174ad5f 100644 --- a/ginac/symbol.h +++ b/ginac/symbol.h @@ -3,7 +3,7 @@ * Interface to GiNaC's symbolic objects. */ /* - * GiNaC Copyright (C) 1999-2006 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2008 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 @@ -31,26 +31,13 @@ namespace GiNaC { /** Basic CAS symbol. It has a name because it must know how to output itself. - * It may be assigned an expression, but this feature is only intended for - * programs like 'ginsh' that want to associate symbols with expressions. - * If you want to replace symbols by expressions in your code, you should - * use ex::subs() or use objects of class ex instead of class symbol in the - * first place. */ + */ class symbol : public basic { GINAC_DECLARE_REGISTERED_CLASS(symbol, basic) friend class realsymbol; - -// types - - /** Symbols as keys to expressions - only for ginsh. */ - class assigned_ex_info : public refcounted { - public: - assigned_ex_info() throw(); ///< Default ctor - bool is_assigned; ///< True if there is an expression assigned - ex assigned_expression; ///< The actual expression - }; + friend class possymbol; // member functions @@ -65,7 +52,7 @@ public: // functions overriding virtual functions from base classes public: bool info(unsigned inf) const; - ex eval(int level = 0) const; + ex eval(int level = 0) const { return *this; } // for performance reasons ex evalf(int level = 0) const { return *this; } // overwrites basic::evalf() for performance reasons ex series(const relational & s, int order, unsigned options = 0) const; ex subs(const exmap & m, unsigned options = 0) const { return subs_one_level(m, options); } // overwrites basic::subs() for performance reasons @@ -75,6 +62,9 @@ public: unsigned return_type() const { return ret_type; } tinfo_t return_type_tinfo() const { return ret_type_tinfo; } ex conjugate() const; + ex real_part() const; + ex imag_part() const; + bool is_polynomial(const ex & var) const; protected: ex derivative(const symbol & s) const; bool is_equal_same_type(const basic & other) const; @@ -82,8 +72,6 @@ protected: // non-virtual functions in this class public: - void assign(const ex & value); - void unassign(); void set_name(const std::string & n) { name = n; } std::string get_name() const { return name; } unsigned get_domain() const { return domain; } @@ -99,7 +87,6 @@ private: // member variables protected: - ptr asexinfop; ///< assigned expression, only for private use by ginsh unsigned serial; ///< unique serial number for comparison std::string name; ///< printname of this symbol std::string TeX_name; ///< LaTeX name of this symbol @@ -124,17 +111,38 @@ public: }; +/** Specialization of symbol to real domain */ +class possymbol : public symbol +{ + // constructors +public: + possymbol(); + explicit possymbol(const std::string & initname, unsigned domain = domain::positive); + possymbol(const std::string & initname, const std::string & texname, unsigned domain = domain::positive); + possymbol(const std::string & initname, unsigned rt, tinfo_t rtt, unsigned domain = domain::positive); + possymbol(const std::string & initname, const std::string & texname, unsigned rt, tinfo_t rtt, unsigned domain = domain::positive); +}; + + // utility functions /** Specialization of is_exactly_a(obj) for realsymbol objects. */ template<> inline bool is_exactly_a(const basic & obj) { - return (obj.tinfo() == &symbol::tinfo_static) && (static_cast(obj).get_domain() == domain::real); + if (obj.tinfo() != &symbol::tinfo_static) + return false; + unsigned domain = static_cast(obj).get_domain(); + return domain==domain::real || domain==domain::positive; } -// wrapper functions around member functions -inline void unassign(symbol & symarg) -{ symarg.unassign(); } +/** Specialization of is_exactly_a(obj) for possymbol objects. */ +template<> inline bool is_exactly_a(const basic & obj) +{ + if (obj.tinfo() != &symbol::tinfo_static) + return false; + unsigned domain = static_cast(obj).get_domain(); + return domain == domain::positive; +} } // namespace GiNaC