dirichlet_values_leftright.m

function Udirichlet = dirichlet_values_leftright(glob, params)
% function Udirichlet = dirichlet_values_leftright(glob, params)
%
% function computing dirichlet-values by pointwise evaluations
%
% required fields of params:
%   dir_middle         : x coordinate where left and right regions are separated
%   c_dir_left         : left dirichlet value
%   c_dir_right        : right dirichlet value
%
% optional fields of params:
%   mu_names : names of fields to be regarded as parameters in vector mu
%
% in 'coefficients' mode, the parameters in brackets are empty

% glob column check
if params.debug
  if ~isempty(glob) && size(glob,1) < size(glob,2)
    warning('coordinates in variable glob are given row-wise, but expected them to be column-wise');
    if params.debug > 2
      keyboard;
    end
  end
end
decomp_mode = params.decomp_mode;


if decomp_mode == 2
  Udirichlet = [params.c_dir_left, params.c_dir_right];
else
  i = (glob(:,1) <= params.dir_middle);
  if decomp_mode == 1
    Udirichlet = {i, (1-i)};
  elseif decomp_mode == 0
    Udirichlet = params.c_dir_left * i + (1-i) * params.c_dir_right;
  end
end

%| \docupdate