如何在Haskell中获取给定Float/Double的相邻浮点数（等效C++ std::nextafter）？

这个问题问得好！要在Haskell里实现和C++ std::nextafter等效的功能，而且只用base库不用第三方工具，核心思路是利用IEEE 754浮点数的二进制编码特性——标准浮点数的整数位模式顺序和数值大小顺序是一致的（除了NaN、正负零这类特殊值）。下面我会一步步讲怎么实现，还有各种特殊情况的处理。

IEEE 754单精度（Float）和双精度（Double）浮点数，在内存中的二进制位可以直接映射成对应的无符号整数：Float对应Word32，Double对应Word64。对于正常的非特殊值：

• 把浮点数转成整数后加1，再转回浮点数，就是原数的下一个更大值
• 把浮点数转成整数后减1，再转回浮点数，就是原数的下一个更小值

首先我们需要导入几个基础模块：

• Data.Word：提供Word32和Word64类型
• Unsafe.Coerce：用来在浮点数和整数之间做位级转换（这里是安全的，因为我们明确知道类型对应关系）
• Data.Float：用来判断NaN、无穷大、负零等特殊值

1. 实现Float版本的nextUp/nextDown #

先写针对Float的“上一个”和“下一个”函数：

import Data.Word (Word32)
import Unsafe.Coerce (unsafeCoerce)
import Data.Float (isNaN, isInfinite, isNegativeZero)

-- 获取比给定Float大一点的数
nextUpFloat :: Float -> Float
nextUpFloat x
  | isNaN x = x  -- NaN返回自身
  | isInfinite x && x > 0 = x  -- 正无穷没有更大的值，返回自身
  | isNegativeZero x = 0.0  -- 负零的下一个值是正零
  | x == 0.0 = unsafeCoerce (1 :: Word32)  -- 正零的下一个是最小的正Float
  | otherwise = unsafeCoerce (unsafeCoerce x + 1 :: Word32)

-- 获取比给定Float小一点的数
nextDownFloat :: Float -> Float
nextDownFloat x
  | isNaN x = x
  | isInfinite x && x < 0 = x  -- 负无穷没有更小的值，返回自身
  | x == 0.0 = unsafeCoerce (maxBound - 1 :: Word32)  -- 正零的前一个是最小的负Float
  | isNegativeZero x = unsafeCoerce (maxBound :: Word32)  -- 负零的前一个是最大的负Float
  | otherwise = unsafeCoerce (unsafeCoerce x - 1 :: Word32)

2. 实现Double版本的nextUp/nextDown #

逻辑和Float完全一致，只是把Word32换成Word64：

import Data.Word (Word64)

nextUpDouble :: Double -> Double
nextUpDouble x
  | isNaN x = x
  | isInfinite x && x > 0 = x
  | isNegativeZero x = 0.0
  | x == 0.0 = unsafeCoerce (1 :: Word64)
  | otherwise = unsafeCoerce (unsafeCoerce x + 1 :: Word64)

nextDownDouble :: Double -> Double
nextDownDouble x
  | isNaN x = x
  | isInfinite x && x < 0 = x
  | x == 0.0 = unsafeCoerce (maxBound - 1 :: Word64)
  | isNegativeZero x = unsafeCoerce (maxBound :: Word64)
  | otherwise = unsafeCoerce (unsafeCoerce x - 1 :: Word64)

3. 实现类似std::nextafter的通用函数 #

std::nextafter的逻辑是：返回从第一个参数向第二个参数方向的下一个浮点数。我们可以基于上面的函数封装：

-- Float版本的nextAfter
nextAfterFloat :: Float -> Float -> Float
nextAfterFloat x y
  | isNaN y = y  -- 符合C++标准：若第二个参数是NaN，返回NaN
  | isNaN x = x
  | x == y = x
  | y > x = nextUpFloat x
  | otherwise = nextDownFloat x

-- Double版本的nextAfter
nextAfterDouble :: Double -> Double -> Double
nextAfterDouble x y
  | isNaN y = y
  | isNaN x = x
  | x == y = x
  | y > x = nextUpDouble x
  | otherwise = nextDownDouble x

• 输入nextUpFloat 0.0会返回你提到的1.40129846432481707e-45（打印时通常显示为1.0e-45）
• 输入nextDownFloat 1.0会返回0.9999999403953552（打印时显示为0.99999994）

你可能会担心unsafeCoerce的安全性，但在这里它是完全安全的：IEEE 754浮点数的内存表示和对应的无符号整数是一一对应的（除了NaN的多种位模式，但我们已经提前处理了NaN的情况）。只要保证Float对应Word32、Double对应Word64，就不会有问题。

内容的提问来源于stack exchange，提问作者Joseph Sible-Reinstate Monica