"valid indices for 'a' are restricted in PARFOR loops" for unindexed struct?

The reason a.b = k does NOT work is because when the parfor finishes, and there are, say N workers creating 1000 parallel versions of a.b = k, which one is going to be the final a.b? If 1000's of parallel universes converge into 1 universe, which universe will we be the "real" universe? Error error.

But, a(k).b = k works because each worker creates its own separate variable (or isolated universe), a(1).b, a(2).b, ..., a(N).b. No conflict issues.

In the example, this works too because temp is defined inside the parfor, which tell matlab "this is a temporary variable that's created in the parallel world - do not bring it to the real world"

parfor i = 1:4
  temp = struct();  %What's created in the parallel world stays in the parallel world
  temp.myfield1 = rand();
  temp.myfield2 = i;
end

7 个评论
显示 5更早的评论隐藏 5更早的评论

Steven Lord 2018-7-16

a.b=k is a sliced variable

No, it's not. The characteristics a variable in a parfor loop must have in order to be a sliced variable, according to the documentation, are:

Type of First-Level Indexing — The first level of indexing is either parentheses, (), or braces, {}.
Fixed Index Listing — Within the first-level parentheses or braces, the list of indices is the same for all occurrences of a given variable.
Form of Indexing — Within the list of indices for the variable, exactly one index involves the loop variable.
Shape of Array — The array maintains a constant shape. In assigning to a sliced variable, the right side of the assignment cannot be [] or '', because these operators attempt to delete elements.

Your code, a.b = k, fails the first characteristic. In fact, the first "Reference for Variable Not Sliced" example in that section of the documentation page is your example, except that the documentation example uses the names A and x instead of a and b respectively.

Note that in the "Structure Arrays in parfor-Loops" section on this documentation page, the "valid" code makes a temporary variable, not a sliced variable. To slice, you need to do what's in the "Slicing Structure Fields" sub-section of that section of the page.

OCDER 2018-7-16

在 MATLAB Online 中打开

More examples for clarity (or confusion).

parfor k = 1:10 %<=== look for the iterator, 'k'
   a(k)      % sliced, I see a direct use of 'k'
   a(k).b    % sliced, I see a direct use of 'k' in the 1st level of a struct
   a.b       % NOT sliced. No use of 'k' anywhere
   a.b(k)    % Tricky. Not sliced at the 1st lv of struct 'a', but is sliced at 2nd lv 'a.b'. 
             % Matlab doesn't know if you want to save 'a' or treat as temporary var
   a(k).b(k) % Tricky. Sliced variable! Sliced at 1st lv of struct a.
end

Hmm, yeah, the explanation isn't obvious as to why a.b = k is not just assumed to be a temporary variable. My guess is that structures are normally assumed to behave like objects that persistent across iteration, but parfor makes this assumption invalid. To prevent confusion, users have to be VERY CLEAR they know this structure will delete after each iteration unless it's sliced.

for j = 1:10
 a.b = j;  %structure "a" is created when j == 1, and it persists from here on out. 
           %a.b = 10 when loop finishes. No ambiguity. 
end
parfor j = 1:10
 a.b = j;  %structure a is created 10 times at random order! a.b = 10, a.b = 2, ... a.b = 5. 
           %How do we assemble this back to a single "a" structure? 
           %Should this even be a single structure as an output?
           %Was this a temporary variable?
           %Let's just make this throw an error. Unclear how to return "a".
end

Evan 2018-7-17

编辑：Evan 2018-7-17

@OCDER I appreciate your perseverance. I stepped away from this for a day to give it a fresh look. I think I have thought of a logical reason why a.b=k doesn't work. If the struct a is defined before the parfor loop, it would seem strange for a.b=k to completely wipe the previous a object, whereas a=struct() does this explicitly. And it is not possible to know by looking at the code whether a will be in the workspace before the code is run (e.g. a load operation).

Or maybe someone from Mathworks will jump in and give the true reasoning. But in any case, I think I would have figured this out more quickly if the error message was more precise about what I was doing wrong. Hopefully the enhancement request goes through.

请先登录，再进行评论。

"valid indices for 'a' are restricted in PARFOR loops" for unindexed struct?

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

回答（1 个）

7 个评论
显示 5更早的评论隐藏 5更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

"valid indices for 'a' are restricted in PARFOR loops" for unindexed struct?

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

回答（1 个）

7 个评论 显示 5更早的评论隐藏 5更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

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

7 个评论
显示 5更早的评论隐藏 5更早的评论