内容来源于 Stack Overflow,并遵循CC BY-SA 3.0许可协议进行翻译与使用
如何构建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对于编译器来说非常重要)。