c – 迭代特征中的对称稀疏矩阵

我按照以下示例循环遍历稀疏矩阵的元素.

SparseMatrix<double> mat(rows,cols);
for (int k=0; k<mat.outerSize(); ++k)
  for (SparseMatrix<double>::InnerIterator it(mat,k); it; ++it)
  {
    it.value();
    it.row();   // row index
    it.col();   // col index (here it is equal to k)
    it.index(); // inner index, here it is equal to it.row()
  }

这个例子我完全没问题.但是,我所拥有的矩阵是对称的,我只想通过下半部分进行迭代.有没有简单的方法来迭代对称矩阵?

最佳答案 您可以在每个循环的开头检查它:

for (int k=0; k<mat.outerSize(); ++k)
  for (SparseMatrix<double>::InnerIterator it(mat,k); it; ++it)
  {
      if(it.row() < it.col())
          continue;
  }
点赞