Question

如何存储和计算超大规模稀疏矩阵

Answer 1

代码如下：
#生成单位矩阵，用普通矩阵无法创建，所以用sparseMatrix来创建。
num_samples=180000
rownumber=c(1:num_samples)
colnumber=c(1:num_samples)
value=c(rep(1,each=num_samples))
M=sparseMatrix(rownumber,colnumber,x=value)
#给单位矩阵赋值，全部赋值完成后是一个4*50*180000个数值的对称稀疏矩阵。
for(i in 1:num_samples){
w=W #有50个浮点型数值。
n=neighbors #有50个整型数值。
M=M-t(w)
M=M-w
M=M+w%*%t(w)
}