ENH: Allow `Splinterpolator` instances to be created from an existing set of spline coefficients, and extend extrapolation options

• The Splinterpolator class can now accept a pointer to an array of pre-calculated coefficients - in this case, a copy of the pointer is used (the data is not copied), and pre-filtering (coefficient calculation) is skipped
• New SoftZeros extrapolation method - out of bounds coefficients will be set to zero, but out-of-bounds coordinates will still be constructed via interpolation. This differs from the existing Zeros extrapolation method, where out-of-bounds coordinates are forced to zero.
• Adjustment to the coefficient initialisation method used when Constant extrapolation is requested; previously, a mirroring boundary condition was used.
• New Splinterpolator::Copy function, which allows the spline coefficients to be copied/extracted.

test/test_splinterpolator.cc

@@ -104,10 +104,10 @@ public:
   }
 
   // compare two float vectors
-  void compare(std::vector<float> a, std::vector<float> b, float tol=1e-6) {
+  void compare(std::vector<float> a, std::vector<float> b, float tol=1e-5) {
     BOOST_CHECK_EQUAL(a.size(), b.size());
     for (auto i = 0; i < a.size(); i++) {
-      BOOST_CHECK_CLOSE(a[i], b[i], tol);
+      BOOST_CHECK_SMALL(a[i] - b[i], tol);
     }
   }