Чтобы получить скорости сравнимые с numpy.cumsum()
для больших массивов, можно Cython использовать:
cpdef void clipped_cumsum(int[:] a) nogil:
cdef Py_ssize_t i, n = a.shape[0]
for i in range(n-1):
a[i+1] += a[i]
if a[i+1] > 0:
a[i+1] = 0
Для примера, можно сохранить этот код в clipped_cumsum.pyx
файл и выполнить pip install cython
, чтобы попробовать:
#!/usr/bin/env python
import numpy as np
import pyximport # $ pip install cython
pyximport.install() # compile on-the-fly
from clipped_cumsum import clipped_cumsum
a = np.array([0, 1, -2, -3, 4, 5, -6], dtype=np.intc)
clipped_cumsum(a)
print(a) # -> [ 0 0 -2 -5 -1 0 -6]
pyximport
модуль из cython
пакета позволяет налету скомпилировать clipped_cumsum.pyx
файл, что делает доступным модуль расширения clipped_cumsum
.
Для крошечных массивов как в вашем примере, подобные изменения кода навряд ли будут полезны.
Для a.repeat(1000_000)
массива, Cython версия с явным типом (int[:]
) ~100 раз быстрее версии, которая работает с Питон int
на моей машине. Можно ещё ускорить, если выключить проверки при индексировании (boundscheck=False
, wraparound=False
), но вероятно уже существующего прироста достаточно и проверки можно оставить для будущих модификаций кода.
int[:]
наfloat[:]
замените)—на моей машине ~0.1 секунды для 25 миллионов float32 (для сравнения np.cumsum(a) в полтора раза медленнее). Если бы было больше возможности для параллелизации, то можно было бы theano попробовать, чтобы на GPU numpy выражения считать). Даже без GPU, так как я использую nogil, то вы можете разные массивы одновременно в разных потоках считать.