Merge pull request #88 from baggepinnen/patch-1

Make cache vector of same type as u0
SciML · Nov 28, 2019 · 8dd5183 · 8dd5183
2 parents 53b396b + 072134d
commit 8dd5183
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/simple_regular_solve.jl b/src/simple_regular_solve.jl
@@ -16,7 +16,7 @@ function DiffEqBase.solve(jump_prob::JumpProblem,alg::SimpleTauLeaping;
   c = rj.c
   dc = zero(rj.dc)
   fill!(dc,0)
-  rate_cache = zeros(size(dc,2))
+  rate_cache = zeros(eltype(prob.u0), size(dc,2))
 
   u0 = copy(prob.u0)
   tspan = prob.tspan