fix: fix initialization_status for SCCNonlinearProblem initprob

Author: AayushSabharwal

src/initialization.jl

@@ -216,8 +216,10 @@ Keyword arguments:
   If the former is `nothing`, this keyword argument will be used. If it is also not provided,
   an error will be thrown.
 
+All additional keyword arguments are forwarded to `solve`.
+
 In case the initialization problem is trivial, `nlsolve_alg`, `abstol` and `reltol` are
-not required.
+not required. `solve` is also not called.
 """
 function get_initial_values(prob, valp, f, alg::OverrideInit,
         iip::Union{Val{true}, Val{false}}; nlsolve_alg = nothing, abstol = nothing, reltol = nothing, kwargs...)
@@ -257,7 +259,7 @@ function get_initial_values(prob, valp, f, alg::OverrideInit,
         else
             throw(OverrideInitNoTolerance(:reltol))
         end
-        nlsol = solve(initprob, nlsolve_alg; abstol = _abstol, reltol = _reltol)
+        nlsol = solve(initprob, nlsolve_alg; abstol = _abstol, reltol = _reltol, kwargs...)
         success = SciMLBase.successful_retcode(nlsol)
     end
 

test/initialization.jl

@@ -270,6 +270,15 @@ end
         @test success
     end
 
+    @testset "Initialization status for `SCCNonlinearProblem`" begin
+        initprob = SCCNonlinearProblem([initprob], [Returns(nothing)])
+        initialization_data = SciMLBase.OverrideInitData(
+            initprob, nothing, nothing, nothing)
+        fn = ODEFunction(rhs2; initialization_data)
+        prob = ODEProblem(fn, [2.0, 0.0], (0.0, 1.0), 0.0)
+        @test SciMLBase.initialization_status(prob) == SciMLBase.FULLY_DETERMINED
+    end
+
     @testset "Trivial initialization" begin
         initprob = NonlinearProblem(Returns(nothing), nothing, [1.0])
         update_initializeprob! = function (iprob, integ)