Merge branch 'master' into myb/unit

Author: ChrisRackauckas

.github/dependabot.yml

@@ -5,3 +5,6 @@ updates:
     directory: "/" # Location of package manifests
     schedule:
       interval: "weekly"
+    ignore:
+      - dependency-name: "crate-ci/typos"
+        update-types: ["version-update:semver-patch"]

.github/workflows/SpellCheck.yml

@@ -0,0 +1,13 @@
+name: Spell Check
+
+on: [pull_request]
+
+jobs:
+  typos-check:
+    name: Spell Check with Typos
+    runs-on: ubuntu-latest
+    steps:
+      - name: Checkout Actions Repository
+        uses: actions/checkout@v4
+      - name: Check spelling
+        uses: crate-ci/[email protected] 

.github/workflows/ci.yml

@@ -21,9 +21,12 @@ jobs:
   test:
     runs-on: ubuntu-latest
     strategy:
+      fail-fast: false
       matrix:
         group:
-          - All
+          - InterfaceI
+          - InterfaceII
+          - Extensions
         version:
           - '1'
     steps:
@@ -43,6 +46,8 @@ jobs:
             ${{ runner.os }}-
       - uses: julia-actions/julia-buildpkg@v1
       - uses: julia-actions/julia-runtest@v1
+        env:
+          GROUP: ${{ matrix.group }}
       - uses: julia-actions/julia-processcoverage@v1
       - uses: codecov/codecov-action@v3
         with:

.typos.toml

@@ -0,0 +1,6 @@
+[default.extend-words]
+nin = "nin"
+nd = "nd"
+Strat = "Strat"
+eles = "eles"
+ser = "ser"

Project.toml

@@ -1,7 +1,7 @@
 name = "ModelingToolkit"
 uuid = "961ee093-0014-501f-94e3-6117800e7a78"
 authors = ["Yingbo Ma <[email protected]>", "Chris Rackauckas <[email protected]> and contributors"]
-version = "8.72.2"
+version = "8.75.0"
 
 [deps]
 AbstractTrees = "1520ce14-60c1-5f80-bbc7-55ef81b5835c"
@@ -69,6 +69,7 @@ DataStructures = "0.17, 0.18"
 DiffEqBase = "6.103.0"
 DiffEqCallbacks = "2.16"
 DiffRules = "0.1, 1.0"
+Distributed = "1"
 Distributions = "0.23, 0.24, 0.25"
 DocStringExtensions = "0.7, 0.8, 0.9"
 DomainSets = "0.6"
@@ -77,26 +78,31 @@ ForwardDiff = "0.10.3"
 FunctionWrappersWrappers = "0.1"
 Graphs = "1.5.2"
 IfElse = "0.1"
+InteractiveUtils = "1"
 JuliaFormatter = "1"
 JumpProcesses = "9.1"
 LabelledArrays = "1.3"
 Latexify = "0.11, 0.12, 0.13, 0.14, 0.15, 0.16"
+Libdl = "1"
+LinearAlgebra = "1"
 MLStyle = "0.4.17"
 MacroTools = "0.5"
 NaNMath = "0.3, 1"
 OrdinaryDiffEq = "6"
 PrecompileTools = "1"
-RecursiveArrayTools = "2.3"
+RecursiveArrayTools = "2.3, 3"
 Reexport = "0.2, 1"
 RuntimeGeneratedFunctions = "0.5.9"
 SciMLBase = "2.0.1"
+Serialization = "1"
 Setfield = "0.7, 0.8, 1"
-SimpleNonlinearSolve = "0.1.0"
+SimpleNonlinearSolve = "0.1.0, 1"
+SparseArrays = "1"
 SpecialFunctions = "0.7, 0.8, 0.9, 0.10, 1.0, 2"
 StaticArrays = "0.10, 0.11, 0.12, 1.0"
-SymbolicIndexingInterface = "0.1, 0.2"
+SymbolicIndexingInterface = "0.3.1"
 SymbolicUtils = "1.0"
-Symbolics = "5.0"
+Symbolics = "5.7"
 URIs = "1"
 UnPack = "0.1, 1.0"
 Unitful = "1.1"
@@ -105,7 +111,6 @@ julia = "1.9"
 [extras]
 AmplNLWriter = "7c4d4715-977e-5154-bfe0-e096adeac482"
 BenchmarkTools = "6e4b80f9-dd63-53aa-95a3-0cdb28fa8baf"
