eigen3，如何用密度和稀疏的矩阵实现多态性？

Question

我想构建一个阵列，其中包含大型稀疏和密集矩阵的指针。然后，我希望使用多态性对它们进行操作，并避免复制/分配大部分内存。例如，

std::vector<UnkownType*> matrices;

matrices.emplace_back(get_large_dense_matrix_ptr());
matrices.emplace_back(get_large_sparse_matrix_ptr());

Eigen::MatrixXd result(known_size, known_size);
for(mat : matrices){
   result+=*mat;
}

unkowntype应该是什么？ 我认为选择使用一系列参考文献或参考包装器，但这是禁止的。

Answer 1

如果您有C ++ 17，则可以使用 std :: variant

using matrix_variant = std::variant<Eigen::MatrixXd, Eigen::SparseMatrix<double>>;

std::vector<matrix_variant> matrices;

matrices.emplace_back(Eigen::MatrixXd::Zeros(10, 10).eval());
matrices.emplace_back(Eigen::SparseMatrix<double>{});

Eigen::MatrixXd result;
for(const matrix_variant& matrix: matrices) {
    if(std::holds_alternative<Eigen::MatrixXd>(matrix)) {
        const Eigen::MatrixXd& dense = std::get<Eigen::MatrixXd>(matrix);
        result += dense;
    } else if(std::holds_alternative<Eigen::SparseMatrix<double>>(matrix)) {
        const Eigen::SparseMatrix<double>& sparse =
              std::get<Eigen::SparseMatrix<double>>(matrix);
        result += sparse;
    }
}