[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