unionjoint.m

Go to the documentation of this file.

function [unionCell,indLeft,indRight,indLocLeftVec,indLocRightVec]=unionjoint(leftCell,rightCell)
import mxberry.core.ismemberjoint;
import mxberry.core.uniquejoint;
catCell=cellfun(@smartcat,leftCell,rightCell,'UniformOutput',false);
%
unionCell=uniquejoint(catCell);
[~,indLeft]=ismemberjoint(leftCell,unionCell);
[~,indRight]=ismemberjoint(rightCell,unionCell);
nUnionElems=max(max(indLeft),max(indRight));
indLocLeftVec=zeros(1,nUnionElems);
indLocRightVec=indLocLeftVec;
%
indLocLeftVec(indLeft)=1:numel(indLeft);
indLocRightVec(indRight)=1:numel(indRight);
%

function catArray=smartcat(xArray,yArray)
import mxberry.core.throwerror;
%
isXOk=numel(xArray)==length(xArray);
isYOk=numel(yArray)==length(yArray);
%
if ~(isXOk||isYOk)
    throwerror('incorrectInput',...
        ['corresponding cell items in input cell arrays should ',...
        'be both either rows or colums']);
end
isXRow=size(xArray,1)==1;
isYRow=size(xArray,1)==1;
catDim=max(isXRow,isYRow)+1;
catArray=cat(catDim,xArray,yArray);