--- a/+parametrization/Ti.m	Mon Sep 05 16:51:22 2016 +0200
+++ b/+parametrization/Ti.m	Tue Sep 06 13:22:28 2016 +0200
@@ -36,6 +36,8 @@
             obj.S = @S_fun;
         end
 
+        % Does this funciton make sense?
+        % Should it always be eval?
         function [X,Y] = map(obj,u,v)
             default_arg('v',u);
 
@@ -65,6 +67,19 @@
             end
         end
 
+        % Evaluate S for each pair of u and v,
+        % Return same shape as u
+        function [x, y] = eval(obj, u, v)
+            x = zeros(size(u));
+            y = zeros(size(u));
+
+            for i = 1:numel(u)
+                p = obj.S(u(i), v(i));
+                x(i) = p(1,:);
+                y(i) = p(2,:);
+            end
+        end
+
         function h = plot(obj,nu,nv)
             S = obj.S;