-BifurcationKit = "0f109fa4-8a5d-4b75-95aa-f515264e7665"
 ControlSystemsMTK = "687d7614-c7e5-45fc-bfc3-9ee385575c88"
 DeepDiffs = "ab62b9b5-e342-54a8-a765-a90f495de1a6"
 ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"
@@ -116,6 +121,7 @@ NonlinearSolve = "8913a72c-1f9b-4ce2-8d82-65094dcecaec"
 Optimization = "7f7a1694-90dd-40f0-9382-eb1efda571ba"
 OptimizationMOI = "fd9f6733-72f4-499f-8506-86b2bdd0dea1"
 OptimizationOptimJL = "36348300-93cb-4f02-beb5-3c3902f8871e"
+Pkg = "44cfe95a-1eb2-52ea-b672-e2afdf69b78f"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 ReferenceTests = "324d217c-45ce-50fc-942e-d289b448e8cf"
 SafeTestsets = "1bc83da4-3b8d-516f-aca4-4fe02f6d838f"
@@ -128,4 +134,4 @@ Sundials = "c3572dad-4567-51f8-b174-8c6c989267f4"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [targets]
-test = ["AmplNLWriter", "BenchmarkTools", "BifurcationKit", "ControlSystemsMTK", "NonlinearSolve", "ForwardDiff", "Ipopt", "Ipopt_jll", "ModelingToolkitStandardLibrary", "Optimization", "OptimizationOptimJL", "OptimizationMOI", "Random", "ReferenceTests", "SafeTestsets", "StableRNGs", "Statistics", "SteadyStateDiffEq", "Test", "StochasticDiffEq", "Sundials", "StochasticDelayDiffEq"]
+test = ["AmplNLWriter", "BenchmarkTools", "ControlSystemsMTK", "NonlinearSolve", "ForwardDiff", "Ipopt", "Ipopt_jll", "ModelingToolkitStandardLibrary", "Optimization", "OptimizationOptimJL", "OptimizationMOI", "Random", "ReferenceTests", "SafeTestsets", "StableRNGs", "Statistics", "SteadyStateDiffEq", "Test", "StochasticDiffEq", "Sundials", "StochasticDelayDiffEq", "Pkg"]

docs/Project.toml

@@ -27,14 +27,14 @@ Distributions = "0.25"
 Documenter = "1"
 ModelingToolkit = "8.33"
 ModelingToolkitDesigner = "1"
-NonlinearSolve = "0.3, 1, 2"
+NonlinearSolve = "0.3, 1, 2, 3"
 Optim = "1.7"
 Optimization = "3.9"
 OptimizationOptimJL = "0.1"
 OrdinaryDiffEq = "6.31"
 Plots = "1.36"
 StochasticDiffEq = "6"
-StructuralIdentifiability = "0.4"
+StructuralIdentifiability = "0.4, 0.5"
 SymbolicUtils = "1"
 Symbolics = "5"
 Unitful = "1.12"

docs/src/basics/AbstractSystem.md

@@ -121,7 +121,7 @@ patterns via an abstract interpretation without requiring differentiation.
 
 At the end, the system types have `DEProblem` constructors, like `ODEProblem`,
 which allow for directly generating the problem types required for numerical
-methods. The first argument is always the `AbstractSystem`, and the proceding
+methods. The first argument is always the `AbstractSystem`, and the next
 arguments match the argument order of their original constructors. Whenever an
 array would normally be provided, such as `u0` the initial condition of an
 `ODEProblem`, it is instead replaced with a variable map, i.e., an array of

docs/src/basics/Linearization.md

@@ -51,7 +51,7 @@ If linearization is to be performed around multiple operating points, the simpli
 
 ## Symbolic linearization
 
-The function [`ModelingToolkit.linearize_symbolic`](@ref) works simiar to [`ModelingToolkit.linearize`](@ref) but returns symbolic rather than numeric Jacobians. Symbolic linearization have several limitations and no all systems that can be linearized numerically can be linearized symbolically.
+The function [`ModelingToolkit.linearize_symbolic`](@ref) works similar to [`ModelingToolkit.linearize`](@ref) but returns symbolic rather than numeric Jacobians. Symbolic linearization have several limitations and no all systems that can be linearized numerically can be linearized symbolically.
 
 ## Input derivatives
 

docs/src/basics/MTKModel_Connector.md

@@ -29,7 +29,7 @@ equations.
   - `@icon` : for embedding the model icon
   - `@parameters`: for specifying the symbolic parameters
   - `@structural_parameters`: for specifying non-symbolic parameters
-  - `@variables`: for specifing the states
+  - `@variables`: for specifying the states
 
 Let's explore these in more detail with the following example:
 
