- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
假设我有一个 Cython 代码,其中的函数计算滚动移动平均值并返回与输入大小相同的数组(该函数为初始部分添加了 nan,但这不是对于当前的问题很重要)。
我编写了三个 Cython 函数(如下所示):
(a) sma_vec
处理一维的 numpy
数组
(b) sma_mat
处理二维 numpy
数组
(c) 第三个 sma
根据大小从 sma_vec
或 sma_mat
返回值。 (我的动机是最终将 sma_vec
或 sma_mat
之前的 cpdef
替换为 cdef
,以便 Python 代码只看到sma
函数)
函数 1 - 处理一维的 numpy
数组
cimport cython
import numpy as np
cimport numpy as np
from numpy cimport ndarray as ar
ctypedef double dtype_t
@cython.boundscheck(False)
@cython.wraparound(False)
cpdef sma_vec(ar[dtype_t, ndim=1] x, int m):
cdef int n
cdef Py_ssize_t i, j
cdef ar[dtype_t, ndim=1] y
if m == 1:
return x.copy()
else:
y = np.zeros_like(x) * np.nan
n = x.shape[0]
if n < m:
return y
else:
for i in range(m-1, n):
for j in range(i-m+1, i+1):
if j == i-m+1:
y[i] = x[j]
else:
y[i] += x[j]
y[i] /= float(m)
return y
函数 2 - 处理二维 numpy
数组(在 ndarray
的每一行上调用函数 1)
@cython.boundscheck(False)
@cython.wraparound(False)
cpdef sma_mat(ar[dtype_t, ndim=2] x, int m):
cdef int n
cdef Py_ssize_t i
cdef ar[dtype_t, ndim=2] y
if m == 1:
return x.copy()
else:
y = np.zeros_like(x) * np.nan
n = x.shape[0]
if n < m:
return y
else:
for i in range(0, x.shape[0]):
y[i] = sma_vec(x[i], m)
return y
函数3-根据维度调用函数1或函数2
@cython.boundscheck(False)
@cython.wraparound(False)
cpdef sma(ar[dtype_t] x, int m):
if x.ndim == 1:
return sma_vec(x, m)
elif x.ndim == 2:
return sma_mat(x, m)
else:
raise ValueError('Cannot handle more than two dimensions')
测试代码
import numpy as np
import common.movavg as mv
x1 = np.array([1.0, 1.4, 1.3, 5.3, 2.3])
y1 = mv.sma_vec(x1, 3)
y1a = mv.sma(x1, 3)
y1
和 y1a
均正确返回 array([nan, nan, 1.233333, 2.666667, 2.966667])
x2 = np.array([[1.0, 1.4, 1.3, 5.3, 2.3], [4.2, 1.3, 2.3, 5.7, -1.3]])
y2 = mv.sma_mat(x2, 2)
y2
正确返回
array([[ nan, 1.2 , 1.35, 3.3 , 3.8 ],
[ nan, 2.75, 1.8 , 4. , 2.2 ]])
但是当我尝试时:
y2a = mv.sma(x2, 2)
我收到错误:
Traceback (most recent call last):
File "C:\PF\WinPython-64bit-3.4.2.4\python-3.4.2.amd64\lib\site-packages\IPython\core\interactiveshell.py", line 2883, in run_code
exec(code_obj, self.user_global_ns, self.user_ns)
File "<ipython-input-4-dc092e343714>", line 3, in <module>
y2a = mv.sma(x2, 2)
File "movavg.pyx", line 54, in movavg.sma (stat\movavg.c:2206)
ValueError: Buffer has wrong number of dimensions (expected 1, got 2)
在sma
函数中,问题似乎是ar[dtype_t] x
(即np.ndarray[double] x
)自动假定 x.ndim
的尺寸应为 1
。
如何重写 sma
函数,以便它可以接受尺寸未知的 np.ndarray
?
最佳答案
找到了答案。
从此链接:numpy_tutorial ,“...“ndim”仅关键字参数,如果未提供,则假定为一维...”
解决方案是将函数3转换为:
@cython.boundscheck(False)
@cython.wraparound(False)
cpdef sma(ar x, int m):
if x.ndim == 1:
return sma_vec(x, m)
elif x.ndim == 2:
return sma_mat(x, m)
else:
raise ValueError('Cannot handle more than two dimensions')
我们需要完全删除[]
中的所有内容。
关于python - 在 cython 中编写未知维度的 numpy 代码,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/30385584/
c 不做边界检查。那么cython是如何检查是否编译成c的呢? %%cython --annotate cimport cython @cython.boundscheck(True) cpdef m
可以直接声明用于 Cython 构造函数? 据我了解,这是可能的: # Cython cdef int[3] li = [1, 2, 3] # C++ int[3] li = {1, 2, 3} 但
所以,如果你有一个头文件。 %%file test.h struct mystruct{ int i; int j; }; 然后你将它包装在 Cython 中: cdef extern fr
我正在构建一个独立于平台的 cython 项目,我想根据正在使用的编译器传递编译器参数。我可以猜测基于平台的编译器,或者假设它与用于 Python 的编译器相同,但不能保证匹配。通常我注入(injec
我使用诗歌构建我的 cython 包。我在所有函数和类中都有 NumPy 风格的文档字符串。我现在要做的是添加 Sphinx 自动文档并发布在 Read the Docs。 我已阅读此主题 How d
赛通 libcpp模块包含 priority_queue 的模板,这很好,除了一件事:我不能通过自定义比较器(或者,至少,我不知道如何)。 我需要这个,因为我需要 priority_queue做一个a
以下代码定义了一个简单的 Cython 函数(为方便起见,使用 Ipython 魔法)。 %load_ext cython %%cython def f(float x, float y=2):
我正在尝试使用 cython 进行复数计算。在示例代码中,我想计算复数的复指数函数。问题是我不知道如何将我的整数乘以虚数单位。python的虚数单位1.0j乘以cython执行时报错。 这是我的代码:
在这里停留在一些基本的 Cython 上 - 在 Cython 中定义字符串数组的规范且有效的方法是什么? 具体来说,我想定义一个定长常量数组char . (请注意,此时我不想引入 NumPy。) 在
是否有可能,如果是,如何确定 Cython 中整数数据类型的大小(以位为单位)? 我正在尝试做这样的事情,以获得整数大小: cdef WORD_BITS = 0 IF sizeof(unsigned
我只是想打印 cython 变量的地址,但我无法绕过错误消息: cdef int myvar print &myvar 抛出 Cannot convert 'int *' to Python obje
我有一个 C 头文件,它在宏中定义了一个函数。我需要从 Cython 调用它。有没有办法在 Cython 中使用宏并使其完全扩展?我已经有了 C 类型的参数。 我尝试像使用函数一样使用 cdef,我认
令人惊讶的是,我似乎找不到通过名称获取结构体元素的单个示例(无论是在网络上还是在 cython 示例中)。 所以我收到了一个指向 C 函数结构体的指针,并且想要一一访问这些元素并将它们重新打包到 py
我尝试围绕 C++ 库编写一个 Cython 包装器 http://primesieve.org/ 它包装了一个函数count。到目前为止,它可以正确安装 python setup.py instal
我正在尝试将 cython 模块 data.pyx 导入另一个 cython 模块 user.pyx。一切都编译得很好,但是当我尝试在 python 模块中调用 user.pyx 时,我收到错误“Im
更新:内存 View 获胜。Cython 使用类型化内存 View :0.0253449 特别感谢 lothario,他指出了几个关键的变化。 荒谬。当然现在的问题是,似乎不能对它们做太多算术(加法和
我有一个使用 memoryview 数组的 cython 模块,即... double[:,:] foo 我想使用多处理并行运行这个模块。但是我得到了错误: PicklingError: Can't
我正在尝试使用 Cython 加速 PEP 484 类型的 python 脚本。我想保持一些语义和可读性。 之前,我有一个 Flags = int def difference(f1: Flags,
这个问题已经有答案了: Collapse multiple submodules to one Cython extension (5 个回答) 已关闭 3 年前。 我在一个包中有多个 .py 文件
我已经能够在我的 .pyx 脚本上使用 cython 在 linux 上创建一个 .so 文件。我也可以成功地在我的 python 解释器上进行导入。 我的问题是如何在不使用 cython 的情况下将
我是一名优秀的程序员,十分优秀!