Add support for Nth order tensor computations

[GiudGiud]

The tests show I have the simple ones down correctly.

My only worry is slice. I think I might be picking a convention that is not what we assume when we use slice in the dphi calculations.