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开发人员。 也许存在一个易于修复的错误。