diff --git a/src/solve.jl b/src/solve.jl index 9d0036aa02..8061d564e6 100644 --- a/src/solve.jl +++ b/src/solve.jl @@ -210,7 +210,7 @@ function DiffEqBase.__init(prob::Union{DiffEqBase.AbstractODEProblem, if prob isa DAEProblem rate_prototype = prob.du0 else - rate_prototype = DiffEqBase._rate_prototype(u, t, one(t)) + rate_prototype = zero.(DiffEqBase._rate_prototype(u, t, one(t))) end rateType = typeof(rate_prototype) ## Can be different if unitful