Mercurial > hg > ltpda
view m-toolbox/classes/@cdata/applyoperator.m @ 15:ce3fbb7ebe71 database-connection-manager
Remove broken functions from utils.jmysql
author | Daniele Nicolodi <nicolodi@science.unitn.it> |
---|---|
date | Mon, 05 Dec 2011 16:20:06 +0100 |
parents | f0afece42f48 |
children |
line wrap: on
line source
% APPLYOPERATOR applys the given operator to the two input data objects. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % DESCRIPTION: APPLYOPERATOR applys the given operator to the two input % cdata objects. % % CALL: d = applyoperator(d1, d2, pl) % % INPUTS: d1 - a cdata object % d2 - a ltpda_data object (cdata, tsdata, fsdata, xydata) % pl - a plist of configuration options % % PARAMETERS: 'op' - the operator to apply, e.g. 'power' % % VERSION: $Id: applyoperator.m,v 1.15 2011/02/18 16:48:51 ingo Exp $ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function varargout = applyoperator(varargin) % Get the objects we have in the correct order objs = varargin(1:2); % Get the operator to apply op = varargin{3}; if numel(objs) ~= 2 error('### cdata/applyoperator requires two input data objects to work on.'); end % Decide the type of the output object if isa(objs{1}, 'data2D') %%% The first object in objs could not be a data2D because %%% then calls MATLAB data2D/applyoperator error('### Could/should not happen. Class obj1 [%s] class obj2 [%s]', class(objs{1}), class(objs{2})) elseif isa(objs{2}, 'data2D') dout = objs{2}; else dout = objs{1}; end % Here we make some checks s1 = size(objs{1}.y); s2 = size(objs{2}.y); trans = false; % The exponent size should be the same OR one of the two should be 1 if ~isequal(s1, s2) if isequal(s1, flipdim(s2, 2)) % We will need to transpose one vector trans = true; elseif ~isequal(s1, [1 1]) && ~isequal(s2, [1 1]) error('### Mismatch between the data and the exponent size') end end % Calculate the uncertainty if ~isempty(objs{1}.dy) switch op case 'power' if trans dout.dy = objs{1}.dy .* abs(objs{2}.y' .* objs{1}.y.^(objs{2}.y'-1)); else dout.dy = objs{1}.dy .* abs(objs{2}.y .* objs{1}.y.^(objs{2}.y-1)); end otherwise end end % Finally, apply the operator to the data if trans dout.y = feval(op, objs{1}.y, objs{2}.y'); else dout.y = feval(op, objs{1}.y, objs{2}.y); end varargout{1} = dout; end