Here are some spline routines I've used. They are
intended to be fast enough for large arrays.
http://www.astro.umd.edu/~jph/splines.ijs
(toward the bottom of http://www.astro.umd.edu/~jph/J_page.html)
Patrick Harrington