Add new feature: power operator for matrix.

The implementation is using Addition-chain exponentiation and Exponentiation by squaring algorithms.
For Addition-chain exponentiation, we need some precomputed sequences that are in additionchains.h, so
exponentiation by squaring will be used when the exponent is greater than 2048.