Merge pull request #2381 from SciML/os/complex-eigen_est

fix Complex `eigen est`

Author: ChrisRackauckas

lib/OrdinaryDiffEqLowOrderRK/src/low_order_rk_perform_step.jl

@@ -762,7 +762,7 @@ end
         g7 = u
         # Hairer II, page 22 modified to use the Inf norm
         integrator.eigen_est = integrator.opts.internalnorm(
-            maximum(abs.(k7 .- k6) ./ (g7 .- g6)), t)
+            maximum(abs.((k7 .- k6) ./ (g7 .- g6))), t)
     end
     if integrator.opts.adaptive
         utilde = dt *

lib/OrdinaryDiffEqTsit5/src/tsit_perform_step.jl

@@ -141,7 +141,7 @@ end
         g7 = u
         # Hairer II, page 22 modified to use the Inf norm
         integrator.eigen_est = integrator.opts.internalnorm(
-            maximum(abs.(k7 .- k6) ./ (g7 .- g6)), t)
+            maximum(abs.((k7 .- k6) ./ (g7 .- g6))), t)
     end
     if integrator.opts.adaptive
         utilde = dt *

test/interface/default_solver_tests.jl

@@ -110,3 +110,11 @@ using StaticArrays
 cb = ContinuousCallback((u, t, integrator) -> t - 1, (integrator) -> nothing)
 SA_ode_problem = ODEProblem((u, p, t) -> zero(u), SA[0], 2)
 @test solve(SA_ode_problem; callback = cb).retcode == ReturnCode.Success
+
+# test Complex numbers
+H(s) = (1-s) * complex([0 1; 1 0]) + s * complex([1 0; 0 -1])
+schrod_eq(state, time, s) = -im * time * H(s) * state
+
+prob_complex = ODEProblem(schrod_eq, complex([1, -1]/sqrt(2)), (0,1), 100)
+complex_sol = solve(prob_complex)
+@test complex_sol.retcode == ReturnCode.Success