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

% % Usage example:
% c = structCellArray({'a','b'}, {
%     1, 2;
%     3, 4;
% });

function c = structCellArray(fields, values)
    assert(length(fields) == size(values, 2), 'Number of fields and number of colums of ''values'' must be equal');
    c = cell(1, size(values, 1));

    for i = 1:size(values, 1)
        for j = 1:length(fields)
            c{i}.(fields{j}) = values{i,j};
        end
    end
end