function ldct = FastLDCTAnalysis(x,w,overlap,pars3)
% FastLDCTAnalysis -- Local transform using a DCT dictionary for 1D signals
%  Usage
%    ldct = FastLDCTAnalysis(x,w,overlap) 
%  Inputs
%    x        	1-d signal:  length(x)=2^J
%    overlap    blocks overlapping, fraction of the window width (>=0 & <1), 
%               eg: 0.75, 0.5, 0.25, 0.125
%    w        	width of window (must be a power of two)
%  Outputs
%    ldct    	Local DCT coefficients
%  Description
%    The ldct contains coefficients of the Local DCT Decomposition.
% See Also
%   FastLDCTSynthesis, dct, idct
%		

    if overlap < 0 | overlap >0.5
        error('The blocks overlapping must be >=0 and <=0.5');
    end
    
	[n,J] = dyadlength(x);
	
	d = floor(n/w);
	m = floor(w*overlap);
	
	ldct = zeros(w+2*m,d);

        xtmp  = [zeros(m,1);ShapeAsRow(x)';zeros(m,1)];
%
	for p=0:d-1
        	xp = xtmp(p*w+1:(p+1)*w+2*m);
		c = dct(xp);                % DCT analysis
		ldct(:,p+1) = c;   	    % store
	end
    
	ldct = ldct(:);