How to convert dependent equation to independent for parallel computing (parfor)
1 view (last 30 days)
Show older comments
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
0 Comments
Answers (1)
Walter Roberson
on 11 Feb 2019
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 Comments
Walter Roberson
on 11 Feb 2019
why bother to take a cumulative sum of snapshot_i but then only extract the last layer? Why not just sum(snapshot_i,3) ?
Walter Roberson
on 11 Feb 2019
Also it looks like
snapshot_i = snapshot_c .* snapshot_d;
with no loop needed ?
See Also
Categories
Find more on Loops and Conditional Statements in Help Center and File Exchange
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!