simple single system reduction test

Author: shashi

test/reduction.jl

@@ -1,6 +1,6 @@
 using ModelingToolkit, OrdinaryDiffEq, Test
 
-@parameters t σ ρ β
+@parameters t σ ρ β F(t)
 @variables x(t) y(t) z(t) a(t)
 @derivatives D'~t
 
@@ -12,14 +12,16 @@ eqs = [D(x) ~ σ*(y-x),
 lorenz1 = ODESystem(eqs,t,[x,y,z,a],[σ,ρ,β],name=:lorenz1)
 
 lorenz1_aliased = alias_elimination(lorenz1)
-length(equations(lorenz1_reduced)) = 3
-length(states(lorenz1_reduced)) = 3
+@test length(equations(lorenz1_aliased)) == 3
+@test length(states(lorenz1_aliased)) == 3
 
 eqs = [D(x) ~ σ*(y-x),
        D(y) ~ x*(ρ-z)-y,
        D(z) ~ x*y - β*z]
-test_lorenz1_aliased = ODESystem(eqs,t,[x,y,z],[σ,ρ,β],observed=[a ~ x],name=:lorenz1)
 
+@test lorenz1_aliased == ODESystem(eqs,t,[x,y,z],[σ,ρ,β],observed=[a ~ x],name=:lorenz1)
+
+#=
 # Multi-System Reduction
 
 eqs1 = [D(x) ~ σ*(y-x) + F,
@@ -83,3 +85,4 @@ observed(reduced_flattened_system) == [
         lorenz2.F ~ lorenz1.u
         a ~ -lorenz1.x + lorenz2.y
 ]
+=#