<br><br><div class="gmail_quote">2011/9/14 Stéfan van der Walt <span dir="ltr"><<a href="mailto:stefan@sun.ac.za">stefan@sun.ac.za</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi all,<br>
<br>
There were some failures in the polynomial tests earlier today, and<br>
while investigating I saw that <a href="http://numpy.ma" target="_blank">numpy.ma</a> implements its own root<br>
finder.  It uses inversion of a Van der Monde matrix, which I believe<br>
may suffer from some numerical instability problems.  Given that<br>
Charles has gone to some length to implement good polynomial root<br>
finders, I think it would be best to employ those instead, and simply<br>
pre-filter the data that comes from the masked array module, if<br>
possible.<br>
<br></blockquote><div><br>The test failure arises because the test compares against np.polyfit where Travis changed the column scaling, which in turn changed the singular values, and it is that comparison that fails.<br>
<br>The Vandermonde matrix needs to be used for the fitting so nothing should be changed there. To fix the test failure, either the test could be modified or the ma.polyfit routine can be changed. Probably the latter is the best option, the new scaling code can be brought over and maybe the weight and covariance options added.Travis' changes needs to be audited in any case to make sure all the corner cases are covered.<br>
<br>Chuck<br></div></div>