c++ - Eigen block 不是左值？

Question

根据 https://eigen.tuxfamily.org/dox/group__TutorialBlockOperations.html，我正在使用 Eigen 库及其 block 操作这些 block 都是左值和右值。但是下面的程序不能编译:

#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

block方法不返回 Matrix对象而是一个 Block<>表达。在这种情况下，最好使用 auto。关键字(但必须是 careful ):

auto column = m.col(0);

或者使用方便的 typedef:

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

甚至是 Ref<> :

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