从MatrixWrapper创建EigenSolver
从MatrixWrapper创建EigenSolver
提问于 2018-12-06 03:53:19
如何从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 {回答 1
Stack Overflow用户
回答已采纳
发布于 2018-12-06 04:59:03
求解器都需要一个简单的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对编译器来说太重了)。
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/53639778
复制相关文章
相似问题