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