- Fix bug when constructing closure for narrow stencils - Update the residual outside of the RK time steps. At least for the inner convergence rate, updating the residual inside does not seem to be required.

% Get the index of the points p within the tall array of points ps
function [I, ok] = pointIndex(p, ps)
    [ok, I] = ismember(p,  ps, 'rows');
end