Блог

Как мне эффективно выполнять итерацию в ячейке?

У меня есть такая ячейка:

 H{i,j} = 

{1} {2} {3} {4} ...
{5} {6} {7} {8} ...
....
  

Фактический размер составляет ~ 300 * 300, всего ~ 1000 непустых элементов.

 for i = 1 : numel(H)
    if isempty(H{i}), continue, end
    for j = 1 : numel(H)
        if i==j || isempty(H{j}), continue, end
            COMPLEX_OPERATION(H{i}, H{J});
        end
    end
end
  

Это двухуровневая (простите за мой английский) итерация для одной ячейки.

Этот код оказался неэффективным.

Слишком много времени тратится впустую на итерации.

Я пытаюсь найти эффективный способ выполнить эту итерацию.

Одним из возможных решений является использование cellfun(@COMPLEX_OPERATION), но я не уверен, как выполнить этот двойной цикл с помощью cellfun.

Другое возможное решение — сначала удалить пустые ячейки. Однако я не уверен, решит ли это проблему в значительной степени.

Это также должно быть интересно. Используйте двойной вызов cellfun, чтобы получить все комбинации. Просто измените C1 C2 на свою complex_function

 %simulation
H=cell(300,300);
indx=randperm(numel(H));
[H{indx(1:1000)}]=deal(1);

%code you want
empty_H=cellfun(@isempty, H);
non_empty_H=H(~empty_H);
all_combinations=cellfun(@(C2) (cellfun(@(C1) C1 C2, non_empty_H, 'UniformOutput', false)), non_empty_H, 'UniformOutput', false);
  

Почему вы не можете просто сделать

 ind = ~cellfun(@isempty, H).* ~eye(size(H))

ind =

   0   1
   0   0

[r,c] = find(ind==1);
cellfun(@(H,r,c) (whatever_your_complex_function_does), H{[find(ind==1)].'}, cell2mat(r,[1 1]),cell2mat(c,[1 1]);
  

Более подробное объяснение…

 H ={[],[1 2]; [] [3 4]}



ind = cellfun(@isempty, H)

ind =

   1   0
   1   0
  

Это ваши пустые индексы. Мы инвертируем их в конце.

Затем найдите те, где i=j . Затем мы исключим их в конце.

Создайте идентификационную матрицу, представляющую i=j .

 I = eye(size(ind))

I =   

   1   0
   0   1
  

Затем умножьте их на ind. Мы используем ~I , чтобы получить те, где i~=j и (~ind) чтобы получить непустые элементы.

 (~ind).*~I

ans =

0   1
0   0
  

Теперь у вас есть индексы, где H не пусто и i~=j .

Получите индексы для них.

 [r,c] = find(ind==1);
  

Затем вы создаете анонимную функцию, которая принимает значения H , соответствующие этим индексам, и применяет ее по элементам к H . H , r и c будет использоваться в качестве входных данных для анонимной функции, где вы можете использовать ее как H{r} или H{c}

 cellfun(@(H,r,c) (whatever_your_complex_function_does), H{[find(ind==1)].'}, cell2mat(r,[1 1]),cell2mat(c,[1 1]);