combvec.m

Go to the documentation of this file.

function indMat = combvec(varargin)
nInputs=nargin;
if nInputs==0
    indMat = [];
else
    indMat = varargin{1};
    for iInp=2:nInputs
        curVec = varargin{iInp};
        indMat = [addBlockedInd(indMat,size(curVec,2));...
            addInterleavedInd(curVec,size(indMat,2))];
    end
end
%
function indResMat = addBlockedInd(indMat,nElems)
[nDims,nCombs] = size(indMat);
indResMat = zeros(nDims,nCombs*nElems);
indVec = 1:nCombs;
for iCombToAdd=(0:(nElems-1))*nCombs
    indResMat(:,indVec+iCombToAdd) = indMat;
end
%
function indResMat = addInterleavedInd(indMat,nElems)
[nDims,nCombs] = size(indMat);
indResMat = zeros(nDims*nElems,nCombs);
ind = 1:nDims;
for iCombToAdd=(0:(nElems-1))*nDims
    indResMat(ind+iCombToAdd,:) = indMat;
end
indResMat = reshape(indResMat,nDims,nElems*nCombs);