如何为待填充的matlab::data::Array构建Eigen::Map？

你的代码崩溃的本质是**getOutDataPtr函数采用值传递参数，生成了原数组的副本，返回的指针指向副本内存；函数返回后副本被销毁，指针变为野指针**，后续Eigen::Map操作野指针触发未定义行为，导致Matlab崩溃。

具体细节：

• 调用getOutDataPtr<double>(outputs[0])时，参数是值传递，arr是outputs[0]的完整副本。
• matlab::data::TypedArray<T> arr_t = arr;又生成了副本的TypedArray包装。
• 你通过迭代器获取的指针指向这个临时副本的内存，函数返回后arr和arr_t被销毁，内存释放，ptrRecon成为野指针，操作它必然引发崩溃。
• 额外问题：用it.operator->()获取基本类型指针的方式不规范，正确做法是用&*it取元素地址。

修改getOutDataPtr为引用传递，直接操作原数组并获取合法指针：

方法1：通过迭代器取地址（兼容全版本API） #

template <typename T>
T* getOutDataPtr(matlab::data::Array& arr) {
    // 将Array转换为TypedArray引用，避免拷贝
    auto& typed_arr = static_cast<matlab::data::TypedArray<T>&>(arr);
    // 获取指向第一个元素的合法指针
    return &*typed_arr.begin();
}

方法2：使用data()方法（Matlab R2020b及以上版本支持） #

template <typename T>
T* getOutDataPtr(matlab::data::Array& arr) {
    auto& typed_arr = static_cast<matlab::data::TypedArray<T>&>(arr);
    // 直接调用API提供的data()方法获取底层指针
    return typed_arr.data();
}

调用逻辑无需大幅修改，此时指针指向outputs[0]的真实内存空间：

// Initialize input parameters
long int numEl = inputs[0][0];
long int na = inputs[0][1];
long int nPoints = inputs[0][2];
        
// Preallocate output and define map to output array
outputs[0] = factory.createArray<double>({static_cast<size_t>(nPoints), static_cast<size_t>(numEl*na)});
        
auto ptrRecon = getOutDataPtr<double>(outputs[0]);
// Eigen::MatrixXd默认是列优先，和Matlab数组存储顺序完全匹配
Eigen::Map<Eigen::MatrixXd> Recon(ptrRecon, nPoints, numEl*na);

for(int i = 0; i < nPoints; i++) {
    for(int n = 0; n < na; n++) {
        for (int j = 0; j < numEl; j++) {
            Recon(i, n*numEl + j) = 5.0;
        }
    }    
}

• Matlab数组和Eigen默认的MatrixXd都采用列优先存储，因此内存布局完全对齐，无需额外指定存储顺序。
• 如果需要使用行优先的Eigen矩阵，构造Map时需显式声明：

  Eigen::Map<Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor>> Recon(ptrRecon, nPoints, numEl*na);
  

  此时要注意数组索引的对应关系，避免数据错位。

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