numpy/numpy — pull request #191 by charris (2012-01-31 03:20:17+00:00)
Title: ENH: Improve the computation of polynomials from roots.

The original method was overly sensitive to roundoff. Of the two
approaches considered, gauss integration or binary subdivision of
the roots, the latter is more compatible with using other number
representations such as mpmath and was used on that account.
No method is going to be suitable for large numbers of arbitrary zeros
but the current method is a significant improvement over the previous one.