1

我下载了 Eigen (3) 库并开始使用它。我写了一个模板函数，并在函数内部声明了一个“模板类型”的局部变量。我收到以下编译错误。

$ g++ EigenTest.cpp

EigenTest.cpp: In instantiation of ‘void myFunc(Eigen::MatrixBase<Derived>&) [with Type1 = Eigen::Matrix<double, -1, -1>]’:
EigenTest.cpp:24:10:   required from here
EigenTest.cpp:16:26: error: conversion from ‘Eigen::DenseCoeffsBase<Eigen::Matrix<double, -1, -1>, 1>::Scalar {aka double}’ to non-scalar type ‘Eigen::Matrix<double, -1, -1>’ requested
   Type1 tmp = matrix(0, 0);

“EigenTest.cpp”如下所示。

#include "Eigen/Dense"

#include <iostream>

template<typename Type1>
void myFunc(Eigen::MatrixBase<Type1>& matrix)
{
int i=matrix.rows();
Type1 tmp = matrix(0, 0);           // getting compiler error here
std::cout<<"tmp is ->"<<tmp<<std::endl;
}

int main()
{
Eigen::MatrixXd m(2,2);
m.setConstant(100); 
myFunc(m);
return 0;
}

我也尝试使用'typename Type1 tmp = matrix(0, 0);'
这也没有用！

如何解决这个问题？在普通的 C++ 模板编程中（没有 Eigen），我可以在模板函数中定义一个局部变量为 'Type1 tmp;"

4

3 回答 3

3

In Eigen::MatrixBase<Type1>,Type1不是标量类型，而是实际表达式的类型。在您的示例中，它将是 MatrixXd，但如果调用了 myFunc，例如 m.block(...)，那么Type1它将是一个 Block<...>。要获取标量类型，可以使用 Type1::Scalar：

template<typename Type1>
void myFunc(Eigen::MatrixBase<Type1>& matrix)
{
  typename Type1::Scalar Scalar;
  Scalar tmp = matrix(0, 0);
}

如果您需要类似于 的矩阵类型Type1，请使用Type1::PlainObject，例如：

typename Type1::PlainObject mat = 2 * matrix * matrix.transpose();
于 2013-10-31T14:43:59.843 回答
1

看起来像MatrixBase使用“CRTP”（见这里），模板参数实际上是从它派生的类型。因此在您使用的方法myFunc()中，Type1实际上是表示Eigen::MatrixXd，而我认为您认为是Type1双重的。所以，这一行：

Type1 tmp = matrix(0, 0);

在这个库的文档中（见这里），typedef forMatrixXd是一个双精度矩阵，所以我猜返回 frommatrix(0, 0)是一个双精度数，并且由于 tmp 是Type1其中的Eigen::MatrixXd，一个不会进入另一个。

扫描文档我认为您的函数将 aMatrix作为参数可能会更好，这样标量类型应该可用。像这样的东西：

template<class T, int rows, int cols, int opts, int maxR, int maxC > 
void myFunc( Eigen::Matrix<T, rows, cols, opts, maxR, maxC>& matrix )
{
    T tmp = matrix(0, 0);
}

（虽然看起来很可怕！！！;-)）

于 2013-10-31T09:20:12.920 回答
0

在您的代码中，Type1推断为double（因为Eigen::MatrixXd以这种方式定义）。

然后你试图做

Type1 tmp = matrix(0, 0);

而且我担心我的 Eigen 知识还不够，所以我通过 Clang 3.3 运行它，并得到了这个错误：

test.cpp:9:7: error: no viable conversion from 'Scalar' (aka 'double') to
      'Eigen::Matrix<double, -1, -1, 0, -1, -1>'
Type1 tmp = matrix(0, 0);           // getting compiler error here
      ^     ~~~~~~~~~~~~
test.cpp:17:1: note: in instantiation of function template specialization
      'myFunc<Eigen::Matrix<double, -1, -1, 0, -1, -1> >' requested here
myFunc(m);
^
/usr/include/eigen3/Eigen/src/Core/Matrix.h:210:5: note: candidate constructor not viable:
      no known conversion from 'Scalar' (aka 'double') to
      'internal::constructor_without_unaligned_array_assert' for 1st argument
    Matrix(internal::constructor_without_unaligned_array_assert)
    ^
/usr/include/eigen3/Eigen/src/Core/Matrix.h:284:25: note: candidate constructor not
      viable: no known conversion from 'Scalar' (aka 'double') to 'const
      Eigen::Matrix<double, -1, -1, 0, -1, -1> &' for 1st argument
    EIGEN_STRONG_INLINE Matrix(const Matrix& other)
                        ^
/usr/include/eigen3/Eigen/src/Core/Matrix.h:272:25: note: candidate template ignored:
      could not match 'MatrixBase<type-parameter-0-0>' against 'double'
    EIGEN_STRONG_INLINE Matrix(const MatrixBase<OtherDerived>& other)
                        ^
/usr/include/eigen3/Eigen/src/Core/Matrix.h:292:25: note: candidate template ignored:
      could not match 'ReturnByValue<type-parameter-0-0>' against 'double'
    EIGEN_STRONG_INLINE Matrix(const ReturnByValue<OtherDerived>& other)
                        ^
/usr/include/eigen3/Eigen/src/Core/Matrix.h:303:25: note: candidate template ignored:
      could not match 'EigenBase<type-parameter-0-0>' against 'double'
    EIGEN_STRONG_INLINE Matrix(const EigenBase<OtherDerived> &other)
                        ^
1 error generated.

这告诉我你不能那样调用matrix，用两个 0 作为参数。这也是一种奇怪的语法，因为MatrixBase该类没有operator()您似乎试图调用的。

于 2013-10-31T09:19:57.103 回答