MATLAB: Parallel for loop: Using broadcasting array variable values as indexes

parallel computingParallel Computing Toolbox

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.

Best Answer

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
Related Question