view +scheme/+bc/verifyFormat.m @ 1031:2ef20d00b386 feature/advectionRV

For easier comparison, return both the first order and residual viscosity when evaluating the residual. Add the first order and residual viscosity to the state of the RungekuttaRV time steppers
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 17 Jan 2019 10:25:06 +0100
parents 51cc7b05b4ab
children
line wrap: on
line source

% Errors with a more or less detailed error message if there is a problem with the bc specification
function verifyBcFormat(bcs, diffOp)
    assertType(bcs, 'cell');
    for i = 1:length(bcs)
        assertType(bcs{i}, 'struct');
        assertStructFields(bcs{i}, {'type', 'boundary'});

        if ~isfield(bcs{i}, 'data') || isempty(bcs{i}.data)
            continue
        end

        if ~isa(bcs{i}.data, 'function_handle')
            error('bcs{%d}.data should be a function of time or a function of time and space',i);
        end

        % Find dimension of boundary
        b = diffOp.grid.getBoundary(bcs{i}.boundary);
        dim = size(b,2);

        % Assert that the data function has a valid number of input arguments
        if ~(nargin(bcs{i}.data) == 1 || nargin(bcs{i}.data) == 1 + dim)
            error('sbplib:scheme:bcSetup:DataWrongNumberOfArguments', 'bcs{%d}.data has the wrong number of input arguments. Must be either only time or time and space.', i);
        end

        if nargin(bcs{i}.data) == 1
            % Grid data (only function of time)
            % Assert that the data has the correct dimension
            assertSize(bcs{i}.data(0), 1, size(b,1));
        end
    end
end