view sequences/@seq/integrate.m @ 61:eff6bddf82e3 tip

Finally implemented perceptual brightness thing.

% integrate - integrate sequence
%
% integrate ::
%    seq [[N1,N2]]    ~'sequence to integrate', 
%    D:natural        ~'dimension to integrate along',
%    [[M1,M2]]        ~'initial value',
%    ( [[M1,M2]] -> [[M1,M2]]) ~'function to filter initial value each block'
% -> seq [[N,M]].
%
% cumulative sum taken along Dth dimension of the array

function y=integrate(x,n,a,ff)
	if nargin<4, ff=@(t)t; end
	y=cumsum(map(ff,x),a,n);
end