ENH: interpolate: add a class for t,c,k B-spline representation

This uses FITPACK for evaluation, and documents the B-spline
representation that is used.
ENH: interpolate: add a Fortran routine for evaluating many B-splines at once
TMP