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

% Replaces string old on the console with new.
% new may be a format string like the one accepted by fprintf
% old has to be the last string printed.
function s = replace_string(old,new,varargin)
    reverseStr = repmat(sprintf('\b'), 1, length(old));
    blankStr   = repmat(sprintf(' ') , 1, length(old));
    fprintf([reverseStr,blankStr,reverseStr, new],varargin{:});
    s = sprintf(new, varargin{:});
end