I’m attempting to use scipy’s solve_ivp in order to perform some numerical analysis on Loewner equations, specifically on the first hitting time between the solution and a certain function. My question is, what happens if part of solve_ivp’s solution becomes complex? For example, if I am calculating the solution for -2 / (h(t)-sqrt(t)) with some h(0), and I know for certainty through complex analysis that the solution contains complex values after time t, how will that output look on the solution? So far the output of my solution after time t when the solution becomes complex still spits out real values. Are these the real parts of the complex values? Similarly, can you tell if the solution provided by solve_ivp contains complex values if you weren’t sure if it did? I’m using the RK45 method which supports complex values. I also find that when I am using solve_ivp and that solution contains real and complex curves, solve_ivp will stall if it tries solving for complex curves that don’t have real parts.

I tried reading the documentation but I couldn’t find anything specific about what happens when a real solution becomes complex.