I want to study effects on a DDE model when I add another component to the system.
Using resize! and callbacks doesn’t work.
Therefore I want to solve the DDE model for n components across some time span, then use this as the initial condition and system history for a second DDE problem with n + 1 components. (The history of the n+1 component would be also defined here, probably something constant).
However it is not clear to me how to use the output of solve(prob, alg) to define the history function h(p, t) to achieve this, as solve returns a vector with (often) variable timesteps between values.
From reading the DifferentialEquations.jl docs, it looks like h is assumed to be a continuous function, so presumably the internal mechanisms in solve create some continuous interpolation between solved values as a proxy for h? If this is true I’d like to replicate this as close as possible.
Does anyone have any suggestions?