- c - 在位数组中找到第一个零
- linux - Unix 显示有关匹配两种模式之一的文件的信息
- 正则表达式替换多个文件
- linux - 隐藏来自 xtrace 的命令
我在 cython
中编写了以下函数来估计对数似然
@cython.boundscheck(False)
@cython.wraparound(False)
def likelihood(double m,
double c,
np.ndarray[np.double_t, ndim=1, mode='c'] r_mpc not None,
np.ndarray[np.double_t, ndim=1, mode='c'] gtan not None,
np.ndarray[np.double_t, ndim=1, mode='c'] gcrs not None,
np.ndarray[np.double_t, ndim=1, mode='c'] shear_err not None,
np.ndarray[np.double_t, ndim=1, mode='c'] beta not None,
double rho_c,
np.ndarray[np.double_t, ndim=1, mode='c'] rho_c_sigma not None):
cdef double rscale = rscaleConstM(m, c,rho_c, 200)
cdef Py_ssize_t ngals = r_mpc.shape[0]
cdef np.ndarray[DTYPE_T, ndim=1, mode='c'] gamma_inf = Sh(r_mpc, c, rscale, rho_c_sigma)
cdef np.ndarray[DTYPE_T, ndim=1, mode='c'] kappa_inf = Kap(r_mpc, c, rscale, rho_c_sigma)
cdef double delta = 0.
cdef double modelg = 0.
cdef double modsig = 0.
cdef Py_ssize_t i
cdef DTYPE_T logProb = 0.
#calculate logprob
for i from ngals > i >= 0:
modelg = (beta[i]*gamma_inf[i] / (1 - beta[i]*kappa_inf[i]))
delta = gtan[i] - modelg
modsig = shear_err[i]
logProb = logProb -.5*(delta/modsig)**2 - logsqrt2pi - log(modsig)
return logProb
但是当我运行这个函数的编译版本时,我得到以下错误信息:
File "Tools.pyx", line 3, in Tools.likelihood
def likelihood(double m,
ValueError: ndarray is not C-contiguous
我不太明白为什么会出现这个问题??!!!我将不胜感激获得任何有用的提示。
最佳答案
就在您收到错误之前,尝试打印您传递给 likelihood
的 numpy 数组的 flags
属性。您可能会看到如下内容:
In [2]: foo.flags
Out[2]:
C_CONTIGUOUS : False
F_CONTIGUOUS : True
OWNDATA : True
WRITEABLE : True
ALIGNED : True
UPDATEIFCOPY : False
请注意它说 C_CONTIGUOUS : False
的地方,因为这就是问题所在。要修复它,只需将其转换为 C 顺序:
In [6]: foo = foo.copy(order='C')
In [7]: foo.flags
Out[7]:
C_CONTIGUOUS : True
F_CONTIGUOUS : False
OWNDATA : True
WRITEABLE : True
ALIGNED : True
UPDATEIFCOPY : False
关于python - 值错误 : ndarray is not C-contiguous in cython,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/26778079/
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 的情况下将
我是一名优秀的程序员,十分优秀!