ensure correct type

byuflowlab · Apr 12, 2023 · 980672a · 980672a
1 parent 42a1e55
commit 980672a
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/unsteady.jl b/src/unsteady.jl
@@ -75,7 +75,8 @@ function odesolve(initialize, onestep!, t, xd, xc, p)
     @views y0 = initialize(t[1], xd, xc[:, 1], p)
 
     # allocate array
-    y = zeros(length(y0), length(t))
+    T = promote_type(eltype(xd), eltype(xc))
+    y = zeros(T, length(y0), length(t))
     @views y[:, 1] .= y0
 
     # step through