fix initialization bug

HISTORY.md

@@ -0,0 +1,9 @@
+# Breaking updates and feature summaries across releases
+
+## JumpProcesses unreleased (master branch)
+- Support for "bounded" `VariableRateJump`s that can be used with the `Coevolve`
+ aggregator for faster simulation of jump processes with time-dependent rates.
+ In particular, if all `VariableRateJump`s in a pure-jump system are bounded one
+ can use `Coevolve` with `SSAStepper` for better performance. See the
+ documentation, particularly the first and second tutorials, for details on
+ defining and using bounded `VariableRateJump`s.

src/aggregators/coevolve.jl

@@ -65,8 +65,8 @@ function aggregate(aggregator::Coevolve, u, p, t, end_time, constant_jumps,
     affects! = Vector{AffectWrapper}(undef, nrjs)
     rates = Vector{RateWrapper}(undef, nvrjs)
     lrates = similar(rates)
-    urates = similar(rates)
     rateintervals = similar(rates)
+    urates = Vector{RateWrapper}(undef, nrjs)
     haslratevec = zeros(Bool, nvrjs)
 
     idx = 1