julia - 如何在稀疏矩阵中找到所有空行的索引，

Question

我有一个大的稀疏矩阵 M。我想找到矩阵中所有空行的索引。你怎么能那样做？

标签： julia

Answer 1

Julia 使用压缩稀疏列 (CSC)格式进行稀疏矩阵存储，这意味着所有存储值的行索引都可用。1:NROWS因此，您可以通过获取行索引集和行索引集之间的集差来找到所有没有存储值的行：

julia> using SparseArrays

julia> A = rand(10, 10); A[3,:] .= 0; A[5,:] .= 0; S = sparse(A);

julia> idx = setdiff(Set(1:size(A, 1)), Set(S.rowval))
Set{Int64} with 2 elements:
  3
  5