More tests

Author: YingboMa

examples/serial_inductor.jl

@@ -14,3 +14,20 @@ eqs = [connect(source.p, resistor.p)
 
 @named ll_model = ODESystem(eqs, t)
 ll_model = compose(ll_model, [source, resistor, inductor1, inductor2, ground])
+
+@named source = ConstantVoltage(V = 10.0)
+@named resistor1 = Resistor(R = 1.0)
+@named resistor2 = Resistor(R = 1.0)
+@named inductor1 = Inductor(L = 1.0e-2)
+@named inductor2 = Inductor(L = 2.0e-2)
+@named ground = Ground()
+
+eqs = [connect(source.p, inductor1.p)
+    connect(inductor1.n, resistor1.p)
+    connect(inductor1.n, resistor2.p)
+    connect(resistor1.n, resistor2.n)
+    connect(resistor2.n, inductor2.p)
+    connect(source.n, inductor2.n)
+    connect(inductor2.n, ground.g)]
+@named ll2_model = ODESystem(eqs, t)
+ll2_model = compose(ll2_model, [source, resistor1, resistor2, inductor1, inductor2, ground])

test/components.jl

@@ -159,6 +159,12 @@ u0 = states(sys) .=> 0
 prob = DAEProblem(sys, Differential(t).(states(sys)) .=> 0, u0, (0, 0.5))
 @test_nowarn sol = solve(prob, DFBDF())
 
+sys2 = structural_simplify(ll2_model)
+@test length(equations(sys2)) == 3
+u0 = states(sys) .=> 0
+prob = ODEProblem(sys, u0, (0, 10.0))
+@test_nowarn sol = solve(prob, FBDF())
+
 @variables t x1(t) x2(t) x3(t) x4(t)
 D = Differential(t)
 @named sys1_inner = ODESystem([D(x1) ~ x1], t)