3
<title>Guru vector and transform sizes - FFTW 3.1.2</title>
4
<meta http-equiv="Content-Type" content="text/html">
5
<meta name="description" content="FFTW 3.1.2">
6
<meta name="generator" content="makeinfo 4.8">
7
<link title="Top" rel="start" href="index.html#Top">
8
<link rel="up" href="Guru-Interface.html#Guru-Interface" title="Guru Interface">
9
<link rel="prev" href="Interleaved-and-split-arrays.html#Interleaved-and-split-arrays" title="Interleaved and split arrays">
10
<link rel="next" href="Guru-Complex-DFTs.html#Guru-Complex-DFTs" title="Guru Complex DFTs">
11
<link href="http://www.gnu.org/software/texinfo/" rel="generator-home" title="Texinfo Homepage">
13
This manual is for FFTW
14
(version 3.1.2, 23 June 2006).
16
Copyright (C) 2003 Matteo Frigo.
18
Copyright (C) 2003 Massachusetts Institute of Technology.
20
Permission is granted to make and distribute verbatim copies of
21
this manual provided the copyright notice and this permission
22
notice are preserved on all copies.
24
Permission is granted to copy and distribute modified versions of
25
this manual under the conditions for verbatim copying, provided
26
that the entire resulting derived work is distributed under the
27
terms of a permission notice identical to this one.
29
Permission is granted to copy and distribute translations of this
30
manual into another language, under the above conditions for
31
modified versions, except that this permission notice may be
32
stated in a translation approved by the Free Software Foundation.
34
<meta http-equiv="Content-Style-Type" content="text/css">
35
<style type="text/css"><!--
36
pre.display { font-family:inherit }
37
pre.format { font-family:inherit }
38
pre.smalldisplay { font-family:inherit; font-size:smaller }
39
pre.smallformat { font-family:inherit; font-size:smaller }
40
pre.smallexample { font-size:smaller }
41
pre.smalllisp { font-size:smaller }
42
span.sc { font-variant:small-caps }
43
span.roman { font-family:serif; font-weight:normal; }
44
span.sansserif { font-family:sans-serif; font-weight:normal; }
50
<a name="Guru-vector-and-transform-sizes"></a>
51
Next: <a rel="next" accesskey="n" href="Guru-Complex-DFTs.html#Guru-Complex-DFTs">Guru Complex DFTs</a>,
52
Previous: <a rel="previous" accesskey="p" href="Interleaved-and-split-arrays.html#Interleaved-and-split-arrays">Interleaved and split arrays</a>,
53
Up: <a rel="up" accesskey="u" href="Guru-Interface.html#Guru-Interface">Guru Interface</a>
57
<h4 class="subsection">4.5.2 Guru vector and transform sizes</h4>
59
<p>The guru interface introduces one basic new data structure,
60
<code>fftw_iodim</code>, that is used to specify sizes and strides for
61
multi-dimensional transforms and vectors:
63
<pre class="example"> typedef struct {
69
<p><a name="index-fftw_005fiodim-234"></a>
70
Here, <code>n</code> is the size of the dimension, and <code>is</code> and <code>os</code>
71
are the strides of that dimension for the input and output arrays. (The
72
stride is the separation of consecutive elements along this dimension.)
74
<p>The meaning of the stride parameter depends on the type of the array
75
that the stride refers to. <em>If the array is interleaved complex,
76
strides are expressed in units of complex numbers
77
(</em><code>fftw_complex</code><em>). If the array is split complex or real, strides
78
are expressed in units of real numbers (</em><code>double</code><em>).</em> This
79
convention is consistent with the usual pointer arithmetic in the C
80
language. An interleaved array is denoted by a pointer <code>p</code> to
81
<code>fftw_complex</code>, so that <code>p+1</code> points to the next complex
82
number. Split arrays are denoted by pointers to <code>double</code>, in
83
which case pointer arithmetic operates in units of
84
<code>sizeof(double)</code>.
85
<a name="index-stride-235"></a>
86
The guru planner interfaces all take a (<code>rank</code>, <code>dims[rank]</code>)
87
pair describing the transform size, and a (<code>howmany_rank</code>,
88
<code>howmany_dims[howmany_rank]</code>) pair describing the “vector” size (a
89
multi-dimensional loop of transforms to perform), where <code>dims</code> and
90
<code>howmany_dims</code> are arrays of <code>fftw_iodim</code>.
92
<p>For example, the <code>howmany</code> parameter in the advanced complex-DFT
93
interface corresponds to <code>howmany_rank</code> = 1,
94
<code>howmany_dims[0].n</code> = <code>howmany</code>, <code>howmany_dims[0].is</code> =
95
<code>idist</code>, and <code>howmany_dims[0].os</code> = <code>odist</code>.
96
<a name="index-howmany-loop-236"></a><a name="index-dist-237"></a>(To compute a single transform, you can just use <code>howmany_rank</code> = 0.)
98
<p>A row-major multidimensional array with dimensions <code>n[rank]</code>
99
(see <a href="Row_002dmajor-Format.html#Row_002dmajor-Format">Row-major Format</a>) corresponds to <code>dims[i].n</code> =
100
<code>n[i]</code> and the recurrence <code>dims[i].is</code> = <code>n[i+1] *
101
dims[i+1].is</code> (similarly for <code>os</code>). The stride of the last
102
(<code>i=rank-1</code>) dimension is the overall stride of the array.
103
e.g. to be equivalent to the advanced complex-DFT interface, you would
104
have <code>dims[rank-1].is</code> = <code>istride</code> and
105
<code>dims[rank-1].os</code> = <code>ostride</code>.
107
<p>In general, we only guarantee FFTW to return a non-<code>NULL</code> plan if
108
the vector and transform dimensions correspond to a set of distinct
109
indices, and for in-place transforms the input/output strides should