Parallel for loop: Using broadcasting array variable values as indexes

Question

Barend 2011-11-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/21148-parallel-for-loop-using-broadcasting-array-variable-values-as-indexes

Good day.

I have recently aquired the parallel computing toolbox to speed up some of our applications. I have run into the following challenge and hope that someone can give me advice on this.

topNodes    = [2 3 4 5 6 7];
bottomNodes = [1 2 3 4 5 6];
polarity    = [1 1 -1 1 1 1];
nodes       = max([topNodes bottomNodes]);
T           = zeros(numel(polarity),nodes);
for i = 1:numel(polarity);
  if polarity(i) < 0
      T(i,bottomNodes(i)) = -1;
      T(i,topNodes(i)) = 1;
  else
      T(i,bottomNodes(i)) = 1;
      T(i,topNodes(i)) = -1;
  end
end

Now, if I use a parfor loop, it complains. Some assistance will be much appreciated.

Kind regards,

Barend.

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

请先登录，再进行评论。

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

Answer 1

Sarah Wait Zaranek 2011-11-14

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/21148-parallel-for-loop-using-broadcasting-array-variable-values-as-indexes#answer_27927

在 MATLAB Online 中打开

For parfor loops, when you index into a sliced variables, restrictions are placed on the first-level variable indices. This allows parfor to easily distribute the right part of the variable to the right workers. One of these first-level indices must be the loop counter variable or the counter variable plus or minus a constant. Every other first-level index must be a constant, a non-loop counter variable, a colon, or an end.

I explain this more in the blog post where you also left this question.

http://blogs.mathworks.com/loren/2009/10/02/using-parfor-loops-getting-up-and-running/

See a possible fix below.

topNodes    = [2 3 4 5 6 7];
bottomNodes = [1 2 3 4 5 6];
polarity    = [1 1 -1 1 1 1];
nodes       = max([topNodes bottomNodes]);
T           = zeros(numel(polarity),nodes);
parfor i = 1:numel(polarity);
  myData = T(i,:);  
  if polarity(i) < 0
      myData(i,bottomNodes(i)) = -1;
      myData(i,topNodes(i)) = 1;
  else
      myData(i,bottomNodes(i)) = 1;
      myData(i,topNodes(i)) = -1;
  end
  T(i,:) = myData;
end

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

Sarah Wait Zaranek 2011-11-15

Yes, there is an overhead to pass the data to and from the workers - so for code that takes seconds to run - there isn't often a speedup. Also, you are correct that vectorization often is a better choice then de-vectorizing and using parfor.

Konrad Malkowski 2011-11-17

You should also try running the script as a function. In general MATLAB performance is better with functions that with scripts, as scripts are interpreted one line at a time.

请先登录，再进行评论。