Merge branch 'SciML:master' into dg/nonlinear

Author: DhairyaLGandhi

.github/workflows/Downstream.yml

@@ -35,14 +35,13 @@ jobs:
           - {user: SciML, repo: StochasticDiffEq.jl, group: AlgConvergence}
           - {user: SciML, repo: Sundials.jl, group: All}
           - {user: SciML, repo: LinearSolve.jl, group: All}
-          - {user: SciML, repo: NonlinearSolve.jl, group: All}
+          - {user: SciML, repo: NonlinearSolve.jl, group: Core}
           - {user: SciML, repo: Optimization.jl, group: All}
           - {user: SciML, repo: Integrals.jl, group: All}
           - {user: SciML, repo: BoundaryValueDiffEq.jl, group: All}
           - {user: SciML, repo: ModelingToolkit.jl, group: All}
           - {user: SciML, repo: ModelingToolkitStandardLibrary.jl, group: Core}
           - {user: SciML, repo: StochasticDelayDiffEq.jl, group: All}
-          - {user: SciML, repo: SimpleNonlinearSolve.jl, group: All}
           - {user: SciML, repo: SimpleDiffEq.jl, group: All}
           - {user: SciML, repo: SciMLSensitivity.jl, group: Core1}
           - {user: SciML, repo: SciMLSensitivity.jl, group: Core2}

Project.toml

@@ -1,7 +1,7 @@
 name = "SciMLBase"
 uuid = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 authors = ["Chris Rackauckas <[email protected]> and contributors"]
-version = "2.75.1"
+version = "2.79.0"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"
@@ -85,7 +85,7 @@ RecipesBase = "1.3.4"
 RecursiveArrayTools = "3.27.2"
 Reexport = "1"
 RuntimeGeneratedFunctions = "0.5.12"
-SciMLOperators = "0.3.7"
+SciMLOperators = "0.3.13"
 SciMLStructures = "1.1"
 StableRNGs = "1.0"
 StaticArrays = "1.7"

src/SciMLBase.jl

@@ -575,21 +575,6 @@ $(TYPEDEF)
 abstract type AbstractSensitivitySolution{T, N, S} <: AbstractTimeseriesSolution{T, N, S} end
 
 # Misc
-# TODO - deprecate AbstractDiffEqOperator family
-"""
-$(TYPEDEF)
-"""
-abstract type AbstractDiffEqOperator{T} <: AbstractSciMLOperator{T} end
-
-"""
-$(TYPEDEF)
-"""
-abstract type AbstractDiffEqLinearOperator{T} <: AbstractDiffEqOperator{T} end
-
-"""
-$(TYPEDEF)
-"""
-abstract type AbstractDiffEqCompositeOperator{T} <: AbstractDiffEqLinearOperator{T} end
 
 """
 $(TYPEDEF)
@@ -720,10 +705,6 @@ $(TYPEDEF)
 abstract type AbstractParameterizedFunction{iip} <: AbstractODEFunction{iip} end
 
 include("retcodes.jl")
-include("operators/operators.jl")
-include("operators/basic_operators.jl")
-include("operators/diffeq_operator.jl")
-include("operators/common_defaults.jl")
 include("symbolic_utils.jl")
 include("performance_warnings.jl")
 
@@ -853,10 +834,6 @@ export EnsembleAnalysis, EnsembleSummary
 
 export tuples, intervals, TimeChoiceIterator
 
-export AffineDiffEqOperator, DiffEqScaledOperator
-
-export DiffEqScalar, DiffEqArrayOperator, DiffEqIdentity
-
 export step!, deleteat!, addat!, get_tmp_cache,
        full_cache, user_cache, u_cache, du_cache,
        rand_cache, ratenoise_cache,

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
 
@@ -304,7 +306,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