--- a/+scheme/Beam.m	Mon Dec 03 16:49:43 2018 -0800
+++ b/+scheme/Beam.m	Mon Jan 07 16:26:00 2019 +0100
@@ -19,7 +19,7 @@
         alphaII
         alphaIII
 
-        opt
+        opt % TODO: Get rid of this and use the interface type instead
     end
 
     methods
@@ -170,7 +170,7 @@
             end
         end
 
-        function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
+        function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary, type)
             % u denotes the solution in the own domain
             % v denotes the solution in the neighbour domain
             [e_u,d1_u,d2_u,d3_u,s_u] = obj.get_boundary_ops(boundary);