Chris@10: Chris@10: Chris@10: Guru Real-data DFTs - FFTW 3.3.3 Chris@10: Chris@10: Chris@10: Chris@10: Chris@10: Chris@10: Chris@10: Chris@10: Chris@10: Chris@10: Chris@10: Chris@10: Chris@10: Chris@10:
Chris@10: Chris@10: Chris@10:

Chris@10: Next: , Chris@10: Previous: Guru Complex DFTs, Chris@10: Up: Guru Interface Chris@10:


Chris@10:
Chris@10: Chris@10:

4.5.4 Guru Real-data DFTs

Chris@10: Chris@10:
     fftw_plan fftw_plan_guru_dft_r2c(
Chris@10:           int rank, const fftw_iodim *dims,
Chris@10:           int howmany_rank, const fftw_iodim *howmany_dims,
Chris@10:           double *in, fftw_complex *out,
Chris@10:           unsigned flags);
Chris@10:      
Chris@10:      fftw_plan fftw_plan_guru_split_dft_r2c(
Chris@10:           int rank, const fftw_iodim *dims,
Chris@10:           int howmany_rank, const fftw_iodim *howmany_dims,
Chris@10:           double *in, double *ro, double *io,
Chris@10:           unsigned flags);
Chris@10:      
Chris@10:      fftw_plan fftw_plan_guru_dft_c2r(
Chris@10:           int rank, const fftw_iodim *dims,
Chris@10:           int howmany_rank, const fftw_iodim *howmany_dims,
Chris@10:           fftw_complex *in, double *out,
Chris@10:           unsigned flags);
Chris@10:      
Chris@10:      fftw_plan fftw_plan_guru_split_dft_c2r(
Chris@10:           int rank, const fftw_iodim *dims,
Chris@10:           int howmany_rank, const fftw_iodim *howmany_dims,
Chris@10:           double *ri, double *ii, double *out,
Chris@10:           unsigned flags);
Chris@10: 
Chris@10:

Chris@10: Plan a real-input (r2c) or real-output (c2r), multi-dimensional DFT with Chris@10: transform dimensions given by (rank, dims) over a Chris@10: multi-dimensional vector (loop) of dimensions (howmany_rank, Chris@10: howmany_dims). dims and howmany_dims should point Chris@10: to fftw_iodim arrays of length rank and Chris@10: howmany_rank, respectively. As for the basic and advanced Chris@10: interfaces, an r2c transform is FFTW_FORWARD and a c2r transform Chris@10: is FFTW_BACKWARD. Chris@10: Chris@10:

The last dimension of dims is interpreted specially: Chris@10: that dimension of the real array has size dims[rank-1].n, but Chris@10: that dimension of the complex array has size dims[rank-1].n/2+1 Chris@10: (division rounded down). The strides, on the other hand, are taken to Chris@10: be exactly as specified. It is up to the user to specify the strides Chris@10: appropriately for the peculiar dimensions of the data, and we do not Chris@10: guarantee that the planner will succeed (return non-NULL) for Chris@10: any dimensions other than those described in Real-data DFT Array Format and generalized in Advanced Real-data DFTs. (That is, Chris@10: for an in-place transform, each individual dimension should be able to Chris@10: operate in place.) Chris@10: Chris@10: Chris@10:

in and out point to the input and output arrays for r2c Chris@10: and c2r transforms, respectively. For split arrays, ri and Chris@10: ii point to the real and imaginary input arrays for a c2r Chris@10: transform, and ro and io point to the real and imaginary Chris@10: output arrays for an r2c transform. in and ro or Chris@10: ri and out may be the same, indicating an in-place Chris@10: transform. (In-place transforms where in and io or Chris@10: ii and out are the same are not currently supported.) Chris@10: Chris@10:

flags is a bitwise OR (‘|’) of zero or more planner flags, Chris@10: as defined in Planner Flags. Chris@10: Chris@10:

In-place transforms of rank greater than 1 are currently only Chris@10: supported for interleaved arrays. For split arrays, the planner will Chris@10: return NULL. Chris@10: Chris@10: Chris@10: Chris@10: Chris@10: