特征块不是左值？

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

该block方法不返回Matrix对象而是返回Block<>表达式.在这种情况下,最好使用auto关键字(但要小心):

auto column = m.col(0);

或者使用方便的typedef:

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

甚至是Ref<>:

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