scipy/scipy — pull request #339 by ewmoore (2012-10-17 01:15:51+00:00)
Title: Make the orthogonal polynomial eval_* functions ufuncs.

As I proposed [here](http://mail.scipy.org/pipermail/scipy-dev/2012-September/017990.html) and branching from Pauli's work in #333 and #334, I've moved the functions for evaluating orthogonal polynomials to cython and made them into ufuncs. At present there are no changes to how things are actually evaluated. Since this is based on #333 and #334, only ed50214 and on are of interest here.

At present none of the polynomials can be evaluated for a complex number. This probably needs to be fixed since it removes an existing feature. I'll add this soon.

Ticket [1435](http://projects.scipy.org/scipy/ticket/1435) still needs to be addressed, but this is easy, too (Need to change orthopoly1d.__call__).

Some other polynomials can also be evaluated using recurrences as long as care is taken, (c.f. http://nalag.cs.kuleuven.be/papers/paul/report74/report74.ps.gz). This ought to be faster for small n at least, so I'd like to add it where it makes sense.

I'm submitting this as something of a work in progress because wanted to get something up before I ended up just letting it rot on my machine.