binaryunionstruct.m

Go to the documentation of this file.

function varargout=binaryunionstruct(leftStruct,rightStruct,fieldOpCell,leftFieldOpCell,rightFieldOpCell)
import mxberry.core.throwerror;
if nargin<2
    throwerror('incorrectInput',...
        'there should be at least 2 arguments');
end
%
sizeStructVec=size(leftStruct);
if ~(isstruct(leftStruct)&&isstruct(rightStruct)&&...
        mxberry.core.checksize(leftStruct,rightStruct,sizeStructVec))
    throwerror('incorrectInput',...
        'first and second inputs should be structure arrays of the same size');
end
%
nDims=ndims(leftStruct);
%
if nargin<=2
    fieldOpCell={@binarynotdefined};
elseif isa(fieldOpCell,'function_handle')
    fieldOpCell={fieldOpCell};
end
%

nOps=length(fieldOpCell);
%
if nargin<=4
    rightFieldOpCell=repmat({@(x)deal(x)},1,nOps);
    if nargin<=3
        leftFieldOpCell=repmat({@(x)deal(x)},1,nOps);
    end
    
end
%
if isa(leftFieldOpCell,'function_handle')
    leftFieldOpCell={leftFieldOpCell};
end
%
if isa(rightFieldOpCell,'function_handle')
    rightFieldOpCell={rightFieldOpCell};
end
%
isValid=auxcheckfuncell(fieldOpCell);
if ~isValid
    throwerror('incorrectInput',...
        ['the third argument(fieldOpCell) should be a function handle',...
        'or cell array of function handles']);
end
%
isValid=auxcheckfuncell(leftFieldOpCell);
if ~isValid
    throwerror('incorrectInput',...
        ['the fourth argument(leftFieldOpCell) should be a function ',...
        'handle or cell array of function handles']);
end
%
isValid=auxcheckfuncell(rightFieldOpCell);
if ~isValid
    throwerror('incorrectInput',...
        ['the fifth argument(rightFieldOpCell) should be a function ',...
        'handle or cell array of function handles']);
end
%
leftFieldList=fieldnames(leftStruct);
rightFieldList=fieldnames(rightStruct);
unitedFieldList=reshape(union(leftFieldList,rightFieldList),[],1);
isLeftField=ismember(unitedFieldList,leftFieldList);
isRightField=ismember(unitedFieldList,rightFieldList);
isCommonField=isLeftField&isRightField;
%
nFields=length(unitedFieldList);
nFieldOps=length(fieldOpCell);
%
nArgout=nargout;
%this cell array will contain result
resCell=cell(nArgout,nFields);
%this field is a temporary variable containing result for a specific field
resTmpCell=cell(nArgout,1);
%
for iField=1:nFields
    fieldName=unitedFieldList{iField};
    if isCommonField(iField)
        for iFieldOp=1:nFieldOps
            %
            fieldOp=fieldOpCell{iFieldOp};
            %
            [resTmpCell{:}]=cellfun(fieldOp,...
                reshape({leftStruct.(fieldName)},sizeStructVec),...
                reshape({rightStruct.(fieldName)},sizeStructVec),...
                'UniformOutput',false);
            %
            updaterescell();
            %
        end
    elseif isLeftField(iField)
        for iFieldOp=1:nFieldOps
            leftFieldOp=leftFieldOpCell{iFieldOp};
            %
            [resTmpCell{:}]=cellfun(leftFieldOp,reshape(...
                {leftStruct.(fieldName)},sizeStructVec),...
                'UniformOutput',false);
            %
            updaterescell();
            %
        end
    else
        for iFieldOp=1:nFieldOps
            rightFieldOp=rightFieldOpCell{iFieldOp};
            %
            [resTmpCell{:}]=cellfun(rightFieldOp,reshape(...
                {rightStruct.(fieldName)},sizeStructVec),...
                'UniformOutput',false);
            %
            updaterescell();
            %
        end
    end
end
%
unitedFieldList=transpose(unitedFieldList);
%
varargout=cell(1,nArgout);
for iOut=1:nArgout
    defCell=[unitedFieldList;resCell(iOut,:)];
    varargout{iOut}=struct(defCell{:});
end
    function updaterescell
        resCell(:,iField)=cellfun(@(x,y)cat(nDims+1,x,y),...
            resCell(:,iField),resTmpCell,'UniformOutput',false);
    end
end

function varargout=binarynotdefined(varargin)
varargout{:}=[];
mxberry.core.throwerror('illegalOperation',...
    ['common fields cannot be processed by default, ',...
    'binary operation required']);
end

function isValid=auxcheckfuncell(fieldOpCell)
isValid=iscell(fieldOpCell);
if ~isValid
    isAllFuncVec=cellfun(@(x)isa(x,'function_handle'),fieldOpCell);
    isValid=all(isAllFuncVec);
end
end