BVP Solvers do not accept certain MTK kwargs (:jac, etc.)

[vyudu]

With codegen to BVProblem being added to MTK soon it'd be nice if the BVP solvers could work with the Jacobian, W matrix, etc. that can be generated from it. Currently doing something like

bvp = SciMLBase.BVProblem(lksys, u0map, tspan; guesses, constraints, jac = true)
solve(bvp, MIRK4(), dt = 0.01)

will result in a Unrecognized keyword arguments: [:jac]

[ErikQQY]

In ModelingToolkit.jl the similar usage for ODESystem(sys::ODESystem, ..., jac=true) is implemented as a dispatch in https://github.com/SciML/ModelingToolkit.jl/blob/2ca9ecfdee95f3026bfbc6b484f9403331615ed7/src/systems/diffeqs/abstractodesystem.jl#L786, so I think this feature need to be implemented in MTK as well.

[vyudu]

Right, that part should be coming with SciML/ModelingToolkit.jl#3323