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

function print_issparse(A)
    warning('Deprecated! Use printExpr() instead!');
    b = issparse(A);
    if b
        s = 'true';
    else
        s = 'false';
    end
    fprintf('%8s is sparse: %s\n',inputname(1),s);
end