- simplify_indexed() simplifies c*(M).i.j -> (M').i.j (c: numeric, M: matrix)

[ginac.git] / ginac / basic.cpp

diff --git a/ginac/basic.cpp b/ginac/basic.cpp
index 6415f33fde5a79456ebb97ba0e81f36139110c97..1c6fe6fa0d1bf5553cafb38e81cc8f1c991de894 100644 (file)
--- a/ginac/basic.cpp
+++ b/ginac/basic.cpp
@@ -310,6 +310,18 @@ ex basic::add_indexed(const ex & self, const ex & other) const
        return self + other;
 }
 
+/** Multiply an indexed expression with a scalar. This function is used
+ *  internally by simplify_indexed().
+ *
+ *  @param self Indexed expression; it's base object is *this
+ *  @param other Numeric value
+ *  @return product of self and other
+ *  @see ex::simplify_indexed() */
+ex basic::scalar_mul_indexed(const ex & self, const numeric & other) const
+{
+       return self * other;
+}
+
 /** Try to contract two indexed expressions that appear in the same product. 
  *  If a contraction exists, the function overwrites one or both of the
  *  expressions and returns true. Otherwise it returns false. It is