How can i find out all indices of nonzero elements in a sets of matrices?

Question

Jiahong Zou 2023-6-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1977229-how-can-i-find-out-all-indices-of-nonzero-elements-in-a-sets-of-matrices

评论： Jiahong Zou 2023-6-2

I have 19 sparse matrices which are 30000 * 30000 big. Each matrix is slightly different from another. How can i use [row, col] = find() function to get row and col vectors, which contain all position of nonzero element?

For example, A = [1 0; 0 0], B = [ 0 1; 0 0], C = [ 0 0 ; 1 0], D = [ 0 1; 1 0]

row = [ 1 1 2], col = [1 2 1]

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

请先登录，再进行评论。

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

Answer 1

Aman 2023-6-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1977229-how-can-i-find-out-all-indices-of-nonzero-elements-in-a-sets-of-matrices#answer_1248969

在 MATLAB Online 中打开

Hi Jiahong,

Assuming that you want to have unique pairs of row and column, the following code should work:

% Create a cell array of example sparse matrices
A = sparse([1 0; 0 0]);
B = sparse([0 1; 0 0]);
C = sparse([0 0; 1 0]);
D = sparse([0 1; 1 0]);
matrices = {A, B, C, D};
% Initialize empty row and column vectors
row = [];
col = [];
% Loop over each matrix in the cell array
for i = 1:numel(matrices)
    % Use find to get the row and column indices of all nonzero elements
    [r, c, ~] = find(matrices{i});
    % Append the row and column indices to the overall row and column vectors
    row = [row; r];
    col = [col; c];
end
% Remove duplicate pairs from row and col
[row_col, idx] = unique([row, col], 'rows');
row = row_col(:,1);
col = row_col(:,2);
% Display the row and column vectors
row
col

If you want all pairs of row and column, use the following code:

% Create a cell array of example sparse matrices
A = sparse([1 0; 0 0]);
B = sparse([0 1; 0 0]);
C = sparse([0 0; 1 0]);
D = sparse([0 1; 1 0]);
matrices = {A, B, C, D};
% Initialize empty row and column vectors
row = [];
col = [];
% Loop over each matrix in the cell array
for i = 1:numel(matrices)
    % Use find to get the row and column indices of all nonzero elements
    [r, c, ~] = find(matrices{i});
    % Append the row and column indices to the overall row and column vectors
    row = [row; r];
    col = [col; c];
end
% Display the row and column vectors
row
col

Hope this helps!

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

Jiahong Zou 2023-6-2

Hi Aman, thx for the answer. There are actually about 180000 nonzero elements in each matrix. I'm afraid it will run out of memory before removing duplicate pairs. Is there a more memory-saving way?

请先登录，再进行评论。