稀疏矩阵上的块运算 - Eigen Toolbox - C++

Question

Block operations for sparse matrices - Eigen Toolbox - C++ 


#include "Eigen/Dense"
#include "Eigen/Sparse"
#include <iostream>
using namespace std;
using namespace Eigen;

    int main()
    {
    MatrixXd silly(6, 3);
    silly << 0, 1, 2,
            0, 3, 0,
            2, 0, 0,
            3, 2, 1,
            0, 1, 0,
            2, 0, 0;

        SparseMatrix<double> sparse_silly,temp;
        sparse_silly= Eigen::SparseMatrix<double>(6, 3);
        temp = Eigen::SparseMatrix<double>(6, 3);
        sparse_silly = silly.sparseView();

        std::cout << "Whole Matrix" << std::endl;
        std::cout << sparse_silly << std::endl;



         temp.block(0, 0, 2, 2)=sparse_silly.block(0, 0, 2, 2);

        std::cout << "block of matrix" << std::endl;
        std::cout << temp.block(0, 0, 2, 2) << std::endl;

        getchar();
        return 0;
    }

在上面的代码中，稀疏矩阵的块运算无法使用 Eigen 工具箱进行。我想要将稀疏愚蠢的块分配给临时矩阵中的块。临时矩阵的打印输出为零。如果我在概念上遗漏了某些内容，任何人都可以帮助我。最近的文档称块操作可用于稀疏矩阵。

Answer 1

Eigen 中的稀疏矩阵块并非全部可写。某些是（例如.col(Index)在列主矩阵中），但通用.block(Index, Index, Index, Index)不是。文档在这方面相当令人困惑，但如果仔细观察，所有示例都是密集矩阵而不是稀疏矩阵。该col()文档也使用了密集矩阵示例，但如果您尝试一下，您会发现它有效。