Adding capabilities of Roots.jl and HomotopyContinuation.jl with NonlinearSolve.jl

[yewalenikhil65]

Multiple root like capabilities from Roots.jl ?.
Many Roots.jl solvers can be included or Roots.jl package can be included in NonlinearSolve.jl for its use as common interface

I think @ChrisRackauckas has plans later to include HomotopyContinuation.jl

Summarising this here as issue

[ChrisRackauckas]

First things first is interface. It probably needs a separate entry point that makes it explicit about giving multiple solutions. Maybe some meta-algorithm wrapper? Deflation(NewtonRhapson()), where then those algorithms always return u as an array of solutions?