#include <Eigen/Dense>
#include <iostream>

int indexing() {
    Eigen::ArrayX2d M(3, 2);
    M << 0, 1, 2, 3, 4, 5;
    std::array<int, 3> i = {2, 0, 1};
    const auto& j = i;
    std::cout << M(j, Eigen::all) << std::endl;
    return 0;
}

int main() {
    Eigen::ArrayXi row(5);
    row << 0, 1, 2, 3, 4;
    Eigen::ArrayX2i zeros(5, 2);
    zeros.fill(0);
    Eigen::ArrayX2i ones(5, 2);
    ones.fill(1);
    auto res = (row > 2).replicate<1,2>().select(zeros, ones);
    std::cout << res << std::endl;
    return 0;
}