Question

为什么我不能在这段代码中使用parfor？

parfor i=1:r

    for j=1:N/r

        xr(j + (N/r) * (i-1)) = x(i + r * (j-1));

    end

end

这是错误：

错误：parfor中的变量xr无法分类。请参见MATLAB中的并行循环，“概述”。

Answer 1

这里的问题是对切片阵列进行不正确的索引。 parfor循环以异步运行，这意味着每次迭代的执行顺序是随机的。来自documentation：

MATLAB工作人员不按特定顺序评估迭代，并且彼此独立。因为每次迭代都是独立的，所以不能保证迭代以任何方式同步，也不需要这样做。

您可以通过在命令行中键入以下内容来轻松验证上述语句：

parfor i=1:100
    i
end

你会看到排序是任意的。因此，如果您在不同的工作者之间拆分并行作业，则一个工作人员无法确定是否已完成不同的迭代。因此，您的变量索引不能依赖于迭代器的过去/未来值。

让我用一个简单的例子来证明这一点。考虑Fibonacci系列1,1,2,3,5,8,...。您可以轻松地（在一个天真的for循环中）生成该系列的前10个术语：

f=zeros(1,10);
f(1:2)=1;
for i=3:10
    f(i)=f(i-1)+f(i-2);
end

现在让我们对parfor循环做同样的事情。

f=zeros(1,10);
f(1:2)=1;
parfor i=3:10
    f(i)=f(i-1)+f(i-2);
end

???错误：parfor中的变量f无法分类。请参见MATLAB中的并行循环，“概述”

但为什么会出错呢？

我已经证明迭代是以任意顺序执行的。因此，假设一个worker得到循环索引i=7和表达式f(i)=f(i-1)+f(i-2);。现在应该执行表达式并将结果返回给主节点。现在迭代i=6完成了吗？ f(6)中存储的值是否可靠？那么f(5)呢？你看到我得到了什么吗？假设f(5)和f(6)未完成，那么您将错误地计算出Fibonnaci系列中的第7项为0！

由于MATLAB无法确定您的计算是否可以保证每次都正确运行并重现相同的结果，因此明确禁止此类模糊分配。

Parfor在MATLAB中的问题

1 个答案: