ODE Event Detection in Parfor loop

2019 8 28
1 个回答
更新时间：2019 8 30
5 次查看（30 天）

0 个投票

I have an ODE45 function that uses an event detection function, and normally it works great! But, I'm trying to parallelize the process and putting the ode45 call inside the parfor loop fails to record the events. My basic code is below. When I use parfor rather than a simple for loop, the .te, .ye, and .ie fields of my struct are all blank (but are correctly populated with "for"), while the .y, and .t fields are fine.
State(1:100) = struct;
parfor i = 1:length(State)
tspan = 1:1000;
Opt = odeset('RelTol',1e-10,'AbsTol',1e-10,'Events',@myEvent);
[t,y,te,ye,ie] = ode45(@myOde, tspan, Y0, Opt);
State(i).t = t;
State(i).y = y;
State(i).te = te;
State(i).ye = ye;
State(i).ie = ie;
end
Is there a way to use Event Detection in a parfor loop? Am I just initializing something wrong, or do I need to try something other than parfor?

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Edric Ellis
Edric Ellis 2019-8-29
Your example is not executable, so I can't work out why this isn't working - please could you set this up as a "Minimal, Reproducible Example".

请先登录，再进行评论。

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

回答（1 个）

James Penrod
James Penrod 2019-8-29

0 个投票

I've managed to resolve this myself, while trying to set up the Mininal, Reproducible Example (Thank you Edric for your time in pointing that out). The example I had come up with is below, and it appears that the issue was trying to use a global variable in the Event Function in the parfor loop. When I remove the "global R" from the event function and change "value = dist-R" to "value = dist-1" my final State struct is properly populated
State(1:10) = struct;
global R
R = 1;
parfor i = 1:length(State)
tspan = 1:10;
Y0 = [0;0;0;10;10;10;2;2;2;0;0;0];
Opt = odeset('RelTol',1e-10,'AbsTol',1e-10,'Events',@myEvent);
[t,y,te,ye,ie] = ode45(@myOde, tspan, Y0, Opt);
State(i).t = t;
State(i).y = y;
State(i).te = te;
State(i).ye = ye;
State(i).ie = ie;
end
function[xdot] = myOde(tsim,x0)
xdot = [x0(7:9);x0(10:12);zeros(6,1)];
end
function[value, isterminal, direction] = myEvent(t,y)
global R
X1 = y(1:3);
X2 = y(4:6);
dist = norm(X1-X2);
value = dist-R;
isterminal = 1;
direction = 0;
end

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Edric Ellis
Edric Ellis 2019-8-30
Glad you got it sorted (ah, the value of the "Minimal, Reproducible Example"!). As you have found out, global variables are not synchronized from client to workers, or across workers.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Parallel and Cloud 的更多信息

产品

版本

R2017b

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by