implemented relational::subs()

[ginac.git] / ginac / relational.cpp

diff --git a/ginac/relational.cpp b/ginac/relational.cpp
index 24fd7a05a0d8c67606bf5475bf871d0706258f5f..e4696b88b4630c3bcc8837d1517fd6c3f616215d 100644 (file)
--- a/ginac/relational.cpp
+++ b/ginac/relational.cpp
@@ -188,6 +188,17 @@ ex relational::eval(int level) const
        return (new relational(lh.eval(level-1),rh.eval(level-1),o))->setflag(status_flags::dynallocated | status_flags::evaluated);
 }
 
+ex relational::subs(const lst & ls, const lst & lr, unsigned options) const
+{
+       const ex & subsed_lh = lh.subs(ls, lr, options);
+       const ex & subsed_rh = rh.subs(ls, lr, options);
+
+       if (!are_ex_trivially_equal(lh, subsed_lh) || !are_ex_trivially_equal(rh, subsed_rh))
+               return relational(subsed_lh, subsed_rh, o).basic::subs(ls, lr, options);
+       else
+               return basic::subs(ls, lr, options);
+}
+
 ex relational::eval_ncmul(const exvector & v) const
 {
        return lh.eval_ncmul(v);