% APPLYOPERATOR applys the given operator to the two input data objects.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% DESCRIPTION: APPLYOPERATOR applys the given operator to the two input
%              2D data objects.
%
% CALL:        d = applyoperator(d1, d2, pl)
%
% INPUTS:      d1 - a data2D object (tsdata, fsdata, xydata)
%              d2 - a ltpda_data object (cdata, tsdata, fsdata, xydata)
%              pl     - a plist of configuration options
%
% PARAMETERS: 'op'     - the operator to apply, e.g. 'power'
%
% VERSION:     $Id: applyoperator.m,v 1.20 2011/02/18 16:48:51 ingo Exp $
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function varargout = applyoperator(varargin)

  % Get the objects we have in the correct order
  objs = varargin(1:2);

  % Get the operator to apply
  op = varargin{3};
  
  if numel(objs) ~= 2
    error('### data2D/applyoperator requires two input data objects to work on.');
  end
  
  % Decide the type of the output object
  if isa(objs{1}, 'data2D') && isa(objs{2}, 'data2D')
    if numel(objs{1}.y) > 1
      dout = objs{1};
    elseif numel(objs{2}.y) > 1
      dout = objs{2};
    else
      dout = objs{1};
    end
  elseif isa(objs{1}, 'data2D') && isa(objs{2}, 'cdata')
    dout = objs{1};
  else
    %%% The first object in objs could not be a cdata because
    %%% then calls MATLAB cdata/applyoperator
    error('### Could/should not happen. Class obj1 [%s] class obj2 [%s]', class(objs{1}), class(objs{2}))
  end
    
  % Here we make some checks  
  s1 = size(objs{1}.y);
  s2 = size(objs{2}.y);
  trans = false;
  
  % The exponent size should be the same OR one of the two should be 1
  if ~isequal(s1, s2)
    if isequal(s1, flipdim(s2, 2))
      % We will need to transpose one vector
      trans = true;
    elseif ~isequal(s1, [1 1]) && ~isequal(s2, [1 1])
      error('### Mismatch between the data and the exponent size')
    end
  end
  
  % Calculate the uncertainty
  if ~isempty(objs{1}.dy)
    switch op
      case 'power'
        if trans
          dout.dy = objs{1}.dy .* abs(objs{2}.y' .* objs{1}.y.^(objs{2}.y'-1));
        else
          dout.dy = objs{1}.dy .* abs(objs{2}.y .* objs{1}.y.^(objs{2}.y-1));
        end
      otherwise
    end
  end
  
  % Finally, apply the operator to the data
  if trans
    dout.y = feval(op, objs{1}.y, objs{2}.y');
  else
    dout.y = feval(op, objs{1}.y, objs{2}.y);
  end

  varargout{1} = dout;
end