change higehr order test

Author: TorkelE

test/reactionsystem_structure/higher_order_reactions.jl

@@ -92,28 +92,28 @@ let
     higher_order_network_alt2 = ConstantRateJump.([rate1, rate2, rate3, rate4, rate5, rate6, rate7, rate8], 
                                 [affect1!, affect2!, affect3!, affect4!, affect5!, affect6!, affect7!, affect8!])
 
-    # For the systems created via Catalyst and manually, compares that they yield identical simulations.
-    for n in [5, 50]   
-        # Prepares JumpProblem via Catalyst.       
-        u0_base = rnd_u0_Int64(base_higher_order_network, rng; n)
-        ps_base = rnd_ps(base_higher_order_network, rng; factor = n/10.0)
-        dprob_base = DiscreteProblem(base_higher_order_network, u0_base, (0.0, 100.0), ps_base)
-        jprob_base = JumpProblem(base_higher_order_network, dprob_base, Direct(); rng = StableRNG(1234))
-
-        # Prepares JumpProblem partially declared manually.   
-        dprob_alt1 = DiscreteProblem(higher_order_network_alt1, u0_base, (0.0, 100.0), ps_base)
-        jprob_alt1 = JumpProblem(higher_order_network_alt1, dprob_alt1, Direct(); rng = StableRNG(1234))
-
-        # Prepares JumpProblem via manually declared system.    
-        u0_alt2 = map_to_vec(u0_base, [:X1, :X2, :X3, :X4, :X5, :X6, :X7, :X8, :X9, :X10])
-        ps_alt2 = map_to_vec(ps_base, [:p, :r1, :r2, :K, :r3, :r4, :r5, :r6, :d])
-        dprob_alt2 = DiscreteProblem(u0_alt2, (0.0, 100.0), ps_alt2)
-        jprob_alt2 = JumpProblem(dprob_alt2, Direct(), higher_order_network_alt2...; rng = StableRNG(1234))
-
-        # Compares the simulations
-        sol_base = solve(jprob_base, SSAStepper(); seed, saveat = 1.0)
-        sol_alt1 = solve(jprob_alt1, SSAStepper(); seed, saveat = 1.0)
-        sol_alt2 = solve(jprob_alt2, SSAStepper(); seed, saveat = 1.0)
-        @test_broken sol_base == sol_alt1 == sol_alt2
-    end
+    # Prepares JumpProblem via Catalyst.       
+    u0_base = rnd_u0_Int64(base_higher_order_network, rng; n)
+    ps_base = rnd_ps(base_higher_order_network, rng; factor = n/10.0)
+    dprob_base = DiscreteProblem(base_higher_order_network, u0_base, (0.0, 100.0), ps_base)
+    jprob_base = JumpProblem(base_higher_order_network, dprob_base, Direct(); rng = StableRNG(1234))
+
+    # Prepares JumpProblem partially declared manually.   
+    dprob_alt1 = DiscreteProblem(higher_order_network_alt1, u0_base, (0.0, 100.0), ps_base)
+    jprob_alt1 = JumpProblem(higher_order_network_alt1, dprob_alt1, Direct(); rng = StableRNG(1234))
+
+    # Prepares JumpProblem via manually declared system.    
+    u0_alt2 = map_to_vec(u0_base, [:X1, :X2, :X3, :X4, :X5, :X6, :X7, :X8, :X9, :X10])
+    ps_alt2 = map_to_vec(ps_base, [:p, :r1, :r2, :K, :r3, :r4, :r5, :r6, :d])
+    dprob_alt2 = DiscreteProblem(u0_alt2, (0.0, 100.0), ps_alt2)
+    jprob_alt2 = JumpProblem(dprob_alt2, Direct(), higher_order_network_alt2...; rng = StableRNG(1234))
+
+    # Simualtes the models.
+    sol_base = solve(jprob_base, SSAStepper(); seed, saveat = 1.0)
+    sol_alt1 = solve(jprob_alt1, SSAStepper(); seed, saveat = 1.0)
+    sol_alt2 = solve(jprob_alt2, SSAStepper(); seed, saveat = 1.0)
+    
+    # Checks that species means in the simulations are similar
+    @test mean(sol_base[:X10]) ≈ mean(sol_alt1[:X10]) atol = 1e-1 rtol = 1e-1
+    @test mean(sol_alt1[:X10]) ≈ mean(sol_alt2[10,:]) atol = 1e-1 rtol = 1e-1
 end