我可以設置float128作爲numpy中的標準浮點數組

Question

所以我的數值程序有問題，我很好奇它是否是精度問題（即舍入誤差）。有沒有一種快速的方法可以將我的程序中的所有浮點數組更改爲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

   和進口它像

   import mynumpy as numpy 
   

2. 重構代碼總是使用dtype=myfloat。這將使這種變化在未來變得容易。您可以將此方法與numpy.empty_like(),numpy.zeros_like()和numpy.ones_like()的使用結合使用，儘可能將硬編碼的實際數據類型放置在儘可能少的地方。

3. 子類numpy.ndarray並且只使用您的自定義構造函數來創建新數組。