@@ -3,7 +3,7 @@ using Catalyst, SBML
33using Test
44
55cd (@__DIR__ )
6- sbmlfile = joinpath (" data" , " reactionsystem_01.xml" )
6+ sbmlfile = joinpath (" test " , " data" , " reactionsystem_01.xml" )
77const IV = default_t ()
88@parameters k1, c1
99@species s1 (IV), s2 (IV), s1s2 (IV)
@@ -83,10 +83,11 @@ trueeqs = Equation[default_time_deriv()(s1) ~ k1]
8383@test isequal (Catalyst. get_eqs (odesys), trueeqs)
8484@test isequal (Catalyst. get_iv (odesys), IV)
8585@test isequal (Catalyst. get_unknowns (odesys), [s1])
86- @test issetequal (Catalyst . get_ps (odesys), [k1, c1])
86+ @test issetequal (ModelingToolkit . parameters_toplevel (odesys), [k1, c1])
8787u0 = [s1 => 1.0 ]
8888par = [k1 => 1.0 , c1 => 2.0 ]
89- @test isequal (ModelingToolkit. defaults (odesys), ModelingToolkit. _merge (u0, par)) # PL: @Anand: for some reason this does not work with `Catalyst.get_default()`
89+ testdef = merge (Dict (u0), Dict (par))
90+ @test issubset (testdef, defaults (odesys))
9091@named odesys = ODESystem (MODEL1)
9192isequal (nameof (odesys), :odesys )
9293@test structural_simplify (odesys) isa ODESystem
@@ -99,10 +100,11 @@ trueeqs = Equation[default_time_deriv()(s1) ~ -((k1 * s1 * s2) / c1),
99100@test isequal (Catalyst. get_eqs (odesys), trueeqs)
100101@test isequal (Catalyst. get_iv (odesys), IV)
101102@test isequal (Catalyst. get_unknowns (odesys), [s1, s1s2, s2])
102- @test issetequal (Catalyst. get_ps (odesys), [k1, c1])
103+ @test issetequal (Catalyst. parameters_toplevel (odesys), [k1, c1])
103104u0 = [s1 => 2 * 1.0 , s2 => 2 * 1.0 , s1s2 => 2 * 1.0 ]
104105par = [k1 => 1.0 , c1 => 2.0 ]
105- @test isequal (ModelingToolkit. defaults (odesys), ModelingToolkit. _merge (u0, par))
106+ testdef = merge (Dict (u0), Dict (par))
107+ @test issubset (testdef, ModelingToolkit. defaults (odesys))
106108@named odesys = ODESystem (MODEL1)
107109isequal (nameof (odesys), :odesys )
108110
0 commit comments