Problem with parfor loop

Question

B.E. 2021-3-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/761721-problem-with-parfor-loop

编辑： B.E. 2021-3-5

Hello,

I need help, I have this problem when I use parfor loop for this code, where x a square matrix of order N = 75692 of size approximately equal to 13 Gb and t a row vector of

N=length(t);
alpha= [2.4e-7, 2.5e-7, 2.5e-7, 2.3e-6];
beta = [0.5,    0.8,    0.6,    0.5 ];
YT=zeros(N,4);
YD=zeros(N,4);
parfor ima1:4
    [YT(:,ima),YD(:,ima)]=f(t,x,alpha(ima),beta(ima));
end

Error using distcomp.remoteparfor/getCompleteIntervals (line 133)

2 个评论
显示无隐藏无

Jan 2021-3-3

This is a part of the error message. Please copy the complete message.

B.E. 2021-3-3

The complete message is

Error using distcomp.remoteparfor/getCompleteIntervals (line 133)

The parallel pool that parfor was using has shut down. To start a new parallel pool, run your parfor

code again or use parpool.

Error in Code_Method_of_Caracteristics (line 83)

parfor ima=2:4

A write error occurred while sending to worker 19.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Edric Ellis 2021-3-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/761721-problem-with-parfor-loop#answer_638781

在 MATLAB Online 中打开

That error basically means a worker crashed while trying to run the parfor loop. You mention that x is large. If you are using the 'local' cluster, then please be aware that x must be copied to each of the worker processes. This can cause a large amount of memory usage, and it's possible (probable?) that this is causing the workers to shut down.

If you're using a recent version of MATLAB, you might be able to use parpool('threads'), which uses multiple computational threads in a single process, and can avoid some memory duplication. (But not all MATLAB functions can operate in this environment).

Otherwise, you are going to be constrained by the memory on your system. Transferring x to the workers from the client incurs additional duplication while the messages are in transit, so if you're only just exceeding the memory, and you can build x directly on the workers by executing a function, then the following pattern might help:

xC = parallel.pool.Constant(@myFunctionThatBuildsX); % build 'x' directly on the worker
parfor ...
    [YT(:,ima),YD(:,ima)]=f(t,xC.Value,alpha(ima),beta(ima));   
end