Python: вычислить точность с плавающей запятой как функцию от значения с плавающей запятой

Question

Я хотел бы воспроизвести следующий сюжет со страницы википедии по стандарту IEEE 754 :

К сожалению, в статье в Википедии не указано, как они создали эти данные, и вычисление не кажется тривиальным.

Answer 1

Я не уверен на 100%, что именно там изображено, но np.nextafter, похоже, дает довольно похожие результаты:

>>> x = 10**np.linspace(-12, 12, 50)
>>> y = np.nextafter(x, 2*x) - x

и

>>> x = 10**np.linspace(-12, 12, 50, dtype=np.float32)
>>> y = np.nextafter(x, 2*x) - x