Replacing a LOT of data by zero (in a matrix)

2011 8 31
2 个回答
回答已采纳
12 次查看（30 天）

0 个投票

I have a matrix M1, below you see the whos information for that matrix:
Name Size Bytes Class Attributes M1 29624x29624 108387640 double sparse
I want to replace quite a bit of data by zeros, using logical indexing: M1(abs(M1)<1e-1) = 0;
This needs more than 16GB memory, there must be a more efficient way. I also tried the other way round, that is M1 = M1(abs(M1)>1e-1);
Same problem. What's the trick?

请先登录，再回答此问题。

 采纳的回答

Herwig Peters
Herwig Peters 2011-8-31

1 个投票

4 way.
Find the indices of the desired matrix elements and the desired matrix elements themselves. With these 3 vectors, rebuild the sparse Matrix M1.
[i1,~,s1] = find(M1(abs(M1)>1e-1));
[i1,j1] = ind2sub(size(M1),i1);
M1 = sparse(i1,j1,s1,m,n);

2 个评论
显示 无 隐藏 无

Herwig Peters
Herwig Peters 2011-8-31
Actually, I just found it can be done even nicer:
M1 = (abs(M1)>1e-1).*M1;

请先登录，再进行评论。

更多回答（1 个）

Dmitry Borovoy
Dmitry Borovoy 2011-8-31

0 个投票

1 way. Use vectorization (your way). fast but requires a lot of RAM 2 way. Use loops. Very slow but no need a lot of RAM 3 way. For perverts. Write your code on C/C++/Fortran with loop and call it from Matlab. If you have restriction on RAM or on calculation speed, I suppose, it's your way.

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Herwig Peters
Herwig Peters 2011-8-31
Yes, that sums it up pretty well. 1 needs too much RAM, 2 needs too much time (CPU) and 3 needs too much time as well, my time.
See my own answer for another possibility.

请先登录，再进行评论。

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by