view spzeros.m @ 1012:1e437c9e5132 feature/advectionRV

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
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 05 Dec 2018 13:44:10 +0100
parents 13d2f20c0c0d
children
line wrap: on
line source

function S = spzeros(varargin)
    switch length(varargin)
        case 2
            S = sparse(varargin{1}, varargin{2});
        case 1
            v = varargin{1};
            switch length(v)
                case 1
                    S = sparse(v,v);
                case 2
                    S = sparse(v(1), v(2));
                otherwise
                    error('Input must be either one integer, two integers or a vector with two integers');
            end
        otherwise
            error('Too many input arguments.');
    end
end