Small performance improvement

I implemented exponentiation by squaring as the new default __pow__ for spmatrix. It should be Python 3 compliant. On a 1000 x 1000 CSC matrix with 1000 entries, I saw about a 20% speedup.