有没有办法加快从上一次迭代更新其值的双循环?
在代码中:
def calc(N, m):
x = 1.0
y = 2.0
container = np.zeros((N, 2))
for i in range(N):
for j in range(m):
x=np.random.gamma(3,1.0/(y*y+4))
y=np.random.normal(1.0/(x+1),1.0/sqrt(x+1))
container[i, 0] = x
container[i, 1] = y
return container
calc(10, 5)
如您所见,内循环正在更新变量 x 和 y,而外循环每次都以不同的 x 值开始。我不认为这是可矢量化的,但也许还有其他可能的改进。
谢谢!