Synced to HEAD

[ginac.git] / ginac / symbol.h

diff --git a/ginac/symbol.h b/ginac/symbol.h
index ab2d52a5477f78ea93d2dd7e2c51f15ff63fb2ac..51164a972de73971d7f51b99d08f5a6800bbb626 100644 (file)
--- a/ginac/symbol.h
+++ b/ginac/symbol.h
@@ -3,7 +3,7 @@
  *  Interface to GiNaC's symbolic objects. */
 
 /*
- *  GiNaC Copyright (C) 1999-2003 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
@@ -40,6 +40,8 @@ class symbol : public basic
 {
        GINAC_DECLARE_REGISTERED_CLASS(symbol, basic)
 
+       friend class realsymbol;
+
 // types
        
        /** Symbols as keys to expressions - only for ginsh. */
@@ -54,10 +56,10 @@ class symbol : public basic
        
        // 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, unsigned rt, unsigned rtt, unsigned domain = domain::complex);
+       symbol(const std::string & initname, const std::string & texname, unsigned rt, unsigned rtt, unsigned domain = domain::complex);
        
        // functions overriding virtual functions from base classes
 public:
@@ -67,10 +69,11 @@ public:
        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; }
+       ex conjugate() const;
 protected:
        ex derivative(const symbol & s) const;
        bool is_equal_same_type(const basic & other) const;
@@ -82,6 +85,7 @@ public:
        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;
@@ -98,6 +102,7 @@ protected:
        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()
 private:
@@ -105,12 +110,31 @@ private:
 };
 
 
+/** 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, unsigned rtt, unsigned domain = domain::real);
+       realsymbol(const std::string & initname, const std::string & texname, unsigned rt, unsigned rtt, unsigned domain = domain::real);
+};
+
+
 // utility functions
 
 /** Specialization of is_exactly_a<symbol>(obj) for symbol objects. */
 template<> inline bool is_exactly_a<symbol>(const basic & obj)
 {
-       return obj.tinfo()==TINFO_symbol;
+       return obj.tinfo() == TINFO_symbol;
+}
+
+/** 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) && (static_cast<const symbol &>(obj).get_domain() == domain::real);
 }
 
 // wrapper functions around member functions