从MatrixWrapper创建EigenSolver

内容来源于 Stack Overflow,并遵循CC BY-SA 3.0许可协议进行翻译与使用

  • 回答 (1)
  • 关注 (0)
  • 查看 (81)

如何构建EigenSolver一个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 {
提问于
用户回答回答于

解算器都需要普通Matrix<...>(或a 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对于编译器来说非常重要)。

所属标签

可能回答问题的人

  • 优惠活动秘书

    0 粉丝2 提问8 回答
  • 天使的炫翼

    17 粉丝531 提问7 回答
  • 最爱开车啦

    8 粉丝503 提问6 回答
  • uncle_light

    5 粉丝518 提问5 回答

扫码关注云+社区

领取腾讯云代金券