I am a bit confused: should we use DAEProblem(sys, du0, u0, tspan) or ODEProblem(sys, u0, tspan) from a differential-algebraic system using ModelingToolkit? Using ODEProblem does work in my case, but then I cannot use DAE solvers like DFBDF.
There is zero benefit to using DFBDF vs FBDF and there are many downsides.
Good to know. Could this be made more clear in the docs?
Which documentation do you think would be most useful to have this note?
It would make sense to have an example DAEProblem, but I can’t find any tutorial on a symbolic DAEProblem, just a modelingtoolkitize-based DAEProblem… It could also make sense to add it directly in the ODESystem docstring, to make it clear that it can be used for DAESystems as well: