stokes_infsup_rbf_interpolant.m

function [infsup, varargout] = stokes_infsup_rbf_interpolant(model, model_data)
%function [infsup, varargout] = stokes_infsup_rbf_interpolant(model, model_data)
% 

varargout = cell(1, nargout-1);

if isfield(model, 'RB_infsup_interpol_method')
    method = model.RB_infsup_interpol_method;
elseif isfield(model, 'descr') && isfield(model.descr, 'RB_infsup_interpol_method')
    method = model.descr.RB_infsup_interpol_method;
else
    % default
    method = 'adaptive_cgrid';
end

switch method
    case 'adaptive_cgrid'
        %infsup = stokes_infsup_adaptive_cgrid(model, model_data);
        params = model;
        params.range = model.mu_ranges;
        for i = 1:length(params.range)
            params.range{i} = [0 1];
        end
        func = @(mu) infsup_func(mu, model, model_data);
        [infsup, varargout{:}] = cubegrid_interp_greedy(func, params);
    case 'greedy_ccriterion'
        infsup = stokes_infsup_greedy_ccriterion(model, model_data);
    otherwise
        error('stokes_infsup_rbf_interpolation:method unknown.');
end

end

function res = infsup_func(mu, model, model_data)
for i = 1:length(mu)
    mu(i) = model.mu_ranges{i}(1) + (model.mu_ranges{i}(2)-model.mu_ranges{i}(1))*mu(i);
end
model = set_mu(model, mu);
res = model.infsup_constant(model, model_data);
end