Re: [Numpy-discussion] Polynomial implementation in numpy.ma

2011-09-17 Thread Charles R Harris
2011/9/17 Stéfan van der Walt > On Sat, Sep 17, 2011 at 9:40 AM, Charles R Harris > wrote: > > The Vandermonde matrix needs to be used for the fitting so nothing should > be > > changed there. > > I remember now where I heard this: > > http://people.maths.ox.ac.uk/trefethen/mythstalk.pdf > > Spe

Re: [Numpy-discussion] Polynomial implementation in numpy.ma

2011-09-17 Thread Stéfan van der Walt
On Sat, Sep 17, 2011 at 9:40 AM, Charles R Harris wrote: > The Vandermonde matrix needs to be used for the fitting so nothing should be > changed there. I remember now where I heard this: http://people.maths.ox.ac.uk/trefethen/mythstalk.pdf Specifically, the statement: "Some well-known algorith

Re: [Numpy-discussion] Polynomial implementation in numpy.ma

2011-09-17 Thread Charles R Harris
On Sat, Sep 17, 2011 at 10:40 AM, Charles R Harris < charlesr.har...@gmail.com> wrote: > > > 2011/9/14 Stéfan van der Walt > >> Hi all, >> >> There were some failures in the polynomial tests earlier today, and >> while investigating I saw that numpy.ma implements its own root >> finder. It uses

Re: [Numpy-discussion] Polynomial implementation in numpy.ma

2011-09-17 Thread Charles R Harris
2011/9/14 Stéfan van der Walt > Hi all, > > There were some failures in the polynomial tests earlier today, and > while investigating I saw that numpy.ma implements its own root > finder. It uses inversion of a Van der Monde matrix, which I believe > may suffer from some numerical instability pr

[Numpy-discussion] Polynomial implementation in numpy.ma

2011-09-14 Thread Stéfan van der Walt
Hi all, There were some failures in the polynomial tests earlier today, and while investigating I saw that numpy.ma implements its own root finder. It uses inversion of a Van der Monde matrix, which I believe may suffer from some numerical instability problems. Given that Charles has gone to som