2016-01-23 21 views
1

У меня есть массив x=np.ones(1000) и хотели бы выполнять операции на этой основе по рекуррентным формулам:рецидивы в NumPy/numexpr питона

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

Какой самый быстрый способ сделать это в Vectorize образом?

Я знаю, что xx= ne.evaluate(xx1*xx2) может управлять уже заполненным массивом.

+0

Не могли бы вы поделиться более подробной информацией о func 'f'? – Divakar

+0

f использует стандартную математическую функцию, такую ​​как sin, cos, exp, .. –

ответ

0

я нашел частичный путь, по крайней мере, для некоторых функций:

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

с F (х) = F (х-1) ехр (а х)

Я считаю родовое cum было бы полезно