从 MatrixWrapper 创建 EigenSolver

Question

如何从 MaxtrixWrapper 构建 EigenSolver？

测试（也在 godbolt.org）

#include <Eigen/Eigen>
using namespace Eigen;

template<typename D>
void f(const Eigen::DenseBase<D>& a) {
    const Eigen::MatrixWrapper<const D> a2(a.derived());
    Eigen::EigenSolver<typename Eigen::MatrixWrapper<const D>> 
    es(a2);
}

int main() {
    ArrayXXf a(3, 3);
    a = 1.0f;
    f(a);
}

第一个错误：

<...>/eigen3/Eigen/src/Eigenvalues/EigenSolver.h:71:10: error: 
‘Options’ is not a member of 
‘Eigen::EigenSolver<
    Eigen::MatrixWrapper<
        const Eigen::Array<float, -1, -1> 
    > 
 >::MatrixType {
 aka Eigen::MatrixWrapper<const Eigen::Array<float, -1, -1> >}’
     enum {

Answer 1

你不知道。求解器都需要一个普通的 Matrix<...>（或 Ref<Matrix<...> >）作为模板参数。您可以使用以下方法获得正确的 Matrix 类型：

template<typename D>
void f(const Eigen::DenseBase<D>& a) {
    Eigen::EigenSolver<typename D::PlainMatrix> es(a.derived().matrix());
}

.derived().matrix() 实际上在这里是可选的，因为 ArrayXXf 隐式转换为 MatrixXf。（godbolt 在此超时——EigenSolver 对编译器来说相当繁重）。

从 MatrixWrapper 创建 EigenSolver

Create EigenSolver from MatrixWrapper

eigen3