diff m-toolbox/classes/@pest/setYunitsForParameter.m @ 0:f0afece42f48

Import.
author Daniele Nicolodi <nicolodi@science.unitn.it>
date Wed, 23 Nov 2011 19:22:13 +0100
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/m-toolbox/classes/@pest/setYunitsForParameter.m	Wed Nov 23 19:22:13 2011 +0100
@@ -0,0 +1,176 @@
+% SETYUNITSFORPARAMETER Sets the according y-unit for the specified parameter.
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% DESCRIPTION: SETYUNITSFORPARAMETER Sets the according y-unit for the
+%              specified parameter.
+%
+% CALL:        obj = obj.setYunitsForParameter('a', 'Hz');
+%              obj = obj.setYunitsForParameter('a', unit('s'), 'b', 'Hz');
+%              obj = obj.setYunitsForParameter(plist('units', [unit('s') unit('Hz')],
+%                                               'params'  {'a', 'b'}))
+%
+% INPUTS:      obj - one pest object.
+%              pl  - parameter list
+%
+% <a href="matlab:utils.helper.displayMethodInfo('pest', 'setYunitsForParameter')">Parameters Description</a>
+%
+% VERSION:     $Id: setYunitsForParameter.m,v 1.8 2011/04/08 08:56:25 hewitson Exp $
+%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+function varargout = setYunitsForParameter(varargin)
+  
+  % Check if this is a call from a class method
+  callerIsMethod = utils.helper.callerIsMethod;
+  
+  if callerIsMethod
+    %%% make sure that the parameters and units are cell-arrays
+    objs   = varargin{1};
+    params = cellstr(varargin{2});
+    units  = varargin{3};
+    if ischar(units), units = cellstr(units); end
+    for uu = 1:numel(units)
+      units{uu} = unit(units{uu});
+    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 = {};
+    units  = {};
+    %%% If pls contains only one plist with the two keys 'params' and
+    %%% 'units' then set the property with a plist.
+    if length(pls) == 1 && isa(pls, 'plist') && nparams(pls) == 2 && isparam(pls, 'units')
+      params = cellstr(find(pls, 'params'));
+      units  = find(pls, 'units');
+      if ischar(units)
+        units = cellstr(units);
+      end
+      for uu = 1:numel(units)
+        units{uu} = unit(units{uu});
+      end
+    end
+    
+    while numel(rest) >= 2
+      params = [params cellstr(rest{1})];
+      if ~iscell(rest{2})
+        units   = [units {unit(rest{2})}];
+      else
+        for uu = 1:numel(rest{2})
+          units   = [units {unit(unit(rest{2}{uu}))}];
+        end
+      end
+      rest(2) = [];
+      rest(1) = [];
+    end
+    
+    %%% Combine plists
+    pls = combine(pls, plist('params', params, 'units', units));
+    
+    if numel(params) ~= numel(units)
+      error('### Please specify one value 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).yunits = [objs(objLoop).yunits; units{ll}];
+      else
+        objs(objLoop).yunits(idx) = units{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: setYunitsForParameter.m,v 1.8 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);
+  
+  % Value
+  p = param({'units', 'A cell-array of units, one for each parameter to set.'}, {1, {'{}'}, paramValue.OPTIONAL});
+  pl.append(p);
+  
+end
+