samer@3: function y=dynfilter(ab,x,z,dim)
samer@3: % dynfilter - filter for sequences of array with dynamically varying filter coeffs
samer@3: %
samer@3: % dynfilter :: seq {[[P]],[[P]]}, seq [Size], [[P]], D:natural -> seq [Size].
samer@3: % dynfilter :: seq {[[P]],[[P]]}, seq [Size], [[P]]  -> seq [Size].
samer@3: % dynfilter :: seq {[[P]],[[P]]}, seq [Size]  -> seq [Size].
samer@3: %
samer@3: % filtering is done along Dth dimension of the array (default=1)
samer@3: 
samer@3: if nargin<5, 
samer@3: 	dim=find(size(x)>1,1);
samer@3: 	if nargin<4, z=[]; end
samer@3: end
samer@3: 
samer@25: y = zipaccum(@(ab1,x1,z1)filter(ab1{1},ab1{2},x1,z1,dim),z,ab,x);