用Cython遍历Numpy数组

Question

设a和b为两个长度为1024的numpy.float数组，定义为

cdef numpy.ndarray a
cdef numpy.ndarray b

我注意到：

cdef int i
for i in range(1024):
    b[i] += a[i]

比：

b += a 

为什么？

我真的需要能够手动循环遍历数组。

Answer 1

如果您告诉Cython a和b的数据类型和维数，则差异会更小：

cdef numpy.ndarray[np.float64_t, ndim=1] a, b

尽管差异会更小，但您不会击败b += a因为这是使用NumPy的SIMD增强功能（该功能将取决于您的CPU是否支持SIMD）。