Use a completely type stable workflow for precompiling

Author: avik-pal

src/BoundaryValueDiffEq.jl

@@ -74,8 +74,13 @@ end
         TwoPointBVProblem(f1, (bc1_a, bc1_b), u0, tspan; bcresid_prototype),
     ]
 
+    jac_alg = BVPJacobianAlgorithm(AutoForwardDiff(; chunksize = 2))
+
     @compile_workload begin
-        for prob in probs, alg in (MIRK2(), MIRK3(), MIRK4(), MIRK5(), MIRK6())
+        for prob in probs,
+            alg in (MIRK2(; jac_alg), MIRK3(; jac_alg), MIRK4(; jac_alg),
+                MIRK5(; jac_alg), MIRK6(; jac_alg))
+
             solve(prob, alg; dt = 0.2)
         end
     end

src/solve/mirk.jl

@@ -139,7 +139,7 @@ function SciMLBase.solve!(cache::MIRKCache)
 
     while SciMLBase.successful_retcode(info) && defect_norm > abstol
         nlprob = __construct_nlproblem(cache, recursive_flatten(y₀))
-        sol_nlprob = solve(nlprob, alg.nlsolve; abstol, kwargs...)
+        sol_nlprob = __solve(nlprob, alg.nlsolve; abstol, kwargs...)
         recursive_unflatten!(cache.y₀, sol_nlprob.u)
 
         info = sol_nlprob.retcode