在 Eigen 2 库中使用 linearRegression 时模板参数 deduction/substitution 失败

Question

我在现有的 C++ 代码库中工作，但在调用函数时遇到问题。

我正在尝试使用代码库中已有的旧 Eigen 库中的线性回归函数，您可以看到 that source code here 但下面是相关声明。

template<typename VectorType>
void linearRegression(int numPoints,
                      VectorType **points,
                      VectorType *result,
                      int funcOfOthers )
{ ... }

下面是我的代码的匿名副本：

// MyClass has member variable int maxIdx = 5
void MyClass::myFunction()
{
    Eigen::Vector2d points[maxIdx];
    Eigen::Vector2d coeffs;
    for(int i = 0; i < maxIdx; i++)
    {
        points[i] = Eigen::Vector2d(data->mydata[i].x, data->mydata[i].y);
    }
    Eigen::linearRegression( maxIdx, &points, &coeffs, 1 );

    // do more stuff with coeffs
}

这是我在尝试编译时收到的错误消息：

myfile.cpp:803:67: error: no matching function for call to 'linearRegression(int, Eigen::Vector2d (*)[((MyClass*)this)->MyClass::maxIdx], Eigen::Vector2d*, int)'
         Eigen::linearRegression( maxIdx, &points, &coeffs, 1 );
                                                                   ^
myfile.cpp:803:67: note: candidate is:
In file included from [redacted]:
lib/Eigen/src/Eigen2Support/LeastSquares.h:85:6: note: template<class VectorType> void Eigen::linearRegression(int, VectorType**, VectorType*, int)
 void linearRegression(int numPoints,
      ^
lib/Eigen/src/Eigen2Support/LeastSquares.h:85:6: note:   template argument deduction/substitution failed:
myfile.cpp:803:67: note:   mismatched types 'VectorType*' and 'Eigen::Vector2d [((MyClass*)this)->MyClass::maxIdx] {aka Eigen::Matrix<double, 2, 1> [((MyClass*)this)->MyClass::maxIdx]}'

这几乎是库源代码中示例代码的精确副本，所以我有点不知道如何解决这个问题。我对模板不是很熟悉，所以类型错误可能与此有关？

Answer 1

根据 that 来源（与他们的示例中所写的不同），points 应该是指向点的指针数组。也就是说，

Eigen::Vector2d *points[maxIdx];

所以，

Eigen::Vector2d points[maxIdx];
Eigen::Vector2d *point_ptrs[maxIdx];
for(int i = 0; i < maxIdx; i++)
    point_ptrs[i] = &points[i];
...
Eigen::linearRegression(maxIdx, point_ptrs, &coeffs, 1);

这对 Vector2d 来说似乎很愚蠢，但可能是为难以复制的巨大向量设计的。

Answer 2

根据@numzero 的回答，该代码已编译，但当我尝试执行它时出现运行时崩溃。看起来第二个参数需要指向该指针数组中的 第一个索引 。解释了发生了什么，以及它在运行时失败，即使它已编译。

以下是最终对我有用的方法：

Eigen::Vector2d points[numPoints];
Eigen::Vector2d* point_ptrs[numPoints];
Eigen::Vector2d coeffs;
for(int i = 0; i < numPoints; i++)
{
    points[i] = Eigen::Vector2d(data->mydata[i].x, data->mydata[i].y);
    point_ptrs[i] = &points[i];
}
Eigen::linearRegression( maxIdx, &(point_ptrs[0]), &coeffs, 1 );

在 Eigen 2 库中使用 linearRegression 时模板参数 deduction/substitution 失败

Template argument deduction/substitution failed when using linearRegression in Eigen 2 library

c++

eigen