PARFOR 循环的 MATLAB 切片变量

Question

我正在尝试在MATLAB 中使以下循环并行友好，以便我可以使用parfor：

for ivert = 1 : nVerts
    b = obj.f( obj.neighIDs{ ivert } ); 
    x = obj.coeffMatrix{ ivert } \ b;
    obj.solution( ivert, : ) = x( 1 : 3 );
end

我尝试根据here 发布的 MATLAB 文档对变量进行切片：

parfor ivert = 1 : nVerts
    i = obj.neighIDs{ ivert };
    b = obj.f( i ); 
    A = obj.coeffMatrix{ ivert }
    x =  A \ b;
    obj.solution( ivert, : ) = x( 1 : 3 );
end

但 MATLAB 抱怨说：

Valid indices for `obj` are restricted in PARFOR loops.

有人可以给我一些提示如何在上述循环中分割变量吗？

Answer 1

这里的问题是 MATLAB 看到您的 parfor 循环的前三行，并将它们视为 obj 上的索引表达式 - 它得出结论 obj 必须是 parfor“广播”变量。 parfor 循环的最后一行被视为obj 的索引赋值（即使它看起来像是obj 字段的索引赋值）。因为obj 已被归类为“广播”，所以您不能分配给它。要解决此问题，我建议您执行以下操作：

tmpSolution = zeros(nVerts, 3);
parfor ivert = 1:nVerts
    ... %# calculate 'x'
    tmpSolution(ivert, :) = x(1:3);
end
obj.solution = tmpSolution;