[SciPy-dev] RFR: fixing scipy.special.orthogonal

Charles R Harris charlesr.harris@gmail....
Sun Oct 25 21:21:29 CDT 2009


On Sun, Oct 25, 2009 at 7:37 PM, Pauli Virtanen <pav@iki.fi> wrote:

> su, 2009-10-25 kello 19:23 -0600, Charles R Harris kirjoitti:
> > On Sun, Oct 25, 2009 at 6:56 PM, Pauli Virtanen <pav@iki.fi> wrote:
> >         http://codereview.appspot.com/140053
> >
> >         The orthogonal polynomial routines in scipy.special.orthogonal
> >         are
> >         cumbersome to use if you only want *values* of the polynomials
> >         at a few
> >         points.
> >
> > I'm not familiar with these routines. In what way is the evaluation at
> > a few points cumbersome?
>
> They IIRC compute polynomial roots via eigenvalues of the companion
> matrix, so chebyt(1000)(x) takes ages, and it's not blazingly fast for
> smaller n either.
>
>
Ah, that's right. Pretty much anything would be an improvement there.

Chuck
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/scipy-dev/attachments/20091025/0740564b/attachment.html 


More information about the Scipy-dev mailing list