add rng in JumpProblems that are simualted

Author: TorkelE

test/simulation_and_solving/hybrid_models.jl

@@ -101,7 +101,7 @@ let
     rn = complete(extend(rn2, rn))
     tspan = (0.0, 200.0)
     jinputs = JumpInputs(rn, u0map, tspan, pmap; save_positions = (false, false))
-    jprob = JumpProblem(jinputs; rng, save_positions = (false, false))
+    jprob = JumpProblem(jinputs; rng, save_positions = (false, false); rng)
     Xsamp = 0.0
     Nsims = 4000
     for n in 1:Nsims
@@ -146,8 +146,8 @@ let
     ps_hybrid = [ps_ode; ps_jump]
     tspan = (0.0, 10000.0)
     ode_prob = ODEProblem(rn_ode, u0_ode, tspan, ps_ode)
-    jump_prob = JumpProblem(JumpInputs(rn_jump, u0_jump, tspan, ps_jump; remake_warn = false); save_positions = (false,false))
-    hybrid_prob = JumpProblem(JumpInputs(rn_hybrid, u0_hybrid, tspan, ps_hybrid; remake_warn = false); save_positions = (false,false))
+    jump_prob = JumpProblem(JumpInputs(rn_jump, u0_jump, tspan, ps_jump; remake_warn = false); save_positions = (false,false); rng)
+    hybrid_prob = JumpProblem(JumpInputs(rn_hybrid, u0_hybrid, tspan, ps_hybrid; remake_warn = false); save_positions = (false,false); rng)
 
     # Performs simulations. Checks that ODE parts are identical. Check that jump parts have similar statistics.
     ode_sol = solve(ode_prob, Tsit5(); saveat = 1.0, abstol = 1e-10, reltol = 1e-10)