--- a/+scheme/Scheme.m	Fri Nov 03 10:43:27 2017 +0100
+++ b/+scheme/Scheme.m	Fri Nov 03 10:53:15 2017 +0100
@@ -25,9 +25,12 @@
         %       neighbour_boundary  is a string specifying which boundary to
         %                           interface to.
         %       penalty  may be a cell array if there are several penalties with different weights
-        [closure, penalty] = boundary_condition(obj,boundary,type)
+        [closure, penalty] = boundary_condition(obj,boundary,type) % TODO: Change name to boundaryCondition
         [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
 
+        % TODO: op = getBoundaryOperator()??
+        %   makes sense to have it available through a method instead of random properties
+
         % Returns the number of degrees of freedom.
         N = size(obj)
     end