Update tests

Author: YingboMa

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, [0, 1, 0, 0, 2, 0, 3, 0, 4])
+@test StructuralTransformations.matching(sss, varassoc .== 0) == map(x -> x == 0 ? StructuralTransformations.UNASSIGNED : x, [1, 2, 3, 4, 0, 0, 0, 0, 0])
 
 sys, assign, eqassoc = StructuralTransformations.pantelides!(pendulum)
 sss = structure(sys)
@@ -48,8 +48,8 @@ scc = StructuralTransformations.find_scc(graph, assign)
                            [6],
                           ]
 
-@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]
+@test graph.fadjlist == [[1, 7], [2, 8], [3, 5, 9], [4, 6, 9], [5, 6], [1, 2, 5, 6], [1, 3, 7, 10], [2, 4, 8, 11], [1, 2, 5, 6, 10, 11]]
+@test varassoc == [10, 11, 0, 0, 1, 2, 3, 4, 0, 0, 0]
 #1: D(x) ~ w
 #2: D(y) ~ z
 #3: D(w) ~ T*x

test/structural_transformation/tearing.jl

@@ -145,7 +145,7 @@ eqs = [
 daesys = ODESystem(eqs, t)
 newdaesys = tearing(daesys)
 @test equations(newdaesys) == [D(x) ~ z; 0 ~ x + sin(z) - p*t]
-@test isequal(states(newdaesys), [z, x])
+@test isequal(states(newdaesys), [x, z])
 prob = ODAEProblem(newdaesys, [x=>1.0], (0, 1.0), [p=>0.2])
 du = [0.0]; u = [1.0]; pr = 0.2; tt = 0.1
 @test (@ballocated $(prob.f)($du, $u, $pr, $tt)) == 0

test/structural_transformation/utils.jl

@@ -20,12 +20,12 @@ sys = initialize_system_structure(pendulum)
 StructuralTransformations.find_solvables!(sys)
 sss = structure(sys)
 @unpack graph, solvable_graph, fullvars, varassoc = sss
-@test isequal(fullvars, [w, D(x), x, z, D(y), y, T, D(w), D(z)])
-@test graph.fadjlist == sort.([[1, 2], [4, 5], [3, 7, 8], [6, 7, 9], [3, 6]])
+@test isequal(fullvars, [D(x), D(y), D(w), D(z), x, y, w, z, T])
+@test graph.fadjlist == [[1, 7], [2, 8], [3, 5, 9], [4, 6, 9], [5, 6]]
 @test graph.badjlist == 9 == length(fullvars)
 @test ne(graph) == nnz(incidence_matrix(graph)) == 12
 @test nv(solvable_graph) == 9 + 5
-@test varassoc == [8, 0, 2, 9, 0, 5, 0, 0, 0]
+@test varassoc == [0, 0, 0, 0, 1, 2, 3, 4, 0]
 
 se = collect(StructuralTransformations.𝑠edges(graph))
 @test se == mapreduce(vcat, enumerate(graph.fadjlist)) do (s, d)