- Move it to its own compilation unit.
- Restore dynamic allocation init function.
- Precompute bounds checking.
- Precompute invalid value checking.
- Apply the same speedups to the 1D cubic interpolator.
- Expand unit tests.
Original: 12bbc4674e410f7cf2eeeb516bfe6eeafb1835cf