addat deleteat

Author: ChrisRackauckas

src/DiffEqBase.jl

@@ -124,7 +124,9 @@ module DiffEqBase
 
   export tuples, intervals
 
-  export resize!,deleteat!,full_cache,user_cache,u_cache,du_cache,terminate!,
+  export resize!,deleteat!,addat!,full_cache,user_cache,u_cache,du_cache,
+         resize_non_user_cache!,deleteat_non_user_cache!,addat_non_user_cache!,
+         terminate!,
          add_tstop!,add_saveat!,set_abstol!,
          set_reltol!,get_du,get_dt,get_proposed_dt,modify_proposed_dt!,u_modified!,
          savevalues!

src/integrator_interface.jl

@@ -1,9 +1,13 @@
 resize!(i::DEIntegrator,ii::Int) = error("This method has not been implemented for the integrator")
-deleteat!(i::DEIntegrator,ii::Int) = error("This method has not been implemented for the integrator")
+deleteat!(i::DEIntegrator,ii) = error("This method has not been implemented for the integrator")
+addat!(i::DEIntegrator,ii,val=zeros(length(idxs))) = error("This method has not been implemented for the integrator")
 user_cache(i::DEIntegrator) = error("This method has not been implemented for the integrator")
 u_cache(i::DEIntegrator) = error("This method has not been implemented for the integrator")
 du_cache(i::DEIntegrator) = error("This method has not been implemented for the integrator")
 full_cache(i::DEIntegrator) = error("This method has not been implemented for the integrator")
+resize_non_user_cache!(i::DEIntegrator,ii::Int) = error("This method has not been implemented for the integrator")
+deleteat_non_user_cache!(i::DEIntegrator,idxs) = error("This method has not been implemented for the integrator")
+addat_non_user_cache!(i::DEIntegrator,idxs) = error("This method has not been implemented for the integrator")
 terminate!(i::DEIntegrator) = error("This method has not been implemented for the integrator")
 get_du(i::DEIntegrator) = error("This method has not been implemented for the integrator")
 get_dt(i::DEIntegrator) = error("This method has not been implemented for the integrator")
@@ -16,6 +20,13 @@ add_saveat!(i::DEIntegrator,t) = error("This method has not been implemented for
 set_abstol!(i::DEIntegrator,t) = error("This method has not been implemented for the integrator")
 set_reltol!(i::DEIntegrator,t) = error("This method has not been implemented for the integrator")
 
+### Addat isn't a real thing. Let's make it a real thing Gretchen
+
+function addat!(a::AbstractArray,idxs,val=zeros(length(idxs)))
+  flip_range = UnitRange(idxs.start,idxs.start-length(idxs))
+  splice!(a,flip_range,val)
+end
+
 ### Abstract Interface
 
 immutable IntegratorTuples{I}