+
−% SETDYFORPARAMETER Sets the according dy-error for the specified parameter.
+
−%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
−%
+
−% DESCRIPTION: SETDYFORPARAMETER Sets the according dy-error for the
+
−%              specified parameter.
+
−%
+
−% CALL:        obj = obj.setDyForParameter('a', 2);
+
−%              obj = obj.setDyForParameter('a', 2, 'b', 2);
+
−%              obj = obj.setDyForParameter(plist('errors', [2 2],
+
−%                                                'params'  {'a', 'b'}))
+
−%
+
−% INPUTS:      obj - one pest model.
+
−%              pl  - parameter list
+
−%
+
−% <a href="matlab:utils.helper.displayMethodInfo('pest', 'setDyForParameter')">Parameters Description</a>
+
−%
+
−% VERSION:     $Id: setDyForParameter.m,v 1.7 2011/04/08 08:56:25 hewitson Exp $
+
−%
+
−%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
−
+
−function varargout = setDyForParameter(varargin)
+
−  
+
−  % Check if this is a call from a class method
+
−  callerIsMethod = utils.helper.callerIsMethod;
+
−  
+
−  if callerIsMethod
+
−    
+
−    %%% make sure that the parameters and errors are cell-arrays
+
−    obj    = varargin{1};
+
−    params = cellstr(varargin{2});
+
−    errs   = varargin{3};
+
−    if ~iscell(errs)
+
−      errs = num2cell(errs);
+
−    end
+
−  else
+
−    % 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(:), '', in_names);
+
−    [pls,  invars, rest]     = utils.helper.collect_objects(rest(:), 'plist');
+
−    
+
−    params = {};
+
−    errs   = {};
+
−    %%% If pls contains only one plist with the two keys 'params' and
+
−    %%% 'errors' then set the property with a plist.
+
−    if length(pls) == 1 && isa(pls, 'plist') && nparams(pls) == 2 && isparam(pls, 'errors')
+
−      params = cellstr(find(pls, 'params'));
+
−      errs   = find(pls, 'errors');
+
−      if ~iscell(errs)
+
−        errs = num2cell(errs);
+
−      end
+
−    end
+
−    
+
−    while numel(rest) >= 2
+
−      params = [params cellstr(rest{1})];
+
−      if ~iscell(rest{2})
+
−        errs   = [errs num2cell(rest{2})];
+
−      else
+
−        errs   = [errs rest{2}];
+
−      end
+
−      rest(2) = [];
+
−      rest(1) = [];
+
−    end
+
−    
+
−    %%% Combine plists
+
−    pls = combine(pls, plist('params', params, 'errors', errs));
+
−    
+
−    if numel(params) ~= numel(errs)
+
−      error('### Please specify one error per parameter.');
+
−    end
+
−    
+
−  end % callerIsMethod
+
−  
+
−  %%% decide whether we modify the input pest object, or create a new one.
+
−  objs = copy(objs, nargout);
+
−  
+
−  for objLoop = 1:numel(objs)
+
−    
+
−    for ll=1:numel(params)
+
−      pname = params{ll};
+
−      % get index of this param
+
−      idx = find(strcmp(objs(objLoop).names, pname));
+
−      if isempty(idx)
+
−        objs(objLoop).names = [objs(objLoop).names {pname}];
+
−        objs(objLoop).dy     = [objs(objLoop).dy;    errs{ll}];
+
−      else
+
−        objs(objLoop).dy(idx) = errs{ll};
+
−      end
+
−    end
+
−    
+
−    if ~callerIsMethod
+
−      % set output history
+
−      objs(objLoop).addHistory(getInfo('None'), pls, obj_invars(objLoop), objs(objLoop).hist);
+
−    end
+
−  end
+
−  
+
−  %%% Set output
+
−  if nargout == numel(objs)
+
−    % List of outputs
+
−    for ii = 1:numel(objs)
+
−      varargout{ii} = objs(ii);
+
−    end
+
−  else
+
−    % Single output
+
−    varargout{1} = objs;
+
−  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: setDyForParameter.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);
+
−  
+
−  % Error
+
−  p = param({'errors', 'A cell-array of errors, one for each parameter to set.'}, {1, {'{}'}, paramValue.OPTIONAL});
+
−  pl.append(p);
+
−  
+
−end
+
−