如何删除Eigen3矩阵中某些不连续的行和列？

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

int main() {
    Eigen::MatrixXd matrix(4, 4);
    matrix << 1, 2, 3, 4,
              5, 6, 7, 8,
              9, 10, 11, 12,
              13, 14, 15, 16;

    Eigen::VectorXi rowsToRemove(2);
    rowsToRemove << 1, 3;

    Eigen::VectorXi colsToRemove(2);
    colsToRemove << 0, 2;

    int numRows = matrix.rows() - rowsToRemove.size();
    int numCols = matrix.cols() - colsToRemove.size();

    Eigen::MatrixXd newMatrix(numRows, numCols);

    int newRow = 0;
    for (int row = 0; row < matrix.rows(); ++row) {
        if (rowsToRemove.contains(row)) {
            continue;
        }

        int newCol = 0;
        for (int col = 0; col < matrix.cols(); ++col) {
            if (colsToRemove.contains(col)) {
                continue;
            }

            newMatrix(newRow, newCol) = matrix(row, col);
            ++newCol;
        }

        ++newRow;
    }

    std::cout << "Original Matrix:" << std::endl << matrix << std::endl;
    std::cout << "New Matrix:" << std::endl << newMatrix << std::endl;

    return 0;
}