外部产品所有维度的总和

Question

如果我要实现此功能：

我知道我可以编写这样的循环：

result = 0
for i in range(len(x)):
    for j in range(len(y)):
        result += x[i] * y[j]

但是，如果我想使用numpy来完成该怎么办？

Answer 1

您可以为此使用广播

np.sum(x * y[:, None])

Answer 2

使用np.einsum-

np.einsum('i,j->',x,y)

或者简单地求和减去然后得到标量的乘积-

x.sum()*y.sum()