|
81 | 81 | D(x) ~ -x + u
|
82 | 82 | y ~ x
|
83 | 83 | z(k) ~ z(k - 2) + yd(k - 2)]
|
84 |
| - @mtkbuild sys = ODESystem(eqs, t; parameter_dependencies = [kq => 2kp]) |
85 |
| - |
86 |
| - Tf = 1.0 |
87 |
| - prob = ODEProblem(sys, [x => 0.0, y => 0.0], (0.0, Tf), |
88 |
| - [kp => 1.0; z(k - 1) => 3.0; yd(k - 1) => 0.0; z(k - 2) => 4.0; yd(k - 2) => 2.0]) |
89 |
| - @test_nowarn solve(prob, Tsit5(); kwargshandle = KeywordArgSilent) |
90 |
| - |
91 |
| - @mtkbuild sys = ODESystem(eqs, t; parameter_dependencies = [kq => 2kp], |
92 |
| - discrete_events = [[0.5] => [kp ~ 2.0]]) |
93 |
| - prob = ODEProblem(sys, [x => 0.0, y => 0.0], (0.0, Tf), |
94 |
| - [kp => 1.0; z(k - 1) => 3.0; yd(k - 1) => 0.0; z(k - 2) => 4.0; yd(k - 2) => 2.0]) |
95 |
| - @test prob.ps[kp] == 1.0 |
96 |
| - @test prob.ps[kq] == 2.0 |
97 |
| - @test_nowarn solve(prob, Tsit5(), kwargshandle = KeywordArgSilent) |
98 |
| - prob = ODEProblem(sys, [x => 0.0, y => 0.0], (0.0, Tf), |
99 |
| - [kp => 1.0; z(k - 1) => 3.0; yd(k - 1) => 0.0; z(k - 2) => 4.0; yd(k - 2) => 2.0]) |
100 |
| - integ = init(prob, Tsit5(), kwargshandle = KeywordArgSilent) |
101 |
| - @test integ.ps[kp] == 1.0 |
102 |
| - @test integ.ps[kq] == 2.0 |
103 |
| - step!(integ, 0.6) |
104 |
| - @test integ.ps[kp] == 2.0 |
105 |
| - @test integ.ps[kq] == 4.0 |
| 84 | + @test_throws ModelingToolkit.HybridSystemNotSupportedExcpetion @mtkbuild sys = ODESystem( |
| 85 | + eqs, t; parameter_dependencies = [kq => 2kp]) |
| 86 | + |
| 87 | + @test_skip begin |
| 88 | + Tf = 1.0 |
| 89 | + prob = ODEProblem(sys, [x => 0.0, y => 0.0], (0.0, Tf), |
| 90 | + [kp => 1.0; z(k - 1) => 3.0; yd(k - 1) => 0.0; z(k - 2) => 4.0; |
| 91 | + yd(k - 2) => 2.0]) |
| 92 | + @test_nowarn solve(prob, Tsit5(); kwargshandle = KeywordArgSilent) |
| 93 | + |
| 94 | + @mtkbuild sys = ODESystem(eqs, t; parameter_dependencies = [kq => 2kp], |
| 95 | + discrete_events = [[0.5] => [kp ~ 2.0]]) |
| 96 | + prob = ODEProblem(sys, [x => 0.0, y => 0.0], (0.0, Tf), |
| 97 | + [kp => 1.0; z(k - 1) => 3.0; yd(k - 1) => 0.0; z(k - 2) => 4.0; |
| 98 | + yd(k - 2) => 2.0]) |
| 99 | + @test prob.ps[kp] == 1.0 |
| 100 | + @test prob.ps[kq] == 2.0 |
| 101 | + @test_nowarn solve(prob, Tsit5(), kwargshandle = KeywordArgSilent) |
| 102 | + prob = ODEProblem(sys, [x => 0.0, y => 0.0], (0.0, Tf), |
| 103 | + [kp => 1.0; z(k - 1) => 3.0; yd(k - 1) => 0.0; z(k - 2) => 4.0; |
| 104 | + yd(k - 2) => 2.0]) |
| 105 | + integ = init(prob, Tsit5(), kwargshandle = KeywordArgSilent) |
| 106 | + @test integ.ps[kp] == 1.0 |
| 107 | + @test integ.ps[kq] == 2.0 |
| 108 | + step!(integ, 0.6) |
| 109 | + @test integ.ps[kp] == 2.0 |
| 110 | + @test integ.ps[kq] == 4.0 |
| 111 | + end |
106 | 112 | end
|
107 | 113 |
|
108 | 114 | @testset "SDESystem" begin
|
|
0 commit comments