Merge pull request #437 from isaacsas/make_switchvals_constants

create variables for the various auto-alg thresholds

src/JumpProcesses.jl

@@ -30,6 +30,12 @@ import Base.Threads
 
 const DEFAULT_RNG = Random.default_rng()
 
+# thresholds for auto-alg below which the listed alg is used
+# see select_aggregator for details
+const USE_DIRECT_THRESHOLD = 20
+const USE_RSSA_THRESHOLD = 100
+const USE_SORTINGDIRECT_THRESHOLD = 200
+
 include("jumps.jl")
 include("massaction_rates.jl")
 include("aggregators/aggregators.jl")

src/aggregators/aggregators.jl

@@ -208,18 +208,18 @@ function select_aggregator(jumps::JumpSet; vartojumps_map = nothing,
     end
 
     # if the number of jumps is small, return the Direct
-    num_jumps(jumps) < 20 && return Direct
+    num_jumps(jumps) < USE_DIRECT_THRESHOLD && return Direct
 
     # if there are only massaction jumps, we can build the species-jumps dependency graphs
     can_build_dgs = num_crjs(jumps) == 0 && num_vrjs(jumps) == 0
     have_species_to_jumps_dgs = !isnothing(vartojumps_map) && !isnothing(jumptovars_map)
 
     # if we have the species-jumps dgs or can build them, use a Rejection-based methods
     if can_build_dgs || have_species_to_jumps_dgs
-        (num_jumps(jumps) < 100) && return RSSA
+        (num_jumps(jumps) < USE_RSSA_THRESHOLD) && return RSSA
         return RSSACR
     elseif !isnothing(dep_graph)   # if only have a normal dg
-        (num_jumps(jumps) < 200) && return SortingDirect
+        (num_jumps(jumps) < USE_SORTINGDIRECT_THRESHOLD) && return SortingDirect
         return DirectCR
     else
         return Direct