At the moment the tensor rotations are done point-wise. This is very slow now that I try with large cases. Make use of vectorized operations.