如何更改特征块？

Question

我无法重新分配Block。在下面的代码我的矩阵A存储在两种不同的方式：

1. 为3 ArrayXd S，一个用于每个行
2. 作为ArrayXXd

。

// data 

ArrayXXd A (3, 3); 
A << 0, 1, 2, 3, 4, 5, 6, 7, 8; 

std::vector<ArrayXd> A_rows = {A.row(0), A.row(1), A.row(2)}; 

// std::vector<ArrayXd> solution 

// first row 
ArrayXd & current_row = A_rows[0]; 
// read it, write it, do stuff 
// start working with the second row 
current_row = std::ref(A_rows[1]); 
cout << current_row << endl << endl; // prints 3 4 5 
cout << A << endl; // A is unchanged 

// Eigen solution 

// first row 
Block<ArrayXXd, 1, -1> && current_row_block = A.row(0); 
// read it, write it, do stuff 
// start working with the second row 
current_row_block = std::ref(A.row(1)); // doesn't compile 
cout << current_row_block << endl; 
cout << A << endl; 

的错误信息是：

error: use of deleted function 'void std::ref(const _Tp&&) [with _Tp = Eigen::Block<Eigen::Array<double, -1, -1>, 1, -1, false>]' 
current_row_block = std::ref(A.row(1)); 
            ^

是否有可能以固定第二方法或我应该移动到存储所述矩阵为std::vector<ArrayXd>？

相关问题：Passing a reference of a vector element to a threaded function

Answer 1

你并不需要一个Block<...>引用一行。你只需要一个索引。

int current_row_id = 0; 
std::out << A.row(current_row_id) << std::end; 
current_row_id = 1; 
std::out << A.row(current_row_id) << std::end; 

为了您std::vector<ArrayXd>的方法，因为你是让行的副本，你不能改变原有的A。