--- a/+scheme/Scheme.m	Mon Sep 11 14:12:54 2017 +0200
+++ b/+scheme/Scheme.m	Mon Sep 11 14:17:15 2017 +0200
@@ -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