function W = ALLHWVF(A,k,t,show)
% ALLHWVF displays sum of 2-dim Haar wavelets of frequency k type t
%
% use for amplitudes (wavelet transform / coefficients) of 2-D-Haar
% wavelets as coded in fasthaaar2.
%
%
%
%   Matt Kawski 8-8-2008
%   http://math.asu.edu/~kawski
%   kawski@asu.edu
%
if nargin < 1; A=fasthaar2(0); end
if nargin < 2; k=3; end
if nargin < 3; t=0; end
if nargin < 4; show=1; else show=0; end

n=log(size(A,1))/log(2);
colormap(gray)

W = zeros(size(A));
if t==0;
    for t=1:4
        W = W+allhwvf(A,k,t,show);
        if show; image(32*(1+W)), pause(1); end
    end
    else
    for i=0:2^(n-k)-1,
        for j=0:2^(n-k)-1
            W=W+A(2^k*i+1,2^k*j+1)*hwv(n,k,i,j,t);
            if show; image(32*(1+W)), pause(0.1); end
        end
    end
end
if show; image(32*(1+W)),  end
end