switch to parameter set where CLE won't risk going negative

Author: TorkelE

test/miscellaneous_tests/symbolic_stoichiometry.jl

@@ -218,8 +218,8 @@ let
     end
 
     u0 = [:X => 8, :Y => 8]
-    ps_int = (:p => 1.0, :k => 2.0, :n1 => 3, :n2 => 4, :d => 0.1)
-    ps_dec = [:p => 1.0, :k => 2.0, :n1 => 0.5, :n2 => 2.5, :d => 0.1]
+    ps_int = (:p => 2.0, :k => 0.01, :n1 => 3, :n2 => 4, :d => 0.2)
+    ps_dec = [:p => 2.0, :k => 0.01, :n1 => 0.5, :n2 => 2.5, :d => 0.2]
     tspan = (0.0, 10.0)
 
     # Test ODE simulations with integer coefficients.
@@ -243,11 +243,13 @@ let
     @test solve(sprob_dec, ImplicitEM(); seed) ≈ solve(sprob_dec_ref, ImplicitEM(); seed)
 
     # Tests jump simulations with integer coefficients.
-    dprob_int = DiscreteProblem(rs_int, u0, tspan, ps_int)
-    dprob_int_ref = DiscreteProblem(rs_ref_int, u0, tspan, ps_int)
+    dprob_int = DiscreteProblem(rs_int, u0, (0.0, 100000.0), ps_int)
+    dprob_int_ref = DiscreteProblem(rs_ref_int, u0, (0.0, 100000.0), ps_int)
     jprob_int = JumpProblem(rs_int, dprob_int, Direct(); rng)
     jprob_int_ref = JumpProblem(rs_ref_int, dprob_int_ref, Direct(); rng)
-    @test solve(jprob_int, SSAStepper(); seed) ≈ solve(jprob_int_ref, SSAStepper(); seed)
+    sol_int = solve(jprob_int, SSAStepper(); seed)
+    sol_int_ref = solve(jprob_int_ref, SSAStepper(); seed)
+    @test mean(sol_int[:Y]) ≈ mean(sol_int_ref[:Y]) atol = 1e-2 rtol = 1e-2
 end
 
 # Check that jump simulations (implemented with and without symbolic stoichiometries) yield simulations