MATLAB: How to stop an ode-solver if integration takes too long

interrupt ode solver

Hi all,

I am working with implicit ode15i solver.

I receive results within seconds for more than 80 parameters to solve. However, sometimes it seems it gets stuck – no progress after hours. I can reproduce that, so I am sure I have chosen bad initial conditions. Is there a way to interrupt the integration after certain time to start new with different initial conditions? Something like a timing option?

I know about the event-function but this is for events given by the system only as far as I understood.

Thanks a lot in advance. Best, Franziska

Best Answer

At the start of your code:

StartTime=clock;

At the end of each loop/the point you might want to exit

TimeElapsed=clock-StartTime;
if TimeElapsed(end)>10 %Set it to a value that you want (I chose 10 seconds)
return
end

Look at the clock function to understand this more.

Related Solutions

MATLAB: Solving ODEs with different sets of initial conditions

function main
  n_times = 1000;
  t_i = 1; t_f = 10; tspan = [t_i,t_f];
  y0 = randn(2, n_times);
  y0 = reshape(y0, 2*n_times, 1);
  [t, y] = ode45(@(t,y) harmosc(t,y,n_times), tspan, y0);
  y = reshape(y,size(y,1),2,n_times)
  plot(t,y(:,1,15),t,y(:,2,15))
end
function dy = harmosc(t,y,n_times)
    y = reshape(y,2,n_times);
    dy = zeros(size(y));
    dy(1,:)=y(2,:);
    dy(2,:)=-y(1,:);  
    dy = reshape(dy,2*n_times,1);
end

MATLAB: How to stop ode integration

function [value,isterminal,direction] = events(t,y)
  miny = [0.0000001, 0.000000025, -83.9];   %for example
  value = y < miny; 
  isterminal = 1;
  direction = 0;

Best Answer

Related Solutions

MATLAB: Solving ODEs with different sets of initial conditions

MATLAB: How to stop ode integration

Related Question