% fusions a block defined matrix stored inside cell array into one matrix
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% DESCRIPTION: blockMatFusion fusions a block defined matrix stored inside cell
% array into one matrix
%
% CALL: [mat3] = ssm.blockMatFusion(cell1, rowsizes, colsizes)
%
% INPUTS:
%       cell1    - cell array of matrices representing a matrix by blocs.
%                  blocs may be empty
%       rowsizes - vector giving block height
%       colsizes - vector giving block width
%
% OUTPUTS:
%       mat3 - double or symbolic array
%
% NOTE : function is private to the ssm class
%
% VERSION: '$Id: blockMatFusion.m,v 1.3 2010/08/27 13:02:45 adrien Exp $'
%
% TO DO :
% check ME in case of mixed symbolic and double
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function a_out = blockMatFusion(a, rowsizes, colsizes)
  % to deal with matrices whose size is not defined
  rowrank = cumsum([1 rowsizes]);
  colrank = cumsum([1 colsizes]);
  Nrow = length(rowsizes);
  Ncol = length(colsizes);
  isEmpty = cellfun(@isempty, a);

  a_out = zeros(rowrank(Nrow+1)-1, colrank(Ncol+1)-1);
  for ii=1:Nrow
    for jj=1:Ncol
      if ~isEmpty(ii,jj)
        rowmin = rowrank(ii);
        rowmax = rowrank(ii+1)-1;
        colmin = colrank(jj);
        colmax = colrank(jj+1)-1;
        try
          a_out(rowmin:rowmax, colmin:colmax) = a{ii,jj};
        catch ME
          a_out = sym(a_out);
          a_out(rowmin:rowmax, colmin:colmax) = sym(a{ii,jj});
        end
      end
    end
  end
end