声明numpy的阵列的列表/矢量/阵列在用Cython

Question

说我有几个阵列，可能大小不同的：

A0 = rand(3,3) 
A1 = rand(4,4) 

在用Cython，我可以声明它们的类型得到快速项目访问：

cdef np.ndarray[double, ndim=2] A0 
cdef np.ndarray[double, ndim=2] A1 

不过，说我要通过索引来访问它们：

A = (A0,A1) 
A[0][2,1] += A[1][1,0] 

不过，现在用Cython不知道的类型和A[1]，这使访问变慢。我不认为Cython具有“类型元组”的概念。那么我怎么能声明A（或类似的对象），以便我仍然可以在上述表达式中快速访问项目？

Answer 1

这是否是有效cython是投机性的，但灵感来自cython例如在结束：

https://docs.scipy.org/doc/numpy/reference/arrays.nditer.html#putting-the-inner-loop-in-cython

我建议：

cdef foo(aTuple): 
    cdef np.ndarray[double, ndim=2] A0 
    cdef np.ndarray[double, ndim=2] A1 
    A0, A1 = aTuple # use unpacking 
    A0[2,1] += A1[1,0]  

https://github.com/cython/cython/blob/master/tests/run/unpack.pyx - 是解压测试pyx。

也许使用A0 = aTuple[0]和A1 = aTuple[1];这个想法是使用数组键入它重要的地方，用于它自己的索引。该元组只保存指针，所以打字无关紧要。

在cython使用数组时，另一个有用的页面是http://cython.readthedocs.io/en/latest/src/userguide/memoryviews.html