有效地更新特定特定稀疏基质的某些块的内置块？

Question

假设我的稀疏矩阵具有以下模式：

• 每列的非齐射率及其位置仅固定，
• 仅矩阵块A和b会改变，其余矩阵保持静态； （A块A和B本身在固定的非零位置也很稀疏）

如，我已经通过

• 的每个列的确切数量
• 保留“列列稀疏矩阵插入列”列
• 最小行索引的列插入列

插入程序后面的 矩阵，仅更新A，B块Inploph。可能的方法是：

1. 通过coeffref访问现有条目，将引入二进制搜索，因此在这里不受欢迎。
2. 如有记录的在这里

似乎有点，它似乎有点 不需要在所有非零条目上进行迭代，因为稀疏矩阵的大部分时间都保持不变。

是否可以在矩阵中的所有非Zeros上进行迭代？

Answer 1

据我所知，Inneriterator可以用于此过程并在恒定时间内运行。

Eigen::Index col = 1;
Eigen::Index offset_in_col = 1;

using SparseMatrixD = Eigen::SparseMatrix<double>;
SparseMatrixD mat = ...;

SparseMatrixD::InnerIterator i =
      SparseMatrixD::InnerIterator(mat, col) + offset_in_col;
assert(i.row() == 1);
assert(i.col() == 1);
assert(i.value() == C);

这应该访问值C。您需要知道的只是每个列的非零元素有多少个（或一般内在维度）。您不需要知道存储多少个非零列（外部维度），因为该数组（sparsematrix.outerindexptr（））每列有一个条目。