feat: forward kwargs to solve in OverrideInit

Author: AayushSabharwal

src/initialization.jl

@@ -304,7 +304,8 @@ function initialization_status(prob::AbstractSciMLProblem)
     iprob = prob.f.initialization_data.initializeprob
     isnothing(prob) && return nothing
 
-    nunknowns = iprob.u0 === nothing ? 0 : length(iprob.u0)
+    iu0 = state_values(iprob)
+    nunknowns = iu0 === nothing ? 0 : length(iu0)
     neqs = if __has_resid_prototype(iprob.f) && iprob.f.resid_prototype !== nothing
         length(iprob.f.resid_prototype)
     else

test/initialization.jl

@@ -203,6 +203,14 @@ end
             @test p ≈ 1.0
             @test success
         end
+        @testset "with kwargs provided to `get_initial_values`" begin
+            u0, p, success = SciMLBase.get_initial_values(
+                prob, integ, fn, SciMLBase.OverrideInit(),
+                Val(false); nlsolve_alg = NewtonRaphson(), abstol, reltol, u0 = [-1.0, 1.0])
+            @test u0 ≈ [2.0, -2.0]
+            @test p ≈ 1.0
+            @test success
+        end
     end
 
     @testset "Solves with non-integrator value provider" begin