Numpy 3D阵列矩阵乘法函数

Question

假设我有一个形状（m，n，n）和维数为（m，n）的矢量C的ndarray。我需要用以下方法乘以这两个：

result = np.empty(m,n) 
for i in range(m): 
    result[i] = W[i] @ C[i] 

如何以矢量化的方式执行此操作，而无需循环和全部？

Answer 1

因为，你需要保持来自W和C第一轴一致，从他们与矩阵乘法失去最后的轴做的，我会建议使用np.einsum一个非常有效的方法，像这样 -

np.einsum('ijk,ik->ij',W,C）

np.tensordot或np.dot没有保持轴对齐的功能，这是np.einsum改进之处。

Answer 2

np.dot应该是你的朋友。 http://docs.scipy.org/doc/numpy-1.10.0/reference/generated/numpy.dot.html

Answer 3

它使用np.tensordot

ans=np.tensordot(W,C,axes=[2,1])[np.arange(m),:,np.arange(m)] 
assert np.all(result==ans)