How to convert dependent equation to independent for parallel computing (parfor)

1 次查看（过去 30 天）
Abhishek Saini
Abhishek Saini 2019-2-11
评论： Walter Roberson 2019-2-11
I am trying to run my code using 'parfor' but getting error due to dependent term. How to resolve that.
The code similar to below. I am getting error for third loop expression
image(: , : , i) = image(: , : , i) + image(: , : , i -1);
parfor ix = 1:nx
% expression to calculate image_c and image_d
for i = 1:nt
image(:,:,i) = image_c(:,:,i) +image_d(:,:,i);
end
for i=2:nt
image(:,:,i) = image(:,:,i) +image(:,:,i-1);
end
end

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

回答（1 个）

Walter Roberson
Walter Roberson 2019-2-11
You cannot do that in parallel.
What you can do is roughly
total = zeros(size(image));
parfor i = 1 : nt
total = total + image(:,:,i);
end
MATLAB recognizes total as being a "reduction variable" and will take care of doing the relevant additions and locking.
  5 个评论
显示 3更早的评论
Walter Roberson
Walter Roberson 2019-2-11
why bother to take a cumulative sum of snapshot_i but then only extract the last layer? Why not just sum(snapshot_i,3) ?

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Images 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by