Merge pull request #475 from isaacsas/initial_ensemble_propensity_test

add test initial variable rate propensities are unique in ensembles
SciML · Feb 11, 2025 · 4d52a24 · 4d52a24
2 parents 5d454dc + 16ac07f
commit 4d52a24
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/test/thread_safety.jl b/test/thread_safety.jl
@@ -24,8 +24,12 @@ let
     ode_prob = ODEProblem(f!, u_0, (0.0, 10))
     rate(u, p, t) = 1.0
     jump!(integrator) = nothing
-    jump_prob = JumpProblem(ode_prob, Direct(), VariableRateJump(rate, jump!); rng)
+    jump_prob = JumpProblem(ode_prob, Direct(), VariableRateJump(rate, jump!))
     prob_func(prob, i, repeat) = deepcopy(prob)
     prob = EnsembleProblem(jump_prob,prob_func = prob_func)
     solve(prob, Tsit5(), EnsembleThreads(), trajectories=10)
+
+    sol = solve(prob, Tsit5(), EnsembleThreads(), trajectories=400)
+    init_props = [sol[i].u[1][2] for i = 1:length(sol)]    
+    @test allunique(init_props)
 end