Add the following methods in LaplaceCurvilinear: getBoundaryOperator, getBoundaryQuadrature, getBoundaryBorrowing. Remove get_boundary_ops. Make interface and boundary_condition methods use the new methods.

function s = stuffStruct(varargin)
    s = struct();

    for i = 1:nargin
        assert(~isempty(inputname(i)), 'All inputs must be variables.');
        s.(inputname(i)) = varargin{i};
    end
end