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 o = labelTask(in)

    switch class(in)
        case 'char'
            fprintf(in);
            o = tic();
        case 'uint64'
            o = toc(in);
            fprintf(' - done %fs\n', o);
        otherwise
            error('Unknow input type: %s', class(in))
    end