% APPLYMETHOD applys the given method to the input 2D data.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% DESCRIPTION: APPLYMETHOD applys the given method to the input 2D data.
%
% CALL:        pl = applymethod(d, pl)
%              pl = applymethod(d, pl, fcns)
%
% INPUTS:      d      - a 2D data object (tsdata, fsdata, xydata)
%              pl     - a plist of configuration options
%              fcns   - function handle(s) for the evaluation of the uncertainty
%                       (alone or in a cell array)
%
% PARAMETERS:
%
%       'method' - the method to apply to the data
%       'axis'   - which axis vector to apply the method to. Possible values
%                  are: 'X', 'Y', 'XY' [default: 'Y']
%       'dim'    - the dimension of the chosen vector to apply the method
%                  to. This is necessary for functions like mean() when
%                  applied to matrices held in cdata objects. For tsdata,
%                  fsdata or xydata, this option has no effect.
%                  [default: 1]
%       'option' - any additional option to pass to the method.
%
% VERSION:     $Id: applymethod.m,v 1.15 2011/04/17 09:13:18 hewitson Exp $
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function pl = applymethod(ds, pl, method, getDefaultPlist, varargin)
  
  % Get function handles
  dxFcn = {};
  for jj = 1:numel(varargin)
    if isa(varargin{jj}, 'function_handle')
      dxFcn = varargin{jj};
    end
    if iscell(varargin{jj})
      list = varargin{jj};
      if ~isempty(list) && isa(list{1}, 'function_handle')
        dxFcn = list{1};
      end
    end
  end
  
  pl = applyDefaults(getDefaultPlist('2D'), pl);
  
  % Get the axis we are dealing with
  axis = find(pl, 'axis');
  % Get the dimension to operate along
  dim = find(pl, 'dim');
  % Get any additional option
  opt = find(pl, 'option');
  
  % Loop over data objects
  for jj=1:numel(ds)
    switch lower(axis)
      case 'x'
        if ~isempty(ds(jj).dx) && ~isempty(dxFcn)
          ds(jj).dx = feval(dxFcn, ds(jj).getX, ds(jj).dx);
        else
          ds(jj).dx = [];
        end
        ds(jj).x  = apply(ds(jj).getX, method, dim, opt);
      case 'y'
        if ~isempty(ds(jj).dy) && ~isempty(dxFcn)
          ds(jj).dy = feval(dxFcn, ds(jj).y, ds(jj).dy);
        else
          ds(jj).dy = [];
        end
        ds(jj).y  = apply(ds(jj).y, method, dim, opt);
      case 'xy'
        if ~isempty(ds(jj).dx) && ~isempty(dxFcn)
          ds(jj).dx = feval(dxFcn, ds(jj).getX, ds(jj).dx);
        else
          ds(jj).dx = [];
        end
        if ~isempty(ds(jj).dy) && ~isempty(dxFcn)
          ds(jj).dy = feval(dxFcn, ds(jj).y, ds(jj).dy);
        else
          ds(jj).dy = [];
        end
        ds(jj).x  = apply(ds(jj).getX, method, dim, opt);
        ds(jj).y  = apply(ds(jj).y, method, dim, opt);
      otherwise
        error('### Unsupported axis ''%s'' to operate on.', axis);
    end
  end
  
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                               Local Functions                               %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%-----------------------------------------------
% Apply method to the vector v
%-----------------------------------------------
function v = apply(v, method, dim, opt)
  if ~isempty(dim) && ~isempty(opt)
    % User supplied a dimension and an option
    v = feval(method, v, dim, opt);
  elseif ~isempty(dim)
    % User supplied only a dimension
    v = feval(method, v, dim);
  elseif ~isempty(opt)
    % User supplied only an option
    v = feval(method, v, opt);
  else
    % User supplied only a method
    v = feval(method, v);
  end
end