我通过双for
循环计算元素如下。
N,l=20,10
a=np.random.rand(N,l)
b=np.random.rand(N,l)
r=np.zeros((N,N,l))
for i in range(N):
for j in range(N):
r[i,j]=a[i]*a[j]*(b[i]-b[j])-a[i]/a[j]
题:
i
not equalsj
,这意味着将对角元素保留为零。我也可以通过矢量化来做到这一点吗?您可以广播所有算术并删除loop.s
r2 = (a[:,None]*a) * (b[:,None]-b) - (a[:,None]/a)
# Verify the correctness
np.array_equal(r, r2)
# True
最后,要将对角线设置为零,请使用就地分配
r2[(np.arange(N),)*2] = 0
或者,numpy.fill_diagonal
,也填补就地。
np.fill_diagonal(r2, 0)