我可以将float128设置为numpy中的标准float数组吗

Question

所以我的数字程序有问题，我很好奇它是否是一个精度问题（即四舍五入误差）。有没有一种快速的方法可以将程序中的所有float数组都更改为float128数组，而无需遍历我的代码并dtype='float128'在各处键入内容。我的数组都是float64，但我从未明确写过dtype='float64'，所以我希望有办法改变这种默认行为。

Answer 1

我不认为您可以更改中央“配置”来实现此目的。您可以执行的一些选择：

1. 如果仅通过NumPy的工厂函数中的很少一个来创建数组，请用您自己的版本替换这些函数。如果您导入这些功能，例如

   from numpy import empty
   

你可以做

    from numpy import float128, empty as _empty
def empty(*args, **kwargs):
    kwargs.update(dtype=float128)
    _empty(*args, **kwargs)

如果你在做

    import numpy

你可以写一个模块 mynumpy.py

    from numpy import *
_empty = empty
def empty(*args, **kwargs):
    kwargs.update(dtype=float128)
    _empty(*args, **kwargs)

并像导入

    import mynumpy as numpy

1. 重构代码以始终使用dtype=myfloat。这将使将来的更改变得容易。您可以将此方法与结合使用numpy.empty_like()，numpy.zeros_like()并numpy.ones_like()在适当的地方将实际数据类型硬编码在尽可能少的位置。

2. 子类numpy.ndarray，仅使用您的自定义构造函数创建新数组。