用Numpy打领带

Question

标准的numpy回合平局打破了IEEE 754约定，将一半舍入到最接近的偶数。有没有一种方法可以指定不同的舍入行为，例如，向零或向-
inf舍入？我说的不是天花板或地板，我只需要换一条领带即可。

Answer 1

NumPy不对内部舍入模式进行任何控制。这有两种选择：

1. 使用此答案中gmpy2概述的内容。这使您可以完全控制舍入模式，但是用于简单的浮点数学运算可能会比NumPy慢。gmpy2

2. 使用fesetroundviactypes手动设置舍入模式。这是特定于系统的，因为常量可能因平台而异。检查fenv.h平台上的常量值。在我的机器上（Mac OS X）：

   import numpy as np
   

   import ctypes
   FE_TONEAREST = 0x0000
   FE_DOWNWARD = 0x0400
   FE_UPWARD = 0x0800
   FE_TOWARDZERO = 0x0c00
   libc = ctypes.CDLL(‘libc.dylib’)
   
   

   v = 1. / (1<<23)
   print repr(np.float32(1+v) - np.float32(v/2)) # prints 1.0
   libc.fesetround(FE_UPWARD)
   print repr(np.float32(1+v) - np.float32(v/2)) # prints 1.0000002