package: pari
version: 2.8.0, [2, 8, 0, 18831, "b283532"]
According to the documentation, polrootsreal includes multiple roots according to their multiplicity. I don't have any strong feelings on whether they should be or not, but it looks like they are suppressed when they are 0 and retained otherwise.
> polrootsreal(x^3 - x^2)
%1 = [0.E-38, 1.0000000000000000000000000000000000000]~
> polrootsreal(x^3 - x^2 - x + 1)
%2 = [-1.0000000000000000000000000000000000000, 1.0000000000000000000000000000000000000, 1.0000000000000000000000000000000000000]~
This should, hopefully, be my last report on polrootsreal -- disregarding this case I can't find any problems in the latest git version.
Charles Greathouse
Analyst/Programmer
Case Western Reserve University