A fix that actually works

SciML · Oct 31, 2018 · 90871fa · ChrisRackauckas · Oct 31, 2018 · 90871fa
1 parent a7feb13
commit 90871fa
Showing 1 changed file with 2 additions and 5 deletions.
diff --git a/src/initdt.jl b/src/initdt.jl
@@ -18,16 +18,13 @@
     f₀ = integrator.fsallast
     f(f₀,u0,p,t)
   else
-    # TODO: use more caches
-    Tf = typeof(oneunit(eltype(u0))/oneunit(t))
-    f₀ = similar(u0,Tf)
+    f₀ = OrdinaryDiffEq.get_du(integrator)
     f(f₀,u0,p,t)
   end
 
   # TODO: use more caches
   #tmp = cache[2]
-  Ttmp = typeof(oneunit(eltype(u0))/oneunit(eltype(sk)))
-  tmp = similar(u0,Ttmp)
+  tmp = similar(u0)
   @. tmp = u0/sk
 
   d₀ = internalnorm(tmp)