add test

Author: oscardssmith

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