Update tests

Author: YingboMa

src/bipartite_graph.jl

@@ -102,11 +102,7 @@ end
 Base.eltype(::Type{<:BipartiteGraph{I}}) where I = I
 function Base.empty!(g::BipartiteGraph)
     foreach(empty!, g.fadjlist)
-    if g.badjlist isa AbstractVector
-        foreach(empty!, g.badjlist)
-    else
-        g.badjlist = 0
-    end
+    g.badjlist isa AbstractVector && foreach(empty!, g.badjlist)
     g.ne = 0
     if g.metadata !== nothing
         foreach(empty!, g.metadata)

test/structural_transformation/index_reduction.jl

@@ -34,7 +34,7 @@ pendulum = ODESystem(eqs, t, [x, y, w, z, T], [L, g], name=:pendulum)
 pendulum = initialize_system_structure(pendulum)
 sss = structure(pendulum)
 @unpack graph, fullvars, varassoc = sss
-@test StructuralTransformations.matching(sss, varassoc .== 0) == map(x -> x == 0 ? StructuralTransformations.UNASSIGNED : x, [1, 0, 2, 0, 3, 4, 0, 0, 0])
+@test StructuralTransformations.matching(sss, varassoc .== 0) == map(x -> x == 0 ? StructuralTransformations.UNASSIGNED : x, [0, 1, 0, 0, 2, 0, 3, 0, 4])
 
 sys, assign, eqassoc = StructuralTransformations.pantelides!(pendulum)
 sss = structure(sys)
@@ -48,18 +48,8 @@ scc = StructuralTransformations.find_scc(graph, assign)
                            [6],
                           ]
 
-@test graph.fadjlist == sort.([
- [1, 2],
- [3, 4],
- [5, 6, 7],
- [6, 8, 9],
- [7, 9],
- [1, 3, 7, 9],
- [1, 2, 5, 10],
- [3, 4, 8, 11],
- [1, 3, 7, 9, 10, 11],
-])
-@test varassoc == [10, 5, 11, 8, 0, 0, 1, 0, 3, 0, 0]
+@test graph.fadjlist == sort.([[1, 2], [4, 5], [3, 7, 8], [6, 7, 9], [3, 6], [2, 3, 5, 6], [1, 2, 8, 10], [4, 5, 9, 11], [2, 3, 5, 6, 10, 11]])
+@test varassoc == [8, 10, 2, 9, 11, 5, 0, 0, 0, 0, 0]
 #1: D(x) ~ w
 #2: D(y) ~ z
 #3: D(w) ~ T*x
@@ -100,22 +90,24 @@ using LinearAlgebra
 prob = ODEProblem(ODEFunction(first_order_idx1_pendulum),
         #  [x, y, w, z, xˍt, yˍt, T]
            [1, 0, 0, 0,   0,   0, 0.0],# 0, 0, 0, 0],
-           (0, 100.0),
+           (0, 10.0),
            [1, 9.8],
            mass_matrix=calculate_massmatrix(first_order_idx1_pendulum))
 sol = solve(prob, Rodas5());
 #plot(sol, vars=(1, 2))
 
 new_sys = dae_index_lowering(ModelingToolkit.ode_order_lowering(pendulum2))
 
-prob_auto = ODEProblem(ODEFunction(new_sys),
-        #  [xˍt, yˍt, x, y,   T]
-           [0,     0, 1, 0, 0.0],# 0, 0, 0, 0],
-           (0, 100.0),
-           [1, 9.8],
-           mass_matrix=calculate_massmatrix(new_sys))
+prob_auto = ODEProblem(new_sys,
+                       [D(x)=>0,
+                        D(y)=>0,
+                        x=>1,
+                        y=>0,
+                        T=>0.0],
+                       (0, 100.0),
+                       [1, 9.8])
 sol = solve(prob_auto, Rodas5());
-#plot(sol, vars=(3, 4))
+#plot(sol, vars=(x, y))
 
 # Define some variables
 @parameters t L g
@@ -146,7 +138,7 @@ p = [
     g => 9.8
 ]
 
-prob_auto = ODEProblem(new_sys,u0,(0.0,100.0),p)
+prob_auto = ODEProblem(new_sys,u0,(0.0,10.0),p)
 sol = solve(prob_auto, Rodas5());
 #plot(sol, vars=(D(x), y))
 

test/structural_transformation/utils.jl

@@ -24,7 +24,7 @@ sss = structure(sys)
 @test graph.fadjlist == sort.([[1, 2], [4, 5], [3, 7, 8], [6, 7, 9], [3, 6]])
 @test graph.badjlist == 9 == length(fullvars)
 @test ne(graph) == nnz(incidence_matrix(graph)) == 12
-@test nv(solvable_graph) == nv(solvable_graph) == 9 + 5
+@test nv(solvable_graph) == 9 + 5
 @test varassoc == [8, 0, 2, 9, 0, 5, 0, 0, 0]
 
 se = collect(StructuralTransformations.𝑠edges(graph))