我做了一些cython编码,但没有更基本的c工作。但既然没有其他人回答,我将抛出一些想法。在
另一个迭代数组的工具是nditer。我熟悉它在python和cython中的用法。有一个很好的介绍页。在
另一件事是查看numpy源代码。我将转到github存储库,对这个变量进行搜索。在
对SHAPE的迭代相当于:for i in range(arr.shape[0]):
print arr[i]
您的arrayData[j]正在使用与Python版本相同的数组索引代码。在
二维等价物(在Python中)是:
^{pr2}$
我怀疑c等价物也有类似的结构。在for i in range (arr.size):
print arr[i]
获取IndexError: index 3 is out of bounds for axis 0 with size 3。我假设没有边界检查,你会得到一个分段错误。在
我怀疑使用SIZE你需要直接通过DATAbufer,而不是使用某种形式的arrayData[j]。换句话说,获取缓冲区指针,并按itemsizeSIZE的次数执行。在
nditer将为您节省大量的详细工作。在
A numpy/core计算中的几个地方。c我看到一个块,如:n = PyArray_SIZE(ap)/m;
rptr = (npy_intp *)PyArray_DATA(rp);
for (ip = PyArray_DATA(ap), i = 0; i < n; i++, ip += elsize*m) {
arg_func(ip, m, rptr, ap);
rptr += 1;
}
如果m=1,这将使ip指针逐级通过数据缓冲区,每个元素(每个elsize字节长),将其传递给arg_func执行操作。通常,PyArray_SIZE用于检查数组是否为空(=0)。在
另一个选择是使用PyArray_Flatten(或更好的PyArray_Ravel)来展平数组,然后使用SHAPE迭代。在