Add utility for constructing the operators used by a discretization emplying RV-stabilization - Add constructDiffOps which creates the operatores used for a scheme with RV-stabilization - Minor clean up in RungekuttaExteriorRV

% 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