comparison src/fftw-3.3.3/doc/html/Basic-distributed_002dtranspose-interface.html @ 95:89f5e221ed7b

Add FFTW3
author Chris Cannam <cannam@all-day-breakfast.com>
date Wed, 20 Mar 2013 15:35:50 +0000
parents
children
comparison
equal deleted inserted replaced
94:d278df1123f9 95:89f5e221ed7b
1 <html lang="en">
2 <head>
3 <title>Basic distributed-transpose interface - FFTW 3.3.3</title>
4 <meta http-equiv="Content-Type" content="text/html">
5 <meta name="description" content="FFTW 3.3.3">
6 <meta name="generator" content="makeinfo 4.13">
7 <link title="Top" rel="start" href="index.html#Top">
8 <link rel="up" href="FFTW-MPI-Transposes.html#FFTW-MPI-Transposes" title="FFTW MPI Transposes">
9 <link rel="prev" href="FFTW-MPI-Transposes.html#FFTW-MPI-Transposes" title="FFTW MPI Transposes">
10 <link rel="next" href="Advanced-distributed_002dtranspose-interface.html#Advanced-distributed_002dtranspose-interface" title="Advanced distributed-transpose interface">
11 <link href="http://www.gnu.org/software/texinfo/" rel="generator-home" title="Texinfo Homepage">
12 <!--
13 This manual is for FFTW
14 (version 3.3.3, 25 November 2012).
15
16 Copyright (C) 2003 Matteo Frigo.
17
18 Copyright (C) 2003 Massachusetts Institute of Technology.
19
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.
23
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.
28
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.
33 -->
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; }
45 --></style>
46 </head>
47 <body>
48 <div class="node">
49 <a name="Basic-distributed-transpose-interface"></a>
50 <a name="Basic-distributed_002dtranspose-interface"></a>
51 <p>
52 Next:&nbsp;<a rel="next" accesskey="n" href="Advanced-distributed_002dtranspose-interface.html#Advanced-distributed_002dtranspose-interface">Advanced distributed-transpose interface</a>,
53 Previous:&nbsp;<a rel="previous" accesskey="p" href="FFTW-MPI-Transposes.html#FFTW-MPI-Transposes">FFTW MPI Transposes</a>,
54 Up:&nbsp;<a rel="up" accesskey="u" href="FFTW-MPI-Transposes.html#FFTW-MPI-Transposes">FFTW MPI Transposes</a>
55 <hr>
56 </div>
57
58 <h4 class="subsection">6.7.1 Basic distributed-transpose interface</h4>
59
60 <p>In particular, suppose that we have an <code>n0</code> by <code>n1</code> array in
61 row-major order, block-distributed across the <code>n0</code> dimension. To
62 transpose this into an <code>n1</code> by <code>n0</code> array block-distributed
63 across the <code>n1</code> dimension, we would create a plan by calling the
64 following function:
65
66 <pre class="example"> fftw_plan fftw_mpi_plan_transpose(ptrdiff_t n0, ptrdiff_t n1,
67 double *in, double *out,
68 MPI_Comm comm, unsigned flags);
69 </pre>
70 <p><a name="index-fftw_005fmpi_005fplan_005ftranspose-397"></a>
71 The input and output arrays (<code>in</code> and <code>out</code>) can be the
72 same. The transpose is actually executed by calling
73 <code>fftw_execute</code> on the plan, as usual.
74 <a name="index-fftw_005fexecute-398"></a>
75
76 <p>The <code>flags</code> are the usual FFTW planner flags, but support
77 two additional flags: <code>FFTW_MPI_TRANSPOSED_OUT</code> and/or
78 <code>FFTW_MPI_TRANSPOSED_IN</code>. What these flags indicate, for
79 transpose plans, is that the output and/or input, respectively, are
80 <em>locally</em> transposed. That is, on each process input data is
81 normally stored as a <code>local_n0</code> by <code>n1</code> array in row-major
82 order, but for an <code>FFTW_MPI_TRANSPOSED_IN</code> plan the input data is
83 stored as <code>n1</code> by <code>local_n0</code> in row-major order. Similarly,
84 <code>FFTW_MPI_TRANSPOSED_OUT</code> means that the output is <code>n0</code> by
85 <code>local_n1</code> instead of <code>local_n1</code> by <code>n0</code>.
86 <a name="index-FFTW_005fMPI_005fTRANSPOSED_005fOUT-399"></a><a name="index-FFTW_005fMPI_005fTRANSPOSED_005fIN-400"></a>
87
88 <p>To determine the local size of the array on each process before and
89 after the transpose, as well as the amount of storage that must be
90 allocated, one should call <code>fftw_mpi_local_size_2d_transposed</code>,
91 just as for a 2d DFT as described in the previous section:
92 <a name="index-data-distribution-401"></a>
93 <pre class="example"> ptrdiff_t fftw_mpi_local_size_2d_transposed
94 (ptrdiff_t n0, ptrdiff_t n1, MPI_Comm comm,
95 ptrdiff_t *local_n0, ptrdiff_t *local_0_start,
96 ptrdiff_t *local_n1, ptrdiff_t *local_1_start);
97 </pre>
98 <p><a name="index-fftw_005fmpi_005flocal_005fsize_005f2d_005ftransposed-402"></a>
99 Again, the return value is the local storage to allocate, which in
100 this case is the number of <em>real</em> (<code>double</code>) values rather
101 than complex numbers as in the previous examples.
102
103 </body></html>
104