Question

我一辈子都想不出为什么这不正确。它似乎没有返回第k个元素。

typedef Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic> matrix;

double test(matrix& D, int k)
{
    au到 d = D.row(1);
    std::nth_element(d.data(),d.data()+k, d.data()+d.size());
    return d(k) ;
}

我也试过

template <typename ScalarType, typename Derived>
void Sort(Eigen::MatrixBase<Derived> &xValues)
{
std::sort(xValues.derived().data(), xValues.derived().data()+xValues.derived().size());
}


double test(matrix& D, int k)
{
    au到 d = D.row(1);
    Sort<double>(d);
    return d(k) ;
}

非常感谢您的帮助。

编辑

我刚试着改变

au到 d = D.row(1);

到

Eigen::Vec到rXd rowD = D.row(1);
....

and it seems 到 work ok.

对此有点困惑。

Answer 1

特征矩阵是column-major默认值。这意味着矩阵的一行不是连续的C数组，并且不能将数据指针用作迭代器。

例如，3x4矩阵将被存储为：

0 3 6  9
1 4 7 10
2 5 8 11

现在，行（1）将是

1 4 7 10

但是您传递给nth_element（）的指针迭代器将访问

1 2 3 4

如果您将矩阵typedef更改为row major，则代码有效：

typedef Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor> matrix;

更新：由于将行复制到向量中，因此您编辑的示例可以正常工作。对于向量（一维矩阵），数据是以行为主还是以列为主存储并不重要。

友情链接