uniquebyfunc.m

Go to the documentation of this file.

function [outUnqVec,indRightToLeftVec,indLeftToRightVec]=uniquebyfunc(inpVec,fCompare,algoName)
import mxberry.core.ismemberbyfunc;
import mxberry.core.throwerror;
%
DEFAULT_ALGO_NAME='mempreserve';
%
if nargin<3
    algoName=DEFAULT_ALGO_NAME;
    if nargin<2
        fCompare=@isequaln;
    end
end
%
switch algoName
    case 'memhungry'
        [isThereVec,indLeftToRightVec]=ismemberbyfunc(inpVec,inpVec,fCompare); %#ok<ASGLU>
        indVec=1:numel(indLeftToRightVec);
        indVec=indVec.';
        if isrow(inpVec)
            indLeftToRightVec=indLeftToRightVec.';
        end
        %
        isInUniq=ismember(indVec,indLeftToRightVec);
        indRightToLeftVec=find(isInUniq);
        [isThereVec,indLeftToRightVec]=ismember(indLeftToRightVec,...
            indRightToLeftVec); %#ok<ASGLU>
        outUnqVec=inpVec(isInUniq);
    case 'mempreserve'
        % more memory preserving
        outUnqVec=[];
        indRightToLeftVec=[];
        indLeftToRightVec=nan(numel(inpVec),1);
        for iObj=1:numel(inpVec)
            [isThereVec,indThereVec]=ismemberbyfunc(inpVec(iObj),outUnqVec,fCompare);
            if ~isThereVec
                outUnqVec=[outUnqVec; inpVec(iObj)]; %#ok<AGROW>
                indRightToLeftVec=[indRightToLeftVec; iObj]; %#ok<AGROW>
                indLeftToRightVec(iObj)=length(outUnqVec);
            else
                indLeftToRightVec(iObj)=indThereVec;
            end
        end
        if isrow(inpVec)
            outUnqVec=outUnqVec.';
        end
    otherwise
        throwerror('wrongInput','algorithm %s is not supported',algoName);
end