from sympy.tensor.tensor import TensorIndexType, tensor_indices
from sympy.tensor.tensor import tensorhead, TensorManager
from sympy.tensor.tensor import get_symmetric_group_sgs
from sympy import LeviCivita
Lorentz = TensorIndexType('Lorentz', dim=3,
eps_dim=3,
metric=False,
dummy_fmt='L')
Lorentz.data = [1,1,1]
Lorentz.epsilon.data = [[[LeviCivita(i,j,k) for k in range(3)]
for j in range(3)]
for i in range(3)]
i,j,k,l,m, n = tensor_indices('i,j,k,l,m,n', Lorentz)
eps = Lorentz.epsilon
delta = Lorentz.delta