Оптимизация итеративного расчета a[i] = a[i-1]*b[i] + c[i] с помощью NumPy ⇐ Python

[Anonymous]

Я хотел бы оптимизировать вычисления с использованием рекурсивной формулы:

Код: Выделить всё

#a, b, c are arrays of shape (N, )

a[0] = c[0]

for i in range(1, N):
a[i] = a[i-1]*b[i] + c[i]

Есть ли способ вернуть результат с помощью комбинации оптимизированных методов NumPy? например cumprod, cumsum и точка?

Подробнее здесь: https://stackoverflow.com/questions/783 ... with-numpy