Mercurial > hg > ltpda
view m-toolbox/classes/@smodel/assignalias.m @ 36:5eb86f6881ef database-connection-manager
Remove commented-out code
author | Daniele Nicolodi <nicolodi@science.unitn.it> |
---|---|
date | Mon, 05 Dec 2011 16:20:06 +0100 |
parents | f0afece42f48 |
children |
line wrap: on
line source
% ASSIGNALIAS assign values to smodel alias %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % DESCRIPTION: ASSIGNALIAS assign numerical values or vectors to smodel % aliases. This is a processing method which should be run before using the % given smodel inside smodel/fftfilt so to gain evaluation time during the % evaluation of the model. Be careful to insert the correct values for the % parameters otherwhise smodel/double will throw an error. % % CALL: mdl = assignalias(mdl) % % <a href="matlab:utils.helper.displayMethodInfo('smodel', 'assignalias')">Parameters Description</a> % % VERSION: $Id: assignalias.m,v 1.6 2011/05/10 11:07:59 luigi Exp $ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function varargout = assignalias(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 smodels and plists [as, smodel_invars, rest] = utils.helper.collect_objects(varargin(:), 'smodel', in_names); pl = utils.helper.collect_objects(varargin(:), 'plist'); % Copy the object(s), so to inherit parameters etc mdls = copy(as, nargout); % combine plists if isempty(pl) setvar = 'fftfilt'; else xvals = find(pl, 'xvals'); if isempty(xvals) setvar = 'fftfilt'; else setvar = 'UserDefinedXvals'; end end pl = parse(pl, getDefaultPlist(setvar)); % get parameters values from plist switch lower(setvar) case 'fftfilt' nsecs = find(pl, 'nsecs'); npad = find(pl, 'npad'); fs = find(pl, 'fs'); % get fft frequancies nfft = nsecs*fs + npad; xvals = utils.math.getfftfreq(nfft,fs,'one'); xvals = xvals.'; case 'userdefinedxvals' xvals = find(pl, 'xvals'); end % run over input objects for ii=1:numel(mdls) % Recover the mapping factor from xvals and xvar % get xvar xxvar = mdls(ii).xvar; % check dimensions if numel(xxvar)>1 error('Multiple xvar are not supported!') else xxvar = xxvar{:}; end trans = mdls(ii).trans; if isempty(trans) scale = 1.0; else if isnumeric(trans) scale = trans; elseif ischar(trans) scale = eval(trans); elseif iscell(trans) % check dimension if numel(trans)>1 error('Multiple trans are not supported!') else switch lower(class(trans{:})) case 'double' scale = trans{:}; case 'char' scale = eval(trans{:}); end end else error('Unknown format for the transformation!'); end end % assign values for the x getSingleVariable(xxvar,scale.*xvals); % assign alias values for jj=1:numel(mdls(ii).aliasNames) switch class(mdls(ii).aliasValues{jj}) case 'char' getSingleVariable('calias',eval(mdls(ii).aliasValues{jj})); mdls(ii).aliasValues{jj} = calias; case 'smodel' tmd = mdls(ii).aliasValues{jj}; tmd.setXvals(xvals); getSingleVariable('calias',tmd.double); mdls(ii).aliasValues{jj} = calias; otherwise % do nothing end end mdls(ii).setXvals(xvals); % set output history mdls(ii).addHistory(getInfo('None'), pl, smodel_invars(ii), mdls(ii).hist); end %%% Set output if nargout == numel(mdls) % List of outputs for ii = 1:numel(mdls) varargout{ii} = mdls(ii); end else % Single output varargout{1} = mdls; end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Local Functions % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % FUNCTION: getCellVariables % % DESCRIPTION: Assign values to variables % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function getSingleVariable(nms,val) assignin('caller',nms,val); end %-------------------------------------------------------------------------- % Get Info Object %-------------------------------------------------------------------------- function ii = getInfo(varargin) if nargin == 1 && strcmpi(varargin{1}, 'None') sets = {}; pl = []; elseif nargin == 1 && ~isempty(varargin{1}) && ischar(varargin{1}) sets{1} = varargin{1}; pl = getDefaultPlist(sets{1}); else sets = SETS(); % get plists pl(size(sets)) = plist; for kk = 1:numel(sets) pl(kk) = getDefaultPlist(sets{kk}); end end % Build info object ii = minfo(mfilename, 'smodel', 'ltpda', utils.const.categories.helper, '$Id: assignalias.m,v 1.6 2011/05/10 11:07:59 luigi Exp $', sets, pl); end %-------------------------------------------------------------------------- % Defintion of Sets %-------------------------------------------------------------------------- function out = SETS() out = {... 'fftfilt', ... 'UserDefinedXvals' ... }; end %-------------------------------------------------------------------------- % Get Default Plist %-------------------------------------------------------------------------- function plout = getDefaultPlist(set) persistent pl; persistent lastset; if exist('pl', 'var')==0 || isempty(pl) || ~strcmp(lastset, set) pl = buildplist(set); lastset = set; end plout = pl; end function pl = buildplist(set) pl = plist.EMPTY_PLIST; switch lower(set) case 'fftfilt' p = param({'nsecs', 'Number of seconds of the time series will be filtered with fftfilt.'}, paramValue.DOUBLE_VALUE(1)); pl.append(p); p = param({'npad', 'Number of sample pad will will be used in fftfilt.'}, paramValue.EMPTY_DOUBLE); pl.append(p); p = param({'fs', 'Sampling frequency of the time series will be filtered with fftfilt.'}, paramValue.DOUBLE_VALUE(1)); pl.append(p); case 'userdefinedxvals' p = param({'xvals', 'A vector of values for the X variable.'... 'If the smodel implement a transformation for x values then you should'... 'input the value before the transformation. E.g. xvar = s -> trans = 2*pi*i'... 'you should input f values so you finally get inside the code xvals = 2*pi*i*f'}, paramValue.EMPTY_DOUBLE); pl.append(p); end end