特征布尔数组切片

Question

在MATLAB中，通常从矩阵/数组(称为逻辑索引)中切片出满足某些条件的值。

vec = [1 2 3 4 5];
condition = vec > 3;
vec(condition) = 3;

我如何在Eigen中做到这一点?到目前为止，我有:

Eigen::Matrix<bool, 1, 5> condition = vec.array() > 3;

Answer 1

试试这个:

#include <iostream>
#include <Eigen/Dense>
int main()
{
    Eigen::MatrixXi m(1, 5);
    m << 1, 2, 3, 4, 5;
    m = (m.array() > 3).select(3, m);
    std::cout << m << std::endl;
    return 0;
}

Answer 2

正如这里一个类似问题的答案所指出的那样:使用Eigen的子矩阵和索引，libigl将此功能添加到Eigen中。

igl::slice(A,indices,B);

等价于

B = A(indices)