大型稀疏方阵对角元素如何赋值

newdiver

A=sparse（10000，10000）vector=[4 6 10 12 13 15];
for i=1:length(vector)
A(vector(i),vector(i))=1;
end

这段代码速度很慢，怎样把这段代码简化一下，去掉循环呢

rocwoods

S = sparse(vector,vector,1,10000,10000);

newdiver

2# rocwoods
问题是A在以前被赋值过一次，里面有值，只是在（vector，vector）处改为1，
应该怎么弄呢

pasuka

2# rocwoods
问题是A在以前被赋值过一次，里面有值，只是在（vector，vector）处改为1，
应该怎么弄呢
newdiver 发表于 2011-1-12 01:44

1、用Fortran或者C/C++写好代码，再用mex；
2、可以尝试用parfor开启并行，不过是否支持稀疏矩阵没有验证过，且属于拿内存换效率的办法，内存小的话，请勿轻易尝试