% FIND Creates analysis objects from the selected parameter(s).
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% DESCRIPTION: FIND Creates analysis objects from the selected
%              parameter(s).
%
% CALL:        obj = obj.find();
%              obj = obj.find('a');
%              obj = obj.find('a', 'b');
%              obj = obj.find(plist('params', {'a', 'b'}))
%
% INPUTS:      obj - one pest model.
%              pl  - parameter list
%
% <a href="matlab:utils.helper.displayMethodInfo('pest', 'find')">Parameters Description</a>
%
% VERSION:     $Id: find.m,v 1.7 2011/04/08 08:56:25 hewitson Exp $
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function varargout = find(varargin)
  
  % Check if this is a call for parameters
  if utils.helper.isinfocall(varargin{:})
    varargout{1} = getInfo(varargin{3});
    return
  end
  
  %%% Collect input variable names
  in_names = cell(size(varargin));
  for ii = 1:nargin,in_names{ii} = inputname(ii);end
  
  [objs, obj_invars, rest] = utils.helper.collect_objects(varargin(:), 'pest', in_names);
  [pls,  invars, rest]     = utils.helper.collect_objects(rest(:), 'plist');
  
  inParams = {};
  out      = [];
  %%% If pls contains only one plist with the key 'params' then set the
  %%% property with a plist.
  if length(pls) == 1 && isa(pls, 'plist') && pls.nparams == 1 && pls.isparam('params')
    inParams = cellstr(find(pls, 'params'));
  end
  
  if ~iscellstr(rest)
    error('### Can not handle the input: %s', utils.helper.val2str(rest));
  end
  
  inParams = [inParams rest];
  
  for objLoop = 1:numel(objs)
    
    %%% decide whether we modify the input pest object, or create a new one.
    objs(objLoop) = copy(objs(objLoop), nargout);
    
    %%% Create from each parameter an AO if the user doesn't specify a
    %%% parameter name.
    if isempty(inParams)
      params = objs(objLoop).names;
    else
      params = inParams;
    end
    
    for ll=1:numel(params)
      pname = params{ll};
      % get index of this param
      idx = find(strcmp(objs(objLoop).names, pname), 1);
      if isempty(idx)
        % Do nothing because it is not possible to create an AO from a not
        % existing parameter name.
      else
        a = ao(objs(objLoop).y(idx));
        a.setName(pname);
        if ~isempty(objs(objLoop).dy)
          a.setDy(objs(objLoop).dy(idx));
        end
        if ~isempty(objs(objLoop).yunits)
          a.setYunits(objs(objLoop).yunits(idx));
        end
        if numel(objs(objLoop).models) > 1
          a.setProcinfo(plist('model', objs(objLoop).models(idx)));
        else
          a.setProcinfo(plist('model', objs(objLoop).models));
        end
        if ~strcmp(varargin{end}, 'internal')
          a.addHistory(pest.getInfo('find', 'None'), plist('params', pname), obj_invars(objLoop), objs(objLoop).hist);
        end
        out = [out a];
      end
    end
    
  end
  
  %%% Set output
  if nargout == numel(out)
    % List of outputs
    for ii = 1:numel(out)
      varargout{ii} = out(ii);
    end
  else
    % Single output
    varargout{1} = out;
  end
  
end

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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% FUNCTION:    getInfo
%
% DESCRIPTION: Get Info Object
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function ii = getInfo(varargin)
  if nargin == 1 && strcmpi(varargin{1}, 'None')
    sets = {};
    pl   = [];
  else
    sets = {'Default'};
    pl   = getDefaultPlist;
  end
  % Build info object
  ii = minfo(mfilename, 'pest', 'ltpda', utils.const.categories.helper, '$Id: find.m,v 1.7 2011/04/08 08:56:25 hewitson Exp $', sets, pl);
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% FUNCTION:    getDefaultPlist
%
% DESCRIPTION: Get Default Plist
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function plout = getDefaultPlist()
  persistent pl;  
  if exist('pl', 'var')==0 || isempty(pl)
    pl = buildplist();
  end
  plout = pl;  
end

function pl = buildplist()
  
  pl = plist();
  
  % Params
  p = param({'params', 'A cell-array of parameter names.'}, {1, {'{}'}, paramValue.OPTIONAL});
  pl.append(p);

end