* Improved lsolve() of systems containing non-numeric coefficients.

[ginac.git] / ginac / symbol.h

diff --git a/ginac/symbol.h b/ginac/symbol.h
index 9c7e526d6dca881b89c35a49e39fff2feb89c27b..c9c7f91a8a0299bfee3395903a4dcd45f42a51a3 100644 (file)
--- 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-2007 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
@@ -41,6 +41,7 @@ class symbol : public basic
        GINAC_DECLARE_REGISTERED_CLASS(symbol, basic)
 
        friend class realsymbol;
+       friend class possymbol;
 
 // types
        
@@ -127,12 +128,37 @@ 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<realsymbol>(obj) for realsymbol objects. */
 template<> inline bool is_exactly_a<realsymbol>(const basic & obj)
 {
-       return (obj.tinfo() == &symbol::tinfo_static) && (static_cast<const symbol &>(obj).get_domain() == domain::real);
+       if (obj.tinfo() != &symbol::tinfo_static)
+               return false;
+       unsigned domain = static_cast<const symbol &>(obj).get_domain();
+       return domain==domain::real || domain==domain::positive;
+}
+
+/** 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;
 }
 
 // wrapper functions around member functions