fix floating point issues in tstops

Author: ChrisRackauckas

src/callbacks.jl

@@ -35,8 +35,10 @@ end
   # Check if the event occured
   if typeof(callback.idxs) <: Void
     previous_condition = callback.condition(integrator.tprev,integrator.uprev,integrator)
-  else
+  elseif typeof(callback.idxs) <: Number
     previous_condition = callback.condition(integrator.tprev,integrator.uprev[callback.idxs],integrator)
+  else
+    previous_condition = callback.condition(integrator.tprev,@view(integrator.uprev[callback.idxs]),integrator)
   end
   if isapprox(previous_condition,0,rtol=callback.reltol,atol=callback.abstol)
     prev_sign = 0.0
@@ -45,9 +47,11 @@ end
   end
   prev_sign_index = 1
   if typeof(callback.idxs) <: Void
-    next_sign = sign(callback.condition(integrator.tprev+integrator.dt,integrator.u,integrator))
+    next_sign = sign(callback.condition(integrator.t,integrator.u,integrator))
+  elseif typeof(callback.idxs) <: Number
+    next_sign = sign(callback.condition(integrator.t,integrator.u[callback.idxs],integrator))
   else
-    next_sign = sign(callback.condition(integrator.tprev+integrator.dt,integrator.u[callback.idxs],integrator))
+    next_sign = sign(callback.condition(integrator.t,@view(integrator.u[callback.idxs]),integrator))
   end
   if ((prev_sign<0 && !(typeof(callback.affect!)<:Void)) || (prev_sign>0 && !(typeof(callback.affect_neg!)<:Void))) && prev_sign*next_sign<0
     event_occurred = true
@@ -183,7 +187,7 @@ end
   end
 
   integrator.u_modified = true
-  if callback.condition(integrator.tprev+integrator.dt,integrator.u,integrator)
+  if callback.condition(integrator.t,integrator.u,integrator)
     callback.affect!(integrator)
     @inbounds if callback.save_positions[2]
       savevalues!(integrator)

src/integrators/integrator_utils.jl

@@ -175,7 +175,12 @@ end
     integrator.accept_step = (!integrator.isout && integrator.EEst <= 1.0) || (integrator.opts.force_dtmin && abs(integrator.dt) <= abs(integrator.opts.dtmin))
     if integrator.accept_step # Accept
       integrator.tprev = integrator.t
-      integrator.t = ttmp
+      if typeof(integrator.t)<:AbstractFloat && !isempty(integrator.opts.tstops)
+        tstop = top(integrator.opts.tstops)
+        abs(ttmp - tstop) < 10eps(integrator.t) ? (integrator.t = tstop) : (integrator.t = ttmp)
+      else
+        integrator.t = ttmp
+      end
       integrator.qold = max(integrator.EEst,integrator.opts.qoldinit)
       calc_dt_propose!(integrator,dtnew)
       handle_callbacks!(integrator)