Dear all,
The idea of renaming dummy indices when expressions are raised to a
power or when this is necessary in the case of a subsitution, is a great
idea. However, it seems rather inconsistently implemented. E.g., if you
calculate x^2 as power(x,2) then dummy indices are renamed but not if
you do it like x*x. The attached patch makes renaming work if
multiplications are constructed using operator* or by the mul::mul
(exvector) constructor. Furthermore, the function
rename_dummy_indices_uniquely was not aware of the existence of spinidx.
The patch also corrects this.
Best wishes,
Chris