Currently one must do the following ``` >>> Ric Riem(L_0, -d1, -L_0, -d3) >>> R = (Ric(d0, d1)*g(-d0, -d1)).contract_metric(g) >>> R Riem(L_0, -L_1, -L_0, L_1) ```
One cannot write ``` >>> Ric(d1, -d1) Riem(L_0, d1, -L_0, -d1) ``` I think this is a bug; the Einstein summation convention is assumed, so if the substitution indices are (d1, -d1) it should give the contracted tensor. -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To post to this group, send email to [email protected]. Visit this group at http://groups.google.com/group/sympy. For more options, visit https://groups.google.com/groups/opt_out.
