| cannam@127 | 1 <!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" "http://www.w3.org/TR/html4/loose.dtd"> | 
| cannam@127 | 2 <html> | 
| cannam@127 | 3 <!-- This manual is for FFTW | 
| cannam@127 | 4 (version 3.3.5, 30 July 2016). | 
| cannam@127 | 5 | 
| cannam@127 | 6 Copyright (C) 2003 Matteo Frigo. | 
| cannam@127 | 7 | 
| cannam@127 | 8 Copyright (C) 2003 Massachusetts Institute of Technology. | 
| cannam@127 | 9 | 
| cannam@127 | 10 Permission is granted to make and distribute verbatim copies of this | 
| cannam@127 | 11 manual provided the copyright notice and this permission notice are | 
| cannam@127 | 12 preserved on all copies. | 
| cannam@127 | 13 | 
| cannam@127 | 14 Permission is granted to copy and distribute modified versions of this | 
| cannam@127 | 15 manual under the conditions for verbatim copying, provided that the | 
| cannam@127 | 16 entire resulting derived work is distributed under the terms of a | 
| cannam@127 | 17 permission notice identical to this one. | 
| cannam@127 | 18 | 
| cannam@127 | 19 Permission is granted to copy and distribute translations of this manual | 
| cannam@127 | 20 into another language, under the above conditions for modified versions, | 
| cannam@127 | 21 except that this permission notice may be stated in a translation | 
| cannam@127 | 22 approved by the Free Software Foundation. --> | 
| cannam@127 | 23 <!-- Created by GNU Texinfo 5.2, http://www.gnu.org/software/texinfo/ --> | 
| cannam@127 | 24 <head> | 
| cannam@127 | 25 <title>FFTW 3.3.5: Complex Multi-Dimensional DFTs</title> | 
| cannam@127 | 26 | 
| cannam@127 | 27 <meta name="description" content="FFTW 3.3.5: Complex Multi-Dimensional DFTs"> | 
| cannam@127 | 28 <meta name="keywords" content="FFTW 3.3.5: Complex Multi-Dimensional DFTs"> | 
| cannam@127 | 29 <meta name="resource-type" content="document"> | 
| cannam@127 | 30 <meta name="distribution" content="global"> | 
| cannam@127 | 31 <meta name="Generator" content="makeinfo"> | 
| cannam@127 | 32 <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> | 
| cannam@127 | 33 <link href="index.html#Top" rel="start" title="Top"> | 
| cannam@127 | 34 <link href="Concept-Index.html#Concept-Index" rel="index" title="Concept Index"> | 
| cannam@127 | 35 <link href="index.html#SEC_Contents" rel="contents" title="Table of Contents"> | 
| cannam@127 | 36 <link href="Tutorial.html#Tutorial" rel="up" title="Tutorial"> | 
| cannam@127 | 37 <link href="One_002dDimensional-DFTs-of-Real-Data.html#One_002dDimensional-DFTs-of-Real-Data" rel="next" title="One-Dimensional DFTs of Real Data"> | 
| cannam@127 | 38 <link href="Complex-One_002dDimensional-DFTs.html#Complex-One_002dDimensional-DFTs" rel="prev" title="Complex One-Dimensional DFTs"> | 
| cannam@127 | 39 <style type="text/css"> | 
| cannam@127 | 40 <!-- | 
| cannam@127 | 41 a.summary-letter {text-decoration: none} | 
| cannam@127 | 42 blockquote.smallquotation {font-size: smaller} | 
| cannam@127 | 43 div.display {margin-left: 3.2em} | 
| cannam@127 | 44 div.example {margin-left: 3.2em} | 
| cannam@127 | 45 div.indentedblock {margin-left: 3.2em} | 
| cannam@127 | 46 div.lisp {margin-left: 3.2em} | 
| cannam@127 | 47 div.smalldisplay {margin-left: 3.2em} | 
| cannam@127 | 48 div.smallexample {margin-left: 3.2em} | 
| cannam@127 | 49 div.smallindentedblock {margin-left: 3.2em; font-size: smaller} | 
| cannam@127 | 50 div.smalllisp {margin-left: 3.2em} | 
| cannam@127 | 51 kbd {font-style:oblique} | 
| cannam@127 | 52 pre.display {font-family: inherit} | 
| cannam@127 | 53 pre.format {font-family: inherit} | 
| cannam@127 | 54 pre.menu-comment {font-family: serif} | 
| cannam@127 | 55 pre.menu-preformatted {font-family: serif} | 
| cannam@127 | 56 pre.smalldisplay {font-family: inherit; font-size: smaller} | 
| cannam@127 | 57 pre.smallexample {font-size: smaller} | 
| cannam@127 | 58 pre.smallformat {font-family: inherit; font-size: smaller} | 
| cannam@127 | 59 pre.smalllisp {font-size: smaller} | 
| cannam@127 | 60 span.nocodebreak {white-space:nowrap} | 
| cannam@127 | 61 span.nolinebreak {white-space:nowrap} | 
| cannam@127 | 62 span.roman {font-family:serif; font-weight:normal} | 
| cannam@127 | 63 span.sansserif {font-family:sans-serif; font-weight:normal} | 
| cannam@127 | 64 ul.no-bullet {list-style: none} | 
| cannam@127 | 65 --> | 
| cannam@127 | 66 </style> | 
| cannam@127 | 67 | 
| cannam@127 | 68 | 
| cannam@127 | 69 </head> | 
| cannam@127 | 70 | 
| cannam@127 | 71 <body lang="en" bgcolor="#FFFFFF" text="#000000" link="#0000FF" vlink="#800080" alink="#FF0000"> | 
| cannam@127 | 72 <a name="Complex-Multi_002dDimensional-DFTs"></a> | 
| cannam@127 | 73 <div class="header"> | 
| cannam@127 | 74 <p> | 
| cannam@127 | 75 Next: <a href="One_002dDimensional-DFTs-of-Real-Data.html#One_002dDimensional-DFTs-of-Real-Data" accesskey="n" rel="next">One-Dimensional DFTs of Real Data</a>, Previous: <a href="Complex-One_002dDimensional-DFTs.html#Complex-One_002dDimensional-DFTs" accesskey="p" rel="prev">Complex One-Dimensional DFTs</a>, Up: <a href="Tutorial.html#Tutorial" accesskey="u" rel="up">Tutorial</a>   [<a href="index.html#SEC_Contents" title="Table of contents" rel="contents">Contents</a>][<a href="Concept-Index.html#Concept-Index" title="Index" rel="index">Index</a>]</p> | 
| cannam@127 | 76 </div> | 
| cannam@127 | 77 <hr> | 
| cannam@127 | 78 <a name="Complex-Multi_002dDimensional-DFTs-1"></a> | 
| cannam@127 | 79 <h3 class="section">2.2 Complex Multi-Dimensional DFTs</h3> | 
| cannam@127 | 80 | 
| cannam@127 | 81 <p>Multi-dimensional transforms work much the same way as one-dimensional | 
| cannam@127 | 82 transforms: you allocate arrays of <code>fftw_complex</code> (preferably | 
| cannam@127 | 83 using <code>fftw_malloc</code>), create an <code>fftw_plan</code>, execute it as | 
| cannam@127 | 84 many times as you want with <code>fftw_execute(plan)</code>, and clean up | 
| cannam@127 | 85 with <code>fftw_destroy_plan(plan)</code> (and <code>fftw_free</code>). | 
| cannam@127 | 86 </p> | 
| cannam@127 | 87 <p>FFTW provides two routines for creating plans for 2d and 3d transforms, | 
| cannam@127 | 88 and one routine for creating plans of arbitrary dimensionality. | 
| cannam@127 | 89 The 2d and 3d routines have the following signature: | 
| cannam@127 | 90 </p><div class="example"> | 
| cannam@127 | 91 <pre class="example">fftw_plan fftw_plan_dft_2d(int n0, int n1, | 
| cannam@127 | 92                            fftw_complex *in, fftw_complex *out, | 
| cannam@127 | 93                            int sign, unsigned flags); | 
| cannam@127 | 94 fftw_plan fftw_plan_dft_3d(int n0, int n1, int n2, | 
| cannam@127 | 95                            fftw_complex *in, fftw_complex *out, | 
| cannam@127 | 96                            int sign, unsigned flags); | 
| cannam@127 | 97 </pre></div> | 
| cannam@127 | 98 <a name="index-fftw_005fplan_005fdft_005f2d"></a> | 
| cannam@127 | 99 <a name="index-fftw_005fplan_005fdft_005f3d"></a> | 
| cannam@127 | 100 | 
| cannam@127 | 101 <p>These routines create plans for <code>n0</code> by <code>n1</code> two-dimensional | 
| cannam@127 | 102 (2d) transforms and <code>n0</code> by <code>n1</code> by <code>n2</code> 3d transforms, | 
| cannam@127 | 103 respectively.  All of these transforms operate on contiguous arrays in | 
| cannam@127 | 104 the C-standard <em>row-major</em> order, so that the last dimension has the | 
| cannam@127 | 105 fastest-varying index in the array.  This layout is described further in | 
| cannam@127 | 106 <a href="Multi_002ddimensional-Array-Format.html#Multi_002ddimensional-Array-Format">Multi-dimensional Array Format</a>. | 
| cannam@127 | 107 </p> | 
| cannam@127 | 108 <p>FFTW can also compute transforms of higher dimensionality.  In order to | 
| cannam@127 | 109 avoid confusion between the various meanings of the the word | 
| cannam@127 | 110 “dimension”, we use the term <em>rank</em> | 
| cannam@127 | 111 <a name="index-rank"></a> | 
| cannam@127 | 112 to denote the number of independent indices in an array.<a name="DOCF2" href="#FOOT2"><sup>2</sup></a>  For | 
| cannam@127 | 113 example, we say that a 2d transform has rank 2, a 3d transform has | 
| cannam@127 | 114 rank 3, and so on.  You can plan transforms of arbitrary rank by | 
| cannam@127 | 115 means of the following function: | 
| cannam@127 | 116 </p> | 
| cannam@127 | 117 <div class="example"> | 
| cannam@127 | 118 <pre class="example">fftw_plan fftw_plan_dft(int rank, const int *n, | 
| cannam@127 | 119                         fftw_complex *in, fftw_complex *out, | 
| cannam@127 | 120                         int sign, unsigned flags); | 
| cannam@127 | 121 </pre></div> | 
| cannam@127 | 122 <a name="index-fftw_005fplan_005fdft"></a> | 
| cannam@127 | 123 | 
| cannam@127 | 124 <p>Here, <code>n</code> is a pointer to an array <code>n[rank]</code> denoting an | 
| cannam@127 | 125 <code>n[0]</code> by <code>n[1]</code> by … by <code>n[rank-1]</code> transform. | 
| cannam@127 | 126 Thus, for example, the call | 
| cannam@127 | 127 </p><div class="example"> | 
| cannam@127 | 128 <pre class="example">fftw_plan_dft_2d(n0, n1, in, out, sign, flags); | 
| cannam@127 | 129 </pre></div> | 
| cannam@127 | 130 <p>is equivalent to the following code fragment: | 
| cannam@127 | 131 </p><div class="example"> | 
| cannam@127 | 132 <pre class="example">int n[2]; | 
| cannam@127 | 133 n[0] = n0; | 
| cannam@127 | 134 n[1] = n1; | 
| cannam@127 | 135 fftw_plan_dft(2, n, in, out, sign, flags); | 
| cannam@127 | 136 </pre></div> | 
| cannam@127 | 137 <p><code>fftw_plan_dft</code> is not restricted to 2d and 3d transforms, | 
| cannam@127 | 138 however, but it can plan transforms of arbitrary rank. | 
| cannam@127 | 139 </p> | 
| cannam@127 | 140 <p>You may have noticed that all the planner routines described so far | 
| cannam@127 | 141 have overlapping functionality.  For example, you can plan a 1d or 2d | 
| cannam@127 | 142 transform by using <code>fftw_plan_dft</code> with a <code>rank</code> of <code>1</code> | 
| cannam@127 | 143 or <code>2</code>, or even by calling <code>fftw_plan_dft_3d</code> with <code>n0</code> | 
| cannam@127 | 144 and/or <code>n1</code> equal to <code>1</code> (with no loss in efficiency).  This | 
| cannam@127 | 145 pattern continues, and FFTW’s planning routines in general form a | 
| cannam@127 | 146 “partial order,” sequences of | 
| cannam@127 | 147 <a name="index-partial-order"></a> | 
| cannam@127 | 148 interfaces with strictly increasing generality but correspondingly | 
| cannam@127 | 149 greater complexity. | 
| cannam@127 | 150 </p> | 
| cannam@127 | 151 <p><code>fftw_plan_dft</code> is the most general complex-DFT routine that we | 
| cannam@127 | 152 describe in this tutorial, but there are also the advanced and guru interfaces, | 
| cannam@127 | 153 <a name="index-advanced-interface-1"></a> | 
| cannam@127 | 154 <a name="index-guru-interface-1"></a> | 
| cannam@127 | 155 which allow one to efficiently combine multiple/strided transforms | 
| cannam@127 | 156 into a single FFTW plan, transform a subset of a larger | 
| cannam@127 | 157 multi-dimensional array, and/or to handle more general complex-number | 
| cannam@127 | 158 formats.  For more information, see <a href="FFTW-Reference.html#FFTW-Reference">FFTW Reference</a>. | 
| cannam@127 | 159 </p> | 
| cannam@127 | 160 <div class="footnote"> | 
| cannam@127 | 161 <hr> | 
| cannam@127 | 162 <h4 class="footnotes-heading">Footnotes</h4> | 
| cannam@127 | 163 | 
| cannam@127 | 164 <h3><a name="FOOT2" href="#DOCF2">(2)</a></h3> | 
| cannam@127 | 165 <p>The | 
| cannam@127 | 166 term “rank” is commonly used in the APL, FORTRAN, and Common Lisp | 
| cannam@127 | 167 traditions, although it is not so common in the C world.</p> | 
| cannam@127 | 168 </div> | 
| cannam@127 | 169 <hr> | 
| cannam@127 | 170 <div class="header"> | 
| cannam@127 | 171 <p> | 
| cannam@127 | 172 Next: <a href="One_002dDimensional-DFTs-of-Real-Data.html#One_002dDimensional-DFTs-of-Real-Data" accesskey="n" rel="next">One-Dimensional DFTs of Real Data</a>, Previous: <a href="Complex-One_002dDimensional-DFTs.html#Complex-One_002dDimensional-DFTs" accesskey="p" rel="prev">Complex One-Dimensional DFTs</a>, Up: <a href="Tutorial.html#Tutorial" accesskey="u" rel="up">Tutorial</a>   [<a href="index.html#SEC_Contents" title="Table of contents" rel="contents">Contents</a>][<a href="Concept-Index.html#Concept-Index" title="Index" rel="index">Index</a>]</p> | 
| cannam@127 | 173 </div> | 
| cannam@127 | 174 | 
| cannam@127 | 175 | 
| cannam@127 | 176 | 
| cannam@127 | 177 </body> | 
| cannam@127 | 178 </html> |