Ubuntu18.04下Eigen::inverse加-ffast-math输出错误矩阵求助

问题场景 #

在Ubuntu 18.04环境中使用Eigen库的Eigen::inverse函数时，发现通过-O3结合-ffast-math编译选项生成的代码，计算得到的逆矩阵部分元素存在符号错误。

重现步骤 #

• 启动容器：docker run -v $PWD:/work -it ubuntu:18.04
• 更新系统源：apt update
• 安装依赖组件：
  • 安装Eigen库：apt install libeigen3-dev -y
  • 安装GCC 10.3.0版本：
    • apt install software-properties-common -y
    • add-apt-repository ppa:ubuntu-toolchain-r/test -y && apt update
    • apt install gcc-10 gcc-10-base gcc-10-doc g++-10 -y
• 创建测试代码文件eigen_test.cpp，内容如下：

#include <iostream>
#include <eigen3/Eigen/Eigen>
int main()
{
    Eigen::Matrix4d t;
    t << 1.0, -0.0, 0.0, -0.0,
         0.0, 1.0, -0.0, -0.0,
         0.0, 0.0, 1.0, -0.0,
         0.0, 0.0, 0.0, 1.0;
    std::cout << "Original matrix: \n" <<t << std::endl;
    std::cout << "Inversed matrix: \n" << t.inverse() << std::endl;
}

• 编译运行验证：
  • 使用g++-10 eigen_test.cpp -O3 -ffast-math -o a.out编译并运行，逆矩阵的(2,2)和(4,4)位置元素为-1，结果错误；
  • 使用g++-10 eigen_test.cpp -O3 -o a.out编译并运行，对应位置元素为1，结果正确。

原因分析 #

1. -ffast-math的浮点数语义修改
   -ffast-math是GCC的激进优化选项，它会打破C++标准对浮点数运算的严格要求：允许编译器对浮点数进行代数重排、忽略NaN/无穷大的特殊处理，并且将+0.0和-0.0视为完全等价的数值。这种优化虽然能提升运算速度，但会破坏依赖精确浮点数符号的算法逻辑。

2. -0.0的特殊符号影响
   测试矩阵中存在多个-0.0元素，虽然从数值大小上它和0.0一致，但Eigen的矩阵求逆实现中，部分内部计算（如伴随矩阵生成、行列式符号推导）会保留并使用浮点数的符号信息。当启用-ffast-math后，编译器在优化过程中篡改了-0.0的符号处理逻辑，导致Eigen内部计算出现符号错误，最终使逆矩阵的对角元素符号反转。

3. Eigen与-ffast-math的兼容性限制
   Eigen官方明确说明，-ffast-math可能导致部分数值算法异常，因为该选项破坏了浮点数运算的严格语义。对于涉及符号零、特殊浮点数状态的场景，-ffast-math的优化会干扰计算逻辑，产生不符合预期的结果。

内容的提问来源于stack exchange，提问作者Sean Xia