Cython会加速erf（）计算吗？

Question

我正在考虑从MATLAB切换到Python。我的MATLAB代码的核心在几千个数字的数组上重复调用erf，如：

r=rand(1,1e5)

erf(r)

这是我在Python中的实现：

import numpy as np
import scipy.special as sps

r=np.random.rand(1e5)

sps.erf(r)

Python版本大约需要三倍的时间。如果我使用Cython来编译程序的这个核心，我会看到一个主要的加速吗？我没有很少的Python经验，也没有C经验，所以我想在尝试找出Cython之前我会在这里查看。

Answer 1

我不希望Cython能够加快速度。所有耗时的工作都在sps.erf和np.random.rand内进行; Python解释器只是传递数组。

在MATLAB术语中，当您无法完全向量化计算时，Cython最有可能提供帮助。

您可能需要考虑向SciPy开发人员报告此性能差异。也许有一个很容易修复的错误。