function X = afdct_usfft(C,isreal)

% afdct_usfft - Adjoint curvelet transform
%
% Input
%     C         Curvelet coefficients
%     isreal    Type of transform
%                   0: complex
%                   1: real
%
% Output
%     X         A double precision matrix
%

  [m,n] = size(C{end}{1});
  nbscales = floor(log2(min(m,n)))-3;
  nbangles_coarse = 16;
  allcurvelets = 0;
  
  if(isreal)
    C = fdct_usfft_c2r(C);
  end
  
  % call mex function
  X = afdct_usfft_mex(m,n,nbscales, nbangles_coarse, allcurvelets, C);