最大稀疏矩阵NumPy

Question

我正在研究一个模型，该模型反映了创新在N个人社会中的传播。我有一个大小为NxN的邻接矩阵A，注意这个邻接矩阵是稀疏的。

我想为NxN=10^7做模拟。我首先尝试了Matlab，但不幸的是Matlab不能处理NxN>10^4。

可以使用NumPy进行模拟吗？

Answer 1

Scipy可以用N=10^7处理NxN稀疏矩阵

import scipy.sparse as sparse
N = 10e7
sparse.bsr_matrix((N, N))

输出：

<100000000x100000000 sparse matrix of type '<class 'numpy.float64'>'
    with 0 stored elements (blocksize = 1x1) in Block Sparse Row format>

它是否适用于您的模拟取决于许多我们不知道的事情。您可能需要使用不同的sparse matrix class。