Scipy中具有很多重复元素的稀疏矩阵

Question

我有一个维数为n的三角稀疏矩阵，其中出现的唯一系数是c_1，c_2，... c_n。矩阵中最多有n个重复的单个系数。有什么方法可以利用矩阵中有很多相同的元素并且使用更少的内存的事实，而不是将相同的值分配给许多不同的（行，列）并且没有利用许多相同的事实。最后，我想将所述矩阵的逆应用于任意向量。

谢谢！

Answer 1

绝对可以使用'scipy.sparse'中的稀疏矩阵数据结构。如果您知道这些值倾向于按特定模式分组，则可以使用针对该模式优化的特定稀疏矩阵结构。如果沿对角线的值组使用“ dia_matrix”，列使用“ csc_matrix”，行使用“ csr_matrix”，块使用“ bsr_matrix”。