Numpy has the einsum function, often very efficient as it avoids allocating memory for intermediate results.

http://docs.scipy.org/doc/numpy/reference/generated/numpy.ei...