Merge pull request #3856 from SciML/fix-structural-simplify-deprecation

Fix deprecated ODESystem and @mtkbuild usage in documentation examples

Author: ChrisRackauckas

README.md

@@ -46,7 +46,7 @@ eqs = [D(D(x)) ~ σ * (y - x),
     D(y) ~ x * (ρ - z) - y,
     D(z) ~ x * y - β * z]
 
-@mtkbuild sys = ODESystem(eqs, t)
+@mtkcompile sys = System(eqs, t)
 
 u0 = [D(x) => 2.0,
     x => 1.0,
@@ -83,13 +83,13 @@ eqs = [D(x) ~ σ * (y - x),
     D(y) ~ x * (ρ - z) - y,
     D(z) ~ x * y - β * z]
 
-@named lorenz1 = ODESystem(eqs, t)
-@named lorenz2 = ODESystem(eqs, t)
+@named lorenz1 = System(eqs, t)
+@named lorenz2 = System(eqs, t)
 
 @variables a(t)
 @parameters γ
 connections = [0 ~ lorenz1.x + lorenz2.y + a * γ]
-@mtkbuild connected = ODESystem(connections, t, systems = [lorenz1, lorenz2])
+@mtkcompile connected = System(connections, t, systems = [lorenz1, lorenz2])
 
 u0 = [lorenz1.x => 1.0,
     lorenz1.y => 0.0,

demo.jl

@@ -20,7 +20,7 @@ using ModelingToolkit: t_nounits as t, D_nounits as D
 eqs = [D(D(x)) ~ λ * x
        D(D(y)) ~ λ * y - g
        x^2 + y^2 ~ 1]
-@mtkbuild pend = System(eqs, t)
+@mtkcompile pend = System(eqs, t)
 prob = ODEProblem(pend, [x => -1, y => 0], (0.0, 10.0), [g => 1], guesses = [λ => 1])
 
 sol = solve(prob, FBDF())
@@ -37,7 +37,7 @@ eqs = [
     D(z) ~ x * y - β * z + 0.1z * a
 ]
 
-@mtkbuild sys1 = System(eqs, t)
+@mtkcompile sys1 = System(eqs, t)
 
 eqs = [
     D(x) ~ σ * (y - x),
@@ -49,7 +49,7 @@ noiseeqs = [0.1*x;
             0.1*y;
             0.1*z;;]
 
-@mtkbuild sys2 = SDESystem(eqs, noiseeqs, t)
+@mtkcompile sys2 = SDESystem(eqs, noiseeqs, t)
 
 u0 = [
     x => 1.0,
@@ -73,7 +73,7 @@ odeprob = ODEProblem(sys1, u0, (0.0, 10.0), p; check_compatibility = false)
 eqs = [0 ~ σ * (y - x),
     y ~ x * (ρ - z),
     β * z ~ x * y]
-@mtkbuild sys = System(eqs)
+@mtkcompile sys = System(eqs)
 
 ## ImplicitDiscrete Affects
 
@@ -83,7 +83,7 @@ eqs = [D(D(x)) ~ λ * x
        D(D(y)) ~ λ * y - g
        x^2 + y^2 ~ 1]
 c_evt = [t ~ 5.0] => [x ~ Pre(x) + 0.1]
-@mtkbuild pend = System(eqs, t, continuous_events = c_evt)
+@mtkcompile pend = System(eqs, t, continuous_events = c_evt)
 prob = ODEProblem(pend, [x => -1, y => 0], (0.0, 10.0), [g => 1], guesses = [λ => 1])
 
 sol = solve(prob, FBDF())
@@ -104,4 +104,4 @@ function SysC(; name)
     @named subsys = SysB(; var1 = x)
     return System([D(x) ~ x], t; systems = [subsys], name)
 end
-@mtkbuild sys = SysC()
+@mtkcompile sys = SysC()