9
$\begingroup$

How to define a Einstein summation convention in Mathematica?

Say, how to let Mathematica know T[i,i] equals T[1,1]+T[2,2]+T[3,3]? How to define the function T[i,j]?

I know there are several Mathematica-abased programs that can do this, but still I want to know how to define my own.

$\endgroup$
1

1 Answer 1

10
$\begingroup$

For a single Function T this is pretty easy:

T[i_, i_] /; ! NumberQ[i] := Sum[T[j, j], {j, 3}]

A quick test gives:

In[7]:= T[i, i]

Out[7]= T[1, 1] + T[2, 2] + T[3, 3]

If you want to always sum over two equal indices for several functions, this might need some extra work though.

$\endgroup$
1
  • 1
    $\begingroup$ Thanks, einbandi. The answer is just what I want. $\endgroup$
    – pengfei_guo
    Nov 15, 2012 at 6:24

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.