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

% Reshapes a matrix as if it was stored in row major order.
function B = reshapeRowMaj(A, m)
    D = length(m);

    if D == 1
        m = [m 1];
        D = 2;
    end

    % Reshape and reverse order of indecies
    B = permute(reshape(permute(A, ndims(A):-1:1), rot90(m,2)), D:-1:1);
end