如何使用c++20的std::lerp进行精确的线性插值？ (数值计算)

std::lerp 并不提升数值精度，其核心价值在于对 NaN、Inf 等特殊浮点值的明确定义行为和边界健壮性；它不重排运算、不提升中间精度，适用于需稳定异常处理的场景，而非高精度计算。

std::lerp 在大多数场景下并不比手动写 a + t * (b - a) 更精确，反而可能因中间计算顺序引入额外误差；它真正价值在于对特殊浮点值（如 NaN、Inf）的明确定义行为。

std::lerp 的设计目标不是提升精度，而是提升健壮性

标准明确要求 std::lerp(a, b, t) 等价于 a + t * (b - a) 仅当所有操作都“按数学意义”进行——但实际中编译器和硬件会做优

化，导致结果差异。尤其当 a 和 b 量级相差极大时，b - a 可能发生有效位丢失，而 std::lerp 并不规避这点。

• 它不重排运算顺序，也不使用更高精度中间类型
• 它对 t 或 t > 1 的外推行为是明确定义的（不是未定义行为）
• 当 a 或 b 是 NaN 时，返回 NaN；当任一为 Inf 且 t 非 0/1 时，按 IEEE 754 规则传播

何时该用 std::lerp，何时该避免

用在需要可预测浮点边界行为的场合，比如图形 API 输入校验、物理引擎状态混合、序列化插值参数传递；不用在对 ULP 误差敏感的科学计算中。

• ✅ 适合：std::lerp(0.0, 1e30, 1e-16) —— 返回接近 0.0，行为稳定
• ✅ 适合：std::lerp(NAN, 42.0, 0.5) —— 明确返回 NAN，不触发浮点异常
• ❌ 不适合：std::lerp(1.0, 1.0 + 0x1p-52, 0.5) —— 和手写公式结果一致，但无精度增益
• ❌ 不适合：需要 long double 中间精度的插值（std::lerp 永远只用参数类型精度）

替代方案：手动控制精度与顺序

若真需更高精度线性插值，应显式选择算法和类型，而非依赖 std::lerp。

• 对 double 输入，用 long double 中间计算：

  long double lerp_highprec(double a, double b, double t) {
      return static_cast(a) + 
             static_cast(t) * 
             (static_cast(b) - static_cast(a));
  }

• 对极小 t，改用 a + fma(t, b - a, 0.0)（若支持 std::fma），减少一次舍入
• 对整数归一化权重（如 0–65535），用定点算术避免浮点累积误差

std::lerp 的“精确”是语义上的精确——它让插值在异常输入下不崩溃、不静默出错；数值层面的精度，它既不承诺也不提供。真正需要精度的地方，得自己选算法、控类型、审汇编。