WIP: add b-splines

This PR implements a base class for evaluation of b-splines without fitpack.
Representation is compatible with `splrep` and `splev` though.

The plan is to have a base class with a sensible API usable by both human users and interpolation / fitting to data (the latter is on an immediate plan).