add_integral_kernels.m

Go to the documentation of this file.

function res = add_integral_kernels(kernels, varargin)
%function res = add_integral_kernels(kernels, varargin)
% Addition of integral kernel evaluations
%
% Integral kernels are evaluated at the same arguments and the
% results are added up. If the results are cell arrays, they are
% concatenated. This function is used in Fem.Assembly.operator().
%
% Parameters:
%  kernels: cell array containing function handles corresponding to
%  Fem integral kernels.
%  varargin: inputs for the integral kernels as variable-length
%  input argument list.
%
% Return values:
%  res: result of addition or cell concatenation.

res_tmp = cell(1, length(kernels));

for q = 1:length(kernels)
    
    res_tmp{q} = kernels{q}(varargin{:});
end

if iscell(res_tmp{1})
    for q = 1:length(res_tmp)
        res_tmp{q} = res_tmp{q}(:)';
    end
    
    res = [res_tmp{:}];
else
    
    res = res_tmp{1};    
    for q = 2:length(res_tmp)
        res = res + res_tmp{q};
    end
end

end