Add integral power routine

A real power routine depends on support for complex matrices which we don't provide yet.
Just noting something while I've thought of it. We should be careful that we use zero correctly with strided matrices. If a matrix is a view of another matrix, zero can unintentionally erase data. It's correct here.