--- a/+scheme/Heat2dCurvilinear.m	Wed Jan 16 11:31:04 2019 -0800
+++ b/+scheme/Heat2dCurvilinear.m	Tue Jan 22 16:47:34 2019 +0100
@@ -321,6 +321,7 @@
         % op        -- string or a cell array of strings
         % boundary  -- string
         function varargout = getBoundaryOperator(obj, op, boundary)
+            assertIsMember(boundary, {'w', 'e', 's', 'n'})
 
             if ~iscell(op)
                 op = {op};
@@ -338,8 +339,6 @@
                         e = obj.e_l{2};
                     case 'n'
                         e = obj.e_r{2};
-                    otherwise
-                        error('No such boundary: boundary = %s',boundary);
                     end
                     varargout{i} = e;
 
@@ -353,8 +352,6 @@
                         d = obj.d1_l{2};
                     case 'n'
                         d = obj.d1_r{2};
-                    otherwise
-                        error('No such boundary: boundary = %s',boundary);
                     end
                     varargout{i} = d;
 
@@ -368,8 +365,6 @@
                         flux = obj.flux_l{2};
                     case 'n'
                         flux = obj.flux_r{2};
-                    otherwise
-                        error('No such boundary: boundary = %s',boundary);
                     end
                     varargout{i} = flux;
                 end
@@ -381,6 +376,7 @@
         %
         % boundary -- string
         function H_b = getBoundaryQuadrature(obj, boundary)
+            assertIsMember(boundary, {'w', 'e', 's', 'n'})
 
             switch boundary
                 case 'w'
@@ -391,34 +387,32 @@
                     H_b = obj.H_boundary_l{2};
                 case 'n'
                     H_b = obj.H_boundary_r{2};
-                otherwise
-                    error('No such boundary: boundary = %s',boundary);
             end
         end
 
         % Returns the boundary sign. The right boundary is considered the positive boundary
         % boundary -- string
         function s = getBoundarySign(obj, boundary)
+            assertIsMember(boundary, {'w', 'e', 's', 'n'})
+
             switch boundary
                 case {'e','n'}
                     s = 1;
                 case {'w','s'}
                     s = -1;
-                otherwise
-                    error('No such boundary: boundary = %s',boundary);
             end
         end
 
         % Returns borrowing constant gamma*h
         % boundary -- string
         function gamm = getBoundaryBorrowing(obj, boundary)
+            assertIsMember(boundary, {'w', 'e', 's', 'n'})
+
             switch boundary
                 case {'w','e'}
                     gamm = obj.h(1)*obj.alpha(1);
                 case {'s','n'}
                     gamm = obj.h(2)*obj.alpha(2);
-                otherwise
-                    error('No such boundary: boundary = %s',boundary);
             end
         end