使用-O2编译时Eigen::Vector2f初始化异常问题咨询

问题现象 #

使用Eigen::Vector2f(x.x(), x.y()) + delta的方式初始化变量并通过auto推导类型时，在-O2编译优化下会出现初始化错误；无优化编译时结果正常。

测试代码（main.cpp） #

#include <Eigen/Core>
#include <iostream>

int main(int argc, char** argv) {
    const auto x = Eigen::Vector2f(1.0, 2.0);
    const auto delta = Eigen::Vector2f(0.1, 0.2);
    const auto add_1 = x + delta;
    const auto add_2 = Eigen::Vector2f(x.x(), x.y()) + delta;

    std::cout << "x: " << x.x() << ", " << x.y() << std::endl;
    std::cout << "delta: " << delta.x() << ", " << delta.y() << std::endl;
    std::cout << "add_1: " << add_1.x() << ", " << add_1.y() << std::endl;
    std::cout << "add_2: " << add_2.x() << ", " << add_2.y() << " <-- " << std::endl;

    return 0;
}

编译输出 #

• 无优化编译结果

> g++ $(pkg-config --cflags eigen3) main.cpp; ./a.out
x: 1, 2
delta: 0.1, 0.2
add_1: 1.1, 2.2
add_2: 1.1, 2.2 <-- 

• -O2优化编译结果

g++ $(pkg-config --cflags eigen3) -O2 main.cpp; ./a.out
x: 1, 2
delta: 0.1, 0.2
add_1: 1.1, 2.2
add_2: 9.97967e+33, 0.2 <-- 

原因分析 #

这个问题并非Eigen::Vector2f(x.x(), x.y())的构造方式不被官方支持，而是Eigen表达式模板机制与auto推导结合产生的未定义行为：

1. Eigen为优化性能采用表达式模板延迟计算——operator+返回的不是直接的Vector2f对象，而是代表加法运算的轻量级表达式对象（如CwiseBinaryOp），该对象会引用参与运算的操作数。
2. 对于add_2，Eigen::Vector2f(x.x(), x.y())是临时对象，语句执行结束后就会被销毁；但auto推导的add_2是表达式模板类型，会持有对该临时对象的引用。
3. 后续访问add_2.x()和add_2.y()时，临时对象的内存已被释放或覆盖，导致读取到非法值。无优化时临时对象内存未被覆盖，看似正常，但本质是未定义行为。
4. add_1的情况中，x和delta是main函数作用域内的变量，生命周期覆盖add_1的使用，表达式模板引用的对象始终有效，因此结果正确。

解决方案 #

可以通过以下方式避免该问题：

• 显式指定变量类型：将auto替换为Eigen::Vector2f，强制表达式立即求值并存储为实际向量对象：

  const Eigen::Vector2f add_2 = Eigen::Vector2f(x.x(), x.y()) + delta;
  

• 使用eval()强制求值：对临时向量调用eval()生成持久化对象，避免表达式模板引用临时对象：

  const auto add_2 = Eigen::Vector2f(x.x(), x.y()).eval() + delta;
  

• 直接复用已有向量：既然x本身就是Vector2f，直接写成x + delta（如add_1的写法），既简洁又避免问题。

官方支持情况 #

Eigen::Vector2f(x.x(), x.y())这种构造方式是官方支持的，问题出在后续表达式模板与临时对象的生命周期冲突上，而非构造函数本身的问题。

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