I see. I noticed that even Mathematica happily returns wrong results if your precision is too low, see e.g.: https://reference.wolfram.com/mathematica/ref/NSolve.html, the section "Possible issues".

So I think we should document in `polyroots` that the user is responsible to do a "convergence with regards to extraprec", otherwise the results might be wrong.