Add operation ^ to Matrix class.

src/Matrix.cpp

@@ -281,6 +281,12 @@ namespace Aurora {
     Matrix operator/(Matrix &&aMatrix, const Matrix &aOther) {
         return operatorMxM_RR(&vdDivI,&vzDivI,aOther,std::forward<Matrix&&>(aMatrix));
     }
+
+    //operator ^ (pow)
+    Matrix Matrix::operator^(int times) const { return operatorMxA(&vdPowI, times, *this);}
+    Matrix operator^( Matrix &&matrix,int times) {
+        return operatorMxA(&vdPowI, times, std::forward<Matrix&&>(matrix));
+    }