--- a/+scheme/Utux.m	Wed Dec 05 16:08:48 2018 +0100
+++ b/+scheme/Utux.m	Wed Dec 05 15:27:44 2018 -0800
@@ -17,27 +17,16 @@
 
 
     methods
-         function obj = Utux(g ,order, operator)
-             default_arg('operator','Standard');
-
-             m = g.size();
-             xl = g.getBoundary('l');
-             xr = g.getBoundary('r');
-             xlim = {xl, xr};
+        function obj = Utux(g, order, opSet)
+            default_arg('opSet',@sbp.D2Standard);
 
-           switch operator
-%                case 'NonEquidistant'
-%               ops = sbp.D1Nonequidistant(m,xlim,order);
-%               obj.D1 = ops.D1;
-               case 'Standard'
-              ops = sbp.D2Standard(m,xlim,order);
-              obj.D1 = ops.D1;
-%                case 'Upwind'
-%               ops = sbp.D1Upwind(m,xlim,order);
-%               obj.D1 = ops.Dm;
-               otherwise
-                   error('Unvalid operator')
-           end
+            m = g.size();
+            xl = g.getBoundary('l');
+            xr = g.getBoundary('r');
+            xlim = {xl, xr};
+
+            ops = opSet(m, xlim, order);
+            obj.D1 = ops.D1;
 
             obj.grid = g;