The companion matrix is already upper Hessenberg, so it should be possible to use a specialized routine to find its eigenvalues, rather than the general non-symmetric eigenvalue routine. In particular,

lib/polynomial.py#L227 should call (d,z)hseqr rather than eigvals, which calls (d,z)geev.