decAddMonomialMatrix

% generate monomial matrices

0001 function addMonomialMatrix(dec, nMax)
0002 
0003 %% generate monomial matrices
0004 mat = cell(nMax+1, 1);
0005 mat{1} = zeros(1, dec.nC);
0006 nMonomial = 1;
0007 
0008 for iExpo = 1:nMax
0009     % obtain indices for exponent
0010     thisMat = obtainOrderedSetMatrix(dec.nC, iExpo, 1);
0011     nRow = size(thisMat, 1);
0012     nMonomial = nMonomial + nRow;
0013     
0014 %     % exponent index matrix
0015 %     mat{iExpo+1} = zeros(nRow, dec.nC);
0016 %     for iRow = 1:nRow
0017 %         mat{iExpo+1}(iRow
0018 %
0019 %     mat = [mat;
0020 end
0021 
0022 dec.monomialMatrix = zeros(nMonomial, nMax);
0023 iRow = 0;
0024 for iExpo = nMax:-1:0
0025     dec.monomialMatrix(iRow + (1:size(mat{iExpo+1}, 1)), 1:iExpo) = mat{iExpo};
0026 end