samer@4: function A=squeeze_sum(A,dims)
samer@4: % squeeze_sum - Squeeze out first N dimensions by summing
samer@4: %
samer@4: % squeeze_sum :: [Size], K:natural -> [Size1] :- Size1=Size(K+1:end).
samer@4: %
samer@4: % Defined as squeeze_sum(A,K)=shiftdim(sumdims(A,1:K),K)
samer@4: if ~isempty(dims),
samer@4: 	for i=dims, A=sum(A,i); end
samer@4: 	A=squeeze(A);
samer@4: end