Node:Real-to-Real Transforms, Next:Real-to-Real Transform Kinds, Previous:Real-data DFT Array Format, Up:Basic Interface

4.3.5 Real-to-Real Transforms

fftw_plan fftw_plan_r2r_1d(int n, double *in, double *out,
                           fftw_r2r_kind kind, unsigned flags);
fftw_plan fftw_plan_r2r_2d(int nx, int ny, double *in, double *out,
                           fftw_r2r_kind kindx, fftw_r2r_kind kindy,
                           unsigned flags);
fftw_plan fftw_plan_r2r_3d(int nx, int ny, int nz,
                           double *in, double *out,
                           fftw_r2r_kind kindx,
                           fftw_r2r_kind kindy,
                           fftw_r2r_kind kindz,
                           unsigned flags);
fftw_plan fftw_plan_r2r(int rank, const int *n, double *in, double *out,
                        const fftw_r2r_kind *kind, unsigned flags);

Plan a real input/output (r2r) transform of various kinds in zero or more dimensions, returning an fftw_plan (see Using Plans).

Once you have created a plan for a certain transform type and parameters, then creating another plan of the same type and parameters, but for different arrays, is fast and shares constant data with the first plan (if it still exists).

The planner returns NULL if the plan cannot be created. A non-NULL plan is always returned by the basic interface unless you are using a customized FFTW configuration supporting a restricted set of transforms, or for size-1 FFTW_REDFT00 kinds (which are not defined).

Arguments

• rank is the dimensionality of the transform (it should be the size of the arrays *n and *kind), and can be any non-negative integer. The _1d, _2d, and _3d planners correspond to a rank of 1, 2, and 3, respectively. A rank of zero is equivalent to a copy of one number from input to output.
• n, or nx/ny/nz, or n[rank], respectively, gives the (physical) size of the transform dimensions. They can be any positive integer.
  • Multi-dimensional arrays are stored in row-major order with dimensions: nx x ny; or nx x ny x nz; or n[0] x n[1] x ... x n[rank-1]. See Multi-dimensional Array Format.
  • FFTW is generally best at handling sizes of the form 2^a 3^b 5^c 7^d 11^e 13^f, where e+f is either 0 or 1, and the other exponents are arbitrary. Other sizes are computed by means of a slow, general-purpose algorithm (which nevertheless retains O(n log n)

    performance even for prime sizes). (It is possible to customize FFTW for different array sizes; see Installation and Customization.) Transforms whose sizes are powers of 2 are especially fast.

  • For a REDFT00 or RODFT00 transform kind in a dimension of size n, it is n-1 or n+1, respectively, that should be factorizable in the above form.
• in and out point to the input and output arrays of the transform, which may be the same (yielding an in-place transform). These arrays are overwritten during planning, unless FFTW_ESTIMATE is used in the flags. (The arrays need not be initialized, but they must be allocated.)
• kind, or kindx/kindy/kindz, or kind[rank], is the kind of r2r transform used for the corresponding dimension. The valid kind constants are described in Real-to-Real Transform Kinds. In a multi-dimensional transform, what is computed is the separable product formed by taking each transform kind along the corresponding dimension, one dimension after another.
• flags is a bitwise OR (|) of zero or more planner flags, as defined in Planner Flags.