function X = ifdct_usfft(C, isreal)

% ifdct_usfft - Inverse curvelet transform
%
% Input
%     C         Curvelet coefficients
%     isreal    Type of transform
%                   0: complex
%                   1: real
%
% Output
%     X         Image
%

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