n = 1500;

W = zeros(n,n);

for j=1:n
    fprintf('processing %d of %d\n', j, n);
    ej = zeros(n,1);
    ej(j) = 1;

    W(:,j) = wavecdf97(ej,2);    
end

fprintf('making inverses and transposes:\n');
Wt = W';
Winv = inv(W);
Winvt = Winv';

fprintf('saving\n');
save('data/wavelet_mats.mat','W','Wt','Winv','Winvt');

