[SciPy-dev] RFR: fixing scipy.special.orthogonal
Pauli Virtanen
pav at iki.fi
Sun Oct 25 20:56:44 EDT 2009
Hi,
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.
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)
0.37728975480815219
instead of the current
>>> special.chebyt(100)(0.98)
-5.2908451051968135e+20
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