如何用广播计算多维数组？
我通过双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，这意味着将对角元素保留为零。我也可以通过矢量化来做到这一点吗？
展开
收起
1
条回答
写回答
写回答
-
您可以广播所有算术并删除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)
2019-07-17 23:29:43赞同 展开评论 打赏
版权声明：本文内容由阿里云实名注册用户自发贡献，版权归原作者所有，阿里云开发者社区不拥有其著作权，亦不承担相应法律责任。具体规则请查看《阿里云开发者社区用户服务协议》和《阿里云开发者社区知识产权保护指引》。如果您发现本社区中有涉嫌抄袭的内容，填写侵权投诉表单进行举报，一经查实，本社区将立刻删除涉嫌侵权内容。
相关问答
问答排行榜
最热
最新
相关课程
更多
推荐问答
