本征块不是左值？

Question

我正在使用 Eigen 库及其块操作，根据 https://eigen.tuxfamily.org/dox/group__TutorialBlockOperations.html 块是左值和右值。但是以下程序无法编译：

#include "Eigen/Dense"
#include <iostream>

int main() {
  Eigen::Matrix<float, Eigen::Dynamic, Eigen::Dynamic> m(4, 4); 
  m << 1, 2, 3, 4,
       5, 6, 7, 8,
       9, 10, 11, 12, 
       13, 14, 15, 16; 
  Eigen::Matrix<float, Eigen::Dynamic, Eigen::Dynamic> &column = m.col(0);
}

在这个程序中，我想使用列作为对矩阵 m 的引用，这样我就可以就地更改元素。错误信息是：

test_eigen_rvalue.cc: In function ‘int main()’:
test_eigen_rvalue.cc:11:73: error: invalid initialization of non-const reference of type ‘Eigen::Matrix<float, -1, -1>&’ from an rvalue of type ‘Eigen::DenseBase<Eigen::Matrix<float, -1, -1> >::ColXpr {aka Eigen::Block<Eigen::Matrix<float, -1, -1>, -1, 1, true>}’
   Eigen::Matrix<float, Eigen::Dynamic, Eigen::Dynamic> &column = m.col(0);

问题似乎是 m.col(0) 不是左值。

Answer 1

您也可以尝试使用块操作，如下所示：https://eigen.tuxfamily.org/dox/group__TutorialBlockOperations.html

例如：&column = m.block<4,1>(0,0);

Answer 2

block 方法不是 return 一个 Matrix 对象而是一个 Block<> 表达式。在这种情况下，最好使用 auto 关键字（但要 careful）：

auto column = m.col(0);

或者使用方便的 typedefs:

MatrixXf::ColXpr column = m.col(0);

甚至 Ref<>:

Ref<VectorXf> column = m.col(0);

本征块不是左值？

Eigen block not lvalue?

c++

eigen