fix callback proposed dt for direction

Author: ChrisRackauckas

src/callbacks.jl

@@ -766,7 +766,7 @@ end
 function apply_callback!(integrator,callback::Union{ContinuousCallback,VectorContinuousCallback},cb_time,prev_sign,event_idx)
 
   if isadaptive(integrator)
-    set_proposed_dt!(integrator, max(nextfloat(integrator.opts.dtmin), callback.dtrelax * integrator.dt))
+    set_proposed_dt!(integrator, integrator.tdir * max(nextfloat(integrator.opts.dtmin), integrator.tdir * callback.dtrelax * integrator.dt))
   end
 
   change_t_via_interpolation!(integrator,integrator.tprev+cb_time)

test/downstream/ode_event_tests.jl

@@ -18,6 +18,21 @@ end
 
 callback = ContinuousCallback(condition,affect!)
 
+sol = solve(prob,Tsit5(),callback=callback)
+@test length(sol) < 20
+
+condition= function (out, u,t,integrator) # Event when event_f(u,t,k) == 0
+  out[1] = - t - 2.95
+end
+
+affect! = function (integrator, idx)
+  if idx == 1
+    integrator.u = integrator.u + 2
+  end
+end
+
+callback = VectorContinuousCallback(condition,affect!,1)
+
 sol = solve(prob,Tsit5(),callback=callback)
 
 f = function (du,u,p,t)
@@ -260,7 +275,7 @@ end
 prob = ODEProblem(model, u0, tspan, perror)
 integrator = init(
     prob,
-    Rosenbrock23(); 
+    Rosenbrock23();
     callback=CallbackSet(
         PositiveDomain(),
         DiscreteCallback(condition, affect!),
@@ -272,4 +287,3 @@ integrator = init(
 )
 
 sol = solve!(integrator)
-