[Octave-bug-tracker] [bug #49890] error if event function stops ode45

[Rik]

Update of bug #49890 (project octave):

                  Status:                    None => Works For Me           

    _______________________________________________________

Follow-up Comment #3:

Without an example from the original reporter, we can close this bug.  It
works for me with a trivial event function


function [value, isterminal, direction] = evfcn (t, y)

  isterminal = 1;
  direction = 0;

  if (t < 1)
    value = 1;
  else
    value = 0;
  endif

end


and the call


opts = odeset ('Events', @evfcn);
[t,y,te,ye,ie] = ode45 (@(t,y) 2, [0, 2], 1, opts)
warning: Solver was stopped by a call of 'break' in the main iteration loop at
time t = 1.146792 before the endpoint at tend = 2.000000 was reached.  This
may happen because the @odeplot function returned 'true' or the @event
function returned 'true'.
warning: called from
    integrate_adaptive at line 319 column 7
    ode45 at line 227 column 12
t =

   0.00000
   0.14679
   0.34679
   0.54679
   0.74679
   0.94679
   1.14679

y =

   1.0000
   1.2936
   1.6936
   2.0936
   2.4936
   2.8936
   3.2936

te =  1.1468
ye =  3.2936
ie =  1




    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?49890>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/