Added all the F4 libraries to the project
[fw/stlink] / exampleF4 / CMSIS / Documentation / DSP_Lib / html / arm__cfft__radix4__f32_8c_source.html
1 <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
2 <html xmlns="http://www.w3.org/1999/xhtml">
3 <head>
4 <meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
5 <title>CMSIS DSP Software Library: arm_cfft_radix4_f32.c Source File</title>
6 <link href="tabs.css" rel="stylesheet" type="text/css"/>
7 <link href="search/search.css" rel="stylesheet" type="text/css"/>
8 <script type="text/javaScript" src="search/search.js"></script>
9 <link href="doxygen.css" rel="stylesheet" type="text/css"/>
10 </head>
11 <body onload='searchBox.OnSelectItem(0);'>
12 <!-- Generated by Doxygen 1.7.2 -->
13 <script type="text/javascript"><!--
14 var searchBox = new SearchBox("searchBox", "search",false,'Search');
15 --></script>
16 <div class="navigation" id="top">
17   <div class="tabs">
18     <ul class="tablist">
19       <li><a href="index.html"><span>Main&#160;Page</span></a></li>
20       <li><a href="modules.html"><span>Modules</span></a></li>
21       <li><a href="annotated.html"><span>Data&#160;Structures</span></a></li>
22       <li class="current"><a href="files.html"><span>Files</span></a></li>
23       <li><a href="examples.html"><span>Examples</span></a></li>
24       <li id="searchli">
25         <div id="MSearchBox" class="MSearchBoxInactive">
26         <span class="left">
27           <img id="MSearchSelect" src="search/mag_sel.png"
28                onmouseover="return searchBox.OnSearchSelectShow()"
29                onmouseout="return searchBox.OnSearchSelectHide()"
30                alt=""/>
31           <input type="text" id="MSearchField" value="Search" accesskey="S"
32                onfocus="searchBox.OnSearchFieldFocus(true)" 
33                onblur="searchBox.OnSearchFieldFocus(false)" 
34                onkeyup="searchBox.OnSearchFieldChange(event)"/>
35           </span><span class="right">
36             <a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
37           </span>
38         </div>
39       </li>
40     </ul>
41   </div>
42   <div class="tabs2">
43     <ul class="tablist">
44       <li><a href="files.html"><span>File&#160;List</span></a></li>
45       <li><a href="globals.html"><span>Globals</span></a></li>
46     </ul>
47   </div>
48 <div class="header">
49   <div class="headertitle">
50 <h1>arm_cfft_radix4_f32.c</h1>  </div>
51 </div>
52 <div class="contents">
53 <a href="arm__cfft__radix4__f32_8c.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/* ----------------------------------------------------------------------   </span>
54 <a name="l00002"></a>00002 <span class="comment">* Copyright (C) 2010 ARM Limited. All rights reserved.   </span>
55 <a name="l00003"></a>00003 <span class="comment">*   </span>
56 <a name="l00004"></a>00004 <span class="comment">* $Date:        15. July 2011  </span>
57 <a name="l00005"></a>00005 <span class="comment">* $Revision:    V1.0.10  </span>
58 <a name="l00006"></a>00006 <span class="comment">*   </span>
59 <a name="l00007"></a>00007 <span class="comment">* Project:      CMSIS DSP Library   </span>
60 <a name="l00008"></a>00008 <span class="comment">* Title:        arm_cfft_radix4_f32.c   </span>
61 <a name="l00009"></a>00009 <span class="comment">*   </span>
62 <a name="l00010"></a>00010 <span class="comment">* Description:  Radix-4 Decimation in Frequency CFFT &amp; CIFFT Floating point processing function   </span>
63 <a name="l00011"></a>00011 <span class="comment">*   </span>
64 <a name="l00012"></a>00012 <span class="comment">*   </span>
65 <a name="l00013"></a>00013 <span class="comment">* Target Processor: Cortex-M4/Cortex-M3/Cortex-M0</span>
66 <a name="l00014"></a>00014 <span class="comment">*  </span>
67 <a name="l00015"></a>00015 <span class="comment">* Version 1.0.10 2011/7/15 </span>
68 <a name="l00016"></a>00016 <span class="comment">*    Big Endian support added and Merged M0 and M3/M4 Source code.  </span>
69 <a name="l00017"></a>00017 <span class="comment">*   </span>
70 <a name="l00018"></a>00018 <span class="comment">* Version 1.0.3 2010/11/29  </span>
71 <a name="l00019"></a>00019 <span class="comment">*    Re-organized the CMSIS folders and updated documentation.   </span>
72 <a name="l00020"></a>00020 <span class="comment">*    </span>
73 <a name="l00021"></a>00021 <span class="comment">* Version 1.0.2 2010/11/11   </span>
74 <a name="l00022"></a>00022 <span class="comment">*    Documentation updated.    </span>
75 <a name="l00023"></a>00023 <span class="comment">*   </span>
76 <a name="l00024"></a>00024 <span class="comment">* Version 1.0.1 2010/10/05    </span>
77 <a name="l00025"></a>00025 <span class="comment">*    Production release and review comments incorporated.   </span>
78 <a name="l00026"></a>00026 <span class="comment">*   </span>
79 <a name="l00027"></a>00027 <span class="comment">* Version 1.0.0 2010/09/20    </span>
80 <a name="l00028"></a>00028 <span class="comment">*    Production release and review comments incorporated.   </span>
81 <a name="l00029"></a>00029 <span class="comment">*   </span>
82 <a name="l00030"></a>00030 <span class="comment">* Version 0.0.5  2010/04/26    </span>
83 <a name="l00031"></a>00031 <span class="comment">*    incorporated review comments and updated with latest CMSIS layer   </span>
84 <a name="l00032"></a>00032 <span class="comment">*   </span>
85 <a name="l00033"></a>00033 <span class="comment">* Version 0.0.3  2010/03/10    </span>
86 <a name="l00034"></a>00034 <span class="comment">*    Initial version   </span>
87 <a name="l00035"></a>00035 <span class="comment">* -------------------------------------------------------------------- */</span>
88 <a name="l00036"></a>00036 
89 <a name="l00037"></a>00037 <span class="preprocessor">#include &quot;<a class="code" href="arm__math_8h.html">arm_math.h</a>&quot;</span>
90 <a name="l00038"></a>00038 
91 <a name="l00177"></a><a class="code" href="group___c_f_f_t___c_i_f_f_t.html#ga521f670cd9c571bc61aff9bec89f4c26">00177</a> <span class="keywordtype">void</span> <a class="code" href="group___c_f_f_t___c_i_f_f_t.html#ga521f670cd9c571bc61aff9bec89f4c26" title="Processing function for the floating-point CFFT/CIFFT.">arm_cfft_radix4_f32</a>(
92 <a name="l00178"></a>00178   <span class="keyword">const</span> <a class="code" href="structarm__cfft__radix4__instance__f32.html" title="Instance structure for the floating-point CFFT/CIFFT function.">arm_cfft_radix4_instance_f32</a> * S,
93 <a name="l00179"></a>00179   <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> * pSrc)
94 <a name="l00180"></a>00180 {
95 <a name="l00181"></a>00181 
96 <a name="l00182"></a>00182   <span class="keywordflow">if</span>(S-&gt;<a class="code" href="structarm__cfft__radix4__instance__f32.html#a25d1da64dd6487c291f04d226f9acc66">ifftFlag</a> == 1u)
97 <a name="l00183"></a>00183   {
98 <a name="l00184"></a>00184     <span class="comment">/*  Complex IFFT radix-4  */</span>
99 <a name="l00185"></a>00185     <a class="code" href="arm__cfft__radix4__f32_8c.html#a2a78df6e4bbf080624f2b6349224ec93" title="Core function for the floating-point CIFFT butterfly process.">arm_radix4_butterfly_inverse_f32</a>(pSrc, S-&gt;<a class="code" href="structarm__cfft__radix4__instance__f32.html#a7e6a6d290ce158ce9a15a45e364b021a">fftLen</a>, S-&gt;<a class="code" href="structarm__cfft__radix4__instance__f32.html#a14860c7544911702ca1fa0bf78204ef3">pTwiddle</a>,
100 <a name="l00186"></a>00186                                      S-&gt;<a class="code" href="structarm__cfft__radix4__instance__f32.html#abe31ea2157dfa233e389cdfd3b9993ee">twidCoefModifier</a>, S-&gt;<a class="code" href="structarm__cfft__radix4__instance__f32.html#ab9eed39e40b8d7c16381fbccf84467cd">onebyfftLen</a>);
101 <a name="l00187"></a>00187   }
102 <a name="l00188"></a>00188   <span class="keywordflow">else</span>
103 <a name="l00189"></a>00189   {
104 <a name="l00190"></a>00190     <span class="comment">/*  Complex FFT radix-4  */</span>
105 <a name="l00191"></a>00191     <a class="code" href="arm__cfft__radix4__f32_8c.html#ae239ddf995d1607115f9e84d5c069b9c" title="Core function for the floating-point CFFT butterfly process.">arm_radix4_butterfly_f32</a>(pSrc, S-&gt;<a class="code" href="structarm__cfft__radix4__instance__f32.html#a7e6a6d290ce158ce9a15a45e364b021a">fftLen</a>, S-&gt;<a class="code" href="structarm__cfft__radix4__instance__f32.html#a14860c7544911702ca1fa0bf78204ef3">pTwiddle</a>,
106 <a name="l00192"></a>00192                              S-&gt;<a class="code" href="structarm__cfft__radix4__instance__f32.html#abe31ea2157dfa233e389cdfd3b9993ee">twidCoefModifier</a>);
107 <a name="l00193"></a>00193   }
108 <a name="l00194"></a>00194 
109 <a name="l00195"></a>00195   <span class="keywordflow">if</span>(S-&gt;<a class="code" href="structarm__cfft__radix4__instance__f32.html#ac10927a1620195a88649ce63dab66120">bitReverseFlag</a> == 1u)
110 <a name="l00196"></a>00196   {
111 <a name="l00197"></a>00197     <span class="comment">/*  Bit Reversal */</span>
112 <a name="l00198"></a>00198     <a class="code" href="arm__cfft__radix4__f32_8c.html#a3d4062fdfa6aaa3f51f41cab868e508b" title="In-place bit reversal function.">arm_bitreversal_f32</a>(pSrc, S-&gt;<a class="code" href="structarm__cfft__radix4__instance__f32.html#a7e6a6d290ce158ce9a15a45e364b021a">fftLen</a>, S-&gt;<a class="code" href="structarm__cfft__radix4__instance__f32.html#acc8cb18a8b901b8321ab9d86491e41a3">bitRevFactor</a>, S-&gt;<a class="code" href="structarm__cfft__radix4__instance__f32.html#a8da0d2ca69749fde8cbb95caeac6fe6a">pBitRevTable</a>);
113 <a name="l00199"></a>00199   }
114 <a name="l00200"></a>00200 
115 <a name="l00201"></a>00201 }
116 <a name="l00202"></a>00202 
117 <a name="l00203"></a>00203 
118 <a name="l00210"></a>00210 <span class="comment">/* ----------------------------------------------------------------------   </span>
119 <a name="l00211"></a>00211 <span class="comment">** Internal helper function used by the FFTs   </span>
120 <a name="l00212"></a>00212 <span class="comment">** ------------------------------------------------------------------- */</span>
121 <a name="l00213"></a>00213 
122 <a name="l00214"></a>00214 <span class="comment">/*   </span>
123 <a name="l00215"></a>00215 <span class="comment"> * @brief  Core function for the floating-point CFFT butterfly process.  </span>
124 <a name="l00216"></a>00216 <span class="comment"> * @param[in, out] *pSrc            points to the in-place buffer of floating-point data type.  </span>
125 <a name="l00217"></a>00217 <span class="comment"> * @param[in]      fftLen           length of the FFT.  </span>
126 <a name="l00218"></a>00218 <span class="comment"> * @param[in]      *pCoef           points to the twiddle coefficient buffer.  </span>
127 <a name="l00219"></a>00219 <span class="comment"> * @param[in]      twidCoefModifier twiddle coefficient modifier that supports different size FFTs with the same twiddle factor table.  </span>
128 <a name="l00220"></a>00220 <span class="comment"> * @return none.  </span>
129 <a name="l00221"></a>00221 <span class="comment"> */</span>
130 <a name="l00222"></a>00222 
131 <a name="l00223"></a><a class="code" href="arm__math_8h.html#ae239ddf995d1607115f9e84d5c069b9c">00223</a> <span class="keywordtype">void</span> <a class="code" href="arm__cfft__radix4__f32_8c.html#ae239ddf995d1607115f9e84d5c069b9c" title="Core function for the floating-point CFFT butterfly process.">arm_radix4_butterfly_f32</a>(
132 <a name="l00224"></a>00224   <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> * pSrc,
133 <a name="l00225"></a>00225   uint16_t fftLen,
134 <a name="l00226"></a>00226   <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> * pCoef,
135 <a name="l00227"></a>00227   uint16_t twidCoefModifier)
136 <a name="l00228"></a>00228 {
137 <a name="l00229"></a>00229 
138 <a name="l00230"></a>00230   <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> co1, co2, co3, si1, si2, si3;
139 <a name="l00231"></a>00231   <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> t1, t2, r1, r2, s1, s2;
140 <a name="l00232"></a>00232   uint32_t ia1, ia2, ia3;
141 <a name="l00233"></a>00233   uint32_t i0, i1, i2, i3;
142 <a name="l00234"></a>00234   uint32_t n1, n2, j, k;
143 <a name="l00235"></a>00235 
144 <a name="l00236"></a>00236 <span class="preprocessor">#ifndef ARM_MATH_CM0</span>
145 <a name="l00237"></a>00237 <span class="preprocessor"></span>
146 <a name="l00238"></a>00238   <span class="comment">/* Run the below code for Cortex-M4 and Cortex-M3 */</span>
147 <a name="l00239"></a>00239 
148 <a name="l00240"></a>00240   <span class="comment">/*  Initializations for the first stage */</span>
149 <a name="l00241"></a>00241   n2 = fftLen;
150 <a name="l00242"></a>00242   n1 = n2;
151 <a name="l00243"></a>00243 
152 <a name="l00244"></a>00244   <span class="comment">/* n2 = fftLen/4 */</span>
153 <a name="l00245"></a>00245   n2 &gt;&gt;= 2u;
154 <a name="l00246"></a>00246   i0 = 0u;
155 <a name="l00247"></a>00247   ia1 = 0u;
156 <a name="l00248"></a>00248 
157 <a name="l00249"></a>00249   j = n2;
158 <a name="l00250"></a>00250 
159 <a name="l00251"></a>00251   <span class="comment">/*  Calculation of first stage */</span>
160 <a name="l00252"></a>00252   <span class="keywordflow">do</span>
161 <a name="l00253"></a>00253   {
162 <a name="l00254"></a>00254     <span class="comment">/*  index calculation for the input as, */</span>
163 <a name="l00255"></a>00255     <span class="comment">/*  pSrc[i0 + 0], pSrc[i0 + fftLen/4], pSrc[i0 + fftLen/2], pSrc[i0 + 3fftLen/4] */</span>
164 <a name="l00256"></a>00256     i1 = i0 + n2;
165 <a name="l00257"></a>00257     i2 = i1 + n2;
166 <a name="l00258"></a>00258     i3 = i2 + n2;
167 <a name="l00259"></a>00259 
168 <a name="l00260"></a>00260     <span class="comment">/*  Butterfly implementation */</span>
169 <a name="l00261"></a>00261 
170 <a name="l00262"></a>00262     <span class="comment">/* xa + xc */</span>
171 <a name="l00263"></a>00263     r1 = pSrc[(2u * i0)] + pSrc[(2u * i2)];
172 <a name="l00264"></a>00264 
173 <a name="l00265"></a>00265     <span class="comment">/* xa - xc */</span>
174 <a name="l00266"></a>00266     r2 = pSrc[2u * i0] - pSrc[2u * i2];
175 <a name="l00267"></a>00267 
176 <a name="l00268"></a>00268     <span class="comment">/* ya + yc */</span>
177 <a name="l00269"></a>00269     s1 = pSrc[(2u * i0) + 1u] + pSrc[(2u * i2) + 1u];
178 <a name="l00270"></a>00270 
179 <a name="l00271"></a>00271     <span class="comment">/* ya - yc */</span>
180 <a name="l00272"></a>00272     s2 = pSrc[(2u * i0) + 1u] - pSrc[(2u * i2) + 1u];
181 <a name="l00273"></a>00273 
182 <a name="l00274"></a>00274     <span class="comment">/* xb + xd */</span>
183 <a name="l00275"></a>00275     t1 = pSrc[2u * i1] + pSrc[2u * i3];
184 <a name="l00276"></a>00276 
185 <a name="l00277"></a>00277     <span class="comment">/* xa&#39; = xa + xb + xc + xd */</span>
186 <a name="l00278"></a>00278     pSrc[2u * i0] = r1 + t1;
187 <a name="l00279"></a>00279 
188 <a name="l00280"></a>00280     <span class="comment">/* (xa + xc) - (xb + xd) */</span>
189 <a name="l00281"></a>00281     r1 = r1 - t1;
190 <a name="l00282"></a>00282 
191 <a name="l00283"></a>00283     <span class="comment">/* yb + yd */</span>
192 <a name="l00284"></a>00284     t2 = pSrc[(2u * i1) + 1u] + pSrc[(2u * i3) + 1u];
193 <a name="l00285"></a>00285 
194 <a name="l00286"></a>00286     <span class="comment">/* ya&#39; = ya + yb + yc + yd */</span>
195 <a name="l00287"></a>00287     pSrc[(2u * i0) + 1u] = s1 + t2;
196 <a name="l00288"></a>00288 
197 <a name="l00289"></a>00289     <span class="comment">/* (ya + yc) - (yb + yd) */</span>
198 <a name="l00290"></a>00290     s1 = s1 - t2;
199 <a name="l00291"></a>00291 
200 <a name="l00292"></a>00292     <span class="comment">/* yb - yd */</span>
201 <a name="l00293"></a>00293     t1 = pSrc[(2u * i1) + 1u] - pSrc[(2u * i3) + 1u];
202 <a name="l00294"></a>00294 
203 <a name="l00295"></a>00295     <span class="comment">/* xb - xd */</span>
204 <a name="l00296"></a>00296     t2 = pSrc[2u * i1] - pSrc[2u * i3];
205 <a name="l00297"></a>00297 
206 <a name="l00298"></a>00298     <span class="comment">/*  index calculation for the coefficients */</span>
207 <a name="l00299"></a>00299     ia2 = ia1 + ia1;
208 <a name="l00300"></a>00300     co2 = pCoef[ia2 * 2u];
209 <a name="l00301"></a>00301     si2 = pCoef[(ia2 * 2u) + 1u];
210 <a name="l00302"></a>00302 
211 <a name="l00303"></a>00303     <span class="comment">/* xc&#39; = (xa-xb+xc-xd)co2 + (ya-yb+yc-yd)(si2) */</span>
212 <a name="l00304"></a>00304     pSrc[2u * i1] = (r1 * co2) + (s1 * si2);
213 <a name="l00305"></a>00305 
214 <a name="l00306"></a>00306     <span class="comment">/* yc&#39; = (ya-yb+yc-yd)co2 - (xa-xb+xc-xd)(si2) */</span>
215 <a name="l00307"></a>00307     pSrc[(2u * i1) + 1u] = (s1 * co2) - (r1 * si2);
216 <a name="l00308"></a>00308 
217 <a name="l00309"></a>00309     <span class="comment">/* (xa - xc) + (yb - yd) */</span>
218 <a name="l00310"></a>00310     r1 = r2 + t1;
219 <a name="l00311"></a>00311 
220 <a name="l00312"></a>00312     <span class="comment">/* (xa - xc) - (yb - yd) */</span>
221 <a name="l00313"></a>00313     r2 = r2 - t1;
222 <a name="l00314"></a>00314 
223 <a name="l00315"></a>00315     <span class="comment">/* (ya - yc) - (xb - xd) */</span>
224 <a name="l00316"></a>00316     s1 = s2 - t2;
225 <a name="l00317"></a>00317 
226 <a name="l00318"></a>00318     <span class="comment">/* (ya - yc) + (xb - xd) */</span>
227 <a name="l00319"></a>00319     s2 = s2 + t2;
228 <a name="l00320"></a>00320 
229 <a name="l00321"></a>00321     co1 = pCoef[ia1 * 2u];
230 <a name="l00322"></a>00322     si1 = pCoef[(ia1 * 2u) + 1u];
231 <a name="l00323"></a>00323 
232 <a name="l00324"></a>00324     <span class="comment">/* xb&#39; = (xa+yb-xc-yd)co1 + (ya-xb-yc+xd)(si1) */</span>
233 <a name="l00325"></a>00325     pSrc[2u * i2] = (r1 * co1) + (s1 * si1);
234 <a name="l00326"></a>00326 
235 <a name="l00327"></a>00327     <span class="comment">/* yb&#39; = (ya-xb-yc+xd)co1 - (xa+yb-xc-yd)(si1) */</span>
236 <a name="l00328"></a>00328     pSrc[(2u * i2) + 1u] = (s1 * co1) - (r1 * si1);
237 <a name="l00329"></a>00329 
238 <a name="l00330"></a>00330     <span class="comment">/*  index calculation for the coefficients */</span>
239 <a name="l00331"></a>00331     ia3 = ia2 + ia1;
240 <a name="l00332"></a>00332     co3 = pCoef[ia3 * 2u];
241 <a name="l00333"></a>00333     si3 = pCoef[(ia3 * 2u) + 1u];
242 <a name="l00334"></a>00334 
243 <a name="l00335"></a>00335 
244 <a name="l00336"></a>00336     <span class="comment">/* xd&#39; = (xa-yb-xc+yd)co3 + (ya+xb-yc-xd)(si3) */</span>
245 <a name="l00337"></a>00337     pSrc[2u * i3] = (r2 * co3) + (s2 * si3);
246 <a name="l00338"></a>00338 
247 <a name="l00339"></a>00339     <span class="comment">/* yd&#39; = (ya+xb-yc-xd)co3 - (xa-yb-xc+yd)(si3) */</span>
248 <a name="l00340"></a>00340     pSrc[(2u * i3) + 1u] = (s2 * co3) - (r2 * si3);
249 <a name="l00341"></a>00341 
250 <a name="l00342"></a>00342     <span class="comment">/*  Twiddle coefficients index modifier */</span>
251 <a name="l00343"></a>00343     ia1 = ia1 + twidCoefModifier;
252 <a name="l00344"></a>00344 
253 <a name="l00345"></a>00345     <span class="comment">/*  Updating input index */</span>
254 <a name="l00346"></a>00346     i0 = i0 + 1u;
255 <a name="l00347"></a>00347 
256 <a name="l00348"></a>00348   }
257 <a name="l00349"></a>00349   <span class="keywordflow">while</span>(--j);
258 <a name="l00350"></a>00350 
259 <a name="l00351"></a>00351   twidCoefModifier &lt;&lt;= 2u;
260 <a name="l00352"></a>00352 
261 <a name="l00353"></a>00353   <span class="comment">/*  Calculation of second stage to excluding last stage */</span>
262 <a name="l00354"></a>00354   <span class="keywordflow">for</span> (k = fftLen / 4; k &gt; 4u; k &gt;&gt;= 2u)
263 <a name="l00355"></a>00355   {
264 <a name="l00356"></a>00356     <span class="comment">/*  Initializations for the first stage */</span>
265 <a name="l00357"></a>00357     n1 = n2;
266 <a name="l00358"></a>00358     n2 &gt;&gt;= 2u;
267 <a name="l00359"></a>00359     ia1 = 0u;
268 <a name="l00360"></a>00360 
269 <a name="l00361"></a>00361     <span class="comment">/*  Calculation of first stage */</span>
270 <a name="l00362"></a>00362     <span class="keywordflow">for</span> (j = 0u; j &lt;= (n2 - 1u); j++)
271 <a name="l00363"></a>00363     {
272 <a name="l00364"></a>00364       <span class="comment">/*  index calculation for the coefficients */</span>
273 <a name="l00365"></a>00365       ia2 = ia1 + ia1;
274 <a name="l00366"></a>00366       ia3 = ia2 + ia1;
275 <a name="l00367"></a>00367       co1 = pCoef[ia1 * 2u];
276 <a name="l00368"></a>00368       si1 = pCoef[(ia1 * 2u) + 1u];
277 <a name="l00369"></a>00369       co2 = pCoef[ia2 * 2u];
278 <a name="l00370"></a>00370       si2 = pCoef[(ia2 * 2u) + 1u];
279 <a name="l00371"></a>00371       co3 = pCoef[ia3 * 2u];
280 <a name="l00372"></a>00372       si3 = pCoef[(ia3 * 2u) + 1u];
281 <a name="l00373"></a>00373 
282 <a name="l00374"></a>00374       <span class="comment">/*  Twiddle coefficients index modifier */</span>
283 <a name="l00375"></a>00375       ia1 = ia1 + twidCoefModifier;
284 <a name="l00376"></a>00376 
285 <a name="l00377"></a>00377       <span class="keywordflow">for</span> (i0 = j; i0 &lt; fftLen; i0 += n1)
286 <a name="l00378"></a>00378       {
287 <a name="l00379"></a>00379         <span class="comment">/*  index calculation for the input as, */</span>
288 <a name="l00380"></a>00380         <span class="comment">/*  pSrc[i0 + 0], pSrc[i0 + fftLen/4], pSrc[i0 + fftLen/2], pSrc[i0 + 3fftLen/4] */</span>
289 <a name="l00381"></a>00381         i1 = i0 + n2;
290 <a name="l00382"></a>00382         i2 = i1 + n2;
291 <a name="l00383"></a>00383         i3 = i2 + n2;
292 <a name="l00384"></a>00384 
293 <a name="l00385"></a>00385         <span class="comment">/* xa + xc */</span>
294 <a name="l00386"></a>00386         r1 = pSrc[(2u * i0)] + pSrc[(2u * i2)];
295 <a name="l00387"></a>00387 
296 <a name="l00388"></a>00388         <span class="comment">/* xa - xc */</span>
297 <a name="l00389"></a>00389         r2 = pSrc[(2u * i0)] - pSrc[(2u * i2)];
298 <a name="l00390"></a>00390 
299 <a name="l00391"></a>00391         <span class="comment">/* ya + yc */</span>
300 <a name="l00392"></a>00392         s1 = pSrc[(2u * i0) + 1u] + pSrc[(2u * i2) + 1u];
301 <a name="l00393"></a>00393 
302 <a name="l00394"></a>00394         <span class="comment">/* ya - yc */</span>
303 <a name="l00395"></a>00395         s2 = pSrc[(2u * i0) + 1u] - pSrc[(2u * i2) + 1u];
304 <a name="l00396"></a>00396 
305 <a name="l00397"></a>00397         <span class="comment">/* xb + xd */</span>
306 <a name="l00398"></a>00398         t1 = pSrc[2u * i1] + pSrc[2u * i3];
307 <a name="l00399"></a>00399 
308 <a name="l00400"></a>00400         <span class="comment">/* xa&#39; = xa + xb + xc + xd */</span>
309 <a name="l00401"></a>00401         pSrc[2u * i0] = r1 + t1;
310 <a name="l00402"></a>00402 
311 <a name="l00403"></a>00403         <span class="comment">/* xa + xc -(xb + xd) */</span>
312 <a name="l00404"></a>00404         r1 = r1 - t1;
313 <a name="l00405"></a>00405 
314 <a name="l00406"></a>00406         <span class="comment">/* yb + yd */</span>
315 <a name="l00407"></a>00407         t2 = pSrc[(2u * i1) + 1u] + pSrc[(2u * i3) + 1u];
316 <a name="l00408"></a>00408 
317 <a name="l00409"></a>00409         <span class="comment">/* ya&#39; = ya + yb + yc + yd */</span>
318 <a name="l00410"></a>00410         pSrc[(2u * i0) + 1u] = s1 + t2;
319 <a name="l00411"></a>00411 
320 <a name="l00412"></a>00412         <span class="comment">/* (ya + yc) - (yb + yd) */</span>
321 <a name="l00413"></a>00413         s1 = s1 - t2;
322 <a name="l00414"></a>00414 
323 <a name="l00415"></a>00415         <span class="comment">/* (yb - yd) */</span>
324 <a name="l00416"></a>00416         t1 = pSrc[(2u * i1) + 1u] - pSrc[(2u * i3) + 1u];
325 <a name="l00417"></a>00417 
326 <a name="l00418"></a>00418         <span class="comment">/* (xb - xd) */</span>
327 <a name="l00419"></a>00419         t2 = pSrc[2u * i1] - pSrc[2u * i3];
328 <a name="l00420"></a>00420 
329 <a name="l00421"></a>00421         <span class="comment">/* xc&#39; = (xa-xb+xc-xd)co2 + (ya-yb+yc-yd)(si2) */</span>
330 <a name="l00422"></a>00422         pSrc[2u * i1] = (r1 * co2) + (s1 * si2);
331 <a name="l00423"></a>00423 
332 <a name="l00424"></a>00424         <span class="comment">/* yc&#39; = (ya-yb+yc-yd)co2 - (xa-xb+xc-xd)(si2) */</span>
333 <a name="l00425"></a>00425         pSrc[(2u * i1) + 1u] = (s1 * co2) - (r1 * si2);
334 <a name="l00426"></a>00426 
335 <a name="l00427"></a>00427         <span class="comment">/* (xa - xc) + (yb - yd) */</span>
336 <a name="l00428"></a>00428         r1 = r2 + t1;
337 <a name="l00429"></a>00429 
338 <a name="l00430"></a>00430         <span class="comment">/* (xa - xc) - (yb - yd) */</span>
339 <a name="l00431"></a>00431         r2 = r2 - t1;
340 <a name="l00432"></a>00432 
341 <a name="l00433"></a>00433         <span class="comment">/* (ya - yc) -  (xb - xd) */</span>
342 <a name="l00434"></a>00434         s1 = s2 - t2;
343 <a name="l00435"></a>00435 
344 <a name="l00436"></a>00436         <span class="comment">/* (ya - yc) +  (xb - xd) */</span>
345 <a name="l00437"></a>00437         s2 = s2 + t2;
346 <a name="l00438"></a>00438 
347 <a name="l00439"></a>00439         <span class="comment">/* xb&#39; = (xa+yb-xc-yd)co1 + (ya-xb-yc+xd)(si1) */</span>
348 <a name="l00440"></a>00440         pSrc[2u * i2] = (r1 * co1) + (s1 * si1);
349 <a name="l00441"></a>00441 
350 <a name="l00442"></a>00442         <span class="comment">/* yb&#39; = (ya-xb-yc+xd)co1 - (xa+yb-xc-yd)(si1) */</span>
351 <a name="l00443"></a>00443         pSrc[(2u * i2) + 1u] = (s1 * co1) - (r1 * si1);
352 <a name="l00444"></a>00444 
353 <a name="l00445"></a>00445         <span class="comment">/* xd&#39; = (xa-yb-xc+yd)co3 + (ya+xb-yc-xd)(si3) */</span>
354 <a name="l00446"></a>00446         pSrc[2u * i3] = (r2 * co3) + (s2 * si3);
355 <a name="l00447"></a>00447 
356 <a name="l00448"></a>00448         <span class="comment">/* yd&#39; = (ya+xb-yc-xd)co3 - (xa-yb-xc+yd)(si3) */</span>
357 <a name="l00449"></a>00449         pSrc[(2u * i3) + 1u] = (s2 * co3) - (r2 * si3);
358 <a name="l00450"></a>00450       }
359 <a name="l00451"></a>00451     }
360 <a name="l00452"></a>00452     twidCoefModifier &lt;&lt;= 2u;
361 <a name="l00453"></a>00453   }
362 <a name="l00454"></a>00454 
363 <a name="l00455"></a>00455   <span class="comment">/*  Initializations of last stage */</span>
364 <a name="l00456"></a>00456   n1 = n2;
365 <a name="l00457"></a>00457   n2 &gt;&gt;= 2u;
366 <a name="l00458"></a>00458 
367 <a name="l00459"></a>00459   <span class="comment">/*  Calculations of last stage */</span>
368 <a name="l00460"></a>00460   <span class="keywordflow">for</span> (i0 = 0u; i0 &lt;= (fftLen - n1); i0 += n1)
369 <a name="l00461"></a>00461   {
370 <a name="l00462"></a>00462     <span class="comment">/*  index calculation for the input as, */</span>
371 <a name="l00463"></a>00463     <span class="comment">/*  pSrc[i0 + 0], pSrc[i0 + fftLen/4], pSrc[i0 + fftLen/2], pSrc[i0 + 3fftLen/4] */</span>
372 <a name="l00464"></a>00464     i1 = i0 + n2;
373 <a name="l00465"></a>00465     i2 = i1 + n2;
374 <a name="l00466"></a>00466     i3 = i2 + n2;
375 <a name="l00467"></a>00467 
376 <a name="l00468"></a>00468     <span class="comment">/*  Butterfly implementation */</span>
377 <a name="l00469"></a>00469 
378 <a name="l00470"></a>00470     <span class="comment">/* xa + xb */</span>
379 <a name="l00471"></a>00471     r1 = pSrc[2u * i0] + pSrc[2u * i2];
380 <a name="l00472"></a>00472 
381 <a name="l00473"></a>00473     <span class="comment">/* xa - xb */</span>
382 <a name="l00474"></a>00474     r2 = pSrc[2u * i0] - pSrc[2u * i2];
383 <a name="l00475"></a>00475 
384 <a name="l00476"></a>00476     <span class="comment">/* ya + yc */</span>
385 <a name="l00477"></a>00477     s1 = pSrc[(2u * i0) + 1u] + pSrc[(2u * i2) + 1u];
386 <a name="l00478"></a>00478 
387 <a name="l00479"></a>00479     <span class="comment">/* ya - yc */</span>
388 <a name="l00480"></a>00480     s2 = pSrc[(2u * i0) + 1u] - pSrc[(2u * i2) + 1u];
389 <a name="l00481"></a>00481 
390 <a name="l00482"></a>00482     <span class="comment">/* xc + xd */</span>
391 <a name="l00483"></a>00483     t1 = pSrc[2u * i1] + pSrc[2u * i3];
392 <a name="l00484"></a>00484 
393 <a name="l00485"></a>00485     <span class="comment">/* xa&#39; = xa + xb + xc + xd */</span>
394 <a name="l00486"></a>00486     pSrc[2u * i0] = r1 + t1;
395 <a name="l00487"></a>00487 
396 <a name="l00488"></a>00488     <span class="comment">/* (xa + xb) - (xc + xd) */</span>
397 <a name="l00489"></a>00489     r1 = r1 - t1;
398 <a name="l00490"></a>00490 
399 <a name="l00491"></a>00491     <span class="comment">/* yb + yd */</span>
400 <a name="l00492"></a>00492     t2 = pSrc[(2u * i1) + 1u] + pSrc[(2u * i3) + 1u];
401 <a name="l00493"></a>00493 
402 <a name="l00494"></a>00494     <span class="comment">/* ya&#39; = ya + yb + yc + yd */</span>
403 <a name="l00495"></a>00495     pSrc[(2u * i0) + 1u] = s1 + t2;
404 <a name="l00496"></a>00496 
405 <a name="l00497"></a>00497     <span class="comment">/* (ya + yc) - (yb + yd) */</span>
406 <a name="l00498"></a>00498     s1 = s1 - t2;
407 <a name="l00499"></a>00499 
408 <a name="l00500"></a>00500     <span class="comment">/* (yb-yd) */</span>
409 <a name="l00501"></a>00501     t1 = pSrc[(2u * i1) + 1u] - pSrc[(2u * i3) + 1u];
410 <a name="l00502"></a>00502 
411 <a name="l00503"></a>00503     <span class="comment">/* (xb-xd) */</span>
412 <a name="l00504"></a>00504     t2 = pSrc[2u * i1] - pSrc[2u * i3];
413 <a name="l00505"></a>00505 
414 <a name="l00506"></a>00506     <span class="comment">/* xc&#39; = (xa-xb+xc-xd)co2 + (ya-yb+yc-yd)(si2) */</span>
415 <a name="l00507"></a>00507     pSrc[2u * i1] = r1;
416 <a name="l00508"></a>00508 
417 <a name="l00509"></a>00509     <span class="comment">/* yc&#39; = (ya-yb+yc-yd)co2 - (xa-xb+xc-xd)(si2) */</span>
418 <a name="l00510"></a>00510     pSrc[(2u * i1) + 1u] = s1;
419 <a name="l00511"></a>00511 
420 <a name="l00512"></a>00512     <span class="comment">/* (xa+yb-xc-yd) */</span>
421 <a name="l00513"></a>00513     r1 = r2 + t1;
422 <a name="l00514"></a>00514 
423 <a name="l00515"></a>00515     <span class="comment">/* (xa-yb-xc+yd) */</span>
424 <a name="l00516"></a>00516     r2 = r2 - t1;
425 <a name="l00517"></a>00517 
426 <a name="l00518"></a>00518     <span class="comment">/* (ya-xb-yc+xd) */</span>
427 <a name="l00519"></a>00519     s1 = s2 - t2;
428 <a name="l00520"></a>00520 
429 <a name="l00521"></a>00521     <span class="comment">/* (ya+xb-yc-xd) */</span>
430 <a name="l00522"></a>00522     s2 = s2 + t2;
431 <a name="l00523"></a>00523 
432 <a name="l00524"></a>00524     <span class="comment">/* xb&#39; = (xa+yb-xc-yd)co1 + (ya-xb-yc+xd)(si1) */</span>
433 <a name="l00525"></a>00525     pSrc[2u * i2] = r1;
434 <a name="l00526"></a>00526 
435 <a name="l00527"></a>00527     <span class="comment">/* yb&#39; = (ya-xb-yc+xd)co1 - (xa+yb-xc-yd)(si1) */</span>
436 <a name="l00528"></a>00528     pSrc[(2u * i2) + 1u] = s1;
437 <a name="l00529"></a>00529 
438 <a name="l00530"></a>00530     <span class="comment">/* xd&#39; = (xa-yb-xc+yd)co3 + (ya+xb-yc-xd)(si3) */</span>
439 <a name="l00531"></a>00531     pSrc[2u * i3] = r2;
440 <a name="l00532"></a>00532 
441 <a name="l00533"></a>00533     <span class="comment">/* yd&#39; = (ya+xb-yc-xd)co3 - (xa-yb-xc+yd)(si3) */</span>
442 <a name="l00534"></a>00534     pSrc[(2u * i3) + 1u] = s2;
443 <a name="l00535"></a>00535   }
444 <a name="l00536"></a>00536 
445 <a name="l00537"></a>00537 
446 <a name="l00538"></a>00538 <span class="preprocessor">#else</span>
447 <a name="l00539"></a>00539 <span class="preprocessor"></span>
448 <a name="l00540"></a>00540   <span class="comment">/* Run the below code for Cortex-M0 */</span>
449 <a name="l00541"></a>00541 
450 <a name="l00542"></a>00542   <span class="comment">/*  Initializations for the fft calculation */</span>
451 <a name="l00543"></a>00543   n2 = fftLen;
452 <a name="l00544"></a>00544   n1 = n2;
453 <a name="l00545"></a>00545   <span class="keywordflow">for</span> (k = fftLen; k &gt; 1u; k &gt;&gt;= 2u)
454 <a name="l00546"></a>00546   {
455 <a name="l00547"></a>00547     <span class="comment">/*  Initializations for the fft calculation */</span>
456 <a name="l00548"></a>00548     n1 = n2;
457 <a name="l00549"></a>00549     n2 &gt;&gt;= 2u;
458 <a name="l00550"></a>00550     ia1 = 0u;
459 <a name="l00551"></a>00551 
460 <a name="l00552"></a>00552     <span class="comment">/*  FFT Calculation */</span>
461 <a name="l00553"></a>00553     <span class="keywordflow">for</span> (j = 0u; j &lt;= (n2 - 1u); j++)
462 <a name="l00554"></a>00554     {
463 <a name="l00555"></a>00555       <span class="comment">/*  index calculation for the coefficients */</span>
464 <a name="l00556"></a>00556       ia2 = ia1 + ia1;
465 <a name="l00557"></a>00557       ia3 = ia2 + ia1;
466 <a name="l00558"></a>00558       co1 = pCoef[ia1 * 2u];
467 <a name="l00559"></a>00559       si1 = pCoef[(ia1 * 2u) + 1u];
468 <a name="l00560"></a>00560       co2 = pCoef[ia2 * 2u];
469 <a name="l00561"></a>00561       si2 = pCoef[(ia2 * 2u) + 1u];
470 <a name="l00562"></a>00562       co3 = pCoef[ia3 * 2u];
471 <a name="l00563"></a>00563       si3 = pCoef[(ia3 * 2u) + 1u];
472 <a name="l00564"></a>00564 
473 <a name="l00565"></a>00565       <span class="comment">/*  Twiddle coefficients index modifier */</span>
474 <a name="l00566"></a>00566       ia1 = ia1 + twidCoefModifier;
475 <a name="l00567"></a>00567 
476 <a name="l00568"></a>00568       <span class="keywordflow">for</span> (i0 = j; i0 &lt; fftLen; i0 += n1)
477 <a name="l00569"></a>00569       {
478 <a name="l00570"></a>00570         <span class="comment">/*  index calculation for the input as, */</span>
479 <a name="l00571"></a>00571         <span class="comment">/*  pSrc[i0 + 0], pSrc[i0 + fftLen/4], pSrc[i0 + fftLen/2], pSrc[i0 + 3fftLen/4] */</span>
480 <a name="l00572"></a>00572         i1 = i0 + n2;
481 <a name="l00573"></a>00573         i2 = i1 + n2;
482 <a name="l00574"></a>00574         i3 = i2 + n2;
483 <a name="l00575"></a>00575 
484 <a name="l00576"></a>00576         <span class="comment">/* xa + xc */</span>
485 <a name="l00577"></a>00577         r1 = pSrc[(2u * i0)] + pSrc[(2u * i2)];
486 <a name="l00578"></a>00578 
487 <a name="l00579"></a>00579         <span class="comment">/* xa - xc */</span>
488 <a name="l00580"></a>00580         r2 = pSrc[(2u * i0)] - pSrc[(2u * i2)];
489 <a name="l00581"></a>00581 
490 <a name="l00582"></a>00582         <span class="comment">/* ya + yc */</span>
491 <a name="l00583"></a>00583         s1 = pSrc[(2u * i0) + 1u] + pSrc[(2u * i2) + 1u];
492 <a name="l00584"></a>00584 
493 <a name="l00585"></a>00585         <span class="comment">/* ya - yc */</span>
494 <a name="l00586"></a>00586         s2 = pSrc[(2u * i0) + 1u] - pSrc[(2u * i2) + 1u];
495 <a name="l00587"></a>00587 
496 <a name="l00588"></a>00588         <span class="comment">/* xb + xd */</span>
497 <a name="l00589"></a>00589         t1 = pSrc[2u * i1] + pSrc[2u * i3];
498 <a name="l00590"></a>00590 
499 <a name="l00591"></a>00591         <span class="comment">/* xa&#39; = xa + xb + xc + xd */</span>
500 <a name="l00592"></a>00592         pSrc[2u * i0] = r1 + t1;
501 <a name="l00593"></a>00593 
502 <a name="l00594"></a>00594         <span class="comment">/* xa + xc -(xb + xd) */</span>
503 <a name="l00595"></a>00595         r1 = r1 - t1;
504 <a name="l00596"></a>00596 
505 <a name="l00597"></a>00597         <span class="comment">/* yb + yd */</span>
506 <a name="l00598"></a>00598         t2 = pSrc[(2u * i1) + 1u] + pSrc[(2u * i3) + 1u];
507 <a name="l00599"></a>00599 
508 <a name="l00600"></a>00600         <span class="comment">/* ya&#39; = ya + yb + yc + yd */</span>
509 <a name="l00601"></a>00601         pSrc[(2u * i0) + 1u] = s1 + t2;
510 <a name="l00602"></a>00602 
511 <a name="l00603"></a>00603         <span class="comment">/* (ya + yc) - (yb + yd) */</span>
512 <a name="l00604"></a>00604         s1 = s1 - t2;
513 <a name="l00605"></a>00605 
514 <a name="l00606"></a>00606         <span class="comment">/* (yb - yd) */</span>
515 <a name="l00607"></a>00607         t1 = pSrc[(2u * i1) + 1u] - pSrc[(2u * i3) + 1u];
516 <a name="l00608"></a>00608 
517 <a name="l00609"></a>00609         <span class="comment">/* (xb - xd) */</span>
518 <a name="l00610"></a>00610         t2 = pSrc[2u * i1] - pSrc[2u * i3];
519 <a name="l00611"></a>00611 
520 <a name="l00612"></a>00612         <span class="comment">/* xc&#39; = (xa-xb+xc-xd)co2 + (ya-yb+yc-yd)(si2) */</span>
521 <a name="l00613"></a>00613         pSrc[2u * i1] = (r1 * co2) + (s1 * si2);
522 <a name="l00614"></a>00614 
523 <a name="l00615"></a>00615         <span class="comment">/* yc&#39; = (ya-yb+yc-yd)co2 - (xa-xb+xc-xd)(si2) */</span>
524 <a name="l00616"></a>00616         pSrc[(2u * i1) + 1u] = (s1 * co2) - (r1 * si2);
525 <a name="l00617"></a>00617 
526 <a name="l00618"></a>00618         <span class="comment">/* (xa - xc) + (yb - yd) */</span>
527 <a name="l00619"></a>00619         r1 = r2 + t1;
528 <a name="l00620"></a>00620 
529 <a name="l00621"></a>00621         <span class="comment">/* (xa - xc) - (yb - yd) */</span>
530 <a name="l00622"></a>00622         r2 = r2 - t1;
531 <a name="l00623"></a>00623 
532 <a name="l00624"></a>00624         <span class="comment">/* (ya - yc) -  (xb - xd) */</span>
533 <a name="l00625"></a>00625         s1 = s2 - t2;
534 <a name="l00626"></a>00626 
535 <a name="l00627"></a>00627         <span class="comment">/* (ya - yc) +  (xb - xd) */</span>
536 <a name="l00628"></a>00628         s2 = s2 + t2;
537 <a name="l00629"></a>00629 
538 <a name="l00630"></a>00630         <span class="comment">/* xb&#39; = (xa+yb-xc-yd)co1 + (ya-xb-yc+xd)(si1) */</span>
539 <a name="l00631"></a>00631         pSrc[2u * i2] = (r1 * co1) + (s1 * si1);
540 <a name="l00632"></a>00632 
541 <a name="l00633"></a>00633         <span class="comment">/* yb&#39; = (ya-xb-yc+xd)co1 - (xa+yb-xc-yd)(si1) */</span>
542 <a name="l00634"></a>00634         pSrc[(2u * i2) + 1u] = (s1 * co1) - (r1 * si1);
543 <a name="l00635"></a>00635 
544 <a name="l00636"></a>00636         <span class="comment">/* xd&#39; = (xa-yb-xc+yd)co3 + (ya+xb-yc-xd)(si3) */</span>
545 <a name="l00637"></a>00637         pSrc[2u * i3] = (r2 * co3) + (s2 * si3);
546 <a name="l00638"></a>00638 
547 <a name="l00639"></a>00639         <span class="comment">/* yd&#39; = (ya+xb-yc-xd)co3 - (xa-yb-xc+yd)(si3) */</span>
548 <a name="l00640"></a>00640         pSrc[(2u * i3) + 1u] = (s2 * co3) - (r2 * si3);
549 <a name="l00641"></a>00641       }
550 <a name="l00642"></a>00642     }
551 <a name="l00643"></a>00643     twidCoefModifier &lt;&lt;= 2u;
552 <a name="l00644"></a>00644   }
553 <a name="l00645"></a>00645 
554 <a name="l00646"></a>00646 <span class="preprocessor">#endif </span><span class="comment">/* #ifndef ARM_MATH_CM0 */</span>
555 <a name="l00647"></a>00647 
556 <a name="l00648"></a>00648 }
557 <a name="l00649"></a>00649 
558 <a name="l00650"></a>00650 <span class="comment">/*   </span>
559 <a name="l00651"></a>00651 <span class="comment"> * @brief  Core function for the floating-point CIFFT butterfly process.  </span>
560 <a name="l00652"></a>00652 <span class="comment"> * @param[in, out] *pSrc            points to the in-place buffer of floating-point data type.  </span>
561 <a name="l00653"></a>00653 <span class="comment"> * @param[in]      fftLen           length of the FFT.  </span>
562 <a name="l00654"></a>00654 <span class="comment"> * @param[in]      *pCoef           points to twiddle coefficient buffer.  </span>
563 <a name="l00655"></a>00655 <span class="comment"> * @param[in]      twidCoefModifier twiddle coefficient modifier that supports different size FFTs with the same twiddle factor table.  </span>
564 <a name="l00656"></a>00656 <span class="comment"> * @param[in]      onebyfftLen      value of 1/fftLen.  </span>
565 <a name="l00657"></a>00657 <span class="comment"> * @return none.  </span>
566 <a name="l00658"></a>00658 <span class="comment"> */</span>
567 <a name="l00659"></a>00659 
568 <a name="l00660"></a><a class="code" href="arm__math_8h.html#a2a78df6e4bbf080624f2b6349224ec93">00660</a> <span class="keywordtype">void</span> <a class="code" href="arm__cfft__radix4__f32_8c.html#a2a78df6e4bbf080624f2b6349224ec93" title="Core function for the floating-point CIFFT butterfly process.">arm_radix4_butterfly_inverse_f32</a>(
569 <a name="l00661"></a>00661   <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> * pSrc,
570 <a name="l00662"></a>00662   uint16_t fftLen,
571 <a name="l00663"></a>00663   <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> * pCoef,
572 <a name="l00664"></a>00664   uint16_t twidCoefModifier,
573 <a name="l00665"></a>00665   <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> onebyfftLen)
574 <a name="l00666"></a>00666 {
575 <a name="l00667"></a>00667   <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> co1, co2, co3, si1, si2, si3;
576 <a name="l00668"></a>00668   <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> t1, t2, r1, r2, s1, s2;
577 <a name="l00669"></a>00669   uint32_t ia1, ia2, ia3;
578 <a name="l00670"></a>00670   uint32_t i0, i1, i2, i3;
579 <a name="l00671"></a>00671   uint32_t n1, n2, j, k;
580 <a name="l00672"></a>00672 
581 <a name="l00673"></a>00673 <span class="preprocessor">#ifndef ARM_MATH_CM0</span>
582 <a name="l00674"></a>00674 <span class="preprocessor"></span>
583 <a name="l00675"></a>00675   <span class="comment">/* Run the below code for Cortex-M4 and Cortex-M3 */</span>
584 <a name="l00676"></a>00676 
585 <a name="l00677"></a>00677   <span class="comment">/*  Initializations for the first stage */</span>
586 <a name="l00678"></a>00678   n2 = fftLen;
587 <a name="l00679"></a>00679   n1 = n2;
588 <a name="l00680"></a>00680 
589 <a name="l00681"></a>00681   <span class="comment">/* n2 = fftLen/4 */</span>
590 <a name="l00682"></a>00682   n2 &gt;&gt;= 2u;
591 <a name="l00683"></a>00683   i0 = 0u;
592 <a name="l00684"></a>00684   ia1 = 0u;
593 <a name="l00685"></a>00685 
594 <a name="l00686"></a>00686   j = n2;
595 <a name="l00687"></a>00687 
596 <a name="l00688"></a>00688   <span class="comment">/*  Calculation of first stage */</span>
597 <a name="l00689"></a>00689   <span class="keywordflow">do</span>
598 <a name="l00690"></a>00690   {
599 <a name="l00691"></a>00691     <span class="comment">/*  index calculation for the input as, */</span>
600 <a name="l00692"></a>00692     <span class="comment">/*  pSrc[i0 + 0], pSrc[i0 + fftLen/4], pSrc[i0 + fftLen/2], pSrc[i0 + 3fftLen/4] */</span>
601 <a name="l00693"></a>00693     i1 = i0 + n2;
602 <a name="l00694"></a>00694     i2 = i1 + n2;
603 <a name="l00695"></a>00695     i3 = i2 + n2;
604 <a name="l00696"></a>00696 
605 <a name="l00697"></a>00697     <span class="comment">/*  Butterfly implementation */</span>
606 <a name="l00698"></a>00698     <span class="comment">/* xa + xc */</span>
607 <a name="l00699"></a>00699     r1 = pSrc[(2u * i0)] + pSrc[(2u * i2)];
608 <a name="l00700"></a>00700 
609 <a name="l00701"></a>00701     <span class="comment">/* xa - xc */</span>
610 <a name="l00702"></a>00702     r2 = pSrc[2u * i0] - pSrc[2u * i2];
611 <a name="l00703"></a>00703 
612 <a name="l00704"></a>00704     <span class="comment">/* ya + yc */</span>
613 <a name="l00705"></a>00705     s1 = pSrc[(2u * i0) + 1u] + pSrc[(2u * i2) + 1u];
614 <a name="l00706"></a>00706 
615 <a name="l00707"></a>00707     <span class="comment">/* ya - yc */</span>
616 <a name="l00708"></a>00708     s2 = pSrc[(2u * i0) + 1u] - pSrc[(2u * i2) + 1u];
617 <a name="l00709"></a>00709 
618 <a name="l00710"></a>00710     <span class="comment">/* xb + xd */</span>
619 <a name="l00711"></a>00711     t1 = pSrc[2u * i1] + pSrc[2u * i3];
620 <a name="l00712"></a>00712 
621 <a name="l00713"></a>00713     <span class="comment">/* xa&#39; = xa + xb + xc + xd */</span>
622 <a name="l00714"></a>00714     pSrc[2u * i0] = r1 + t1;
623 <a name="l00715"></a>00715 
624 <a name="l00716"></a>00716     <span class="comment">/* (xa + xc) - (xb + xd) */</span>
625 <a name="l00717"></a>00717     r1 = r1 - t1;
626 <a name="l00718"></a>00718 
627 <a name="l00719"></a>00719     <span class="comment">/* yb + yd */</span>
628 <a name="l00720"></a>00720     t2 = pSrc[(2u * i1) + 1u] + pSrc[(2u * i3) + 1u];
629 <a name="l00721"></a>00721 
630 <a name="l00722"></a>00722     <span class="comment">/* ya&#39; = ya + yb + yc + yd */</span>
631 <a name="l00723"></a>00723     pSrc[(2u * i0) + 1u] = s1 + t2;
632 <a name="l00724"></a>00724 
633 <a name="l00725"></a>00725     <span class="comment">/* (ya + yc) - (yb + yd) */</span>
634 <a name="l00726"></a>00726     s1 = s1 - t2;
635 <a name="l00727"></a>00727 
636 <a name="l00728"></a>00728     <span class="comment">/* yb - yd */</span>
637 <a name="l00729"></a>00729     t1 = pSrc[(2u * i1) + 1u] - pSrc[(2u * i3) + 1u];
638 <a name="l00730"></a>00730 
639 <a name="l00731"></a>00731     <span class="comment">/* xb - xd */</span>
640 <a name="l00732"></a>00732     t2 = pSrc[2u * i1] - pSrc[2u * i3];
641 <a name="l00733"></a>00733 
642 <a name="l00734"></a>00734     <span class="comment">/*  index calculation for the coefficients */</span>
643 <a name="l00735"></a>00735     ia2 = ia1 + ia1;
644 <a name="l00736"></a>00736     co2 = pCoef[ia2 * 2u];
645 <a name="l00737"></a>00737     si2 = pCoef[(ia2 * 2u) + 1u];
646 <a name="l00738"></a>00738 
647 <a name="l00739"></a>00739     <span class="comment">/* xc&#39; = (xa-xb+xc-xd)co2 - (ya-yb+yc-yd)(si2) */</span>
648 <a name="l00740"></a>00740     pSrc[2u * i1] = (r1 * co2) - (s1 * si2);
649 <a name="l00741"></a>00741 
650 <a name="l00742"></a>00742     <span class="comment">/* yc&#39; = (ya-yb+yc-yd)co2 + (xa-xb+xc-xd)(si2) */</span>
651 <a name="l00743"></a>00743     pSrc[(2u * i1) + 1u] = (s1 * co2) + (r1 * si2);
652 <a name="l00744"></a>00744 
653 <a name="l00745"></a>00745     <span class="comment">/* (xa - xc) - (yb - yd) */</span>
654 <a name="l00746"></a>00746     r1 = r2 - t1;
655 <a name="l00747"></a>00747 
656 <a name="l00748"></a>00748     <span class="comment">/* (xa - xc) + (yb - yd) */</span>
657 <a name="l00749"></a>00749     r2 = r2 + t1;
658 <a name="l00750"></a>00750 
659 <a name="l00751"></a>00751     <span class="comment">/* (ya - yc) + (xb - xd) */</span>
660 <a name="l00752"></a>00752     s1 = s2 + t2;
661 <a name="l00753"></a>00753 
662 <a name="l00754"></a>00754     <span class="comment">/* (ya - yc) - (xb - xd) */</span>
663 <a name="l00755"></a>00755     s2 = s2 - t2;
664 <a name="l00756"></a>00756 
665 <a name="l00757"></a>00757     co1 = pCoef[ia1 * 2u];
666 <a name="l00758"></a>00758     si1 = pCoef[(ia1 * 2u) + 1u];
667 <a name="l00759"></a>00759 
668 <a name="l00760"></a>00760     <span class="comment">/* xb&#39; = (xa+yb-xc-yd)co1 - (ya-xb-yc+xd)(si1) */</span>
669 <a name="l00761"></a>00761     pSrc[2u * i2] = (r1 * co1) - (s1 * si1);
670 <a name="l00762"></a>00762 
671 <a name="l00763"></a>00763     <span class="comment">/* yb&#39; = (ya-xb-yc+xd)co1 + (xa+yb-xc-yd)(si1) */</span>
672 <a name="l00764"></a>00764     pSrc[(2u * i2) + 1u] = (s1 * co1) + (r1 * si1);
673 <a name="l00765"></a>00765 
674 <a name="l00766"></a>00766     <span class="comment">/*  index calculation for the coefficients */</span>
675 <a name="l00767"></a>00767     ia3 = ia2 + ia1;
676 <a name="l00768"></a>00768     co3 = pCoef[ia3 * 2u];
677 <a name="l00769"></a>00769     si3 = pCoef[(ia3 * 2u) + 1u];
678 <a name="l00770"></a>00770 
679 <a name="l00771"></a>00771     <span class="comment">/* xd&#39; = (xa-yb-xc+yd)co3 - (ya+xb-yc-xd)(si3) */</span>
680 <a name="l00772"></a>00772     pSrc[2u * i3] = (r2 * co3) - (s2 * si3);
681 <a name="l00773"></a>00773 
682 <a name="l00774"></a>00774     <span class="comment">/* yd&#39; = (ya+xb-yc-xd)co3 + (xa-yb-xc+yd)(si3) */</span>
683 <a name="l00775"></a>00775     pSrc[(2u * i3) + 1u] = (s2 * co3) + (r2 * si3);
684 <a name="l00776"></a>00776 
685 <a name="l00777"></a>00777     <span class="comment">/*  Twiddle coefficients index modifier */</span>
686 <a name="l00778"></a>00778     ia1 = ia1 + twidCoefModifier;
687 <a name="l00779"></a>00779 
688 <a name="l00780"></a>00780     <span class="comment">/*  Updating input index */</span>
689 <a name="l00781"></a>00781     i0 = i0 + 1u;
690 <a name="l00782"></a>00782 
691 <a name="l00783"></a>00783   }
692 <a name="l00784"></a>00784   <span class="keywordflow">while</span>(--j);
693 <a name="l00785"></a>00785 
694 <a name="l00786"></a>00786   twidCoefModifier &lt;&lt;= 2u;
695 <a name="l00787"></a>00787 
696 <a name="l00788"></a>00788   <span class="comment">/*  Calculation of second stage to excluding last stage */</span>
697 <a name="l00789"></a>00789   <span class="keywordflow">for</span> (k = fftLen / 4; k &gt; 4u; k &gt;&gt;= 2u)
698 <a name="l00790"></a>00790   {
699 <a name="l00791"></a>00791     <span class="comment">/*  Initializations for the first stage */</span>
700 <a name="l00792"></a>00792     n1 = n2;
701 <a name="l00793"></a>00793     n2 &gt;&gt;= 2u;
702 <a name="l00794"></a>00794     ia1 = 0u;
703 <a name="l00795"></a>00795 
704 <a name="l00796"></a>00796     <span class="comment">/*  Calculation of first stage */</span>
705 <a name="l00797"></a>00797     <span class="keywordflow">for</span> (j = 0u; j &lt;= (n2 - 1u); j++)
706 <a name="l00798"></a>00798     {
707 <a name="l00799"></a>00799       <span class="comment">/*  index calculation for the coefficients */</span>
708 <a name="l00800"></a>00800       ia2 = ia1 + ia1;
709 <a name="l00801"></a>00801       ia3 = ia2 + ia1;
710 <a name="l00802"></a>00802       co1 = pCoef[ia1 * 2u];
711 <a name="l00803"></a>00803       si1 = pCoef[(ia1 * 2u) + 1u];
712 <a name="l00804"></a>00804       co2 = pCoef[ia2 * 2u];
713 <a name="l00805"></a>00805       si2 = pCoef[(ia2 * 2u) + 1u];
714 <a name="l00806"></a>00806       co3 = pCoef[ia3 * 2u];
715 <a name="l00807"></a>00807       si3 = pCoef[(ia3 * 2u) + 1u];
716 <a name="l00808"></a>00808 
717 <a name="l00809"></a>00809       <span class="comment">/*  Twiddle coefficients index modifier */</span>
718 <a name="l00810"></a>00810       ia1 = ia1 + twidCoefModifier;
719 <a name="l00811"></a>00811 
720 <a name="l00812"></a>00812       <span class="keywordflow">for</span> (i0 = j; i0 &lt; fftLen; i0 += n1)
721 <a name="l00813"></a>00813       {
722 <a name="l00814"></a>00814         <span class="comment">/*  index calculation for the input as, */</span>
723 <a name="l00815"></a>00815         <span class="comment">/*  pSrc[i0 + 0], pSrc[i0 + fftLen/4], pSrc[i0 + fftLen/2], pSrc[i0 + 3fftLen/4] */</span>
724 <a name="l00816"></a>00816         i1 = i0 + n2;
725 <a name="l00817"></a>00817         i2 = i1 + n2;
726 <a name="l00818"></a>00818         i3 = i2 + n2;
727 <a name="l00819"></a>00819 
728 <a name="l00820"></a>00820         <span class="comment">/* xa + xc */</span>
729 <a name="l00821"></a>00821         r1 = pSrc[(2u * i0)] + pSrc[(2u * i2)];
730 <a name="l00822"></a>00822 
731 <a name="l00823"></a>00823         <span class="comment">/* xa - xc */</span>
732 <a name="l00824"></a>00824         r2 = pSrc[(2u * i0)] - pSrc[(2u * i2)];
733 <a name="l00825"></a>00825 
734 <a name="l00826"></a>00826         <span class="comment">/* ya + yc */</span>
735 <a name="l00827"></a>00827         s1 = pSrc[(2u * i0) + 1u] + pSrc[(2u * i2) + 1u];
736 <a name="l00828"></a>00828 
737 <a name="l00829"></a>00829         <span class="comment">/* ya - yc */</span>
738 <a name="l00830"></a>00830         s2 = pSrc[(2u * i0) + 1u] - pSrc[(2u * i2) + 1u];
739 <a name="l00831"></a>00831 
740 <a name="l00832"></a>00832         <span class="comment">/* xb + xd */</span>
741 <a name="l00833"></a>00833         t1 = pSrc[2u * i1] + pSrc[2u * i3];
742 <a name="l00834"></a>00834 
743 <a name="l00835"></a>00835         <span class="comment">/* xa&#39; = xa + xb + xc + xd */</span>
744 <a name="l00836"></a>00836         pSrc[2u * i0] = r1 + t1;
745 <a name="l00837"></a>00837 
746 <a name="l00838"></a>00838         <span class="comment">/* xa + xc -(xb + xd) */</span>
747 <a name="l00839"></a>00839         r1 = r1 - t1;
748 <a name="l00840"></a>00840 
749 <a name="l00841"></a>00841         <span class="comment">/* yb + yd */</span>
750 <a name="l00842"></a>00842         t2 = pSrc[(2u * i1) + 1u] + pSrc[(2u * i3) + 1u];
751 <a name="l00843"></a>00843 
752 <a name="l00844"></a>00844         <span class="comment">/* ya&#39; = ya + yb + yc + yd */</span>
753 <a name="l00845"></a>00845         pSrc[(2u * i0) + 1u] = s1 + t2;
754 <a name="l00846"></a>00846 
755 <a name="l00847"></a>00847         <span class="comment">/* (ya + yc) - (yb + yd) */</span>
756 <a name="l00848"></a>00848         s1 = s1 - t2;
757 <a name="l00849"></a>00849 
758 <a name="l00850"></a>00850         <span class="comment">/* (yb - yd) */</span>
759 <a name="l00851"></a>00851         t1 = pSrc[(2u * i1) + 1u] - pSrc[(2u * i3) + 1u];
760 <a name="l00852"></a>00852 
761 <a name="l00853"></a>00853         <span class="comment">/* (xb - xd) */</span>
762 <a name="l00854"></a>00854         t2 = pSrc[2u * i1] - pSrc[2u * i3];
763 <a name="l00855"></a>00855 
764 <a name="l00856"></a>00856         <span class="comment">/* xc&#39; = (xa-xb+xc-xd)co2 - (ya-yb+yc-yd)(si2) */</span>
765 <a name="l00857"></a>00857         pSrc[2u * i1] = (r1 * co2) - (s1 * si2);
766 <a name="l00858"></a>00858 
767 <a name="l00859"></a>00859         <span class="comment">/* yc&#39; = (ya-yb+yc-yd)co2 + (xa-xb+xc-xd)(si2) */</span>
768 <a name="l00860"></a>00860         pSrc[(2u * i1) + 1u] = (s1 * co2) + (r1 * si2);
769 <a name="l00861"></a>00861 
770 <a name="l00862"></a>00862         <span class="comment">/* (xa - xc) - (yb - yd) */</span>
771 <a name="l00863"></a>00863         r1 = r2 - t1;
772 <a name="l00864"></a>00864 
773 <a name="l00865"></a>00865         <span class="comment">/* (xa - xc) + (yb - yd) */</span>
774 <a name="l00866"></a>00866         r2 = r2 + t1;
775 <a name="l00867"></a>00867 
776 <a name="l00868"></a>00868         <span class="comment">/* (ya - yc) +  (xb - xd) */</span>
777 <a name="l00869"></a>00869         s1 = s2 + t2;
778 <a name="l00870"></a>00870 
779 <a name="l00871"></a>00871         <span class="comment">/* (ya - yc) -  (xb - xd) */</span>
780 <a name="l00872"></a>00872         s2 = s2 - t2;
781 <a name="l00873"></a>00873 
782 <a name="l00874"></a>00874         <span class="comment">/* xb&#39; = (xa+yb-xc-yd)co1 - (ya-xb-yc+xd)(si1) */</span>
783 <a name="l00875"></a>00875         pSrc[2u * i2] = (r1 * co1) - (s1 * si1);
784 <a name="l00876"></a>00876 
785 <a name="l00877"></a>00877         <span class="comment">/* yb&#39; = (ya-xb-yc+xd)co1 + (xa+yb-xc-yd)(si1) */</span>
786 <a name="l00878"></a>00878         pSrc[(2u * i2) + 1u] = (s1 * co1) + (r1 * si1);
787 <a name="l00879"></a>00879 
788 <a name="l00880"></a>00880         <span class="comment">/* xd&#39; = (xa-yb-xc+yd)co3 - (ya+xb-yc-xd)(si3) */</span>
789 <a name="l00881"></a>00881         pSrc[2u * i3] = (r2 * co3) - (s2 * si3);
790 <a name="l00882"></a>00882 
791 <a name="l00883"></a>00883         <span class="comment">/* yd&#39; = (ya+xb-yc-xd)co3 + (xa-yb-xc+yd)(si3) */</span>
792 <a name="l00884"></a>00884         pSrc[(2u * i3) + 1u] = (s2 * co3) + (r2 * si3);
793 <a name="l00885"></a>00885       }
794 <a name="l00886"></a>00886     }
795 <a name="l00887"></a>00887     twidCoefModifier &lt;&lt;= 2u;
796 <a name="l00888"></a>00888   }
797 <a name="l00889"></a>00889 
798 <a name="l00890"></a>00890   <span class="comment">/*  Initializations of last stage */</span>
799 <a name="l00891"></a>00891   n1 = n2;
800 <a name="l00892"></a>00892   n2 &gt;&gt;= 2u;
801 <a name="l00893"></a>00893 
802 <a name="l00894"></a>00894   <span class="comment">/*  Calculations of last stage */</span>
803 <a name="l00895"></a>00895   <span class="keywordflow">for</span> (i0 = 0u; i0 &lt;= (fftLen - n1); i0 += n1)
804 <a name="l00896"></a>00896   {
805 <a name="l00897"></a>00897     <span class="comment">/*  index calculation for the input as, */</span>
806 <a name="l00898"></a>00898     <span class="comment">/*  pSrc[i0 + 0], pSrc[i0 + fftLen/4], pSrc[i0 + fftLen/2], pSrc[i0 + 3fftLen/4] */</span>
807 <a name="l00899"></a>00899     i1 = i0 + n2;
808 <a name="l00900"></a>00900     i2 = i1 + n2;
809 <a name="l00901"></a>00901     i3 = i2 + n2;
810 <a name="l00902"></a>00902 
811 <a name="l00903"></a>00903     <span class="comment">/*  Butterfly implementation */</span>
812 <a name="l00904"></a>00904     <span class="comment">/* xa + xc */</span>
813 <a name="l00905"></a>00905     r1 = pSrc[2u * i0] + pSrc[2u * i2];
814 <a name="l00906"></a>00906 
815 <a name="l00907"></a>00907     <span class="comment">/* xa - xc */</span>
816 <a name="l00908"></a>00908     r2 = pSrc[2u * i0] - pSrc[2u * i2];
817 <a name="l00909"></a>00909 
818 <a name="l00910"></a>00910     <span class="comment">/* ya + yc */</span>
819 <a name="l00911"></a>00911     s1 = pSrc[(2u * i0) + 1u] + pSrc[(2u * i2) + 1u];
820 <a name="l00912"></a>00912 
821 <a name="l00913"></a>00913     <span class="comment">/* ya - yc */</span>
822 <a name="l00914"></a>00914     s2 = pSrc[(2u * i0) + 1u] - pSrc[(2u * i2) + 1u];
823 <a name="l00915"></a>00915 
824 <a name="l00916"></a>00916     <span class="comment">/* xc + xd */</span>
825 <a name="l00917"></a>00917     t1 = pSrc[2u * i1] + pSrc[2u * i3];
826 <a name="l00918"></a>00918 
827 <a name="l00919"></a>00919     <span class="comment">/* xa&#39; = xa + xb + xc + xd */</span>
828 <a name="l00920"></a>00920     pSrc[2u * i0] = (r1 + t1) * onebyfftLen;
829 <a name="l00921"></a>00921 
830 <a name="l00922"></a>00922     <span class="comment">/* (xa + xb) - (xc + xd) */</span>
831 <a name="l00923"></a>00923     r1 = r1 - t1;
832 <a name="l00924"></a>00924 
833 <a name="l00925"></a>00925     <span class="comment">/* yb + yd */</span>
834 <a name="l00926"></a>00926     t2 = pSrc[(2u * i1) + 1u] + pSrc[(2u * i3) + 1u];
835 <a name="l00927"></a>00927 
836 <a name="l00928"></a>00928     <span class="comment">/* ya&#39; = ya + yb + yc + yd */</span>
837 <a name="l00929"></a>00929     pSrc[(2u * i0) + 1u] = (s1 + t2) * onebyfftLen;
838 <a name="l00930"></a>00930 
839 <a name="l00931"></a>00931     <span class="comment">/* (ya + yc) - (yb + yd) */</span>
840 <a name="l00932"></a>00932     s1 = s1 - t2;
841 <a name="l00933"></a>00933 
842 <a name="l00934"></a>00934     <span class="comment">/* (yb-yd) */</span>
843 <a name="l00935"></a>00935     t1 = pSrc[(2u * i1) + 1u] - pSrc[(2u * i3) + 1u];
844 <a name="l00936"></a>00936 
845 <a name="l00937"></a>00937     <span class="comment">/* (xb-xd) */</span>
846 <a name="l00938"></a>00938     t2 = pSrc[2u * i1] - pSrc[2u * i3];
847 <a name="l00939"></a>00939 
848 <a name="l00940"></a>00940     <span class="comment">/* xc&#39; = (xa-xb+xc-xd)co2 - (ya-yb+yc-yd)(si2) */</span>
849 <a name="l00941"></a>00941     pSrc[2u * i1] = r1 * onebyfftLen;
850 <a name="l00942"></a>00942 
851 <a name="l00943"></a>00943     <span class="comment">/* yc&#39; = (ya-yb+yc-yd)co2 + (xa-xb+xc-xd)(si2) */</span>
852 <a name="l00944"></a>00944     pSrc[(2u * i1) + 1u] = s1 * onebyfftLen;
853 <a name="l00945"></a>00945 
854 <a name="l00946"></a>00946 
855 <a name="l00947"></a>00947     <span class="comment">/* (xa - xc) - (yb-yd) */</span>
856 <a name="l00948"></a>00948     r1 = r2 - t1;
857 <a name="l00949"></a>00949 
858 <a name="l00950"></a>00950     <span class="comment">/* (xa - xc) + (yb-yd) */</span>
859 <a name="l00951"></a>00951     r2 = r2 + t1;
860 <a name="l00952"></a>00952 
861 <a name="l00953"></a>00953     <span class="comment">/* (ya - yc) + (xb-xd) */</span>
862 <a name="l00954"></a>00954     s1 = s2 + t2;
863 <a name="l00955"></a>00955 
864 <a name="l00956"></a>00956     <span class="comment">/* (ya - yc) - (xb-xd) */</span>
865 <a name="l00957"></a>00957     s2 = s2 - t2;
866 <a name="l00958"></a>00958 
867 <a name="l00959"></a>00959     <span class="comment">/* xb&#39; = (xa+yb-xc-yd)co1 - (ya-xb-yc+xd)(si1) */</span>
868 <a name="l00960"></a>00960     pSrc[2u * i2] = r1 * onebyfftLen;
869 <a name="l00961"></a>00961 
870 <a name="l00962"></a>00962     <span class="comment">/* yb&#39; = (ya-xb-yc+xd)co1 + (xa+yb-xc-yd)(si1) */</span>
871 <a name="l00963"></a>00963     pSrc[(2u * i2) + 1u] = s1 * onebyfftLen;
872 <a name="l00964"></a>00964 
873 <a name="l00965"></a>00965     <span class="comment">/* xd&#39; = (xa-yb-xc+yd)co3 - (ya+xb-yc-xd)(si3) */</span>
874 <a name="l00966"></a>00966     pSrc[2u * i3] = r2 * onebyfftLen;
875 <a name="l00967"></a>00967 
876 <a name="l00968"></a>00968     <span class="comment">/* yd&#39; = (ya+xb-yc-xd)co3 + (xa-yb-xc+yd)(si3) */</span>
877 <a name="l00969"></a>00969     pSrc[(2u * i3) + 1u] = s2 * onebyfftLen;
878 <a name="l00970"></a>00970   }
879 <a name="l00971"></a>00971 
880 <a name="l00972"></a>00972 
881 <a name="l00973"></a>00973 <span class="preprocessor">#else</span>
882 <a name="l00974"></a>00974 <span class="preprocessor"></span>
883 <a name="l00975"></a>00975   <span class="comment">/* Run the below code for Cortex-M0 */</span>
884 <a name="l00976"></a>00976 
885 <a name="l00977"></a>00977   <span class="comment">/*  Initializations for the first stage */</span>
886 <a name="l00978"></a>00978   n2 = fftLen;
887 <a name="l00979"></a>00979   n1 = n2;
888 <a name="l00980"></a>00980 
889 <a name="l00981"></a>00981   <span class="comment">/*  Calculation of first stage */</span>
890 <a name="l00982"></a>00982   <span class="keywordflow">for</span> (k = fftLen; k &gt; 4u; k &gt;&gt;= 2u)
891 <a name="l00983"></a>00983   {
892 <a name="l00984"></a>00984     <span class="comment">/*  Initializations for the first stage */</span>
893 <a name="l00985"></a>00985     n1 = n2;
894 <a name="l00986"></a>00986     n2 &gt;&gt;= 2u;
895 <a name="l00987"></a>00987     ia1 = 0u;
896 <a name="l00988"></a>00988 
897 <a name="l00989"></a>00989     <span class="comment">/*  Calculation of first stage */</span>
898 <a name="l00990"></a>00990     <span class="keywordflow">for</span> (j = 0u; j &lt;= (n2 - 1u); j++)
899 <a name="l00991"></a>00991     {
900 <a name="l00992"></a>00992       <span class="comment">/*  index calculation for the coefficients */</span>
901 <a name="l00993"></a>00993       ia2 = ia1 + ia1;
902 <a name="l00994"></a>00994       ia3 = ia2 + ia1;
903 <a name="l00995"></a>00995       co1 = pCoef[ia1 * 2u];
904 <a name="l00996"></a>00996       si1 = pCoef[(ia1 * 2u) + 1u];
905 <a name="l00997"></a>00997       co2 = pCoef[ia2 * 2u];
906 <a name="l00998"></a>00998       si2 = pCoef[(ia2 * 2u) + 1u];
907 <a name="l00999"></a>00999       co3 = pCoef[ia3 * 2u];
908 <a name="l01000"></a>01000       si3 = pCoef[(ia3 * 2u) + 1u];
909 <a name="l01001"></a>01001 
910 <a name="l01002"></a>01002       <span class="comment">/*  Twiddle coefficients index modifier */</span>
911 <a name="l01003"></a>01003       ia1 = ia1 + twidCoefModifier;
912 <a name="l01004"></a>01004 
913 <a name="l01005"></a>01005       <span class="keywordflow">for</span> (i0 = j; i0 &lt; fftLen; i0 += n1)
914 <a name="l01006"></a>01006       {
915 <a name="l01007"></a>01007         <span class="comment">/*  index calculation for the input as, */</span>
916 <a name="l01008"></a>01008         <span class="comment">/*  pSrc[i0 + 0], pSrc[i0 + fftLen/4], pSrc[i0 + fftLen/2], pSrc[i0 + 3fftLen/4] */</span>
917 <a name="l01009"></a>01009         i1 = i0 + n2;
918 <a name="l01010"></a>01010         i2 = i1 + n2;
919 <a name="l01011"></a>01011         i3 = i2 + n2;
920 <a name="l01012"></a>01012 
921 <a name="l01013"></a>01013         <span class="comment">/* xa + xc */</span>
922 <a name="l01014"></a>01014         r1 = pSrc[(2u * i0)] + pSrc[(2u * i2)];
923 <a name="l01015"></a>01015 
924 <a name="l01016"></a>01016         <span class="comment">/* xa - xc */</span>
925 <a name="l01017"></a>01017         r2 = pSrc[(2u * i0)] - pSrc[(2u * i2)];
926 <a name="l01018"></a>01018 
927 <a name="l01019"></a>01019         <span class="comment">/* ya + yc */</span>
928 <a name="l01020"></a>01020         s1 = pSrc[(2u * i0) + 1u] + pSrc[(2u * i2) + 1u];
929 <a name="l01021"></a>01021 
930 <a name="l01022"></a>01022         <span class="comment">/* ya - yc */</span>
931 <a name="l01023"></a>01023         s2 = pSrc[(2u * i0) + 1u] - pSrc[(2u * i2) + 1u];
932 <a name="l01024"></a>01024 
933 <a name="l01025"></a>01025         <span class="comment">/* xb + xd */</span>
934 <a name="l01026"></a>01026         t1 = pSrc[2u * i1] + pSrc[2u * i3];
935 <a name="l01027"></a>01027 
936 <a name="l01028"></a>01028         <span class="comment">/* xa&#39; = xa + xb + xc + xd */</span>
937 <a name="l01029"></a>01029         pSrc[2u * i0] = r1 + t1;
938 <a name="l01030"></a>01030 
939 <a name="l01031"></a>01031         <span class="comment">/* xa + xc -(xb + xd) */</span>
940 <a name="l01032"></a>01032         r1 = r1 - t1;
941 <a name="l01033"></a>01033 
942 <a name="l01034"></a>01034         <span class="comment">/* yb + yd */</span>
943 <a name="l01035"></a>01035         t2 = pSrc[(2u * i1) + 1u] + pSrc[(2u * i3) + 1u];
944 <a name="l01036"></a>01036 
945 <a name="l01037"></a>01037         <span class="comment">/* ya&#39; = ya + yb + yc + yd */</span>
946 <a name="l01038"></a>01038         pSrc[(2u * i0) + 1u] = s1 + t2;
947 <a name="l01039"></a>01039 
948 <a name="l01040"></a>01040         <span class="comment">/* (ya + yc) - (yb + yd) */</span>
949 <a name="l01041"></a>01041         s1 = s1 - t2;
950 <a name="l01042"></a>01042 
951 <a name="l01043"></a>01043         <span class="comment">/* (yb - yd) */</span>
952 <a name="l01044"></a>01044         t1 = pSrc[(2u * i1) + 1u] - pSrc[(2u * i3) + 1u];
953 <a name="l01045"></a>01045 
954 <a name="l01046"></a>01046         <span class="comment">/* (xb - xd) */</span>
955 <a name="l01047"></a>01047         t2 = pSrc[2u * i1] - pSrc[2u * i3];
956 <a name="l01048"></a>01048 
957 <a name="l01049"></a>01049         <span class="comment">/* xc&#39; = (xa-xb+xc-xd)co2 - (ya-yb+yc-yd)(si2) */</span>
958 <a name="l01050"></a>01050         pSrc[2u * i1] = (r1 * co2) - (s1 * si2);
959 <a name="l01051"></a>01051 
960 <a name="l01052"></a>01052         <span class="comment">/* yc&#39; = (ya-yb+yc-yd)co2 + (xa-xb+xc-xd)(si2) */</span>
961 <a name="l01053"></a>01053         pSrc[(2u * i1) + 1u] = (s1 * co2) + (r1 * si2);
962 <a name="l01054"></a>01054 
963 <a name="l01055"></a>01055         <span class="comment">/* (xa - xc) - (yb - yd) */</span>
964 <a name="l01056"></a>01056         r1 = r2 - t1;
965 <a name="l01057"></a>01057 
966 <a name="l01058"></a>01058         <span class="comment">/* (xa - xc) + (yb - yd) */</span>
967 <a name="l01059"></a>01059         r2 = r2 + t1;
968 <a name="l01060"></a>01060 
969 <a name="l01061"></a>01061         <span class="comment">/* (ya - yc) +  (xb - xd) */</span>
970 <a name="l01062"></a>01062         s1 = s2 + t2;
971 <a name="l01063"></a>01063 
972 <a name="l01064"></a>01064         <span class="comment">/* (ya - yc) -  (xb - xd) */</span>
973 <a name="l01065"></a>01065         s2 = s2 - t2;
974 <a name="l01066"></a>01066 
975 <a name="l01067"></a>01067         <span class="comment">/* xb&#39; = (xa+yb-xc-yd)co1 - (ya-xb-yc+xd)(si1) */</span>
976 <a name="l01068"></a>01068         pSrc[2u * i2] = (r1 * co1) - (s1 * si1);
977 <a name="l01069"></a>01069 
978 <a name="l01070"></a>01070         <span class="comment">/* yb&#39; = (ya-xb-yc+xd)co1 + (xa+yb-xc-yd)(si1) */</span>
979 <a name="l01071"></a>01071         pSrc[(2u * i2) + 1u] = (s1 * co1) + (r1 * si1);
980 <a name="l01072"></a>01072 
981 <a name="l01073"></a>01073         <span class="comment">/* xd&#39; = (xa-yb-xc+yd)co3 - (ya+xb-yc-xd)(si3) */</span>
982 <a name="l01074"></a>01074         pSrc[2u * i3] = (r2 * co3) - (s2 * si3);
983 <a name="l01075"></a>01075 
984 <a name="l01076"></a>01076         <span class="comment">/* yd&#39; = (ya+xb-yc-xd)co3 + (xa-yb-xc+yd)(si3) */</span>
985 <a name="l01077"></a>01077         pSrc[(2u * i3) + 1u] = (s2 * co3) + (r2 * si3);
986 <a name="l01078"></a>01078       }
987 <a name="l01079"></a>01079     }
988 <a name="l01080"></a>01080     twidCoefModifier &lt;&lt;= 2u;
989 <a name="l01081"></a>01081   }
990 <a name="l01082"></a>01082   <span class="comment">/*  Initializations of last stage */</span>
991 <a name="l01083"></a>01083   n1 = n2;
992 <a name="l01084"></a>01084   n2 &gt;&gt;= 2u;
993 <a name="l01085"></a>01085 
994 <a name="l01086"></a>01086   <span class="comment">/*  Calculations of last stage */</span>
995 <a name="l01087"></a>01087   <span class="keywordflow">for</span> (i0 = 0u; i0 &lt;= (fftLen - n1); i0 += n1)
996 <a name="l01088"></a>01088   {
997 <a name="l01089"></a>01089     <span class="comment">/*  index calculation for the input as, */</span>
998 <a name="l01090"></a>01090     <span class="comment">/*  pSrc[i0 + 0], pSrc[i0 + fftLen/4], pSrc[i0 + fftLen/2], pSrc[i0 + 3fftLen/4] */</span>
999 <a name="l01091"></a>01091     i1 = i0 + n2;
1000 <a name="l01092"></a>01092     i2 = i1 + n2;
1001 <a name="l01093"></a>01093     i3 = i2 + n2;
1002 <a name="l01094"></a>01094 
1003 <a name="l01095"></a>01095     <span class="comment">/*  Butterfly implementation */</span>
1004 <a name="l01096"></a>01096     <span class="comment">/* xa + xc */</span>
1005 <a name="l01097"></a>01097     r1 = pSrc[2u * i0] + pSrc[2u * i2];
1006 <a name="l01098"></a>01098 
1007 <a name="l01099"></a>01099     <span class="comment">/* xa - xc */</span>
1008 <a name="l01100"></a>01100     r2 = pSrc[2u * i0] - pSrc[2u * i2];
1009 <a name="l01101"></a>01101 
1010 <a name="l01102"></a>01102     <span class="comment">/* ya + yc */</span>
1011 <a name="l01103"></a>01103     s1 = pSrc[(2u * i0) + 1u] + pSrc[(2u * i2) + 1u];
1012 <a name="l01104"></a>01104 
1013 <a name="l01105"></a>01105     <span class="comment">/* ya - yc */</span>
1014 <a name="l01106"></a>01106     s2 = pSrc[(2u * i0) + 1u] - pSrc[(2u * i2) + 1u];
1015 <a name="l01107"></a>01107 
1016 <a name="l01108"></a>01108     <span class="comment">/* xc + xd */</span>
1017 <a name="l01109"></a>01109     t1 = pSrc[2u * i1] + pSrc[2u * i3];
1018 <a name="l01110"></a>01110 
1019 <a name="l01111"></a>01111     <span class="comment">/* xa&#39; = xa + xb + xc + xd */</span>
1020 <a name="l01112"></a>01112     pSrc[2u * i0] = (r1 + t1) * onebyfftLen;
1021 <a name="l01113"></a>01113 
1022 <a name="l01114"></a>01114     <span class="comment">/* (xa + xb) - (xc + xd) */</span>
1023 <a name="l01115"></a>01115     r1 = r1 - t1;
1024 <a name="l01116"></a>01116 
1025 <a name="l01117"></a>01117     <span class="comment">/* yb + yd */</span>
1026 <a name="l01118"></a>01118     t2 = pSrc[(2u * i1) + 1u] + pSrc[(2u * i3) + 1u];
1027 <a name="l01119"></a>01119 
1028 <a name="l01120"></a>01120     <span class="comment">/* ya&#39; = ya + yb + yc + yd */</span>
1029 <a name="l01121"></a>01121     pSrc[(2u * i0) + 1u] = (s1 + t2) * onebyfftLen;
1030 <a name="l01122"></a>01122 
1031 <a name="l01123"></a>01123     <span class="comment">/* (ya + yc) - (yb + yd) */</span>
1032 <a name="l01124"></a>01124     s1 = s1 - t2;
1033 <a name="l01125"></a>01125 
1034 <a name="l01126"></a>01126     <span class="comment">/* (yb-yd) */</span>
1035 <a name="l01127"></a>01127     t1 = pSrc[(2u * i1) + 1u] - pSrc[(2u * i3) + 1u];
1036 <a name="l01128"></a>01128 
1037 <a name="l01129"></a>01129     <span class="comment">/* (xb-xd) */</span>
1038 <a name="l01130"></a>01130     t2 = pSrc[2u * i1] - pSrc[2u * i3];
1039 <a name="l01131"></a>01131 
1040 <a name="l01132"></a>01132     <span class="comment">/* xc&#39; = (xa-xb+xc-xd)co2 - (ya-yb+yc-yd)(si2) */</span>
1041 <a name="l01133"></a>01133     pSrc[2u * i1] = r1 * onebyfftLen;
1042 <a name="l01134"></a>01134 
1043 <a name="l01135"></a>01135     <span class="comment">/* yc&#39; = (ya-yb+yc-yd)co2 + (xa-xb+xc-xd)(si2) */</span>
1044 <a name="l01136"></a>01136     pSrc[(2u * i1) + 1u] = s1 * onebyfftLen;
1045 <a name="l01137"></a>01137 
1046 <a name="l01138"></a>01138 
1047 <a name="l01139"></a>01139     <span class="comment">/* (xa - xc) - (yb-yd) */</span>
1048 <a name="l01140"></a>01140     r1 = r2 - t1;
1049 <a name="l01141"></a>01141 
1050 <a name="l01142"></a>01142     <span class="comment">/* (xa - xc) + (yb-yd) */</span>
1051 <a name="l01143"></a>01143     r2 = r2 + t1;
1052 <a name="l01144"></a>01144 
1053 <a name="l01145"></a>01145     <span class="comment">/* (ya - yc) + (xb-xd) */</span>
1054 <a name="l01146"></a>01146     s1 = s2 + t2;
1055 <a name="l01147"></a>01147 
1056 <a name="l01148"></a>01148     <span class="comment">/* (ya - yc) - (xb-xd) */</span>
1057 <a name="l01149"></a>01149     s2 = s2 - t2;
1058 <a name="l01150"></a>01150 
1059 <a name="l01151"></a>01151     <span class="comment">/* xb&#39; = (xa+yb-xc-yd)co1 - (ya-xb-yc+xd)(si1) */</span>
1060 <a name="l01152"></a>01152     pSrc[2u * i2] = r1 * onebyfftLen;
1061 <a name="l01153"></a>01153 
1062 <a name="l01154"></a>01154     <span class="comment">/* yb&#39; = (ya-xb-yc+xd)co1 + (xa+yb-xc-yd)(si1) */</span>
1063 <a name="l01155"></a>01155     pSrc[(2u * i2) + 1u] = s1 * onebyfftLen;
1064 <a name="l01156"></a>01156 
1065 <a name="l01157"></a>01157     <span class="comment">/* xd&#39; = (xa-yb-xc+yd)co3 - (ya+xb-yc-xd)(si3) */</span>
1066 <a name="l01158"></a>01158     pSrc[2u * i3] = r2 * onebyfftLen;
1067 <a name="l01159"></a>01159 
1068 <a name="l01160"></a>01160     <span class="comment">/* yd&#39; = (ya+xb-yc-xd)co3 + (xa-yb-xc+yd)(si3) */</span>
1069 <a name="l01161"></a>01161     pSrc[(2u * i3) + 1u] = s2 * onebyfftLen;
1070 <a name="l01162"></a>01162   }
1071 <a name="l01163"></a>01163 
1072 <a name="l01164"></a>01164 <span class="preprocessor">#endif </span><span class="comment">/* #ifndef ARM_MATH_CM0 */</span>
1073 <a name="l01165"></a>01165 
1074 <a name="l01166"></a>01166 }
1075 <a name="l01167"></a>01167 
1076 <a name="l01168"></a>01168 <span class="comment">/*   </span>
1077 <a name="l01169"></a>01169 <span class="comment"> * @brief  In-place bit reversal function.  </span>
1078 <a name="l01170"></a>01170 <span class="comment"> * @param[in, out] *pSrc        points to the in-place buffer of floating-point data type.  </span>
1079 <a name="l01171"></a>01171 <span class="comment"> * @param[in]      fftSize      length of the FFT.  </span>
1080 <a name="l01172"></a>01172 <span class="comment"> * @param[in]      bitRevFactor bit reversal modifier that supports different size FFTs with the same bit reversal table.  </span>
1081 <a name="l01173"></a>01173 <span class="comment"> * @param[in]      *pBitRevTab  points to the bit reversal table.  </span>
1082 <a name="l01174"></a>01174 <span class="comment"> * @return none.  </span>
1083 <a name="l01175"></a>01175 <span class="comment"> */</span>
1084 <a name="l01176"></a>01176 
1085 <a name="l01177"></a><a class="code" href="arm__math_8h.html#a3d4062fdfa6aaa3f51f41cab868e508b">01177</a> <span class="keywordtype">void</span> <a class="code" href="arm__cfft__radix4__f32_8c.html#a3d4062fdfa6aaa3f51f41cab868e508b" title="In-place bit reversal function.">arm_bitreversal_f32</a>(
1086 <a name="l01178"></a>01178   <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> * pSrc,
1087 <a name="l01179"></a>01179   uint16_t <a class="code" href="arm__fft__bin__example__f32_8c.html#a9b500899c581f6df3ffc0a9f3a9ef6aa">fftSize</a>,
1088 <a name="l01180"></a>01180   uint16_t bitRevFactor,
1089 <a name="l01181"></a>01181   uint16_t * pBitRevTab)
1090 <a name="l01182"></a>01182 {
1091 <a name="l01183"></a>01183   uint16_t fftLenBy2, fftLenBy2p1;
1092 <a name="l01184"></a>01184   uint16_t i, j;
1093 <a name="l01185"></a>01185   <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> in;
1094 <a name="l01186"></a>01186 
1095 <a name="l01187"></a>01187   <span class="comment">/*  Initializations */</span>
1096 <a name="l01188"></a>01188   j = 0u;
1097 <a name="l01189"></a>01189   fftLenBy2 = fftSize &gt;&gt; 1u;
1098 <a name="l01190"></a>01190   fftLenBy2p1 = (fftSize &gt;&gt; 1u) + 1u;
1099 <a name="l01191"></a>01191 
1100 <a name="l01192"></a>01192   <span class="comment">/* Bit Reversal Implementation */</span>
1101 <a name="l01193"></a>01193   <span class="keywordflow">for</span> (i = 0u; i &lt;= (fftLenBy2 - 2u); i += 2u)
1102 <a name="l01194"></a>01194   {
1103 <a name="l01195"></a>01195     <span class="keywordflow">if</span>(i &lt; j)
1104 <a name="l01196"></a>01196     {
1105 <a name="l01197"></a>01197       <span class="comment">/*  pSrc[i] &lt;-&gt; pSrc[j]; */</span>
1106 <a name="l01198"></a>01198       in = pSrc[2u * i];
1107 <a name="l01199"></a>01199       pSrc[2u * i] = pSrc[2u * j];
1108 <a name="l01200"></a>01200       pSrc[2u * j] = in;
1109 <a name="l01201"></a>01201 
1110 <a name="l01202"></a>01202       <span class="comment">/*  pSrc[i+1u] &lt;-&gt; pSrc[j+1u] */</span>
1111 <a name="l01203"></a>01203       in = pSrc[(2u * i) + 1u];
1112 <a name="l01204"></a>01204       pSrc[(2u * i) + 1u] = pSrc[(2u * j) + 1u];
1113 <a name="l01205"></a>01205       pSrc[(2u * j) + 1u] = in;
1114 <a name="l01206"></a>01206 
1115 <a name="l01207"></a>01207       <span class="comment">/*  pSrc[i+fftLenBy2p1] &lt;-&gt; pSrc[j+fftLenBy2p1] */</span>
1116 <a name="l01208"></a>01208       in = pSrc[2u * (i + fftLenBy2p1)];
1117 <a name="l01209"></a>01209       pSrc[2u * (i + fftLenBy2p1)] = pSrc[2u * (j + fftLenBy2p1)];
1118 <a name="l01210"></a>01210       pSrc[2u * (j + fftLenBy2p1)] = in;
1119 <a name="l01211"></a>01211 
1120 <a name="l01212"></a>01212       <span class="comment">/*  pSrc[i+fftLenBy2p1+1u] &lt;-&gt; pSrc[j+fftLenBy2p1+1u] */</span>
1121 <a name="l01213"></a>01213       in = pSrc[(2u * (i + fftLenBy2p1)) + 1u];
1122 <a name="l01214"></a>01214       pSrc[(2u * (i + fftLenBy2p1)) + 1u] =
1123 <a name="l01215"></a>01215         pSrc[(2u * (j + fftLenBy2p1)) + 1u];
1124 <a name="l01216"></a>01216       pSrc[(2u * (j + fftLenBy2p1)) + 1u] = in;
1125 <a name="l01217"></a>01217 
1126 <a name="l01218"></a>01218     }
1127 <a name="l01219"></a>01219 
1128 <a name="l01220"></a>01220     <span class="comment">/*  pSrc[i+1u] &lt;-&gt; pSrc[j+1u] */</span>
1129 <a name="l01221"></a>01221     in = pSrc[2u * (i + 1u)];
1130 <a name="l01222"></a>01222     pSrc[2u * (i + 1u)] = pSrc[2u * (j + fftLenBy2)];
1131 <a name="l01223"></a>01223     pSrc[2u * (j + fftLenBy2)] = in;
1132 <a name="l01224"></a>01224 
1133 <a name="l01225"></a>01225     <span class="comment">/*  pSrc[i+2u] &lt;-&gt; pSrc[j+2u] */</span>
1134 <a name="l01226"></a>01226     in = pSrc[(2u * (i + 1u)) + 1u];
1135 <a name="l01227"></a>01227     pSrc[(2u * (i + 1u)) + 1u] = pSrc[(2u * (j + fftLenBy2)) + 1u];
1136 <a name="l01228"></a>01228     pSrc[(2u * (j + fftLenBy2)) + 1u] = in;
1137 <a name="l01229"></a>01229 
1138 <a name="l01230"></a>01230     <span class="comment">/*  Reading the index for the bit reversal */</span>
1139 <a name="l01231"></a>01231     j = *pBitRevTab;
1140 <a name="l01232"></a>01232 
1141 <a name="l01233"></a>01233     <span class="comment">/*  Updating the bit reversal index depending on the fft length  */</span>
1142 <a name="l01234"></a>01234     pBitRevTab += bitRevFactor;
1143 <a name="l01235"></a>01235   }
1144 <a name="l01236"></a>01236 }
1145 </pre></div></div>
1146 </div>
1147 <!--- window showing the filter options -->
1148 <div id="MSearchSelectWindow"
1149      onmouseover="return searchBox.OnSearchSelectShow()"
1150      onmouseout="return searchBox.OnSearchSelectHide()"
1151      onkeydown="return searchBox.OnSearchSelectKey(event)">
1152 <a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(0)"><span class="SelectionMark">&#160;</span>All</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(1)"><span class="SelectionMark">&#160;</span>Data Structures</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(2)"><span class="SelectionMark">&#160;</span>Files</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(3)"><span class="SelectionMark">&#160;</span>Functions</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(4)"><span class="SelectionMark">&#160;</span>Variables</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(5)"><span class="SelectionMark">&#160;</span>Typedefs</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(6)"><span class="SelectionMark">&#160;</span>Enumerations</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(7)"><span class="SelectionMark">&#160;</span>Enumerator</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(8)"><span class="SelectionMark">&#160;</span>Defines</a></div>
1153
1154 <!-- iframe showing the search results (closed by default) -->
1155 <div id="MSearchResultsWindow">
1156 <iframe src="" frameborder="0" 
1157         name="MSearchResults" id="MSearchResults">
1158 </iframe>
1159 </div>
1160
1161 <hr class="footer"/><address class="footer"><small>Generated on Fri Jul 15 2011 13:16:15 for CMSIS DSP Software Library by&#160;
1162 <a href="http://www.doxygen.org/index.html">
1163 <img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.7.2 </small></address>
1164 </body>
1165 </html>