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.
Edited by Gordian Edenhofer