Added spline interpolation with pre-defined knot parameters.

1 files changed, 21 insertions, 5 deletions

diff --git a/unsupported/test/splines.cpp b/unsupported/test/splines.cpp
index 02e170c44..fe98bf183 100644
--- a/unsupported/test/splines.cpp
+++ b/unsupported/test/splines.cpp
@@ -216,16 +216,32 @@ void check_global_interpolation2d()
   typedef Spline2d::ControlPointVectorType ControlPointVectorType;
 
   ControlPointVectorType points = ControlPointVectorType::Random(2,100);
-  const Spline2d spline = SplineFitting<Spline2d>::Interpolate(points,3);
 
   KnotVectorType chord_lengths; // knot parameters
   Eigen::ChordLengths(points, chord_lengths);
 
-  for (Eigen::DenseIndex i=0; i<points.cols(); ++i)
+  // interpolation without knot parameters
   {
-    PointType pt = spline( chord_lengths(i) );
-    PointType ref = points.col(i);
-    VERIFY( (pt - ref).matrix().norm() < 1e-14 );
+    const Spline2d spline = SplineFitting<Spline2d>::Interpolate(points,3);  
+
+    for (Eigen::DenseIndex i=0; i<points.cols(); ++i)
+    {
+      PointType pt = spline( chord_lengths(i) );
+      PointType ref = points.col(i);
+      VERIFY( (pt - ref).matrix().norm() < 1e-14 );
+    }
+  }
+
+  // interpolation with given knot parameters
+  {
+    const Spline2d spline = SplineFitting<Spline2d>::Interpolate(points,3,chord_lengths);  
+
+    for (Eigen::DenseIndex i=0; i<points.cols(); ++i)
+    {
+      PointType pt = spline( chord_lengths(i) );
+      PointType ref = points.col(i);
+      VERIFY( (pt - ref).matrix().norm() < 1e-14 );
+    }
   }
 }