c++ - 指向特征矩阵的指针数组

Question

我在我的代码中使用来自 Eigen 的 MatrixXd 矩阵，在某个时候我需要一个 3D 矩阵。由于 Eigen 没有三维矩阵类型，因为它仅针对线性代数进行了优化，因此我创建了一个 MatrixXd 类型的指针数组：

Eigen::MatrixXd *CVM =new Eigen::MatrixXd[100];
for (int i = 0; i < 100; i++){
   CVM[i]= Eigen::MatrixXd::Zero(5,5);
}

但是，稍后我需要访问这个数组上的值，为此我正在做类似的事情：

for (int k = 0; k < 100; k++){
   Eigen::MatrixXd* b=&CVM[k];

   for (int i = 0; i < 5; i++){
      for (int j = 0; j < 5; j++){
         b->coeff(i,j)=47;
      }      
   }
}

由于b是指针而不是MatrixXd本身，b(i,j)显然不起作用，所以我使用的是该coeff()方法，但是，我收到以下错误：

error: assignment of read-only location ‘b->Eigen::Matrix<double, -1, -1>::<anonymous>.Eigen::PlainObjectBase<Derived>::coeff<Eigen::Matrix<double, -1, -1> >(((Eigen::PlainObjectBase<Eigen::Matrix<double, -1, -1> >::Index)i), ((Eigen::PlainObjectBase<Eigen::Matrix<double, -1, -1> >::Index)j))’

编辑：添加输出

cout << b << endl;
cout << CVM[47] << endl;


0x1c34b00
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0

Answer 1

要么使用b->coeffRef(i,j)获取读/写引用，要么取消引用 b: (*b)(i,j)，或者使用 b: 的引用

MatrixXd& b = CVM[k];
b(i,j) = 47;

Answer 2

只需使用 operator()(int,int)

CVM[k].operator()(i,j) = 47;

或者

CVM[k](i,j) = 47;

或者

Eigen::MatrixXd* b = &CVM[k];
b->operator()(i,j) = 47;

这里k是矩阵，i是行，j是列。