python中大矩阵到稀疏矩阵的转换

Question

我需要做以下事情

b = numpy.random.randn(50001,2)

cof = numpy.corrcoef(b)

c= b>=0.3 

return np.dot(c, np.ones([50001,1]))

这是在给我一个分割错误。在

如果我尝试使用稀疏矩阵，例如：

^{pr2}$

我有一个分割错误

如果矩阵的大小很小，则转换是有效的。在

有什么建议吗？在

Answer 1

你是想计算50001维空间的两个样本之间的相关性还是二维空间的50001个样本之间的相关性？在

在您当前的情况下，您正在创建一个隐式50001 x 50001协方差矩阵（这会导致segfault）。做以下事情不会破坏记忆

b = numpy.random.randn(2,50001)
cof = numpy.corrcoef(b)

希望这正是您需要的（另一种方法产生的协方差矩阵和segfault非常差的近似值）。在