Merge pull request #913 from AayushSabharwal/as/no-init

feat: support `NoInit` in `get_initial_values`

Author: ChrisRackauckas

src/initialization.jl

@@ -271,6 +271,16 @@ function get_initial_values(prob, valp, f, alg::OverrideInit,
     return u0, p, success
 end
 
+"""
+    $(TYPEDSIGNATURES)
+
+Do not run any initialization routine, and simply return the state and parameter values
+stored in `integrator` successfully.
+"""
+function get_initial_values(prob, integrator, f, ::NoInit, iip; kwargs...)
+    return state_values(integrator), parameter_values(integrator), true
+end
+
 is_trivial_initialization(::Nothing) = true
 
 function is_trivial_initialization(initdata::OverrideInitData)

test/initialization.jl

@@ -286,3 +286,13 @@ end
         @test success
     end
 end
+
+@testset "NoInit" begin
+    prob = ODEProblem(ones(2), (0.0, 1.0), ones(2)) do u, p, t
+        return u
+    end
+    u, p, success = SciMLBase.get_initial_values(prob, prob, prob.f, SciMLBase.NoInit(), Val(true))
+    @test u == ones(2)
+    @test p == ones(2)
+    @test success
+end