view +scheme/+bc/verifyFormat.m @ 1037:2d7ba44340d0 feature/burgers1d

Pass scheme specific parameters as cell array. This will enabale constructDiffOps to be more general. In addition, allow for schemes returning function handles as diffOps, which is currently how non-linear schemes such as Burgers1d are implemented.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 18 Jan 2019 09:02:02 +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