Dear all, in expression { indexed(A, i, j) * indexed(B, j, k) } GiNaC is saying that the dummy index is "j" and "i, k" are free. Is there a way to determine dummy indices based on Einstein summation notation; namely all the indices in the latter be free. All the bests, Esa