awkward that you overwrite existing Polynomials.roots(); prevents comparison

[ahbarnett]

Thanks for this exciting package, which I am benchmarking against AMRVW.jl, Polynomials.jl, PolynomialRoots.jl, and my own companion matrix eigval implementation.

However, because your package chooses to redefine the roots() method from the Polynomials module, one cannot actually run benchmarks, unless they are done as in your README.md where you can only test Polynomials.roots() before you load in FPR. I want to be able to load both modules, then call one or the other roots() function but cannot because either using or import on your pkg appears to redefine Polynomials.roots(). I think this is not good behavior; FPR should keep its own namespace and not change another pkg.

The offending line is

FastPolynomialRoots.jl/src/FastPolynomialRoots.jl

Line 5 in 17f0ab4

Polynomials.roots(p::Union{Polynomial{Float64},Polynomial{Complex{Float64}}}) = rootsFastPolynomialRoots(coeffs(p))

Is there something I'm missing here (relatively new to Julia)? Thanks, Alex

[andreasnoack]

The original idea here was to speed up the existing implementation. However, this kind of piracy can indeed be confusing and cause issues so I can probably be convinced to change the current behavior and instead introduce a new function.