feat: propagate algebraic_only kwarg in InitializationProblem

Author: AayushSabharwal

src/systems/diffeqs/abstractodesystem.jl

@@ -1296,6 +1296,7 @@ function InitializationProblem{iip, specialize}(sys::AbstractSystem,
         use_scc = true,
         allow_incomplete = false,
         force_time_independent = false,
+        algebraic_only = false,
         kwargs...) where {iip, specialize}
     if !iscomplete(sys)
         error("A completed system is required. Call `complete` or `structural_simplify` on the system before creating an `ODEProblem`")
@@ -1306,12 +1307,12 @@ function InitializationProblem{iip, specialize}(sys::AbstractSystem,
     elseif isempty(u0map) && get_initializesystem(sys) === nothing
         isys = generate_initializesystem(
                 sys; initialization_eqs, check_units, pmap = parammap,
-                guesses, extra_metadata = (; use_scc))
+                guesses, extra_metadata = (; use_scc), algebraic_only)
         simplify_system = true
     else
         isys = generate_initializesystem(
                 sys; u0map, initialization_eqs, check_units,
-                pmap = parammap, guesses, extra_metadata = (; use_scc))
+                pmap = parammap, guesses, extra_metadata = (; use_scc), algebraic_only)
         simplify_system = true
     end