最大尺寸稀疏矩陣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>

它是否適合您的模擬取決於許多我們不知道的事情。 您可能需要使用其他稀疏矩陣類 。