@@ -104,7 +104,7 @@ end
 
 #### `@structural_parameters` begin block
 
-  - This block is for non symbolic input arguements. These are for inputs that usually are not meant to be part of components; but influence how they are defined. One can list inputs like boolean flags, functions etc... here.
+  - This block is for non symbolic input arguments. These are for inputs that usually are not meant to be part of components; but influence how they are defined. One can list inputs like boolean flags, functions etc... here.
   - Whenever default values are specified, unlike parameters/variables, they are reflected in the keyword argument list.
 
 #### `@parameters` and `@variables` begin block
@@ -251,3 +251,104 @@ Dict{Symbol, Any} with 7 entries:
   :extend               => Any[[:p2, :p1], Symbol("#mtkmodel__anonymous__ModelB"), :ModelB]
   :equations            => ["model_a.k ~ f(v)"]
 ```
+
+### Using conditional statements
+
+#### Conditional elements of the system
+
+Both `@mtkmodel` and `@connector` support conditionally defining parameters,
+variables, equations, and components.
+
+The if-elseif-else statements can be used inside `@equations`, `@parameters`,
+`@variables`, `@components`.
+
+```@example branches-in-components
+using ModelingToolkit
+
+@mtkmodel C begin end
+
+@mtkmodel BranchInsideTheBlock begin
+    @structural_parameters begin
+        flag = true
+    end
+    @parameters begin
+        if flag
+            a1
+        else
+            a2
+        end
+    end
+    @components begin
+        if flag
+            sys1 = C()
+        else
+            sys2 = C()
+        end
+    end
+end
+```
+
+Alternatively, the `@equations`, `@parameters`, `@variables`, `@components` can be
+used inside the if-elseif-else statements.
+
+```@example branches-in-components
+@mtkmodel BranchOutsideTheBlock begin
+    @structural_parameters begin
+        flag = true
+    end
+    if flag
+        @parameters begin
+            a1
+        end
+        @components begin
+            sys1 = C()
+        end
+        @equations begin
+            a1 ~ 0
+        end
+    else
+        @parameters begin
+            a2
+        end
+        @equations begin
+            a2 ~ 0
+        end
+    end
+end
+```
+
+The conditional parts are reflected in the `structure`. For `BranchOutsideTheBlock`, the metadata is:
+
+```julia
+julia> BranchOutsideTheBlock.structure
+Dict{Symbol, Any} with 5 entries:
+  :components           => Any[(:if, :flag, [[:sys1, :C]], Any[])]
+  :kwargs               => Dict{Symbol, Any}(:flag=>true)
+  :independent_variable => t
+  :parameters           => Dict{Symbol, Dict{Symbol, Any}}(:a1=>Dict(:condition=>(:if, :flag, Dict{Symbol, Any}(:kwargs => Dict{Any, Any}(:a1 => nothing), :parameters => Any[Dict{Symbol, Dict{Symbol, Any}}(:a1 => Dict())]), Dict{Symbol, Any}(:kwargs => Dict{Any, Any}(:a2 => nothing), :parameters => Any[Dict{Symbol, Dict{Symbol, Any}}(:a2 => Dict())]))
+  :equations            => Any[(:if, :flag, ["a1 ~ 0"], ["a2 ~ 0"])]
+```
+
+Conditional entries are entered in the format of `(branch, condition, [case when it is true], [case when it is false])`;
+where `branch` is either `:if` or `:elseif`.<br>
+The `[case when it is false]` is either an empty vector or `nothing` when only if branch is
+present; it is a vector or dictionary whenever else branch is present; it is a conditional tuple
+whenever elseif branches are present.
+
+For the conditional components and equations these condition tuples are added
+directly, while for parameters and variables these are added as `:condition` metadata.
+
+#### Conditional initial guess of symbolic variables
+
+Using ternary operator or if-elseif-else statement, conditional initial guesses can be assigned to parameters and variables.
+
+```@example branches-in-components
+@mtkmodel DefaultValues begin
+    @structural_parameters begin
+        flag = true
+    end
+    @parameters begin
+        p = flag ? 1 : 2
+    end
+end
+```

docs/src/basics/Variable_metadata.md

@@ -83,6 +83,19 @@ hasbounds(u)
 getbounds(u)
 ```
 
+## Guess
+
+Specify an initial guess for custom initial conditions of an `ODESystem`.
+
+```@example metadata
+@variables u [guess = 1]
+hasguess(u)
+```
+
+```@example metadata
+getguess(u)
+```
+
 ## Mark input as a disturbance
 
 Indicate that an input is not available for control, i.e., it's a disturbance input.