有没有办法修改特征切片？

问题描述

我有一个Eigen::Matrix<double,M,N> matrix;，我有一个只需要3x3切片并对其进行修改的函数

void ModifyBlock(Eigen::Matrix<double,3,3> & block){
   // A stupid example
   block += Eigen::Matrix<double,3>::Identity();
}

无法传递原始矩阵的一部分，因为以下代码无法编译：

int main(int argc,char *argv[]){
    const int M = 10;
    const int N = 20;

    Eigen::Matrix<double,N> matrix;
    // Initialize matrix
    ModifyBlock(matrix(vector<int>{1,7},vector<int>{0,2,4});
    return 0;    
}

cannot bind non-const lvalue reference of type ‘Eigen::Matrix<double,3>&’ to an rvalue of type ‘Eigen::Matrix<double,3>

我想知道是否有办法做到这一点。

在documentation中，我找不到提到切片必须是常量引用的情况。

解决方法

问题不是切片不可写，但您不能将它们转换为Eigen::Matrix的左值引用。使用C ++ 11，您可以简单地将它们作为右值参考模板传递：

template<class Derived>
void ModifyBlock(Eigen::MatrixBase<Derived> && block){
   // A stupid example
   block += Eigen::Matrix<double,3,3>::Identity();
}

c++11 eigen3