--- a/+parametrization/Curve.m	Wed Aug 31 17:31:25 2016 +0200
+++ b/+parametrization/Curve.m	Mon Sep 05 16:34:22 2016 +0200
@@ -24,7 +24,7 @@
             if ~isempty(transformation)
                 transformation.base_g = g;
                 transformation.base_gp = gp;
-                [g,gp] = grid.Curve.transform_g(g,gp,transformation);
+                [g,gp] = parametrization.Curve.transform_g(g,gp,transformation);
             end
 
             obj.g = g;
@@ -106,7 +106,7 @@
             % New function and derivative
             g_new = @(t)obj.g(arcPar(t));
             gp_new = @(t) normalize(obj.gp(arcPar(t)));
-            curve = grid.Curve(g_new,gp_new);
+            curve = parametrization.Curve(g_new,gp_new);
 
         end
 
@@ -131,7 +131,7 @@
         function D = reverse(C)
             % g = C.g;
             % gp = C.gp;
-            % D = grid.Curve(@(t)g(1-t),@(t)-gp(1-t));
+            % D = parametrization.Curve(@(t)g(1-t),@(t)-gp(1-t));
             D = C.transform([],[],-1);
         end
 
@@ -157,7 +157,7 @@
                 transformation.flip = flip*flip_old;
             end
 
-            D = grid.Curve(g,gp,transformation);
+            D = parametrization.Curve(g,gp,transformation);
 
         end
 
@@ -171,7 +171,7 @@
             %     v(2,:) = x(2,:)+a(2);
             % end
 
-            % D = grid.Curve(@g_fun,gp);
+            % D = parametrization.Curve(@g_fun,gp);
 
             D = C.transform([],a);
         end
@@ -206,7 +206,7 @@
             %     v = A*gp(t);
             % end
 
-            % D = grid.Curve(@g_fun,@gp_fun);
+            % D = parametrization.Curve(@g_fun,@gp_fun);
 
             % g = A(g-a)+a = Ag - Aa + a;
             b = - A*a + a;
@@ -237,7 +237,7 @@
             %     v = A*gp(t);
             % end
 
-            % D = grid.Curve(@g_fun,@gp_fun);
+            % D = parametrization.Curve(@g_fun,@gp_fun);
 
 
              % g = A(g-a)+a = Ag - Aa + a;
@@ -268,7 +268,7 @@
             end
             g = @g_fun;
 
-            obj = grid.Curve(g);
+            obj = parametrization.Curve(g);
         end
 
         function obj = circle(c,r,phi)
@@ -288,7 +288,7 @@
                 v(2,:) =  (phi(2)-phi(1))*r*cos(w);
             end
 
-            obj = grid.Curve(@g_fun,@g_fun_deriv);
+            obj = parametrization.Curve(@g_fun,@g_fun_deriv);
         end
 
         function obj = bezier(p0, p1, p2, p3)
@@ -302,7 +302,7 @@
                 v(2,:) = 3*(1-t).^2*(p1(2)-p0(2)) + 6*(1-t).*t*(p2(2)-p1(2)) + 3*t.^2*(p3(2)-p2(2));
             end
 
-            obj = grid.Curve(@g_fun,@g_fun_deriv);
+            obj = parametrization.Curve(@g_fun,@g_fun_deriv);
         end