Python是否具有與java.lang.Math.nextUp等效的功能？

Question

我有一個Python float，我想要float這是1 ULP s越來越小。

在Java中，我會用Math.nextUp(x)和Math.nextAfter(x, Double.NEGATIVE_INFINITY)來做到這一點。

有沒有辦法在Python中做到這一點？我想用math.frexp and math.ldexp自己來實現它，但據我所知Python沒有指定浮點類型的大小。

Answer 1

我不確定這是不是你想要的，但是sys.float_info.epsilon是「1和大於1的最小值之間的差值」，你可以做x * (1 + sys.float_info.epsilon)。

http://docs.python.org/library/sys.html#sys.float_info

Answer 2

你可以看看如何Decimal.next_plus()/Decimal.next_minus()實現：

>>> from decimal import Decimal as D 
>>> d = D.from_float(123456.789) 
>>> d 
Decimal('123456.789') 
>>> d.next_plus() 
Decimal('123456.789') 
>>> d.next_minus() 
Decimal('123456.789') 
>>> d.next_toward(D('-inf')) 
Decimal('123456.789') 

確保decimal context有你需要的值：

>>> from decimal import getcontext 
>>> getcontext() 
Context(prec=28, rounding=ROUND_HALF_EVEN, Emin=-999999999, Emax=999999999, 
capitals=1, flags=[], traps=[InvalidOperation, DivisionByZero, Overflow]) 

---

的替代品：

• 呼叫C99 nextafter()使用​​：

  >>> import ctypes 
  >>> nextafter = ctypes.CDLL(None).nextafter 
  >>> nextafter.argtypes = ctypes.c_double, ctypes.c_double 
  >>> nextafter.restype = ctypes.c_double 
  >>> nextafter(4, float('+inf')) 
  4.000000000000001 
  >>> _.as_integer_ratio() 
  (4503599627370497, 1125899906842624) 
  

  使用numpy：

  >>> import numpy 
  >>> numpy.nextafter(4, float('+inf')) 
  4.0000000000000009 
  >>> _.as_integer_ratio() 
  (4503599627370497, 1125899906842624) 
  

  儘管不同repr()，結果是一樣的。

• 如果我們忽略的邊緣情況，然後從@S.Lott answer作品簡單frexp/ldexp解決方案：

  >>> import math, sys 
  >>> m, e = math.frexp(4.0) 
  >>> math.ldexp(2 * m + sys.float_info.epsilon, e - 1) 
  4.000000000000001 
  >>> _.as_integer_ratio() 
  (4503599627370497, 1125899906842624) 
  

• pure Python next_after(x, y) implementation by @Mark Dickinson，考慮到邊緣情況。結果在這種情況下是相同的。