嵌套的parfor循环

你不能使用parfor循环在另一个parfor循环。作为一个例子,下面的嵌套parfor循环是不允许的:

parfori = 1:10parforj = 1:5…结束结束

在MATLAB代码分析器^®编辑器标记的使用parfor在另一个parfor循环:

你不能嵌套parfor循环,因为并行可以执行在只有一个水平。因此,选择要并行运行的循环,把另一个循环为循环。

考虑以下性能问题在处理嵌套循环:

您还可以使用一个函数,使用parfor并将它嵌入parfor循环。并行化只发生在外层的水平。在下列的示例中,调用一个函数MyFun.m在外面parfor循环。内parfor循环嵌入MyFun.m按顺序运行,而不是平行的。

parfor我= 1:10 MyFun(我)结束函数MyFun(我)parforj = 1:5…结束结束

将嵌套为循环,parfor循环

典型的使用嵌套循环使用一个循环遍历一个数组变量指数一维,和一个内嵌循环变量指数另一个维度。基本形式是:

X = 0 (n, m);为一个= 1:n为b = 1: m X (a, b) =乐趣(a, b)结束结束

下面的代码显示了一个简单的例子。使用抽搐和toc测量所需的计算时间。

一个= 100;抽搐为i = 1:10 0为j = 1:10 0 (i, j) = max (abs (eig(兰德(a))));结束结束toc

运行时间是49.376732秒。

您可以嵌套循环的并行化,但你不能并行运行。原因是,工人们在一个平行池不能启动或访问进一步平行池。

如果循环数我转换为parfor池中循环,然后每个工人执行嵌套循环使用j循环计数器。的j循环本身不能运行parfor在每一个工人。

因为并行处理会增加开销,你必须仔细选择你是否想把内部或外部为循环,parfor循环。下面的例子展示了如何衡量并行开销。

首先将只外为循环,parfor循环。使用抽搐和toc测量所需的计算时间。使用ticBytes和tocBytes测量多少数据转移的工人在并行池。

运行新代码,并再次运行它。第一次运行比后续运行慢,因为平行池需要一些时间开始,使可用的代码工人。

一个= 100;抽搐ticBytes (gcp);parfori = 1:10 0为j = 1:10 0 (i, j) = max (abs (eig(兰德(a))));结束结束tocBytes toc (gcp)

BytesSentToWorkers BytesReceivedFromWorkers _____________ ________________________ 1 32984 24512 33784 25312 3 33784 25312 4 34584 26112 1.3514 1.0125 e + e + 05年05年运行时间是14.130674秒。

下一只转换内心的循环parfor循环。测量所需的时间和数据传输与前面的情况。

一个= 100;抽搐ticBytes (gcp);为i = 1:10 0parforj = 1:10 0 (i, j) = max (abs (eig(兰德(a))));结束结束tocBytes toc (gcp)

BytesSentToWorkers BytesReceivedFromWorkers _____________ ________________________ 1 1.3496 1.3496 5.487 e + e + 06 05年2 5.4858 e + e + 06 05年3 1.3476 1.3677 5.6034 e + e + 06 05年4 e + 06 5.4717总5.4144 e + e + 05 06 2.2048 e + 06年运行时间是48.631737秒。

如果你把内心的循环parfor循环,时间和传输的数据量远远大于在平行的外循环。在这种情况下,运行时间几乎是一样的嵌套为循环的例子。运行的速度小于外循环并行,因为你有更多的数据传输,因此并行开销。因此如果你执行内心的并行循环,得到比运行串行计算没有好处为循环。

如果你想减少并行开销和加速你的计算,并行运行外循环。

如果你把内心的外循环的每次迭代循环,然后启动一个单独的parfor循环。内循环转换创建100年parfor循环。每一个多个parfor死刑会增加开销。如果你想减少并行开销,你应该并行运行外循环,因为开销只发生一次。

嵌套的为循环:需求和局限性

如果你想转换一个嵌套为循环,parfor循环,您必须确保您的循环变量分类是适当的,看到的排除变量parfor-Loops。如果您的代码不遵守指导方针和限制贴上要求,你会得到一个错误。MATLAB捕获其中的一些错误时它读取代码。这些错误都贴上要求(静态)。

在接下来的例子中,左边的代码不工作因为你定义的上限为由一个函数调用循环。右边的代码提供了一个解决方法首先定义一个广播或常量变量之外的parfor循环:

一个= 0 (100、200);parfori = 1:尺寸(1)为j = 1:尺寸(2)(i, j) =我+ j;结束结束

一个= 0 (100、200);n =大小(2);parfori = 1:尺寸(1)为j = 1: n (i, j) =我+ j;结束结束

下面这个限制是必需的。如果嵌套为循环变量是改变了在一个地方parfor循环的为声明中,索引的为循环变量是不能保证可以在每个工人。

左边的代码是无效的,因为它试图修改嵌套的价值为循环变量j在体内的循环。右边的代码提供了一个通过指定嵌套解决方案为循环变量来一个临时变量t,然后更新t。

一个= 0 (10);parfori = 1:10为j = 1:10 (i, j) = 1;j = + 1;结束结束

一个= 0 (10);parfori = 1:10为j = 1:10 (i, j) = 1;t = j;t = t + 1;结束结束

下面这个限制是必需的。如果一个嵌套的为循环变量索引、迭代不能保证独立。

左边的例子是无效的,因为它试图指数嵌套为循环变量j。右边的例子会删除这个索引。

一个= 0 (10);parfori = 1:10为j = 1:10 j (1);结束结束

一个= 0 (10);parfori = 1:10为j = 1:10;结束结束

例如,下面的代码左边不工作,但是右边的代码:

一个= 0 (11);parfori = 1:4为j = 1:10 (i, j + 1) = i + j;结束结束

一个= 0 (11);parfori = 1:4为j = 2:11 (i, j) =我+ j - 1;结束结束

在接下来的例子中,因为左边的代码不工作一个嵌套的内切和索引吗为循环。右边的代码,因为工作v被分配给一个外的嵌套循环:

A = 0 (10);parfori = 1:4为j = 1:10 (i, j) = i + j;结束disp ((i, j))结束

A = 0 (10);parfor我= 1:4 v = 0 (10);为j = 1:10 v (j) = + j;结束disp (v (j))(我:)= v;结束

parfor循环的限制

函数一个= pfeg函数= nfcn(在)= 1 +;结束fcn = @nfcn;parforidx = 1:10 (idx) =函数宏指令(fcn idx);结束结束

> > pfeg开始平行池(parpool)使用“本地”概要文件…连接到4工人。ans = 2 3 4 5 6 7 8 9 10 11

x = 1:10;parfor我= 1:10 temp = x(我);anonymousFunction = @() 2 *温度;x (i) = anonymousFunction () + i;结束disp (x);

一个=“一个”;myFunction (a)函数name = inputname X = myFunction (a) (1);parfor我= 1:2 X (i)。(名)=我;结束结束

myFunction (“一个”,“b”)函数X = myFunction (a, b) nin =输入参数个数;parfor我= 1:2 X (i) = * nin;结束结束