将特征矩阵映射到 C 数组

Question

我最近开始使用Eigen 库。我有一个将特征矩阵映射到 C/C++ 数组的问题。默认情况下，特征矩阵是列专业的。因此，如果我使用以下代码将矩阵映射到 C/C++ 数组，

double a[10];
double *p = &a[0];
MatrixXd(2,5) m;
for (int i=0; i<2;i++)
    for (int j=0; j<5;j++)
        m(i,j) = i+j;
cout<<m<<endl;
Eigen::Map<MatrixXd>(p,2,5) = m;
for (int i=0; i<10; i++)
    cout<<a[i]<<" ";
cout<<endl;

输出是：

0 1 2 3 4
1 2 3 4 5
0 1 1 2 2 3 3 4 4 5

如果我将 m 的定义更改为 row majored：

Matrix <double,2,5,RowMajor> m;

我希望输出如下所示：

0 1 2 3 4
1 2 3 4 5
0 1 2 3 4 1 2 3 4 5

但实际上结果还是和第一个一样。我的问题是，有没有办法将 Eigen 矩阵映射到 C/C++ 数组，以便数组的数据是基于行的？

我发现我可以使用 matrix.data() memember 函数来获得所需的结果，但我想知道我是否可以使用 map 来做到这一点：

使用 matrix.data() 作品：

double a[10];
double *p = &a[0];
Matrix <double,2,5,RowMajor> m;
for (int i=0; i<2;i++)
    for (int j=0; j<5;j++)
        m(i,j) = i+j;
double *p1 = m.data();
for (int i=0; i<10; i++)
    cout<<p1[i]<<" ";
cout<<endl;

Answer 1

重要的不是矩阵m 的类型，而是Map 模板中使用的类型。您必须将 Map 模板中使用的类型更改为行专业。

Eigen::Map<Matrix<double,2,5,RowMajor> >(p,2,5) = m;