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

Pauli Virtanen pav@iki...
Sun Oct 25 20:37:19 CDT 2009


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.

	Pauli





More information about the Scipy-dev mailing list