Parallel nested for-Loop with a compounding variable

Question

Jeongseop Lee 2015-10-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/246902-parallel-nested-for-loop-with-a-compounding-variable

回答： Jeongseop Lee 2015-10-6

Hi I'm trying to execute a code that contains three for-loops with the innermost loop containing a variable whose value gets compounded at every iteration.

Here is an example:

A = [];
for i = 1:10
  for j = 1:10
    B = f(i,j);
    for k = 1:10
      A = A + g(B);
    end
  end
end

The nested for-loop itself cannot be worked around unfortunately but this is taking too long computation-wise I am trying to incorporate "parfor", the matlab does not allow variable "A". Is there a workaround to this?

Thanks!

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

请先登录，再进行评论。

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

Answer 1

Kirby Fears 2015-10-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/246902-parallel-nested-for-loop-with-a-compounding-variable#answer_194790

在 MATLAB Online 中打开

Each iteration of a parfor loop will be performed in an unpredictable order. Therefore Matlab prevents you from changing the value of a variable in one iteration that will be referenced in another iteration (such as A in your code). This general rule is enforced to ensure coherent calculations are performed.

You should store the values in an array or a cell, then sum the terms after the parfor loop has ended. Here's a simple example that should be adaptable to your problem.

% Initialize an array to store all terms
A=NaN(1,10);
parfor j=1:10,
    % Fill each of the j positions independently
    A(j)=j;
end,
% Sum A now that parfor is over
A=sum(A);

Hope this helps.

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

请先登录，再进行评论。

Answer 2

Walter Roberson 2015-10-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/246902-parallel-nested-for-loop-with-a-compounding-variable#answer_194825

If you initialize A as 0 then you can do it with parfor, as A would then be recognized as a "reduction variable". The calculation might internally be done with subtotals that are then totaled, so if your calculation depends critically on the order in which the totals are done then parfor is not appropriate. For example, ((2 + (-2)) + eps(1)) would have a different result than (2 + ((-2) + eps(1)) so if you have carefully programmed your loop to avoid those kinds of round-off problems then you would want to be very careful how you programmed your parfor.

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

请先登录，再进行评论。

Answer 3

Jeongseop Lee 2015-10-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/246902-parallel-nested-for-loop-with-a-compounding-variable#answer_194985

My initial take on this problem based on your answer was to pre-allocate an array, i.e., "Array", to store all the intermediate calculation results under different indices "Array{i,j,k}". The parfor command this did not allow the use of this type of array (it did not recognize the variable type according to the error message). So I had to reduce it to a two-index array instead of the three-index and had a summation along one index take place outside of the nested loop and it works. Thanks.