[SciPy-dev] Generic polynomials class (was Re: Volunteer for Scipy Project)

Fernando Perez fperez.net@gmail....
Thu Oct 15 17:35:23 CDT 2009

Hi all,

On Wed, Oct 14, 2009 at 8:08 AM, Charles R Harris
<charlesr.harris@gmail.com> wrote:
> One reason for the low level functions is that they are short and can be
> easily converted to cython functions or be library functions written in
> fortran or whatever. The z series could be done that way, but at the cython
> level they could probably be dispensed with and the underlying Cheybyshev
> series used instead, that's one reason they are private functions. What they
> brought to the python case was some vectorization of the algorithms. Another
> way to look at them is as symbolic Fourier transforms.

I don't know if this will genuinely add anything to this (very
interesting discussion), but Ariel, Mike and I thought a lot about
fairly similar problems in the nitime package, and went for an
approach with:

- a purely functional library (even if the functions end up having
long signatures and don't cache things they could),

- very simple data cointainers whose design evolves very conservatively

- analysis objects that expose a nice, high-level interface and manage
the containers and library functions.  These are what users should use
90% of the time, but if we get their design wrong, little is lost
because most of the complexities are in the functional library.

This is all summarized in our preprint for the scipy proceedings
(currently under review):


and the nitime docs can be found here:


while the master branch lives at github:


So far we're quite happy with this approach, and I hope this example
may be of some use to you guys as you sort out the poly design.  BTW,
I'm very excited to see this coming along, many thanks to Anne and
Chuck for so much hard work!



More information about the Scipy-dev mailing list