矩阵幂和矩阵点的numpy矢量化

Question

我想计算许多二维矩阵的幂，并希望尽快完成。我找不到一个矢量化的方法。在

例如：

import numpy as np
mat_list = [np.mat(np.random.randn(100,100)) for i in range(1000)]
output = [np.linalg.matrix_power(mat_list[i], 100) for i in xrange(len(mat_list))]

我想对点乘法也一样：

^{pr2}$