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)

我认为通用cum会很有用。