Python：make scipy使用numpy.float128而不是numpy.float64？

Question

是否可能，如果是的话怎么做：默认情况下如何使用scumpy numpy.float128。例如

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

我希望它是

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

如果没有，我需要自己实现norm.pdf功能，这很容易，但不能解决我的问题。

Answer 1

一般情况下，你不能。一些scipy例程是用C或Fortran编写的代码包装器，只能以双精度方式提供。即使你弄清楚哪些是纯python + numpy，并设法确保在计算中执行的操作保留数据类型，你会发现许多函数使用硬编码的64位常量，如numpy.pi