+
−% SUBS substitutes symbolic parameters with the given values.
+
−%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
−%
+
−% DESCRIPTION: SUBS substitutes symbolic parameters with the given values.
+
−%
+
−% CALL:        mdl = mdl.subs(pl)
+
−%              mdl = mdl.subs('all')       % Substitutes all parameters
+
−%              mdl = mdl.subs('a', 'b')    % Substitutes the parameters 'a' and
+
−%                                            'b' with their default values
+
−%              mdl = mdl.subs({'a', 'b'})  % Substitutes the parameters 'a' and
+
−%                                            'b' with their default values
+
−%              mdl = mdl.subs('a', 33)     % Substitutes the parameters 'a'
+
−%                                            with 33
+
−%              mdl = mdl.subs({'a'}, {33}) % Substitutes the parameters 'a'
+
−%                                            with 33
+
−%
+
−% Examples
+
−% --------
+
−%
+
−% 1)   m = subs(m, plist('Params', 'all')) % substitute all default values
+
−% 2)   m = subs(m, plist('Params', {'a', 'b'}, 'Values',{1, 1:10})) % substitute
+
−%                                                            values
+
−% 3)   m = subs(m, plist('Params', {'a', 'b'})) % substitute default values
+
−%
+
−% <a href="matlab:utils.helper.displayMethodInfo('smodel', 'subs')">Parameters Description</a>
+
−%
+
−% VERSION:     $Id: subs.m,v 1.11 2011/04/08 08:56:28 hewitson Exp $
+
−%
+
−%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
−
+
−function varargout = subs(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
+
−  
+
−  % Collect all AOs and plists
+
−  [mdls, mdl_invars, rest] = utils.helper.collect_objects(varargin(:), 'smodel', in_names);
+
−  [pls,  pl_invars,  rest] = utils.helper.collect_objects(rest, 'plist', in_names);
+
−  
+
−  % combine plists
+
−  pl = combine(pls, getDefaultPlist());
+
−  
+
−  % Check if the user uses a plist or not
+
−  if isempty(pls) && ~isempty(rest)
+
−    % The user doesn't use a input plist
+
−    iparams     = {};
+
−    ivals       = {};
+
−    iexceptions = {};
+
−    % Collect each character or cellstr in 'rest' as a parameter
+
−    % Collect each integer of numeric cell 
+
−    for ii=1:numel(rest)
+
−      if ischar(rest{ii}) && strcmp(rest{ii}, 'all') && ii == 1
+
−        iparams = 'all';
+
−        break;
+
−      elseif ischar(rest{ii})
+
−        iparams = [iparams rest{ii}];
+
−      elseif iscellstr(rest{ii})
+
−        iparams = [iparams rest{ii}];
+
−      elseif isnumeric(rest{ii})
+
−        ivals = [ivals rest{ii}];
+
−      elseif iscell(rest{ii})
+
−        c = rest{ii};
+
−        for kk=1:numel(c)
+
−          if isnumeric(c{kk}), ivals = [ivals c{kk}]; end
+
−        end
+
−      end
+
−    end
+
−    pl.pset('params', iparams);
+
−    pl.pset('values', ivals);
+
−    
+
−  else
+
−    % The user uses a input plist
+
−    iparams     = find(pl, 'params');
+
−    ivals       = find(pl, 'values');
+
−    iexceptions = find(pl, 'exceptions');
+
−  end
+
−  
+
−  %
+
−  
+
−  % Loop over input models
+
−  bs = copy(mdls, nargout);
+
−  for j=1:numel(mdls)
+
−    
+
−    mdl = bs(j);
+
−    
+
−    if ischar(iparams)
+
−      if strcmpi(iparams, 'all')
+
−        iparams = mdl.params;
+
−        ivals   = mdl.values;
+
−      else
+
−        iparams = {iparams};
+
−      end
+
−    end
+
−    
+
−    if isempty(ivals)
+
−      ivals = {};
+
−      % get values from the model
+
−      for ll=1:numel(iparams)
+
−        idx = find(strcmp(iparams{ll}, mdl.params));
+
−        ivals = [ivals mdl.values(idx)];
+
−      end
+
−    end
+
−    
+
−    % remove the exceptions from the params list
+
−    if ~isempty(iexceptions)
+
−      idx = ismember(iparams, iexceptions);
+
−      iparams = iparams(~idx);
+
−      ivals   = ivals(~idx);
+
−    end
+
−    
+
−    if isempty(ivals) && ~strcmpi(iparams, 'all')
+
−      error('### Please specify one value per parameter to substitute.');
+
−    end
+
−    if ~ischar(iparams) && numel(iparams) ~= numel(ivals)
+
−      error('### Please specify one value per parameter to substitute.');
+
−    end
+
−    
+
−    if ~iscell(ivals)
+
−      ivals = {ivals};
+
−    end
+
−    
+
−    if numel(iparams) ~= numel(ivals)
+
−      error('### The number of parameters and values doesn''t match for model %s', mdl.name);
+
−    end
+
−    
+
−    % Get remaining parameters and values for the output model
+
−    oparams = {};
+
−    ovals   = {};
+
−    for kk=1:numel(mdl.params)
+
−      if ~utils.helper.ismember(mdl.params{kk}, iparams)
+
−        oparams = [oparams mdl.params(kk)];
+
−        if ~isempty(mdl.values)
+
−          ovals   = [ovals mdl.values(kk)];
+
−        end
+
−      end
+
−    end
+
−    
+
−    if ~isempty(iparams) && ~isnumeric(mdl.expr)
+
−      mdl.expr = subs(mdl.expr, iparams, ivals);
+
−      mdl.values = ovals;
+
−      mdl.params = oparams;
+
−      % add history
+
−      mdl.addHistory(getInfo('None'), pl, mdl_invars(j), mdls(j).hist);
+
−    end
+
−    
+
−  end
+
−  
+
−  % Set outputs
+
−  varargout{1} = bs;
+
−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, 'smodel', 'ltpda', utils.const.categories.helper, '$Id: subs.m,v 1.11 2011/04/08 08:56:28 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', 'The parameters to substitute for.<br>Specify ''all'' to substitute all.'}, {1, {'all'}, paramValue.OPTIONAL});
+
−  pl.append(p);
+
−  
+
−  % Values
+
−  p = param({'values', 'A cell-array of values to set that overide the defaults.'}, {1, {{}}, paramValue.OPTIONAL});
+
−  pl.append(p);
+
−  
+
−  % Exceptions
+
−  p = param({'exceptions', 'A cell-array of parameters which will not be substitute.'}, paramValue.EMPTY_CELL);
+
−  pl.append(p);
+
−  
+
−end
+
−
+
−%
+
−% PARAMETERS:  'Params' - A cell array of the parameter names to
+
−%                         substitute. If you specify 'all' for this
+
−%                         parameter, then the current values of mdl.params
+
−%                         will be substituted.
+
−%              'Values' - The values to substitute.
+
−%