【eigen】解决报错 return type of RCmpEQReturnType Eigen::operator== is not ‘bool’

文章目录

- 1. 发现问题
- 2. GPT 分析问题
- 3. 解决问题

1. 发现问题

在 c++ 项目中使用了 C++ 20 编译标准（GCC10），结果 eigen 中出现报错

/usr/include/eigen3/Eigen/src/Cholesky/LDLT.h:372:58: error: return type of ‘const RCmpEQReturnType Eigen::operator==(const Scalar&, const Eigen::ArrayWrapper<Eigen::Block<Eigen::Block<Eigen::Matrix<double, -1, -1>, -1, 1, true>, -1, 1, false> >&)’ is not ‘bool’
372 |           ret = ret && (mat.col(j).tail(size-j-1).array()==Scalar(0)).all();

2. GPT 分析问题

错误提示指出在 ret = ret && (mat.col(j).tail(size-j-1).array()Scalar(0)).all(); 这行代码中，操作符对 Eigen::ArrayWrapper 和 Scalar 类型不可用。你可以修改这行代码以适应 Eigen 库的比较操作。

首先，确保 mat 是 Eigen::Matrix 类型，Scalar 是正确的类型，然后使用 .isApprox() 函数来比较两者是否近似相等。

3. 解决问题

Eigen/src/Cholesky/LDLT.h:517

ret = ret && (mat.col(j).tail(size-j-1).array()==Scalar(0)).all();

改为

// ret = ret && (mat.col(j).tail(size-j-1).array()==Scalar(0)).all();
ret = ret && (mat.col(j).tail(size-j-1).isApprox(Eigen::Matrix<Scalar, Eigen::Dynamic, 1>::Zero(size-j-1), 1e-6));