scipy/scipy — pull request #221 by pv (2012-05-19 23:25:11+00:00)
Title: ENH: interpolate: allow polyint evaluation on arbitrary shape arrays

Directly from ticket http://projects.scipy.org/scipy/ticket/1021

I didn't read this through or really test it, but the tests seem to pass, and it's from the original author of the module.