function varargout = adjsiza(varargin)% DESCRIPTION [a,b,...] = adjsiza(a,b,...) %  Will extend singleton dimensions on each input argument to match%  the size of the other arguments. If dimensions are not singletons%  nor of the same size an error will occur. It is also an error%  if any argument is empty. The result is matrices of the same size.% INPUT%  any number of matrices% OUTPUT%  one matrix for each input all of the same size% TRY %  [a, b] = adjsiza(ones(1,1,2,2), ones(1,2,2,1)) => size(a)==size(b)% SEE ALSO%  adjsiz, mplus, mprod% by Magnus Almgren 990714% Take chance on size otherwise a reallocation will take place.dim = zeros(nargin,4); % the dimensions are hopefully not more than 4 for j = 1:nargin siz = size(varargin{j})-1; dim(j,1:length(siz)) = siz; % store row by rowenddim = dim+1;% Make a matrix out of max size vector.mdim = repmat(max(dim,[],1), [size(dim,1) 1]); % error if zero or (not 1 and not max)errv = (dim==0) | ((dim~=1) & (dim~=mdim));if any(errv(:)) error('Input arguments mismatch or empty')endrepdim = mdim./dim; % values to expand inputs with  for j = 1:max(1,nargout)  % go through all arguments  varargout{j} = repmat(varargin{j}, repdim(j,:));end