@@ -4,13 +4,12 @@ using ModelingToolkit.BipartiteGraphs
4
4
using ModelingToolkit. StructuralTransformations
5
5
6
6
function check_contract (sys)
7
- graph = ModelingToolkit. get_tearing_state (sys). structure. graph
8
- sys = tearing_substitution (sys)
9
- state = TearingState (sys)
7
+ state = ModelingToolkit. get_tearing_state (sys)
8
+ graph = state. structure. graph
10
9
fullvars = state. fullvars
10
+ sys = tearing_substitution (sys)
11
11
12
12
eqs = equations (sys)
13
- var2idx = Dict (enumerate (fullvars))
14
13
for (i, eq) in enumerate (eqs)
15
14
actual = union (ModelingToolkit. vars (eq. lhs), ModelingToolkit. vars (eq. rhs))
16
15
actual = filter (! ModelingToolkit. isparameter, collect (actual))
@@ -143,9 +142,7 @@ sol = solve(prob, Tsit5())
143
142
include (" ../examples/serial_inductor.jl" )
144
143
sys = structural_simplify (ll_model)
145
144
check_contract (sys)
146
- u0 = [inductor1. i => 0.0
147
- inductor2. i => 0.0
148
- inductor2. v => 0.0 ]
145
+ u0 = states (sys) .=> 0
149
146
@test_throws Any ODEProblem (sys, u0, (0 , 10.0 ))
150
147
@test_throws Any ODAEProblem (sys, u0, (0 , 10.0 ))
151
148
prob = DAEProblem (sys, Differential (t).(states (sys)) .=> 0 , u0, (0 , 0.5 ))
0 commit comments