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

Pauli Virtanen pav at iki.fi
Sun Oct 25 21:37:19 EDT 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 at 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