How to correct the error “Valid indices for 'variable' are restricted in PARFOR loops”

Question

Nathaniel Werner 2018-6-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/404809-how-to-correct-the-error-valid-indices-for-variable-are-restricted-in-parfor-loops

评论： Walter Roberson 2018-6-13

I am trying to set up a parfor nested loop in MatLab R2016a as below.

N = size(A,1);
M = size(v,1);
in = zeros(N*M,1);
parfor i=1:N
   for j=1:M
      k = (i-1)*M-j;
      if sqrt(sum((A(i,:)-v(j,:)).^2))<=tol
        in(k) = i;          
      end
   end
end

However, I am getting the following error Valid indices for 'in' are restricted in PARFOR loops. Is there some way I can correct this because both arrays A and v are considerably large, over 40,000 rows for A. Thanks very much.

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

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2018-6-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/404809-how-to-correct-the-error-valid-indices-for-variable-are-restricted-in-parfor-loops#answer_323913

在 MATLAB Online 中打开

Step 1: go back and de-vectorize your in, making it a 2D array again.

Step 2: use temporary variables to process one row at a time, and afterwards

in(i, :) = temporary

It looks to me as if for each row, the columns that are within the tolerance distance should be replaced with the row number, and the ones not within tolerance should stay 0. If you had a logical vector for the row, that would be equivalent to multiplying the vector by the row number.

To get the logical vector you could use pdist2(), or you could perform the equivalent code using bsxfun or the new implicit expansion facilities of R2016b or later.

2 个评论
显示无隐藏无

Nathaniel Werner 2018-6-10

在 MATLAB Online 中打开

Hello, Thank you for responding and sorry for not replying sooner.

I have some questions for clarification.

When you say de-vectorize in and make it a 2D array again, the variable is already a 2D array. Can you explain what you mean?
The reason I am using in(k) not in(i,:) is because k will count from 1:1:N*M. I don't want to replace the row numbers, I just want to save them, but I'm not too picky as to how that is done. If the way you described will accomplish the job then I can do that. And I don't think I need logicals, simply doubles representing the integer row numbers should be fine. I plan on just removing all rows that are zero once the for loop is finished. Unfortunately, I am limited by the tools in R2016a.

in(in==0)=[]

Walter Roberson 2018-6-13

在 MATLAB Online 中打开

You have

in = zeros(N*M,1);

that is a vector. Instead use

in = zeros(M, N);

Normally I would have said zeros(N, M) given your N*M, but later in the code when you calculated k, you were doing so consistent with zeros(M,N)

"k will count from 1:1:N*M."

Don't do that for parallel work.

in = zeros(N, M);
parfor i=1:N
   Ai = A(i,:);
   ini = zeros(1,M);
   for j=1:M
      if sqrt(sum((Ai-v(j,:)).^2))<=tol
        ini(j) = i;          
      end
      in(i, ;) = ini;
   end
end