Introduce a numpy.einsum wrapper with derivatives

Test said operator for a couple of cases inspired by real-world scenarios.

I am pretty certain that the operator does not work for complex numbers yet but it should work just fine for real numbers.