如何将行主序Eigen稠密矩阵映射到预分配C++ vector？

你的代码有两个核心问题导致段错误：

1. vector未预分配空间：默认构造的vector<double> array是空的，没有任何元素，直接通过array[i]访问下标会越界，触发内存访问错误。
2. 元素访问方式错误：*mat.data() + i*sizeof(double)的逻辑完全错误——mat.data()是指向矩阵第一个元素的指针，指针算术会自动根据double类型的字节数步进，手动乘sizeof(double)会导致访问地址偏移过大，触及非法内存；而且这个表达式是取第一个元素的值加上字节数，并非访问第i个元素。

正确实现方式 #

因为要求预分配vector，先确保vector有足够存储空间，再正确访问矩阵元素：

方式1：循环赋值（手动控制）

// 假设你的行主序矩阵有6个元素（比如2x3）
vector<double> array(6); // 预分配6个元素的空间，完成初始化
for (int i = 0; i < 6; ++i) {
    // 两种等价的正确访问方式
    array[i] = mat.data()[i];
    // 或者 array[i] = *(mat.data() + i);
}

方式2：用std::copy简化代码

利用Eigen矩阵行主序时元素连续存储的特性，直接用标准库的std::copy完成复制：

vector<double> array(6);
std::copy(mat.data(), mat.data() + 6, array.begin());

方式3：直接通过构造函数预分配并初始化

如果不需要先预分配空空间再写入，也可以直接用矩阵的内存区间构造vector，这同样是预分配的方式：

// mat.size()返回矩阵总元素数，自动匹配大小
vector<double> array(mat.data(), mat.data() + mat.size());

额外注意 #

确保你的Eigen矩阵确实是行主序：默认MatrixXd是列主序，若要行主序需要显式声明，比如：

// 声明行主序的动态矩阵
Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor> mat(2, 3);

或者对已有列主序矩阵转换为行主序存储：

auto row_major_mat = mat.rowMajor();

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