feat: MTK Jacobian for CoupledDEs

[oameye]

resolves #228

[codecov-commenter]

Codecov Report

Attention: Patch coverage is 91.66667% with 1 line in your changes missing coverage. Please review.

Project coverage is 86.43%. Comparing base (6c98239) to head (95e9738).
Report is 36 commits behind head on main.

Files with missing lines	Patch %	Lines
ext/src/CoupledSDEs.jl	75.00%	1 Missing ⚠️

Additional details and impacted files

@@            Coverage Diff             @@
##             main     #229      +/-   ##
==========================================
+ Coverage   82.00%   86.43%   +4.42%     
==========================================
  Files          15       17       +2     
  Lines         717      936     +219     
==========================================
+ Hits          588      809     +221     
+ Misses        129      127       -2     

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

[Datseris]

sorry for the late reply!

[Datseris]

This will be nothing if the system is not constructed through MTK. Doesn't this break normal systems? You cna't do prob.f then.

[oameye]

It think you are confused with referrenced_sciml_model. Every ContineousTimeDynamicalSystem is made through a DEProblem. See

DynamicalSystemsBase.jl/src/core_systems/continuous_time_ode.jl

Lines 82 to 88 in 7b7ac7a

	function CoupledODEs(f, u0, p = SciMLBase.NullParameters(); t0 = 0, diffeq = DEFAULT_DIFFEQ)
	IIP = isinplace(f, 4) # from SciMLBase
	s = correct_state(Val{IIP}(), u0)
	T = eltype(s)
	prob = ODEProblem{IIP}(f, s, (T(t0), T(Inf)), p)
	return CoupledODEs(prob, diffeq)
	end