From 84c93b048ecab7c22291503e68c47b6e3264fbb0 Mon Sep 17 00:00:00 2001
From: Hauke Heibel <hauke.heibel@gmail.com>
Date: Fri, 13 Apr 2012 12:50:05 +0200
Subject: Added spline interpolation with pre-defined knot parameters.

---
 unsupported/test/splines.cpp | 26 +++++++++++++++++++++-----
 1 file changed, 21 insertions(+), 5 deletions(-)

(limited to 'unsupported/test/splines.cpp')

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 );
+    }
   }
 }
 
-- 
cgit v1.2.3