a = single(1);

[Achirp,AchirpOut,H,scale] = fracF_init(a);

X = complex(randn(1024,1024,'single'), ...
            randn(1024,1024,'single'));

Fdpw = fracF_dpw(X,...
                 Achirp,...
                 AchirpOut,...
                 H,...
                 scale);

Fref = complex(zeros(512,64,'single'));

for k = 1:1024
    Fref(:,k) = fracF_cg(X(:,k),a);
end

maxErr = max(abs(Fdpw(:)-Fref(:)))