Python: 如何让scipy使用numpy.float128而不是numpy.float64？

Question

是否可能，并且如果可以，如何做到：如何使scipy默认使用numpy.float128。例如

>>> from scipy.stats import norm
>>> type(norm.pdf(10, 10, 1))
<class 'numpy.float64'>

and I want it to be

>>> from scipy.stats import norm
>>> type(norm.pdf(10, 10, 1))
<class 'numpy.float128'>

如果没有的话，我就需要自己实现norm.pdf函数，这个很简单，但是不能解决我的问题。

- Michał Piotr Stankiewicz

---

也许这可以帮到你？http://mail.scipy.org/pipermail/scipy-user/2013-February/034063.html - J'e

Answer 1

一般情况下是不行的。一些scipy例程是用C或Fortran编写的包装器，仅限于双精度。即使你确定哪些是纯python+numpy，并确保计算中执行的操作保留数据类型，你会发现许多函数使用硬编码的64位常数，例如numpy.pi。

- Warren Weckesser