Is your feature request related to a problem? Please describe.

I want to find all the roots of a polynomial function

I find very convenient your wrapper around HomotopyContinuation, but at the moment only returns real solutions

even when Homotopy Continuation has found the complex solutions already.

Right now you even have a test that checks for real solutions only

Describe the solution you’d like

I would like to pass a kwarg like allow_complex=true and get complex solutions if wanted. I think it will be relatively simple since the complex roots are calculated by Homotopy Continuation either way.

Describe alternatives you’ve considered

Right now I know you have the Muller method that finds complex roots:

but I want to find all the roots. That is why at the moment I am using your homotopy_continuation_preprocessing function

and then I use HomotopyContinuation.solve

Additional context

NonlinearSolve v4.10.0
NonlinearSolveHomotopyContinuation v0.1.5