I had inadvertently used the complex-valued Catmull-Rom interpolation
function for real-valued interpolation. The real-valued version is about
twice as fast and includes some special cases for infinite or NaN data
points.
I am also reinstating the Mathematica unit tests.
Original: 6ea910f0b7c187d4672955acf09f66ab6b61ceb6