When solving an ODE with the DifferentialEquations package, is there a way to automatically retrieve the state transition matrix, i.e. the sensitivity of the solution w.r.t. the initial conditions, along the entire trajectory?
I would need something similar to ODEForwardSensitivityProblem in Direct Forward Sensitivity Analysis of ODEs, but capable of computing the sensitivities w.r.t. u0 rather than p.
Right now there isn’t a nice high level utility for this. I’ve been asked to extend ODEForwardSensitivityProblem for this but just haven’t gotten around to it yet. As Steven suggests, just using dual numbers on u0 and extracting it manually is the current technique, though I do hope to just make ODEForwardSensitivityProblem handle this sometime in the future.