@@ -18,7 +18,7 @@ eqs = [
18
18
0 ~ u4 - hypot (u2, u3),
19
19
0 ~ u5 - hypot (u4, u1),
20
20
]
21
- sys = NonlinearSystem (eqs, [u1, u2, u3, u4, u5], [])
21
+ @named sys = NonlinearSystem (eqs, [u1, u2, u3, u4, u5], [])
22
22
sys = initialize_system_structure (sys)
23
23
StructuralTransformations. find_solvables! (sys)
24
24
sss = structure (sys)
@@ -125,7 +125,7 @@ eqs = [
125
125
0 ~ z + y,
126
126
0 ~ x + z,
127
127
]
128
- nlsys = NonlinearSystem (eqs, [x, y, z], [])
128
+ @named nlsys = NonlinearSystem (eqs, [x, y, z], [])
129
129
newsys = tearing (nlsys)
130
130
@test equations (newsys) == [0 ~ z]
131
131
@test isequal (states (newsys), [z])
@@ -142,7 +142,7 @@ eqs = [
142
142
0 ~ x - y
143
143
0 ~ sin (z) + y - p* t
144
144
]
145
- daesys = ODESystem (eqs, t)
145
+ @named daesys = ODESystem (eqs, t)
146
146
newdaesys = tearing (daesys)
147
147
@test equations (newdaesys) == [D (x) ~ z; 0 ~ x + sin (z) - p* t]
148
148
@test isequal (states (newdaesys), [x, z])
@@ -152,7 +152,8 @@ du = [0.0]; u = [1.0]; pr = 0.2; tt = 0.1
152
152
@test du ≈ [- asin (u[1 ] - pr * tt)] atol= 1e-5
153
153
154
154
# test the initial guess is respected
155
- infprob = ODAEProblem (tearing (ODESystem (eqs, t, defaults= Dict (z=> Inf ))), [x=> 1.0 ], (0 , 1.0 ), [p=> 0.2 ])
155
+ @named sys = ODESystem (eqs, t, defaults= Dict (z=> Inf ))
156
+ infprob = ODAEProblem (tearing (sys), [x=> 1.0 ], (0 , 1.0 ), [p=> 0.2 ])
156
157
@test_throws DomainError infprob. f (du, u, pr, tt)
157
158
158
159
sol1 = solve (prob, Tsit5 ())
0 commit comments