Understanding jump model preparation cost

I am using Catalyst.jl to make JumpProblems of around 160k species and a fairly sparse interaction structure that includes what JumpProcesses.jl calls “constant-rate jump” (CRJ) reactions. This works, but JumpProblem construction is quite slow (and so is integrator initialization). I would like to understand why that is and what I could do about it. As an (admittedly contrived) example, consider: using Catalyst using Graphs using JumpProcesses t0 = @elapsed begin g = static_scale_free(n, 2n,...