我有以下想要加速的代码:
import numpy as np
def do_calc(a, b, d, sum_)->None:
idx = np.ndindex(a.shape[:1])
sum_[:] = 0
for i in idx:
sum_[0:d[i[0]]] -= a[i, 0:d[i[0]]][0]
sum_[d[i[0]]] += b[i, d[i[0]]][0]
return sum_
a=np.ones((5,10), order='F')
b=np.ones((5,10), order='F')*10
sum_ = np.empty(a.shape[1])
d = [8,5,2,7,1]
do_calc(a, b, d, sum_)
因此 a 和 b 是静态的 5x10 数组,我想将它们加起来形成sum_
一个长度为 1d 的数组,如下所示:
对于每一行
i
,b
我只想包含索引d[i]
,也就是b[i, d[i]]
将进入位置sum_[d[i]]
对于每一行
i
,a
我想包含所有到索引的值d[i]
,也就是说a[i, 0:d[i]]
,将进入位置sum_[0:d[i]]
在这个例子中,我们有:
a= array([[1., 1., 1., 1., 1., 1., 1., 1., 1., 1.],
[1., 1., 1., 1., 1., 1., 1., 1., 1., 1.],
[1., 1., 1., 1., 1., 1., 1., 1., 1., 1.],
[1., 1., 1., 1., 1., 1., 1., 1., 1., 1.],
[1., 1., 1., 1., 1., 1., 1., 1., 1., 1.]])
b= array([[10., 10., 10., 10., 10., 10., 10., 10., 10., 10.],
[10., 10., 10., 10., 10., 10., 10., 10., 10., 10.],
[10., 10., 10., 10., 10., 10., 10., 10., 10., 10.],
[10., 10., 10., 10., 10., 10., 10., 10., 10., 10.],
[10., 10., 10., 10., 10., 10., 10., 10., 10., 10.]])
d= [8, 5, 2, 7, 1]
sum_= array([-5., 6., 7., -3., -3., 8., -2., 9., 10., 0.])
我尝试了以下方法:
- 使用 Fortan 顺序,因为我们主要复制最后一个索引
sum_
传递给 do_calc 函数,这样就不必每次都启动它。最后两行有一个大循环,我们每次都会为 d 生成新值并运行 do_calc。- 使用 ndindex 来迭代行。
我想知道是否可以更有效地完成此操作,因为代码一直在处理这一部分。而实际上 a 和 b 不是常量,因此我们无法在此示例中利用这一点。