@@ -50,7 +50,7 @@ rs = ReactionSystem(MODEL1)
5050 Catalyst. Reaction[Catalyst. Reaction (k1, nothing , [s1], nothing , [1.0 ])])
5151@test isequal (Catalyst. get_iv (rs), IV)
5252@test isequal (Catalyst. get_species (rs), [s1])
53- @test isequal (Catalyst. get_ps (rs), [k1, c1])
53+ @test issetequal (Catalyst. get_ps (rs), [k1, c1])
5454@named rs = ReactionSystem (MODEL1)
5555isequal (nameof (rs), :rs )
5656
@@ -60,7 +60,7 @@ rs = ReactionSystem(readSBML(sbmlfile))
6060 [1.0 ])])
6161@test isequal (Catalyst. get_iv (rs), IV)
6262@test isequal (Catalyst. get_species (rs), [s1, s1s2, s2])
63- @test isequal (Catalyst. get_ps (rs), [k1, c1])
63+ @test issetequal (Catalyst. get_ps (rs), [k1, c1])
6464@named rs = ReactionSystem (MODEL1)
6565isequal (nameof (rs), :rs )
6666
@@ -72,7 +72,7 @@ rs = ReactionSystem(MODEL2) # Contains reversible reaction
7272 nothing , [1 ])])
7373@test isequal (Catalyst. get_iv (rs), IV)
7474@test isequal (Catalyst. get_species (rs), [s1])
75- @test isequal (Catalyst. get_ps (rs), [k1, c1])
75+ @test issetequal (Catalyst. get_ps (rs), [k1, c1])
7676
7777@test convert (ModelingToolkit. ODESystem, rs) isa ODESystem
7878@test structural_simplify (convert (ModelingToolkit. ODESystem, rs)) isa ODESystem
@@ -83,7 +83,7 @@ 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 isequal (Catalyst. get_ps (odesys), [k1, c1])
86+ @test issetequal (Catalyst. get_ps (odesys), [k1, c1])
8787u0 = [s1 => 1.0 ]
8888par = [k1 => 1.0 , c1 => 2.0 ]
8989@test isequal (ModelingToolkit. defaults (odesys), ModelingToolkit. _merge (u0, par)) # PL: @Anand: for some reason this does not work with `Catalyst.get_default()`
@@ -99,7 +99,7 @@ trueeqs = Equation[default_time_deriv()(s1) ~ -((k1 * s1 * s2) / c1),
9999@test isequal (Catalyst. get_eqs (odesys), trueeqs)
100100@test isequal (Catalyst. get_iv (odesys), IV)
101101@test isequal (Catalyst. get_unknowns (odesys), [s1, s1s2, s2])
102- @test isequal (Catalyst. get_ps (odesys), [k1, c1])
102+ @test issetequal (Catalyst. get_ps (odesys), [k1, c1])
103103u0 = [s1 => 2 * 1.0 , s2 => 2 * 1.0 , s1s2 => 2 * 1.0 ]
104104par = [k1 => 1.0 , c1 => 2.0 ]
105105@test isequal (ModelingToolkit. defaults (odesys), ModelingToolkit. _merge (u0, par))
0 commit comments