Как исправить Допустимые индексы для «переменной» ограничены ошибкой циклов PARFOR в Matlab

Проблема в том, что MATLAB не распознает, что переменная k правильно разрезает матрицу in. Решение должно состоять в том, чтобы индексировать in, используя i и j отдельно:

N = size(A,1);
M = size(v,1);
in = zeros(M,N);
parfor i=1:N
   for j=1:M
      if sqrt(sum((A(i,:)-v(j,:)).^2))<=tol
        in(j,i) = i;          
      end
   end
end
in = in(:); % reshape to a column vector, as the output in the question's code

Другая альтернатива, но она требует больше промежуточной памяти, состоит в том, чтобы вычислить это вообще без цикла:

A = reshape(A,1,N,[]);
v = reshape(v,M,1,[]);
in = sum(bsxfun(@minus,A,v).^2,3) < tol*tol;
in = in(:);

(Или что-то подобное, я не запускал этот код... Пожалуйста, дайте мне знать или исправьте пост, если есть опечатка или другая ошибка.)

N = size(A,1);
M = size(v,1);
in = cell(N,1);
parfor i=1:N
    s=v;
    p=zeros(1:M,1);
   for j=1:M
      k = (i-1)*M+j;
      if sqrt(sum((A(i,:)-s(j,:)).^2))<=tol
        p(k) = i;          
      end
   end
   in{i}=single(p);
end
in=cell2mat(in);
in=reshape(in,[N*M,1]);

иногда Matlab не распознает переменную в цикле parfor как нарезанную переменную, нарезанная переменная — это переменная, которая имеет ссылку вне цикла parfor, и каждый из ее элементов доступен только одному рабочему (в parfor параллельных рабочих)
поэтому вы можно использовать временную переменную и собирать результаты после цикла parfor,

ПРИМЕЧАНИЕ 1. В более старых версиях код лучше векторизовать, потому что циклы не были такими хорошими, как сейчас в R2017, ссылаясь на (это).

ПРИМЕЧАНИЕ 2: если большинство элементов in равно нулю, попробуйте использовать разреженную матрицу, которая может сэкономить много памяти;