numpy - numpy / numexpr python中的重复

Question

我有一个数组x=np.ones(1000)，想根据递归公式对此执行操作：

    x[i]= f(x[i-1], x[i-2])

以矢量化方式最快的方法是什么？

我知道xx= ne.evaluate(xx1*xx2)只能管理已经填充的数组。

Answer 1

我找到了一种部分方法，至少对于某些功能：

list[:,1:]= np.exp(hlist)  
list = np.cumprod(list, axis = 1)  

其中 f(x) = f(x-1) exp(a x)

我认为通用暨会很有用