* Interface to GiNaC's symbolic objects. */
/*
- * GiNaC Copyright (C) 1999-2003 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
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, write to the Free Software
- * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
*/
#ifndef __GINAC_SYMBOL_H__
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)
-// types
-
- /** Symbols as keys to expressions - only for ginsh. */
- class assigned_ex_info {
- friend class ptr<assigned_ex_info>;
- public:
- assigned_ex_info(); ///< Default ctor
- bool is_assigned; ///< True if there is an expression assigned
- ex assigned_expression; ///< The actual expression
- private:
- size_t refcount; ///< Reference counter, managed by ptr<assigned_ex_info>
- };
+ friend class realsymbol;
+ friend class possymbol;
// member functions
// other constructors
public:
- explicit symbol(const std::string & initname);
- symbol(const std::string & initname, const std::string & texname);
- symbol(const std::string & initname, unsigned rt, unsigned rtt);
- symbol(const std::string & initname, const std::string & texname, unsigned rt, unsigned rtt);
+ explicit symbol(const std::string & initname, unsigned domain = domain::complex);
+ symbol(const std::string & initname, const std::string & texname, unsigned domain = domain::complex);
+ symbol(const std::string & initname, const std::string & texname, unsigned rt, unsigned domain);
+ symbol(const std::string & initname, unsigned rt, tinfo_t rtt, unsigned domain = domain::complex);
+ symbol(const std::string & initname, const std::string & texname, unsigned rt, tinfo_t rtt, unsigned domain = domain::complex);
// 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
ex normal(exmap & repl, exmap & rev_lookup, int level = 0) const;
- ex to_rational(lst &repl_lst) const;
- ex to_polynomial(lst &repl_lst) const;
+ ex to_rational(exmap & repl) const;
+ ex to_polynomial(exmap & repl) const;
unsigned return_type() const { return ret_type; }
- unsigned return_type_tinfo() const { return ret_type_tinfo; }
+ 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;
// 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; }
protected:
void do_print(const print_context & c, unsigned level) const;
void do_print_latex(const print_latex & c, unsigned level) const;
// member variables
protected:
- ptr<assigned_ex_info> 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
+ unsigned domain; ///< domain of symbol, complex (default) or real
unsigned ret_type; ///< value returned by return_type()
- unsigned ret_type_tinfo; ///< value returned by return_type_tinfo()
+ tinfo_t ret_type_tinfo; ///< value returned by return_type_tinfo()
private:
static unsigned next_serial;
};
+/** Specialization of symbol to real domain */
+class realsymbol : public symbol
+{
+ // constructors
+public:
+ realsymbol();
+ explicit realsymbol(const std::string & initname, unsigned domain = domain::real);
+ realsymbol(const std::string & initname, const std::string & texname, unsigned domain = domain::real);
+ realsymbol(const std::string & initname, unsigned rt, tinfo_t rtt, unsigned domain = domain::real);
+ realsymbol(const std::string & initname, const std::string & texname, unsigned rt, tinfo_t rtt, unsigned domain = domain::real);
+};
+
+
+/** 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<symbol>(obj) for symbol objects. */
-template<> inline bool is_exactly_a<symbol>(const basic & obj)
+/** Specialization of is_exactly_a<realsymbol>(obj) for realsymbol objects. */
+template<> inline bool is_exactly_a<realsymbol>(const basic & obj)
{
- return obj.tinfo()==TINFO_symbol;
+ if (obj.tinfo() != &symbol::tinfo_static)
+ return false;
+ unsigned domain = static_cast<const symbol &>(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<possymbol>(obj) for possymbol objects. */
+template<> inline bool is_exactly_a<possymbol>(const basic & obj)
+{
+ if (obj.tinfo() != &symbol::tinfo_static)
+ return false;
+ unsigned domain = static_cast<const symbol &>(obj).get_domain();
+ return domain == domain::positive;
+}
} // namespace GiNaC
git://www.ginac.de / ginac.git / blobdiff