Matlab工人怎么工作?

3视图(30天)

显示旧的评论

托拜厄斯Brambier 2022年5月10日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/1715865-how-do-matlab-workers-work

评论道: Edric埃利斯 2022年5月23日

所以我一直在努力优化/ parallelising现有的代码作为项目的一部分,在我的研究。这样做的时候,我遇到了一个非常奇怪的问题,至少给我。

“问题”时我让代码运行在一个正常的循环,而parfor-loop。根据parfor-loop tic-toc命令的运行时使用一个职工大约一半的运行时,它将使用一个标准的for循环。

我的问题是,根据我的理解独立的Matlab是单线程的。但也使用一个线程使用一个职工仍然是更快,甚至更快。这正是我对形势的理解让我困惑。

是的,我有检查,实际上只有一个工人使用单个线程。

最重要的是我需要一个解释,我为这个项目最终的论文。

我真的希望找到任何解释或甚至一个提示,对我来说,奇怪的行为。

4评论
显示3年长的评论隐藏3年长的评论

Edric埃利斯 2022年5月11日

你能缩小这一个简单的复制,你可以分享吗?如果没有,我们只能推测原因可能是什么。提供你使用 parpool(“本地”) (而不是 “线程” ),那么是的,MATLAB工人线程(默认情况下)。桌面MATLAB可以使用多个线程进行一些操作(比如大型矩阵计算),这通常会导致计算速度比在一个工人在桌面上。

一个可能性这是使用的代码转换吗 parfor 可能为职工得到更高效的代码执行。因为 parfor 循环体执行某些限制,转换后的代码可以被更有效率。这里有几件事情我试试(如果你不能轻易地提出一个简洁的复制):

如果你跑parfor没有平行池?(被相同的代码转换,要做到这一点,要么关掉并行优先自动创建一个池塘,否则使用特殊的parfor (i = 1: N, 0)语法)
是你的为循环在脚本或一个函数?如果它是一个脚本,使它成为一个功能——改变了MATLAB语言的优化可以使用

托拜厄斯Brambier 2022年5月11日

因为我不知道如何报价,或者甚至可能:

Edric埃利斯的问题:“如果你运行 parfor 没有平行池?”

答:没有并行运行parfor循环池的运行时接近我得到当运行一个标准的for循环。虽然快一点,仍然离*我使用parpool尺寸1。

从沃尔特·罗伯森推荐:“也试着尝试设置numcompthread 1。”

答:我认为你的意思 NumThreads 在选择菜单选项,然后是的我也试着这样做。实际上这是我的意思,当我写道:“是的,我已经检查,实际上只有一个工人使用单个线程。”I now realise, I could have clarified that with a bit more detail.

从Edric埃利斯的假设:“提供你使用 parpool(“本地”) (而不是 “线程” )……”

答:是的,因为当地集群是默认我的测试。使用线程而导致运行时非常相似,但仍然较慢,使用当地的虽然。

Edric埃利斯的问题:“是你为循环在脚本或一个函数?”

答:目前在脚本中使用。因为我不知道你什么意思,我只是试着使用函数关键词,导致运行时更慢就在脚本中使用它。

Edric埃利斯的问题:“你能够缩小这一个简单的复制,你可以分享吗?”

答:我将问我的上司如果我允许仅仅分享这篇文章的示例代码造成的混乱。虽然我不知道是有帮助的。我也可以试着降低到一个更简单的版本的代码共享,看看情况保持不变。

到目前为止,感谢你的建议!我也忘了说的是用于组装的代码仅仅是stiffnes矩阵中使用“有限元”计划称为“DAEdalon”由指导教授。

沃尔特·罗伯森 2022年5月11日

//www.tatmou.com/help/matlab/ref/maxnumcompthreads.html

登录置评。

在回答这个问题。

答案(1)

托拜厄斯Brambier 2022年5月19日

0
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/1715865-how-do-matlab-workers-work answer_967055

编辑:托拜厄斯Brambier 2022年5月19日

好了,之后我发现以前实现waitbar创建了一些巨大的开销是一个简化版本的代码仍然引起同样的困惑:

                              %这通常只是用于生成外部提供输入数据
                             
                              %这个数字(深圳)是一个增加改变的工作量!
                             
                              深圳= 120;
                             
                              为我= 1:深圳
                             
                              为j = 1:深圳
                             
                              x ((i - 1) *深圳+ j:) = (j - 1, - 1);
                             
                              结束
                             
                              结束
                             
                              为j = 1: sz - 1
                             
                              为我= 1:深圳
                             
                              公司= (j - 1) *深圳+ i;
                             
                              :康涅狄格州(inc) =(公司,公司+深圳);
                             
                              结束
                             
                              结束
                             
                              为j = 1:深圳
                             
                              为i = 1: sz - 1
                             
                              公司= (j - 1) *深圳+ i;
                             
                              康涅狄格州(+ 1,:)=(公司,公司+ 1);
                             
                              结束
                             
                              结束
                             
                              numnp =大小(x, 1);
                             
                              的ndm =大小(x, 2);
                             
                              ndf = ndm;
                             
                              numelem =大小(康涅狄格州,1);
                             
                              k_parallel =稀疏(numnp * ndf, numnp * ndf);
                             
                              %这是重要组成部分
                             
                              parpool (1)
                             
                              抽搐
                             
                              parfore = 1: numelem
                             
                              k_elem_gs =稀疏(numnp * ndf, numnp * ndf);
                             
                              I =康涅狄格州(e, 1);
                             
                              J =康涅狄格州(e, 2);
                             
                              t = x (J:)——(我:);
                             
                              L =规范(t);
                             
                              t = t / L;
                             
                              B = (- t, t) / L;
                             
                              k_elem = B * B * L;
                             
                              为i = 1:2
                             
                              我=康涅狄格州(e,我);
                             
                              为j = 1:2
                             
                              J =康涅狄格州(e, J);
                             
                              k_elem_gs (ndf *(张)+ 1:我* ndf, ndf * (J - 1) + 1: J * ndf) =…
                             
                              k_elem_gs (ndf *(张)+ 1:我* ndf, ndf * (J - 1) + 1: J * ndf)…
                             
                              + k_elem (ndf *(张)+ 1:我* ndf, ndf * (j - 1) + 1: j * ndf);
                             
                              结束
                             
                              结束
                             
                              k_parallel = k_parallel + k_elem_gs;
                             
                              结束
                             
                              toc

我没有一个明确的答案,但这就是我认为是怎么回事。串行MATLAB分析器显示,昂贵的计算的增量 k_elem_gs 的更新 k_parallel 。特别是,我认为的更新 k_parallel 得到更多的非零元素数量的增加。这是重要的,因为一个 parfor 这些添加循环重新排列。即使当使用一个工人, parfor 的工作循环分为“子范围”,分别执行。这些附属的区域将开始增加 k_parallel 从头开始,即从“廉价”添加。因此而客户端为这是循环:

k = k + k + k + k + k + k + k;

的一名员工 parfor 更喜欢这个:

k = {k + k + k} + {k + k} + {k + k};

(每 k 右边当然是不同的)。仍然有相同数量的增加,但并不是所有的昂贵的增加。

登录置评。

在回答这个问题。

类别

并行计算并行计算工具箱并行for循环(parfor)

找到更多的在并行for循环(parfor)在帮助中心和文件交换

下载188bet金宝搏

并行计算工具箱

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!