epsilon tensor contractions evaluate to metric tensors instead of deltas,

[ginac.git] / ginac / relational.h

diff --git a/ginac/relational.h b/ginac/relational.h
index c3b223b5eb0f494f722ffc25f20a10ce72de5431..40c799bf5ce580247a78949a6fd9bf6b557b51cc 100644 (file)
--- a/ginac/relational.h
+++ b/ginac/relational.h
@@ -49,7 +49,7 @@ public:
 public:
        relational(const ex & lhs, const ex & rhs, operators oper=equal);
        
-       // functions overriding virtual functions from bases classes
+       // functions overriding virtual functions from base classes
 public:
        void print(const print_context & c, unsigned level = 0) const;
        unsigned precedence(void) const {return 20;}
@@ -57,10 +57,10 @@ public:
        unsigned nops() const;
        ex & let_op(int i);
        ex eval(int level=0) const;
-       ex evalf(int level=0) const;
-       ex normal(lst &sym_lst, lst &repl_lst, int level=0) const;
        ex simplify_ncmul(const exvector & v) const;
+
 protected:
+       bool match_same_type(const basic & other) const;
        unsigned return_type(void) const;
        unsigned return_type_tinfo(void) const;
 
@@ -82,11 +82,20 @@ protected:
 };
 
 // utility functions
+
+/** Return the relational object handled by an ex.  Deprecated: use ex_to<relational>().
+ *  This is unsafe: you need to check the type first. */
 inline const relational &ex_to_relational(const ex &e)
 {
        return static_cast<const relational &>(*e.bp);
 }
 
+/** Specialization of is_exactly_a<relational>(obj) for relational objects. */
+template<> inline bool is_exactly_a<relational>(const basic & obj)
+{
+       return obj.tinfo()==TINFO_relational;
+}
+
 } // namespace GiNaC
 
 #endif // ndef __GINAC_RELATIONAL_H__