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

% Calculates the maximum field width needed width fprintf for a given format fmt for the values in A
%     A = [1.13232 10.233 1.1];
%     width = findFieldWidth('%.2f',A);
% Gives wdith = 5
function width = findFieldWidth(fmt, A)
    width = 0;
    for i = 1:numel(A)
        elem = A(i);
        if iscell(elem)
            elem = elem{1};
        end
        str = sprintf(fmt,elem);
        width = max(width,length(str));
    end
end