Sum over cell array of sparse matrices

2017 4 2

1 个回答

9 次查看（30 天）

0 个投票

I have 100 sparse matrices of size N by N, where N=65536.

They have about 5500000 nonzero elements(0.13% of N^2) close to the diagonal.

They are stored in a cell array, S{1}, ... , S{100}, and I want to compute the sum S{1}+...+S{100}.

Sum=sparse(N,N);
for i=1:100
    Sum=Sum+S{i};
end

The above for loop code took about 25 seconds. Is there any way to optimize this code?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

Stalin Samuel 2017-4-2

0 个投票

v1=[S{1:100}];
Sum=sum(v1);

Gobi 2017-4-2

This code gives a 1 by N*100 vector. v1 is a matrix stacking S{:}, and just summing up is not what I want.

在帮助中心和 File Exchange 中查找有关 Sparse Matrices 的更多信息

Find the treasures in MATLAB Central and discover how the community can help you!

Translated by