view general/numerical/array/squeeze_prod.m @ 61:eff6bddf82e3 tip

Finally implemented perceptual brightness thing.

function A=squeeze_prod(A,dims)
% squeeze_prod - Squeeze out first N dimensions by multiplying
%
% squeeze_prod :: [Size], K:natural -> [Size1] :- Size1=Size(K+1:end).
%
% Defined as squeeze_prod(A,K)=shiftdim(proddims(A,1:K),K)
if ~isempty(dims)
	for i=dims, A=prod(A,i); end
	A=squeeze(A);
end