Linux下浮点异常处理：如何将其转化为可用的错误提示机制

我有一段用于计算自变量向量（及参数）函数向量的代码，结构如下：

struct instance
{    
    double m_adParam1, m_adParam2;
    std::array<double, OUTSIZE> calculate(const std::array<double, INSIZE>& x) const;
};

要求是：只要计算过程中出现除零（或接近零）、定义域错误这类执行失败的情况，不管是核心计算还是仅用于比较的中间值，整个计算都要终止并告知调用者无法处理该输入。但逐操作加检查太麻烦，只检查结果的NaN/无穷大又不够——毕竟中间异常值可能在后续操作中被掩盖，没法被结果检查捕捉到。

我试过了解fetestexcept()，但搞不清它的状态标志会被什么操作重置；SIGFPE除了调试基本没用，要是能像Windows那样直接转成C++异常就最好了，退一步说，在计算结束后统一检查也能接受。另外，给后续代码加无效输入检查绝对不行，会把计算速度拖垮。

1. 用fetestexcept()做事后统一检查 #

这是最稳妥且低开销的方案，步骤清晰：

• 在calculate函数最开头，先调用feclearexcept(FE_ALL_EXCEPT)清除所有浮点异常标志，避免之前的计算残留影响本次检测。
• 正常执行所有计算逻辑。
• 计算结束后，调用fetestexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW | FE_UNDERFLOW)，这个函数会返回非零值如果指定的异常被触发过：
  • 比如FE_DIVBYZERO对应除零，FE_INVALID对应定义域错误（比如对负数开根号、0取对数）。
  • 如果返回非零，直接抛出自定义C++异常（比如struct FloatCalcError : public std::runtime_error {};），或者返回一个带错误标记的可选值（比如std::optional<std::array<double, OUTSIZE>>）通知调用者。
• 注意：fetestexcept()不会自动重置异常标志，每次检测完或者下次计算前必须手动调用feclearexcept清除，否则下一次检测会误判。另外部分编译器需要开编译选项，比如GCC要加-fnon-call-exceptions才能让浮点异常状态被正确捕获。

2. 把SIGFPE转成C++异常（类Windows模式） #

类Unix系统下可以通过信号处理实现类似Windows的浮点异常转C++异常：

• 首先在程序初始化阶段注册SIGFPE的信号处理函数，在处理函数里直接抛出自定义异常：

  #include <signal.h>
  #include <stdexcept>
  
  struct FloatCalcError : public std::runtime_error {
      using std::runtime_error::runtime_error;
  };
  
  void sigfpe_handler(int) {
      throw FloatCalcError("浮点计算异常：除零或定义域错误");
  }
  
  // 程序启动时调用一次
  void setup_float_exception() {
      signal(SIGFPE, sigfpe_handler);
  }
  

• 编译时要加GCC的-fnon-call-exceptions选项，否则信号处理函数里抛异常会触发未定义行为。
• 这种方式的好处是异常触发时会立刻中断计算，不需要等到最后检查，但要注意多线程环境下的信号处理绑定问题——信号是进程级的，多线程中需要额外处理才能确保每个线程的异常都能正确捕获。

3. 编译选项强制浮点异常陷阱 #

可以通过编译器选项让浮点异常直接触发陷阱，比如GCC的-ffpe-trap=divzero,invalid，开启后只要遇到除零、定义域错误，立刻触发SIGFPE，再结合上面的信号转异常方案，就能实现和Windows几乎一致的体验。

• 但要注意，这个选项是全局的，会影响整个程序的浮点计算行为，如果其他代码依赖“静默产生NaN/无穷大”的行为，就不能用这个方案。

关键提醒 #

• 绝对不要只检查结果的NaN/无穷大：中间计算的异常值可能在后续操作中被“消化”（比如NaN参与比较后结果为假，但不会留下可见痕迹），只有检测浮点异常状态才能确保所有错误都被捕捉。
• 性能问题不用担心：feclearexcept和fetestexcept都是极轻量的操作，开销微乎其微，完全不会像逐操作检查那样拖慢速度。

内容的提问来源于stack exchange，提问作者Frank Puck