%MIRT_DCTN N-D (multidimensional) discrete cosine transform.
% B = MIRT_DCTN(A) returns the discrete cosine transform of A.
% The array B is the same size as A and contains the
% discrete cosine transform coefficients.
%
% This function works much faster than Matlab standard
% dct (for 1D case) or dct2 (for 2D case), and also allows the ND input.
% The function takes advantage of 1) FFT 2) fast permutation through indicies
% 3) persistent precomputation (coefficients and indicies are precomputed during
% the first run).
%
% This transform can be inverted using MIRT_IDCTN.
% Author: Andriy Myronenko, see www.bme.ogi.edu/~myron
% revised on May 06, 2010
%
% See also MIRT_IDCT.