view +scheme/+bc/verifyFormat.m @ 1134:7259b3f3e227 feature/laplace_curvilinear_test

Make LaplaceCurvilinear work for coeff a and b that are scalar, grid functions, or function handles.
author Martin Almquist <malmquist@stanford.edu>
date Mon, 27 May 2019 16:53:32 -0700
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