为何我的Eigen-MEX代码输出的矩阵元素不正确？

从你给出的代码片段来看，大概率是类型不匹配或者Eigen SparseMatrix Map的构造参数错误导致的元素位置或数值异常，下面是几个最常见的坑和对应的解决方法：

1. 最可能的元凶：索引类型不匹配 #

Matlab的MEX API中，mxGetIr和mxGetJc返回的是mwIndex*类型（本质是size_t，64位平台上是无符号64位整数），但你直接强制转换为int*，这会引发两个问题：

• 64位Matlab环境下，当索引值超过32位int的范围时会被截断，直接破坏行/列索引的正确性
• 即使数值没溢出，类型不匹配会触发未定义行为，导致内存访问错误

解决方法：
在Eigen的Map<SparseMatrix>中指定正确的索引类型，直接使用mwIndex而非int，同时不要随意转换指针类型：

// 正确接收Matlab的稀疏矩阵索引
mwIndex *ir = mxGetIr(prhs[0]);
mwIndex *jc = mxGetJc(prhs[0]);

// 构造Eigen的稀疏矩阵Map，指定索引类型为mwIndex
Map<SparseMatrix<double, ColMajor, mwIndex>> eigen_mat(
    m,          // 矩阵行数
    n,          // 矩阵列数
    jc[n],      // 实际非零元素数（而非nzmax）
    jc,         // 列指针数组
    ir,         // 行索引数组
    pr          // 数值数组
);

2. 别把nzmax当成实际非零元素数 #

mxGetNzmax返回的是Matlab为稀疏矩阵分配的存储空间大小，而非实际存在的非零元素数！实际非零元素数应该是jc[n]（Matlab的jc数组最后一个元素就是总非零数）。如果误用nzmax，会读取到未初始化的内存，产生垃圾值。

3. 确认存储顺序和Matlab一致 #

Matlab的稀疏矩阵是**列优先（CSC）**存储，Eigen的SparseMatrix默认也是ColMajor，但如果你不小心指定了RowMajor，就会完全打乱元素的位置。一定要确保Map的第二个模板参数是ColMajor（或省略，使用默认值）。

4. 输出回Matlab时的注意事项 #

如果后续要把Eigen的稀疏矩阵转回Matlab格式，也要保持类型匹配：

// 创建Matlab稀疏矩阵
plhs[0] = mxCreateSparse(m, n, jc[n], mxREAL);
double* out_pr = mxGetPr(plhs[0]);
mwIndex* out_ir = mxGetIr(plhs[0]);
mwIndex* out_jc = mxGetJc(plhs[0]);

// 拷贝数据时确保类型一致
memcpy(out_pr, eigen_mat.valuePtr(), jc[n] * sizeof(double));
memcpy(out_ir, eigen_mat.innerIndexPtr(), jc[n] * sizeof(mwIndex));
memcpy(out_jc, eigen_mat.outerIndexPtr(), (n+1) * sizeof(mwIndex));

快速排查小技巧 #

可以先打印几个关键值验证：

• 对比jc[n]和nz（mxGetNzmax的结果），如果不一致，说明误用了存储空间大小而非实际非零数
• 打印前几个ir和jc的值，和Matlab中原始稀疏矩阵的索引对比，看是否匹配

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