如何在大稀疏矩阵中找到非零元素的索引？

Question

我有两个sq矩阵(a,b)的大小为100000 X 100000.我必须区分这两个矩阵(c = ab).结果矩阵'c'是稀疏矩阵.我想找到所有非零元素的索引.我必须多次做这个操作(> 100).

最简单的方法是使用两个for循环.但这是计算密集型的.你能告诉我任何算法或包/库最好在R/python/c中尽快做到这一点吗？

Answer 1

由于您有两个密集矩阵，因此双 for 循环是您唯一的选择。您根本不需要稀疏矩阵类，因为您只想知道 的索引(i,j)列表a[i,j] != b[i,j]。

在 R 和 Python 等语言中，双 for 循环的性能会很差。我可能会在双 for 循环的本机代码中编写此代码，并将索引添加到列表对象中。但毫无疑问，解释代码的向导（即 R、Python 等）知道无需诉诸本机编码即可实现此目的的有效方法。