Make cache vector of same type as u0

SciML · Nov 28, 2019 · 072134d · 072134d
1 parent 53b396b
commit 072134d
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