我想把两个大的稀疏矩阵相乘。第一个是150,000 x 300,000,第二个是300,000 x 300,000。第一个矩阵具有大约1,000,000个非零项,而第二个矩阵具有大约20,000,000个非零项。有没有一种直接的方法来得到这些矩阵的乘积?我目前正在以csr或csc格式存储矩阵,并尝试使用matrix_a * matrix_b。这将产生错误ValueError: array is too big。我猜我可以用pytable将独立的矩阵存储在磁盘上,将它们分
假设我们有A(NxN)矩阵和列向量(B(Nx1))。我们还有一个函数f(i,j),它在i行和列j处返回矩阵A的元素。如果我们想做一些矩阵运算,例如A和B的矩阵乘积,我们可以使用以下方法(下面,C是矩阵乘积的结果):
returnfor j in range(N): for k in range(N):