Linear Operator: add pow(), to allow `A**(-1)`

* also allows other integer powers
1 job for feature/operator_overloading_sugar in 5 minutes and 35 seconds