I have a DAE built with ModelingToolkit. Between solves I re-initialize from a new operating point. For variables that have an initial constraint*I can update them in place via the Initial parameter:
prob.ps[Initial(x)] = 1.0 # or setp(prob, Initial(x))(prob, 1.0)
For the algebraic unknowns I only provide guesses (they’re solved by the initialization, not constrained). I’d like to update those starting guesses between solves — but there doesn’t seem to be a settable equivalent to Initial(x) for guesses.
As far as I can tell I have to remake the entire problem, parameters used in guesses get pruned by mtkcompile and cannot be used.
Is there a lighter-weight, supported way to mutate the guess values on an existing problem (something like a Guess(x) parameter), or is remaking the problem the intended path?