下面的代码片段来自一个Eigen断言:
MatrixXd L;
VectorXd x, b;
...
ASSERT_MATRIX_EQ(L*x, b);
与,
template <typename DerivedL, typename DerivedR>
void ASSERT_MATRIX_EQ(const Eigen::DenseBase<DerivedL>& A, const Eigen::DenseBase<DerivedR>& B, double tol=1e-7) {
ASSERT_EQ(A.rows(), B.rows());
ASSERT_EQ(A.cols(), B.cols());
for(int i=0; i < A.rows(); i++) {
for(int j=0; j < A.cols(); j++) {
ASSERT_NEAR(A(i,j), B(i,j), tol);
}
}
}
它死于错误:
test_leq: /usr/include/eigen3/Eigen/src/Core/ProductBase.h:154: typename Base::CoeffReturnType Eigen::ProductBase<Eigen::GeneralProduct<Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix<double, -1, 1, 0, -1, 1>, 4>, Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix<double, -1, 1, 0, -1, 1> >::coeff(Index, Index) const: Assertion `this->rows() == 1 && this->cols() == 1' failed.
在对
A(i,j)
的调用中。 (不过,我可以称呼cout << A << endl;
很好。)在第154行,
ProductBase.h
奇怪地拥有断言 // restrict coeff accessors to 1x1 expressions. No need to care about mutators here since this isnt a Lvalue expression
typename Base::CoeffReturnType coeff(Index row, Index col) const
{
#ifdef EIGEN2_SUPPORT
return lhs().row(row).cwiseProduct(rhs().col(col).transpose()).sum();
#else
EIGEN_STATIC_ASSERT_SIZE_1x1(Derived)
eigen_assert(this->rows() == 1 && this->cols() == 1);
return derived().coeff(row,col);
#endif
}
我遵循Eigen's guide编写通用矩阵函数。如何正确编写此泛型函数?
编辑:很高兴知道
ProductBase
为什么期望一个1x1矩阵。 最佳答案
本征邮件列表上的thread表示有意禁用ProductBase
上的系数访问。现在的解决方案是避免使用类似foobar(A*x)
的表达式。
关于c++ - Eigen :访问ProductBase系数,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/18169460/