我想将double数组映射到现有的MatrixXd结构.到目前为止,我已经设法将Eigen矩阵映射到一个简单的数组,但我找不到回来的方法.
- void foo(MatrixXd matrix,int n){
- double arrayd = new double[n*n];
- // map the input matrix to an array
- Map<MatrixXd>(arrayd,n,n) = matrix;
- //do something with the array
- .......
- // map array back to the existing matrix
- }
解决方法
我不确定你想要什么,但我会试着解释一下.
你在代码中混合了double和float(MatrixXf是一个矩阵,每个条目都是一个浮点数).我现在假设这是无意的,你想在任何地方使用双倍;如果这真的是你的意图,请参见下文.
指令Map< MatrixXd>(arrayd,n)=矩阵将矩阵的条目复制到arrayd中.它相当于循环
- for (int i = 0; i < n; ++i)
- for (int j = 0; j < n; ++j)
- arrayd[i + j*n] = matrix(i,j);
要将arrayd的条目复制到矩阵中,可以使用反向赋值：matrix = Map< MatrixXd>(arrayd,n).
但是,通常以下技术更有用：
- void foo(MatrixXd matrix,int n) {
- double* arrayd = matrix.data();
- // do something with the array
- }
现在arrayd指向矩阵中的条目,您可以将其作为任何C数组进行处理.矩阵和arrayd之间共享数据,因此您不必在最后复制任何内容.顺便说一下,你不需要将n传递给函数foo(),因为它存储在矩阵中;使用matrix.rows()和matrix.cols()来查询其值.
如果您确实要将MatrixXf复制到双精度数组,则需要显式包含强制转换. Eigen中的语法是：Map< MatrixXd>(arrayd,n)= matrix.cast< double>().
