feat: add supports_initialization trait

Author: AayushSabharwal

src/systems/abstractsystem.jl

@@ -676,6 +676,8 @@ function SymbolicUtils.maketerm(::Type{<:BasicSymbolic}, ::Initial, args, meta)
     return metadata(val, meta)
 end
 
+supports_initialization(sys::AbstractSystem) = true
+
 function add_initialization_parameters(sys::AbstractSystem)
     @assert !has_systems(sys) || isempty(get_systems(sys))
     is_initializesystem(sys) && return sys

src/systems/discrete_system/discrete_system.jl

@@ -421,3 +421,5 @@ end
 function DiscreteFunctionExpr(sys::DiscreteSystem, args...; kwargs...)
     DiscreteFunctionExpr{true}(sys, args...; kwargs...)
 end
+
+supports_initialization(::DiscreteSystem) = false

src/systems/jumps/jumpsystem.jl

@@ -685,3 +685,5 @@ function (ratemap::JumpSysMajParamMapper{U, V, W})(maj::MassActionJump, newparam
     scale_rates && JumpProcesses.scalerates!(maj.scaled_rates, maj.reactant_stoch)
     nothing
 end
+
+supports_initialization(::JumpSystem) = false

src/systems/optimization/constraints_system.jl

@@ -246,3 +246,5 @@ function get_cmap(sys::ConstraintsSystem, exprs = nothing)
     cmap = map(x -> x ~ getdefault(x), cs)
     return cmap, cs
 end
+
+supports_initialization(::ConstraintsSystem) = false

src/systems/optimization/optimizationsystem.jl

@@ -751,3 +751,5 @@ function structural_simplify(sys::OptimizationSystem; split = true, kwargs...)
     sys = complete(sys; split)
     return sys
 end
+
+supports_initialization(::OptimizationSystem) = false