下面的代码完全符合我的要求(这是克里格方法的一部分)。但问题是它太慢了,我想知道是否有任何选择把for-循环向下推到numpy?如果我推开numpy.sum,并在那里使用axis参数,它会加速一些,但显然这不是瓶颈。我有什么想法可以把how循环推到numpy来加速它,或者其他加速它的方法?)# n = 2116print Gin
在我代码中,我让for循环遍历多维numpy数组,并使用在每次迭代中获得的子数组执行一些操作。它看起来像这样 #do stuff using sub
现在使用sub完成的工作是完全矢量化的,所以它应该是有效的。另一方面,这个循环迭代大约~10^5次,这是一个瓶颈。我不太愿意这样做,因为do stuff using sub使用广播,切片,智能索引技巧,这些技巧用普通的C编写会很繁琐。