Is there a way to avoid repeating time-consuming calculations from odefun in the events function?

Question

Are Mjaavatten 2019-10-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/483972-is-there-a-way-to-avoid-repeating-time-consuming-calculations-from-odefun-in-the-events-function

编辑： Are Mjaavatten 2019-10-8

The right-hand side in my ODE is quite time-consuming. In the course of calculating dy/dt I calculate a variable u that I also use in my events function. Is there a way to do the time-consuming calculations only once at each time step?

Example, where I have replaced the time-consuming calculation with a simple one:

function event_example
   opt = odeset('events',@stopfun);
   y0 = 1;
   [t,y] = ode45(@rhs,[0,3],y0,opt);
   plot(t,y)
end
function [dydt,eventvalue] = time_consuming_fun(y)
% a lot of calculations to find variable u, here simplified to:
 u = y;
 eventvalue = u;
 dydt = -u*y;
end
function dydt = rhs(~,y)
  dydt = time_consuming_fun(y);
end
function [value,isterminal,direction] = stopfun(~,y)
  [~,eventvalue] = time_consuming_fun(y);
  value = eventvalue-0.4;
  isterminal = 1;
  direction = 0;
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Stephen23 2019-10-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/483972-is-there-a-way-to-avoid-repeating-time-consuming-calculations-from-odefun-in-the-events-function#answer_395161

编辑：Stephen23 2019-10-7

"Is there a way to do the time-consuming calculations only once at each time step?"

Not really, because numeric ODE solvers do not call the function at specific time points, as the ode45 documentation makes clear: "However, the solver does not step precisely to each point specified in tspan. Instead, the solver uses its own internal steps to compute the solution, then evaluates the solution at the requested points in tspan". So even if you "precalculated" the values at particular timepoints, they would not be used by the solver, because the solver picks its own arbitrary timepoints to use... which there is no trivial way to predict.

You might be able to use interpolation though, depending on the behavior of the function and if it really is slow enough to justify doing: evaluate the actual ode function at the timepoints before calling ode45, pass that data as extra parameter/s, then interpolate it within an anonymous/nested function.

I doubt that memoize would make much difference, as very few function evaluations are likely to repeat exactly the same input values.

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

Are Mjaavatten 2019-10-8

编辑：Are Mjaavatten 2019-10-8

I was relatively sure that there was no easy solution. I submitted the question in the hope that some clever person out there might prove me wrong. Your answer strengthens my suspicion that no such solution exists. Precalculation and interpolation is an interesting idea that may work in simple cases. However, it is unlikely to work in my real problem, where y is a five component vector.

I guess I just have to accept that my problem takes about twice the time it would "ideally" take.

请先登录，再进行评论。