Eigen库结合MSVC编译出现意外内存行为问题求助

使用MSVC 19.34.31944编译器、Eigen 3.4.0线性代数库，以/Oxt /std:c++20 /EHs选项编译C++代码时出现两处异常行为，已将问题简化为两个可复现的示例。Eigen官方初步判断为MSVC编译器bug，以下是详细分析：

示例一：临时对象成员引用后的内存失效 #

测试代码 #

#include "Eigen/Dense"
#include <iostream>

std::pair<Eigen::Vector2d, double> func() {
    return { Eigen::Vector2d{11, 0}, 0 };
};

int main(int argc, char *argv[]) {
    const Eigen::Vector2d &result = func().first;

    std::cout << *result.data() << "  " << *(result.data() + 1) << std::endl;
    std::cout << result.transpose() << std::endl;
    std::cout << func().first.transpose() << std::endl;

    return 0;
}

预期输出 #

11  0
11  0
11  0

实际输出 #

11  0
6.95158e-310 6.95145e-310
11  0

现象分析 #

• 仅在MSVC Release模式（/Oxt优化）下复现，Debug模式或Linux G++编译时结果正常。
• 直接通过data()访问向量元素时结果正确，但调用transpose()输出时出现野值；若添加无transpose()的std::cout输出，异常会消失。
• 根据C++标准，临时对象成员绑定到const引用的操作合规，且改用值拷贝（非引用）方式可得到正确结果，排除代码逻辑问题。

示例二：临时对象成员与局部变量内存地址复用 #

测试代码 #

#include "Eigen/Dense"
#include <iostream>
#include <numbers>

// Computes the angle between two 2D vectors in [-pi, pi]
inline double computeAngle(const Eigen::Vector2d &vec1, const Eigen::Vector2d &vec2) {
    double angle = std::atan2(vec2(1), vec2(0)) - std::atan2(vec1(1), vec1(0));
    angle = angle > std::numbers::pi ? angle - 2.0 * std::numbers::pi : angle;
    angle = angle < -std::numbers::pi ? angle + 2.0 * std::numbers::pi : angle;
    return angle;
}

std::pair<Eigen::Vector2d, double> func() {
    double coeffMid = std::sin(2 * std::abs(computeAngle(Eigen::Vector2d{ 5, -5 }, Eigen::Vector2d{ -5, -5 })));
    std::cout << coeffMid << std::endl;

    Eigen::Vector2d center = (Eigen::Vector2d{ 10, 0 } + Eigen::Vector2d{ 0, 0 }) / 2;

    return { Eigen::Vector2d{5, 0}, 0 };
}

int main(int argc, char *argv[]) {
    const Eigen::Vector2d &result = func().first;
    std::cout << &result << std::endl;

    for (size_t i = 0; i < 5; i++) {
        const Eigen::Vector2d foo = result + result;
        std::cout << &foo << " " << foo[0] << " " << foo[1] << std::endl;
    }
    return 0;
}

预期输出 #

1.22465e-16
0000001954B3F9B0
0000001954B3F990 10 0
0000001954B3F990 10 0
0000001954B3F990 10 0
0000001954B3F990 10 0
0000001954B3F990 10 0

实际输出 #

1.22465e-16
00000025FB4FFEA0
00000025FB4FFEA0 10 0
00000025FB4FFEA0 20 0
00000025FB4FFEA0 40 0
00000025FB4FFEA0 80 0
00000025FB4FFEA0 160 0

现象分析 #

• 打印内存地址后发现，result（临时对象的Vector2d成员）与局部变量foo复用了同一块内存地址。
• 每次循环中foo = result + result的赋值操作会覆盖内存中的数据，导致后续读取到的是累加后的结果，而非预期的固定值10。

技术解析 #

1. 临时对象生命周期延长规则的违反：
   根据C++标准，当const引用绑定到临时对象的成员时，整个临时对象（此处为std::pair<Eigen::Vector2d, double>）的生命周期应被延长至该引用的生命周期结束。但MSVC在/Oxt优化下，未正确识别这种绑定关系，导致pair临时对象被提前销毁，其内部的Vector2d成员内存被释放，后续通过引用访问时触发未定义行为（如示例一中的野值）。

2. 内存分配优化的误判：
   示例二中，编译器错误地将局部变量foo的内存地址分配给已失效的临时对象成员result的地址，违反了对象内存独立性的规则。这导致每次foo的构造都会覆盖同一块内存，出现数值累加的异常结果，属于优化器在内存布局处理上的逻辑错误。

3. 触发条件的特殊性：
   异常仅在MSVC Release优化模式下出现，Debug模式或G++编译器中无此问题，说明是MSVC优化器在处理临时对象成员引用场景时的特定bug，可能与临时对象生命周期跟踪、内存复用策略的实现缺陷有关。

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