代码之家 › 专栏 › 技术社区 › R zu

从MatrixWrapper创建Eigensolver

eigen3

R zu · 技术社区 · 5 年前

如何构造 EigenSolver 从A MaxtrixWrapper ?

#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 {

1 回复 | 直到 5 年前

chtz 5 年前

你不想。解决者都想要一个简单的 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对于编译器来说相当重)。

推荐文章

avocado · 特征稀疏矩阵的性能调整

7 年前

Daiver · 特征:沿一维复制项,而不进行无用的分配

7 年前

zer0hedge · RcppEigen。h: 安装期间没有此类文件或目录。。。Centos 7上的RcppEigen

7 年前

Daiver · 特征:反转矩形矩阵(最小范数解)

7 年前

German Capuano · 以特征类型为参数的函数的返回类型不明确

7 年前

Daiver · 在没有内存开销的情况下设置特征::稀疏矩阵的稀疏模式

7 年前

Olivier · 使用MKL编译时,Eigen C++运行速度较慢

7 年前

Aditya Kashi · 使用valuePtr()等为特征3稀疏矩阵分配空间

8 年前