与cython的numpy数组

时间:2013-01-06 15:19:50

标签: python numpy cython

我正在尝试将一些python代码移植到cython中,我遇到了一些小问题。

下面您将看到代码的代码段(简化示例)。

cimport numpy as np
cimport cython
@cython.boundscheck(False) # turn of bounds-checking for entire function
@cython.wraparound(False)
@cython.nonecheck(False)
def Interpolation(cells, int nmbcellsx):
    cdef np.ndarray[float,ndim=1] celle
    cdef int cellnonzero
    cdef int i,l
    for i in range(nmbcellsx):
          celle = cells[i].e
          cellnonzero = cells[i].nonzero
          for l in range(cellnonzero):
               celle[l] = celle[l] * celle[l]

我不明白为什么最内层循环没有完全转换为C代码(即最后一行,celle [l] = ...),请参阅cython -a feedback的输出:

enter image description here

我在这里缺少什么?

非常感谢。

1 个答案:

答案 0 :(得分:1)

我终于意识到函数最后的一个简单的“返回0”解决了这个问题。但是,这种行为对我来说似乎很奇怪。这实际上是一个错误吗?