Type stable initialization for eigen_est in DP5

SciML · Mar 19, 2019 · 354ea3c · 354ea3c
1 parent e167b26
commit 354ea3c
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/perform_step/low_order_rk_perform_step.jl b/src/perform_step/low_order_rk_perform_step.jl
@@ -809,7 +809,7 @@ end
     g6 = tmp
     g7 = u
     # Hairer II, page 22
-    ϱu, ϱd = zero(eltype(k7)), zero(eltype(g7))
+    ϱu, ϱd = zero(eltype(k7))^2, zero(eltype(g7))^2
     @inbounds for i in eachindex(k7)
       ϱu += (k7[i] - k6[i])^2
       ϱd += (g7[i] - g6[i])^2