parfor variable classification issue

1 次查看（过去 30 天）
Leos Pohl
Leos Pohl 2021-5-24
评论： Raymond Norris 2021-5-25
I had the code as such (note the size of the map array is (1001,1001,3):
map = zeros(i1max - i1min + 1,j1max - j1min + 1,3);
parfor (i1 = i1min:i1max)
for j1 = j1min:j1max
[a,b] = some_function(1,i1,j1);
%some work
map(i1,j1,1)=r1;
map(i1,j1,2)=phi1;
map(i1,j1,3)=alpha;
end
end
But matlab complained that Parfor cannot run due to the way the variable map is being used. So after reading some online documentation, i changed it to :
map = zeros(i1max - i1min + 1,j1max - j1min + 1,3);
parfor (i1 = i1min:i1max)
map_temp = zeros(j1max - j1min + 1,3);
for j1 = j1min:j1max
[a,b] = some_function(1,i1,j1);
%some work not using i1 or j1 variables but a,b
jdx = j1 - j1min + 1;
map_temp(jdx,1)=r1;
map_temp(jdx,2)=phi1;
map_temp(jdx,3)=alpha;
end
idx = i1 - i1min + 1;
map(idx,:,:) = map_temp;
end
But that does not work either because apparently as matlab complains "Valid indices are restricted", it has something to do with the idx variable. Why? And how to make this work?

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

采纳的回答

Matt J
Matt J 2021-5-25
编辑：Matt J 2021-5-25
map is a sliced variable, so it can only be indexed using the loop variable.
map = zeros(i1max - i1min + 1,j1max - j1min + 1,3);
parfor (idx = 1:i1max-i1min+1)
i1=idx+i1min-1;
map_temp = zeros(j1max - j1min + 1,3);
for j1 = j1min:j1max
[a,b] = some_function(1,i1,j1);
%some work not using i1 or j1 variables but a,b
jdx = j1 - j1min + 1;
map_temp(jdx,1)=r1;
map_temp(jdx,2)=phi1;
map_temp(jdx,3)=alpha;
end
map(idx,:,:) = map_temp;
end
  1 个评论
Raymond Norris
Raymond Norris 2021-5-25
I read the posting, but went back to it later without refreshing the page first. Matt's got the idea.

请先登录，再进行评论。

更多回答（1 个）

Raymond Norris
Raymond Norris 2021-5-25
As a side note, the difference of min and max could be 1000, but there are infinity cases where j1 will start outside the boundaries of map. For instance,
j1min = 200000;
j1max = 201000;
In this case, j1 = 200000, which would far exceed the allocation of map (1001x1001x3). This would still work, but it's probably not what you want.
With that said, try the following. You're using i1 as a parameter for some_function but you really need i1 to slice map (which is what you're defining idx for), so let's swap it. Assign idx before parfor as such
idx = i1min:i1max;
and then use the loop variable, i1, to pass it to some_function. Now use i1 to index into map at the bottom of the loop.
map = zeros(i1max - i1min + 1,j1max - j1min + 1,3);
idx = i1min:i1max; % used for some_function
parfor (i1 = i1min:i1max - i1min + 1)
map_temp = zeros(j1max - j1min + 1,3);
for j1 = j1min:j1max
[a,b] = some_function(1,idx(i1),j1);
%some work not using i1 or j1 variables but a,b
jdx = j1 - j1min + 1;
map_temp(jdx,1)=r1;
map_temp(jdx,2)=phi1;
map_temp(jdx,3)=alpha;
end
% idx = i1 - i1min + 1; % not needed
map(i1,:,:) = map_temp; % use i1, not idx
end

类别

在 Help Center 和 File Exchange 中查找有关 2-D and 3-D Plots 的更多信息

标签

产品


版本

R2018b

Community Treasure Hunt

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

Start Hunting!

Translated by