The following code indicates that the Einstein sum of two 3D (2x2x2) matrices is a 4D (2x2x2x2) matrix.
$ c_{ijlm} = \\Sigma_k a_{i,j,k}b_{k,l,m} $ $ c_{1,1,