Hi all, What is the purpose of the bounds symbolic metadata in ModelingToolkit? When I add bounds to x in the following code, I would at first expect the solution for x to be limited between 0 and 2, but this doesn’t seem to be the case. @variables t, x(t) [bounds=(0, 2)] D = Differential(t) eqs…

You can already do that, see Modeling Optimization Problems · ModelingToolkit.jl