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

Pauli Virtanen pav@iki...
Sun Oct 25 19:56:44 CDT 2009



The orthogonal polynomial routines in scipy.special.orthogonal are
cumbersome to use if you only want *values* of the polynomials at a few

This patch adds separate vectorized eval_* routines for them.

What is more controversial, it also hooks these routines in orthopoly1d,
so that

	>>> special.chebyt(100)(0.98)

instead of the current

	>>> special.chebyt(100)(0.98)

Any arithmetic on the returned polynomial objects, however, discards the
stabler evaluation routines. (Maybe when the scipy/scikits.polynomial
stuff gets ready, we can address also this.)

Comments? Good to go in?

Pauli Virtanen

More information about the Scipy-dev mailing list