--- a/+scheme/Schrodinger1dCurve.m	Wed Feb 15 11:21:04 2017 +0100
+++ b/+scheme/Schrodinger1dCurve.m	Thu Feb 23 09:31:56 2017 +0100
@@ -107,17 +107,17 @@
         function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
             % u denotes the solution in the own domain
             % v denotes the solution in the neighbour domain
-            %             [e_u,d_u,s_u] = obj.get_boundary_ops(boundary);
-            %             [e_v,d_v,s_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary);
+                         [e_u,d_u,s_u] = obj.get_boundary_ops(boundary);
+                         [e_v,d_v,s_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary);
             
-            %             a =  -s_u* 1/2 * 1i ;
-            %             b =  a';
+                         a =  -s_u* 1/2 * 1i ;
+                         b =  a';
             
-            %             tau = b*d_u;
-            %             sig = -a*e_u;
+                         tau = b*d_u;
+                         sig = -a*e_u;
             
-            %             closure = obj.Hi * (tau*e_u' + sig*d_u');
-            %             penalty = obj.Hi * (-tau*e_v' - sig*d_v');
+                         closure = obj.Hi * (tau*e_u' + sig*d_u');
+                         penalty = obj.Hi * (-tau*e_v' - sig*d_v');
         end
         
         % Ruturns the boundary ops and sign for the boundary specified by the string boundary.