I am trying to efficiently reinitialize a SteadyStateProblem, built from a ModelingToolkit System
.
steady_prob = SteadyStateProblem(prob)
steady_prob = remake(steady_prob; u0=new_u0)
But this doesn’t actually change the u0 of the problem. What would be the correct way to update the state? It has to be efficient, so either using setsym
or preferably by providing u0 directly.