Create residual viscosity package +rv and generalize the ResidualViscosity class - Generalize residual viscosity, by passing user-defined flux and calculating the time derivative outside of the update. - Create separate RungekuttaRV specifically using interior RV updates - Separate the artifical dissipation operator from the scheme AdvectionRV1D so that the same scheme can be reused for creating the diff op used by the ResidualViscosity class

% Discretises a shape such that points on the curves are no further than h appart.
function p = shape_discretise(s,h)
    p = [];
    for i = 1:length(s)
        [~,pt] = parametrization.curve_discretise(s{i},h);
        p = [p, pt];
    end
end