Return vector of kept waypoints from method simplify().

Author: gnestlinger

PolygonPath.m

@@ -879,22 +879,25 @@
 
         end%fcn
 
-        function obj = simplify(obj)
+        function [obj,keep] = simplify(obj)
         %SIMPLIFY   Simplify path.
         %   OBJ = SIMPLIFY(OBJ) removes intermediate points from line
         %   segments of the path OBJ.
-        
+        %
+        %   [OBJ,KEEP] = SIMPLIFY(OBJ) returns a logical vector KEEP
+        %   indicating which waypoints are kept.
+            
             h = atan2(diff(obj.y), diff(obj.x));
             dh = [true; diff(h); true];
             keep = (dh ~= 0);
-
+            
             obj.x = obj.x(keep);
             obj.y = obj.y(keep);
             obj.head = obj.head(keep);
             obj.curv = obj.curv(keep);
             obj.ArcLengths = obj.ArcLengths(keep(2:end));
         end%fcn
-
+        
         function [tau,idx] = s2tau(obj, s)
 
             sObj = obj.arcLengths0();

xUnitTests/PolygonPath/SimplifyTest.m

@@ -7,7 +7,7 @@ function testLine(testCase)
             % Create a line with intermediate points
             obj0 = PolygonPath.xy2Path(0:10, 0:10);
 
-            obj1 = obj0.simplify();
+            [obj1,keep] = obj0.simplify();
 
             % The terminal points must match
             [exp0,exp1] = obj0.termPoints();
@@ -20,6 +20,8 @@ function testLine(testCase)
 
             % Exactly two points, i.e. one segment, must survive
             verifyEqual(testCase, obj1.numel(), 1)
+            
+            verifySize(testCase, keep, [obj0.numel()+1 1])
         end%fcn
 
         function testNothingToRemove(testCase)
@@ -31,7 +33,10 @@ function testNothingToRemove(testCase)
                 [0 1 1 2 2 3 3 4 4], ...
                 [0 0 1 1 2 2 3 3 4]);
 
-            verifyEqual(testCase, obj0, obj0.simplify())
+            [obj1,keep] = obj0.simplify();
+            
+            verifyEqual(testCase, obj0, obj1)
+            verifySize(testCase, keep, [obj0.numel()+1 1])
         end%fcn
     end