c++ - vector 元素类型的 Eigen 稀疏矩阵CoeffRef的误差

Question

我将稀疏矩阵定义为Eigen::SparseMatrix< Eigen::Matrix<float, 3, 1> >，这意味着矩阵的每个元素都是3x1 vector 。
但是，当我调用函数CoeffRef以便将 vector 分配给元素时，在SparseMatrix.h中出现以下错误:

no operator "=" matches these operands. 

并且错误来自函数insert，同时它为Eigen::Matrix< float, 3, 1>分配了一个int到m_data.value(p) = 0(考虑到m_data.value(p)是 vector 3x1，而0是一个int)。

似乎在这段代码中(SparseMatrix.h的1235行)，他们没有考虑矩阵的模板类型进行比较。

我想知道您是否有解决这种错误的想法？

typedef Eigen::Matrix< float, 3, 1> Vec3f;

Eigen::SparseMatrix< Vec3f > lA( m, n);
lA.reserve( Eigen::VectorXi::Constant(m, 4) );

for( unsigned int i = 0; i < m; i++)
{

   Vec3f lVec( 0.0, 0.0, 1.0);

   lA.coeffRef(i, i) = lVec; // got the error here!
}

Answer 1

这是因为coeffRef尝试将新创建的元素初始化为0，但不能将0分配给Vector3f。因此解决方案是改用Array3f:

typedef Eigen::Array<float, 3, 1> Vec3f;

当然，请注意，operator *在Array上的行为不同于在 vector 和矩阵上的行为。