function [af, sf] = doubledualfilt

%% USAGE:
%   [af, sf] = doubledualfilt
% OUTPUT:
%   af - cell array containing analysis filters
%       af{1} - analysis filters for tree 1
%       af{2} - analysis filters for tree 2
%   sf - cell array containing synthesis filters
%       sf{1} - synthesis filters for tree 1
%       sf{2} - synthesis filters for tree 2
% See also: afb3, doubledualtree_f1D
% 
% FRAMELET SOFTWARE AT POLYTECHNIC UNIVERSITY, BROOKLYN, NY
% 

% ANALYSIS FILTERS
af{1} = [
   0.00017870679071  -0.00012344587034   0.00001437252392
  -0.02488304194507   0.01718853971559  -0.00200122286479
   0.00737700819766  -0.00675291099550   0.00027261232228
   0.29533805776119   0.02671809818132   0.06840460220387
   0.59529279993637  -0.64763513288874   0.01936710587994
   0.45630440337458   0.47089724990858  -0.68031992557818
   0.11239376309619   0.16040017815754   0.42976785708978
  -0.01971220693439  -0.01484700537727   0.11428688385011
  -0.00813549683439  -0.00588868840296   0.05057805218407
   0.00005956893024   0.00004311757177  -0.00037033761102
];

af{2} = [
   0.00005956893024   0.00004311757177  -0.00037033761102
  -0.00813549683439  -0.00588868840296   0.05057805218407
  -0.01971220693439  -0.01484700537727   0.11428688385011
   0.11239376309619   0.16040017815754   0.42976785708978
   0.45630440337458   0.47089724990858  -0.68031992557818
   0.59529279993637  -0.64763513288874   0.01936710587994
   0.29533805776119   0.02671809818132   0.06840460220387
   0.00737700819766  -0.00675291099550   0.00027261232228
  -0.02488304194507   0.01718853971559  -0.00200122286479
   0.00017870679071  -0.00012344587034   0.00001437252392
];

% Swap the two highpass filters for af{2}
tmp = af{2}(:,2);
af{2}(:,2) = af{2}(:,3);
af{2}(:,3) = tmp;

% SYNTHESIS FILTERS
sf{1} = af{1}(end:-1:1, :);
sf{2} = af{2}(end:-1:1, :);