我通过双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不等于j,这意味着将对角元素保留为零。我也可以通过矢量化来做到这一点吗?
您可以广播所有算术并删除loop.s
r2 = (a[:,None]a) (b[:,None]-b) - (a[:,None]/a)
np.array_equal(r, r2)
最后,要将对角线设置为零,请使用就地分配
r2[(np.arange(N),)*2] = 0
或者,numpy.fill_diagonal也填补了原地。
np.fill_diagonal(r2, 0)
版权声明:本文内容由阿里云实名注册用户自发贡献,版权归原作者所有,阿里云开发者社区不拥有其著作权,亦不承担相应法律责任。具体规则请查看《阿里云开发者社区用户服务协议》和《阿里云开发者社区知识产权保护指引》。如果您发现本社区中有涉嫌抄袭的内容,填写侵权投诉表单进行举报,一经查实,本社区将立刻删除涉嫌侵权内容。