ismembersortableobj.m

Go to the documentation of this file.

function [isThereVec,indThereVec] = ismembersortableobj(firstVec,secVec)
import mxberry.core.uniquesortableobj;
%
if isscalar(firstVec) || isscalar(secVec)
    unitedVec = [firstVec(:);secVec(:)];
    nFirstElems = numel(firstVec);
    isThereVec = unitedVec(1:nFirstElems)==unitedVec(1+nFirstElems:end);
    if any(isThereVec)
        if isscalar(secVec)
            indThereVec = double(isThereVec);
        else
            indThereVec = find(isThereVec);
            indThereVec = indThereVec(1);
            isThereVec = any(isThereVec);
        end
    else
        isThereVec  = false(size(firstVec));
        indThereVec = zeros(size(firstVec));
    end
else
    [unqFirstVec,~,indFirstLeft2RightVec] = uniquesortableobj(firstVec(:));
    if nargout <= 1
        unqSecVec = uniquesortableobj(secVec(:));
    else
        [unqSecVec,indSecRight2LeftVec] = uniquesortableobj(secVec(:));
    end
    %
    [unionSortedVec,indUnionSortedVec] = sort([unqFirstVec;unqSecVec]);
    %
    isSortedUnionNeighborEqVec = unionSortedVec(1:end-1)==...
        unionSortedVec(2:end);
    indSortedUnionNeighborEqVec = indUnionSortedVec(...
        isSortedUnionNeighborEqVec);
    %
    if nargout <= 1
        isThereVec = ismember(indFirstLeft2RightVec,...
            indSortedUnionNeighborEqVec);
    else
        nFirstUnqElems = size(unqFirstVec,1);
        isSortedUnionNeighborEqVec = find(isSortedUnionNeighborEqVec);
        [isThereVec,indThereVec] = ismember(indFirstLeft2RightVec,...
            indSortedUnionNeighborEqVec);
        isFirstAmonthEqSortedUnionNeighborsVec =...
            isSortedUnionNeighborEqVec(indThereVec(isThereVec));
        indWhereFirstAmongEqNeighborsVec =...
            indSecRight2LeftVec(indUnionSortedVec(...
            isFirstAmonthEqSortedUnionNeighborsVec+1)-nFirstUnqElems);
        indThereVec(isThereVec) = indWhereFirstAmongEqNeighborsVec;
    end
end
isThereVec = reshape(isThereVec,size(firstVec));
if nargout > 1
    indThereVec = reshape(indThereVec,size(firstVec));
end