开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > 流媒体/Mpeg4/MP4 > 查看源码
fastfft.c
资源名称：NETVIDEO.rar　[点击查看]
上传用户：sun1608
上传日期：2007-02-02
资源大小：6116k
文件大小：82k
源码类别：

流媒体/Mpeg4/MP4

开发平台：

Visual C++

fastfft.c：源码内容
  1. /*
  2.  * FAAD - Freeware Advanced Audio Decoder
  3.  * Copyright (C) 2001 Menno Bakker
  4.  *
  5.  * This library is free software; you can redistribute it and/or
  6.  * modify it under the terms of the GNU Lesser General Public
  7.  * License as published by the Free Software Foundation; either
  8.  * version 2.1 of the License, or (at your option) any later version.
  9.  *
  10.  * This library is distributed in the hope that it will be useful,
  11.  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  12.  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13.  * Lesser General Public License for more details.
  15.  * You should have received a copy of the GNU Lesser General Public
  16.  * License along with this library; if not, write to the Free Software
  17.  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  18.  *
  19.  * $Id: fastfft.c,v 1.4 2002/01/09 22:25:41 wmay Exp $
  20.  */
  22. #include "fastfft.h"
  23. #include "all.h"
  25. #ifdef _MSC_VER
  26. #pragma warning(disable:4305)
  27. #endif
  30. static fftw_real K980785280[1] =
  31. {FFTW_KONST(+0.980785280403230449126182236134239036973933731)};
  32. static fftw_real K195090322[1] =
  33. {FFTW_KONST(+0.195090322016128267848284868477022240927691618)};
  34. static fftw_real K555570233[1] =
  35. {FFTW_KONST(+0.555570233019602224742830813948532874374937191)};
  36. static fftw_real K831469612[1] =
  37. {FFTW_KONST(+0.831469612302545237078788377617905756738560812)};
  38. static fftw_real K923879532[1] =
  39. {FFTW_KONST(+0.923879532511286756128183189396788286822416626)};
  40. static fftw_real K382683432[1] =
  41. {FFTW_KONST(+0.382683432365089771728459984030398866761344562)};
  42. static fftw_real K707106781[1] =
  43. {FFTW_KONST(+0.707106781186547524400844362104849039284835938)};
  45. static fftw_complex PFFTW(W_64)[30] = {
  46. { 0.995184726672197, 0.0980171403295606 },
  47. { 0.98078528040323, 0.195090322016128 },
  48. { 0.98078528040323, 0.195090322016128 },
  49. { 0.923879532511287, 0.38268343236509 },
  50. { 0.956940335732209, 0.290284677254462 },
  51. { 0.831469612302545, 0.555570233019602 },
  52. { 0.923879532511287, 0.38268343236509 },
  53. { 0.707106781186548, 0.707106781186547 },
  54. { 0.881921264348355, 0.471396736825998 },
  55. { 0.555570233019602, 0.831469612302545 },
  56. { 0.831469612302545, 0.555570233019602 },
  57. { 0.38268343236509, 0.923879532511287 },
  58. { 0.773010453362737, 0.634393284163645 },
  59. { 0.195090322016128, 0.98078528040323 },
  60. { 0.707106781186548, 0.707106781186547 },
  61. { 6.12303176911189e-17, 1 },
  62. { 0.634393284163645, 0.773010453362737 },
  63. { -0.195090322016128, 0.98078528040323 },
  64. { 0.555570233019602, 0.831469612302545 },
  65. { -0.38268343236509, 0.923879532511287 },
  66. { 0.471396736825998, 0.881921264348355 },
  67. { -0.555570233019602, 0.831469612302545 },
  68. { 0.38268343236509, 0.923879532511287 },
  69. { -0.707106781186547, 0.707106781186548 },
  70. { 0.290284677254462, 0.956940335732209 },
  71. { -0.831469612302545, 0.555570233019602 },
  72. { 0.195090322016128, 0.98078528040323 },
  73. { -0.923879532511287, 0.38268343236509 },
  74. { 0.0980171403295608, 0.995184726672197 },
  75. { -0.98078528040323, 0.195090322016129 },
  76. };
  78. static fftw_complex PFFTW(W_128)[62] = {
  79. { 0.998795456205172, 0.049067674327418 },
  80. { 0.995184726672197, 0.0980171403295606 },
  81. { 0.995184726672197, 0.0980171403295606 },
  82. { 0.98078528040323, 0.195090322016128 },
  83. { 0.989176509964781, 0.146730474455362 },
  84. { 0.956940335732209, 0.290284677254462 },
  85. { 0.98078528040323, 0.195090322016128 },
  86. { 0.923879532511287, 0.38268343236509 },
  87. { 0.970031253194544, 0.242980179903264 },
  88. { 0.881921264348355, 0.471396736825998 },
  89. { 0.956940335732209, 0.290284677254462 },
  90. { 0.831469612302545, 0.555570233019602 },
  91. { 0.941544065183021, 0.33688985339222 },
  92. { 0.773010453362737, 0.634393284163645 },
  93. { 0.923879532511287, 0.38268343236509 },
  94. { 0.707106781186548, 0.707106781186547 },
  95. { 0.903989293123443, 0.427555093430282 },
  96. { 0.634393284163645, 0.773010453362737 },
  97. { 0.881921264348355, 0.471396736825998 },
  98. { 0.555570233019602, 0.831469612302545 },
  99. { 0.857728610000272, 0.514102744193222 },
  100. { 0.471396736825998, 0.881921264348355 },
  101. { 0.831469612302545, 0.555570233019602 },
  102. { 0.38268343236509, 0.923879532511287 },
  103. { 0.803207531480645, 0.595699304492433 },
  104. { 0.290284677254462, 0.956940335732209 },
  105. { 0.773010453362737, 0.634393284163645 },
  106. { 0.195090322016128, 0.98078528040323 },
  107. { 0.740951125354959, 0.671558954847018 },
  108. { 0.0980171403295608, 0.995184726672197 },
  109. { 0.707106781186548, 0.707106781186547 },
  110. { 6.12303176911189e-17, 1 },
  111. { 0.671558954847018, 0.740951125354959 },
  112. { -0.0980171403295606, 0.995184726672197 },
  113. { 0.634393284163645, 0.773010453362737 },
  114. { -0.195090322016128, 0.98078528040323 },
  115. { 0.595699304492433, 0.803207531480645 },
  116. { -0.290284677254462, 0.956940335732209 },
  117. { 0.555570233019602, 0.831469612302545 },
  118. { -0.38268343236509, 0.923879532511287 },
  119. { 0.514102744193222, 0.857728610000272 },
  120. { -0.471396736825998, 0.881921264348355 },
  121. { 0.471396736825998, 0.881921264348355 },
  122. { -0.555570233019602, 0.831469612302545 },
  123. { 0.427555093430282, 0.903989293123443 },
  124. { -0.634393284163645, 0.773010453362737 },
  125. { 0.38268343236509, 0.923879532511287 },
  126. { -0.707106781186547, 0.707106781186548 },
  127. { 0.33688985339222, 0.941544065183021 },
  128. { -0.773010453362737, 0.634393284163645 },
  129. { 0.290284677254462, 0.956940335732209 },
  130. { -0.831469612302545, 0.555570233019602 },
  131. { 0.242980179903264, 0.970031253194544 },
  132. { -0.881921264348355, 0.471396736825998 },
  133. { 0.195090322016128, 0.98078528040323 },
  134. { -0.923879532511287, 0.38268343236509 },
  135. { 0.146730474455362, 0.989176509964781 },
  136. { -0.956940335732209, 0.290284677254462 },
  137. { 0.0980171403295608, 0.995184726672197 },
  138. { -0.98078528040323, 0.195090322016129 },
  139. { 0.0490676743274181, 0.998795456205172 },
  140. { -0.995184726672197, 0.0980171403295608 },
  141. };
  143. static fftw_complex PFFTW(W_512)[254] = {
  144. { 0.999924701839145, 0.0122715382857199 },
  145. { 0.999698818696204, 0.0245412285229123 },
  146. { 0.999698818696204, 0.0245412285229123 },
  147. { 0.998795456205172, 0.049067674327418 },
  148. { 0.99932238458835, 0.0368072229413588 },
  149. { 0.99729045667869, 0.0735645635996674 },
  150. { 0.998795456205172, 0.049067674327418 },
  151. { 0.995184726672197, 0.0980171403295606 },
  152. { 0.998118112900149, 0.0613207363022086 },
  153. { 0.99247953459871, 0.122410675199216 },
  154. { 0.99729045667869, 0.0735645635996674 },
  155. { 0.989176509964781, 0.146730474455362 },
  156. { 0.996312612182778, 0.0857973123444399 },
  157. { 0.985277642388941, 0.170961888760301 },
  158. { 0.995184726672197, 0.0980171403295606 },
  159. { 0.98078528040323, 0.195090322016128 },
  160. { 0.993906970002356, 0.110222207293883 },
  161. { 0.975702130038529, 0.21910124015687 },
  162. { 0.99247953459871, 0.122410675199216 },
  163. { 0.970031253194544, 0.242980179903264 },
  164. { 0.99090263542778, 0.134580708507126 },
  165. { 0.96377606579544, 0.266712757474898 },
  166. { 0.989176509964781, 0.146730474455362 },
  167. { 0.956940335732209, 0.290284677254462 },
  168. { 0.987301418157858, 0.158858143333861 },
  169. { 0.949528180593037, 0.313681740398892 },
  170. { 0.985277642388941, 0.170961888760301 },
  171. { 0.941544065183021, 0.33688985339222 },
  172. { 0.983105487431216, 0.183039887955141 },
  173. { 0.932992798834739, 0.359895036534988 },
  174. { 0.98078528040323, 0.195090322016128 },
  175. { 0.923879532511287, 0.38268343236509 },
  176. { 0.978317370719628, 0.207111376192219 },
  177. { 0.914209755703531, 0.40524131400499 },
  178. { 0.975702130038529, 0.21910124015687 },
  179. { 0.903989293123443, 0.427555093430282 },
  180. { 0.97293995220556, 0.231058108280671 },
  181. { 0.893224301195515, 0.449611329654607 },
  182. { 0.970031253194544, 0.242980179903264 },
  183. { 0.881921264348355, 0.471396736825998 },
  184. { 0.966976471044852, 0.254865659604515 },
  185. { 0.870086991108711, 0.492898192229784 },
  186. { 0.96377606579544, 0.266712757474898 },
  187. { 0.857728610000272, 0.514102744193222 },
  188. { 0.960430519415566, 0.278519689385053 },
  189. { 0.844853565249707, 0.534997619887097 },
  190. { 0.956940335732209, 0.290284677254462 },
  191. { 0.831469612302545, 0.555570233019602 },
  192. { 0.953306040354194, 0.302005949319228 },
  193. { 0.817584813151584, 0.575808191417845 },
  194. { 0.949528180593037, 0.313681740398892 },
  195. { 0.803207531480645, 0.595699304492433 },
  196. { 0.945607325380521, 0.325310292162263 },
  197. { 0.788346427626606, 0.615231590580627 },
  198. { 0.941544065183021, 0.33688985339222 },
  199. { 0.773010453362737, 0.634393284163645 },
  200. { 0.937339011912575, 0.348418680249435 },
  201. { 0.757208846506485, 0.653172842953777 },
  202. { 0.932992798834739, 0.359895036534988 },
  203. { 0.740951125354959, 0.671558954847018 },
  204. { 0.928506080473216, 0.371317193951838 },
  205. { 0.724247082951467, 0.689540544737067 },
  206. { 0.923879532511287, 0.38268343236509 },
  207. { 0.707106781186548, 0.707106781186547 },
  208. { 0.919113851690058, 0.393992040061048 },
  209. { 0.689540544737067, 0.724247082951467 },
  210. { 0.914209755703531, 0.40524131400499 },
  211. { 0.671558954847018, 0.740951125354959 },
  212. { 0.909167983090522, 0.416429560097637 },
  213. { 0.653172842953777, 0.757208846506484 },
  214. { 0.903989293123443, 0.427555093430282 },
  215. { 0.634393284163645, 0.773010453362737 },
  216. { 0.898674465693954, 0.438616238538528 },
  217. { 0.615231590580627, 0.788346427626606 },
  218. { 0.893224301195515, 0.449611329654607 },
  219. { 0.595699304492433, 0.803207531480645 },
  220. { 0.887639620402854, 0.46053871095824 },
  221. { 0.575808191417845, 0.817584813151584 },
  222. { 0.881921264348355, 0.471396736825998 },
  223. { 0.555570233019602, 0.831469612302545 },
  224. { 0.876070094195407, 0.482183772079123 },
  225. { 0.534997619887097, 0.844853565249707 },
  226. { 0.870086991108711, 0.492898192229784 },
  227. { 0.514102744193222, 0.857728610000272 },
  228. { 0.863972856121587, 0.503538383725718 },
  229. { 0.492898192229784, 0.870086991108711 },
  230. { 0.857728610000272, 0.514102744193222 },
  231. { 0.471396736825998, 0.881921264348355 },
  232. { 0.851355193105265, 0.524589682678469 },
  233. { 0.449611329654607, 0.893224301195515 },
  234. { 0.844853565249707, 0.534997619887097 },
  235. { 0.427555093430282, 0.903989293123443 },
  236. { 0.838224705554838, 0.545324988422046 },
  237. { 0.40524131400499, 0.914209755703531 },
  238. { 0.831469612302545, 0.555570233019602 },
  239. { 0.38268343236509, 0.923879532511287 },
  240. { 0.824589302785025, 0.565731810783613 },
  241. { 0.359895036534988, 0.932992798834739 },
  242. { 0.817584813151584, 0.575808191417845 },
  243. { 0.33688985339222, 0.941544065183021 },
  244. { 0.810457198252595, 0.585797857456439 },
  245. { 0.313681740398892, 0.949528180593037 },
  246. { 0.803207531480645, 0.595699304492433 },
  247. { 0.290284677254462, 0.956940335732209 },
  248. { 0.795836904608884, 0.605511041404326 },
  249. { 0.266712757474898, 0.96377606579544 },
  250. { 0.788346427626606, 0.615231590580627 },
  251. { 0.242980179903264, 0.970031253194544 },
  252. { 0.780737228572094, 0.624859488142386 },
  253. { 0.21910124015687, 0.975702130038529 },
  254. { 0.773010453362737, 0.634393284163645 },
  255. { 0.195090322016128, 0.98078528040323 },
  256. { 0.765167265622459, 0.643831542889791 },
  257. { 0.170961888760301, 0.985277642388941 },
  258. { 0.757208846506485, 0.653172842953777 },
  259. { 0.146730474455362, 0.989176509964781 },
  260. { 0.749136394523459, 0.662415777590172 },
  261. { 0.122410675199216, 0.99247953459871 },
  262. { 0.740951125354959, 0.671558954847018 },
  263. { 0.0980171403295608, 0.995184726672197 },
  264. { 0.732654271672413, 0.680600997795453 },
  265. { 0.0735645635996675, 0.99729045667869 },
  266. { 0.724247082951467, 0.689540544737067 },
  267. { 0.0490676743274181, 0.998795456205172 },
  268. { 0.715730825283819, 0.698376249408973 },
  269. { 0.0245412285229123, 0.999698818696204 },
  270. { 0.707106781186548, 0.707106781186547 },
  271. { 6.12303176911189e-17, 1 },
  272. { 0.698376249408973, 0.715730825283819 },
  273. { -0.0245412285229121, 0.999698818696204 },
  274. { 0.689540544737067, 0.724247082951467 },
  275. { -0.049067674327418, 0.998795456205172 },
  276. { 0.680600997795453, 0.732654271672413 },
  277. { -0.0735645635996673, 0.99729045667869 },
  278. { 0.671558954847018, 0.740951125354959 },
  279. { -0.0980171403295606, 0.995184726672197 },
  280. { 0.662415777590172, 0.749136394523459 },
  281. { -0.122410675199216, 0.99247953459871 },
  282. { 0.653172842953777, 0.757208846506484 },
  283. { -0.146730474455362, 0.989176509964781 },
  284. { 0.643831542889791, 0.765167265622459 },
  285. { -0.170961888760301, 0.985277642388941 },
  286. { 0.634393284163645, 0.773010453362737 },
  287. { -0.195090322016128, 0.98078528040323 },
  288. { 0.624859488142386, 0.780737228572094 },
  289. { -0.21910124015687, 0.975702130038529 },
  290. { 0.615231590580627, 0.788346427626606 },
  291. { -0.242980179903264, 0.970031253194544 },
  292. { 0.605511041404326, 0.795836904608883 },
  293. { -0.266712757474898, 0.96377606579544 },
  294. { 0.595699304492433, 0.803207531480645 },
  295. { -0.290284677254462, 0.956940335732209 },
  296. { 0.585797857456439, 0.810457198252595 },
  297. { -0.313681740398891, 0.949528180593037 },
  298. { 0.575808191417845, 0.817584813151584 },
  299. { -0.33688985339222, 0.941544065183021 },
  300. { 0.565731810783613, 0.824589302785025 },
  301. { -0.359895036534988, 0.932992798834739 },
  302. { 0.555570233019602, 0.831469612302545 },
  303. { -0.38268343236509, 0.923879532511287 },
  304. { 0.545324988422046, 0.838224705554838 },
  305. { -0.40524131400499, 0.914209755703531 },
  306. { 0.534997619887097, 0.844853565249707 },
  307. { -0.427555093430282, 0.903989293123443 },
  308. { 0.524589682678469, 0.851355193105265 },
  309. { -0.449611329654607, 0.893224301195515 },
  310. { 0.514102744193222, 0.857728610000272 },
  311. { -0.471396736825998, 0.881921264348355 },
  312. { 0.503538383725718, 0.863972856121587 },
  313. { -0.492898192229784, 0.870086991108711 },
  314. { 0.492898192229784, 0.870086991108711 },
  315. { -0.514102744193222, 0.857728610000272 },
  316. { 0.482183772079123, 0.876070094195407 },
  317. { -0.534997619887097, 0.844853565249707 },
  318. { 0.471396736825998, 0.881921264348355 },
  319. { -0.555570233019602, 0.831469612302545 },
  320. { 0.46053871095824, 0.887639620402854 },
  321. { -0.575808191417845, 0.817584813151584 },
  322. { 0.449611329654607, 0.893224301195515 },
  323. { -0.595699304492433, 0.803207531480645 },
  324. { 0.438616238538528, 0.898674465693954 },
  325. { -0.615231590580627, 0.788346427626606 },
  326. { 0.427555093430282, 0.903989293123443 },
  327. { -0.634393284163645, 0.773010453362737 },
  328. { 0.416429560097637, 0.909167983090522 },
  329. { -0.653172842953777, 0.757208846506485 },
  330. { 0.40524131400499, 0.914209755703531 },
  331. { -0.671558954847018, 0.740951125354959 },
  332. { 0.393992040061048, 0.919113851690058 },
  333. { -0.689540544737067, 0.724247082951467 },
  334. { 0.38268343236509, 0.923879532511287 },
  335. { -0.707106781186547, 0.707106781186548 },
  336. { 0.371317193951838, 0.928506080473215 },
  337. { -0.724247082951467, 0.689540544737067 },
  338. { 0.359895036534988, 0.932992798834739 },
  339. { -0.740951125354959, 0.671558954847019 },
  340. { 0.348418680249435, 0.937339011912575 },
  341. { -0.757208846506485, 0.653172842953777 },
  342. { 0.33688985339222, 0.941544065183021 },
  343. { -0.773010453362737, 0.634393284163645 },
  344. { 0.325310292162263, 0.945607325380521 },
  345. { -0.788346427626606, 0.615231590580627 },
  346. { 0.313681740398892, 0.949528180593037 },
  347. { -0.803207531480645, 0.595699304492433 },
  348. { 0.302005949319228, 0.953306040354194 },
  349. { -0.817584813151584, 0.575808191417845 },
  350. { 0.290284677254462, 0.956940335732209 },
  351. { -0.831469612302545, 0.555570233019602 },
  352. { 0.278519689385053, 0.960430519415566 },
  353. { -0.844853565249707, 0.534997619887097 },
  354. { 0.266712757474898, 0.96377606579544 },
  355. { -0.857728610000272, 0.514102744193222 },
  356. { 0.254865659604515, 0.966976471044852 },
  357. { -0.870086991108711, 0.492898192229784 },
  358. { 0.242980179903264, 0.970031253194544 },
  359. { -0.881921264348355, 0.471396736825998 },
  360. { 0.231058108280671, 0.97293995220556 },
  361. { -0.893224301195515, 0.449611329654607 },
  362. { 0.21910124015687, 0.975702130038529 },
  363. { -0.903989293123443, 0.427555093430282 },
  364. { 0.207111376192219, 0.978317370719628 },
  365. { -0.914209755703531, 0.40524131400499 },
  366. { 0.195090322016128, 0.98078528040323 },
  367. { -0.923879532511287, 0.38268343236509 },
  368. { 0.183039887955141, 0.983105487431216 },
  369. { -0.932992798834739, 0.359895036534988 },
  370. { 0.170961888760301, 0.985277642388941 },
  371. { -0.941544065183021, 0.33688985339222 },
  372. { 0.158858143333861, 0.987301418157858 },
  373. { -0.949528180593037, 0.313681740398891 },
  374. { 0.146730474455362, 0.989176509964781 },
  375. { -0.956940335732209, 0.290284677254462 },
  376. { 0.134580708507126, 0.99090263542778 },
  377. { -0.96377606579544, 0.266712757474898 },
  378. { 0.122410675199216, 0.99247953459871 },
  379. { -0.970031253194544, 0.242980179903264 },
  380. { 0.110222207293883, 0.993906970002356 },
  381. { -0.975702130038528, 0.21910124015687 },
  382. { 0.0980171403295608, 0.995184726672197 },
  383. { -0.98078528040323, 0.195090322016129 },
  384. { 0.0857973123444399, 0.996312612182778 },
  385. { -0.985277642388941, 0.170961888760301 },
  386. { 0.0735645635996675, 0.99729045667869 },
  387. { -0.989176509964781, 0.146730474455362 },
  388. { 0.0613207363022086, 0.998118112900149 },
  389. { -0.99247953459871, 0.122410675199216 },
  390. { 0.0490676743274181, 0.998795456205172 },
  391. { -0.995184726672197, 0.0980171403295608 },
  392. { 0.036807222941359, 0.99932238458835 },
  393. { -0.99729045667869, 0.0735645635996677 },
  394. { 0.0245412285229123, 0.999698818696204 },
  395. { -0.998795456205172, 0.049067674327418 },
  396. { 0.0122715382857199, 0.999924701839145 },
  397. { -0.999698818696204, 0.0245412285229123 },
  398. };
  399. ///////////////////////////////////////////////////////////////
  401. void PFFTW(16) (fftw_complex * input) {
  402.      fftw_real tmp332;
  403.      fftw_real tmp331;
  404.      fftw_real tmp330;
  405.      fftw_real tmp329;
  406.      fftw_real tmp328;
  407.      fftw_real tmp327;
  408.      fftw_real tmp326;
  409.      fftw_real tmp325;
  410.      fftw_real tmp324;
  411.      fftw_real tmp323;
  412.      fftw_real tmp322;
  413.      fftw_real tmp321;
  414.      fftw_real tmp320;
  415.      fftw_real tmp319;
  416.      fftw_real tmp318;
  417.      fftw_real tmp317;
  418.      fftw_real tmp316;
  419.      fftw_real tmp315;
  420.      fftw_real tmp314;
  421.      fftw_real tmp313;
  422.      fftw_real tmp312;
  423.      fftw_real tmp311;
  424.      fftw_real tmp310;
  425.      fftw_real tmp309;
  426.      fftw_real tmp308;
  427.      fftw_real tmp307;
  428.      fftw_real tmp306;
  429.      fftw_real tmp305;
  430.      fftw_real tmp304;
  431.      fftw_real tmp303;
  432.      fftw_real tmp302;
  433.      fftw_real tmp301;
  434.      fftw_real st1;
  435.      fftw_real st2;
  436.      fftw_real st3;
  437.      fftw_real st4;
  438.      fftw_real st5;
  439.      fftw_real st6;
  440.      fftw_real st7;
  441.      fftw_real st8;
  442.      st8 = c_re(input[0]);
  443.      st8 = st8 - c_re(input[8]);
  444.      st7 = c_im(input[4]);
  445.      st7 = st7 - c_im(input[12]);
  446.      st6 = c_re(input[4]);
  447.      st5 = st8 - st7;
  448.      st8 = st8 + st7;
  449.      st6 = st6 - c_re(input[12]);
  450.      st4 = c_im(input[0]);
  451.      st4 = st4 - c_im(input[8]);
  452.      st3 = c_re(input[0]);
  453.      st2 = st6 + st4;
  454.      st4 = st4 - st6;
  455.      st3 = st3 + c_re(input[8]);
  456.      st1 = c_re(input[4]);
  457.      st1 = st1 + c_re(input[12]);
  458.      st7 = c_im(input[0]);
  459.      st6 = st3 + st1;
  460.      st3 = st3 - st1;
  461.      st7 = st7 + c_im(input[8]);
  462.      st1 = c_im(input[4]);
  463.      st1 = st1 + c_im(input[12]);
  464.      tmp301 = st4;
  465.      st4 = c_re(input[2]);
  466.      tmp302 = st8;
  467.      st8 = st7 + st1;
  468.      st7 = st7 - st1;
  469.      st4 = st4 + c_re(input[10]);
  470.      st1 = c_re(input[6]);
  471.      st1 = st1 + c_re(input[14]);
  472.      tmp303 = st2;
  473.      st2 = c_im(input[2]);
  474.      tmp304 = st5;
  475.      st5 = st4 + st1;
  476.      st4 = st4 - st1;
  477.      st2 = st2 + c_im(input[10]);
  478.      st1 = st6 + st5;
  479.      st6 = st6 - st5;
  480.      st5 = st4 + st7;
  481.      st7 = st7 - st4;
  482.      st4 = c_im(input[6]);
  483.      st4 = st4 + c_im(input[14]);
  484.      tmp305 = st5;
  485.      st5 = c_re(input[6]);
  486.      tmp306 = st7;
  487.      st7 = st2 + st4;
  488.      st2 = st2 - st4;
  489.      st4 = st8 - st7;
  490.      st8 = st8 + st7;
  491.      st7 = st3 - st2;
  492.      st3 = st3 + st2;
  493.      st5 = st5 - c_re(input[14]);
  494.      st2 = c_im(input[2]);
  495.      st2 = st2 - c_im(input[10]);
  496.      tmp307 = st3;
  497.      st3 = c_re(input[2]);
  498.      tmp308 = st6;
  499.      st6 = st5 + st2;
  500.      st2 = st2 - st5;
  501.      st3 = st3 - c_re(input[10]);
  502.      st5 = c_im(input[6]);
  503.      st5 = st5 - c_im(input[14]);
  504.      tmp309 = st7;
  505.      st7 = c_re(input[5]);
  506.      tmp310 = st8;
  507.      st8 = st3 - st5;
  508.      st3 = st3 + st5;
  509.      st5 = st6 - st8;
  510.      st6 = st6 + st8;
  511.      st5 = st5 * K707106781[0];
  512.      st8 = st2 + st3;
  513.      st6 = st6 * K707106781[0];
  514.      st2 = st2 - st3;
  515.      st8 = st8 * K707106781[0];
  516.      st7 = st7 - c_re(input[13]);
  517.      st2 = st2 * K707106781[0];
  518.      st3 = tmp304 + st5;
  519.      tmp311 = st4;
  520.      st4 = tmp303 + st6;
  521.      st6 = tmp303 - st6;
  522.      st5 = tmp304 - st5;
  523.      tmp312 = st1;
  524.      st1 = tmp302 - st8;
  525.      st8 = tmp302 + st8;
  526.      tmp313 = st8;
  527.      st8 = tmp301 + st2;
  528.      st2 = tmp301 - st2;
  529.      tmp314 = st2;
  530.      st2 = c_im(input[1]);
  531.      st2 = st2 - c_im(input[9]);
  532.      tmp315 = st8;
  533.      st8 = c_re(input[1]);
  534.      tmp316 = st1;
  535.      st1 = st7 + st2;
  536.      st2 = st2 - st7;
  537.      st7 = st1 * K923879532[0];
  538.      st8 = st8 - c_re(input[9]);
  539.      st1 = st1 * K382683432[0];
  540.      tmp317 = st5;
  541.      st5 = c_im(input[5]);
  542.      tmp318 = st6;
  543.      st6 = st2 * K923879532[0];
  544.      st5 = st5 - c_im(input[13]);
  545.      st2 = st2 * K382683432[0];
  546.      tmp319 = st4;
  547.      st4 = st8 - st5;
  548.      st8 = st8 + st5;
  549.      st5 = st4 * K382683432[0];
  550.      tmp320 = st3;
  551.      st3 = c_re(input[7]);
  552.      st4 = st4 * K923879532[0];
  553.      st7 = st7 + st5;
  554.      st5 = st8 * K382683432[0];
  555.      st1 = st1 - st4;
  556.      st8 = st8 * K923879532[0];
  557.      st6 = st6 - st5;
  558.      st2 = st2 + st8;
  559.      st3 = st3 - c_re(input[15]);
  560.      st4 = c_im(input[3]);
  561.      st4 = st4 - c_im(input[11]);
  562.      st5 = c_re(input[3]);
  563.      st8 = st3 + st4;
  564.      st4 = st4 - st3;
  565.      st3 = st8 * K382683432[0];
  566.      st5 = st5 - c_re(input[11]);
  567.      st8 = st8 * K923879532[0];
  568.      tmp321 = st2;
  569.      st2 = c_im(input[7]);
  570.      tmp322 = st6;
  571.      st6 = st4 * K382683432[0];
  572.      st2 = st2 - c_im(input[15]);
  573.      st4 = st4 * K923879532[0];
  574.      tmp323 = st1;
  575.      st1 = st5 - st2;
  576.      st5 = st5 + st2;
  577.      st2 = st1 * K923879532[0];
  578.      tmp324 = st7;
  579.      st7 = c_re(input[1]);
  580.      st1 = st1 * K382683432[0];
  581.      st3 = st3 + st2;
  582.      st2 = st5 * K923879532[0];
  583.      st1 = st1 - st8;
  584.      st5 = st5 * K382683432[0];
  585.      st6 = st6 - st2;
  586.      st4 = st4 + st5;
  587.      st7 = st7 + c_re(input[9]);
  588.      st8 = tmp324 - st3;
  589.      st3 = tmp324 + st3;
  590.      st2 = tmp320 - st8;
  591.      st8 = tmp320 + st8;
  592.      st5 = tmp319 - st3;
  593.      st3 = tmp319 + st3;
  594.      tmp325 = st3;
  595.      st3 = tmp323 + st1;
  596.      st1 = tmp323 - st1;
  597.      tmp326 = st5;
  598.      st5 = tmp318 - st3;
  599.      st3 = tmp318 + st3;
  600.      tmp327 = st5;
  601.      st5 = tmp317 - st1;
  602.      st1 = tmp317 + st1;
  603.      tmp328 = st3;
  604.      st3 = tmp322 - st6;
  605.      tmp329 = st5;
  606.      st5 = tmp321 + st4;
  607.      tmp330 = st1;
  608.      st1 = tmp316 - st3;
  609.      st3 = tmp316 + st3;
  610.      tmp331 = st2;
  611.      st2 = tmp313 - st5;
  612.      st5 = tmp313 + st5;
  613.      st6 = tmp322 + st6;
  614.      c_re(input[9]) = st2;
  615.      c_re(input[1]) = st5;
  616.      st2 = tmp315 - st6;
  617.      st6 = tmp315 + st6;
  618.      st4 = tmp321 - st4;
  619.      st5 = c_re(input[5]);
  620.      c_re(input[5]) = st3;
  621.      st5 = st5 + c_re(input[13]);
  622.      c_re(input[13]) = st1;
  623.      st1 = st7 + st5;
  624.      st7 = st7 - st5;
  625.      st3 = tmp314 - st4;
  626.      st4 = tmp314 + st4;
  627.      st5 = c_im(input[1]);
  628.      c_im(input[1]) = st6;
  629.      st5 = st5 + c_im(input[9]);
  630.      c_im(input[9]) = st2;
  631.      st2 = c_im(input[5]);
  632.      c_im(input[5]) = st3;
  633.      st2 = st2 + c_im(input[13]);
  634.      c_im(input[13]) = st4;
  635.      st6 = st5 - st2;
  636.      st5 = st5 + st2;
  637.      st3 = c_re(input[3]);
  638.      c_re(input[3]) = st8;
  639.      st3 = st3 + c_re(input[11]);
  640.      c_re(input[11]) = tmp331;
  641.      st8 = c_re(input[7]);
  642.      c_re(input[7]) = tmp330;
  643.      st8 = st8 + c_re(input[15]);
  644.      c_re(input[15]) = tmp329;
  645.      st4 = st3 + st8;
  646.      st3 = st3 - st8;
  647.      st2 = st1 + st4;
  648.      st1 = st1 - st4;
  649.      st8 = st3 + st6;
  650.      st6 = st6 - st3;
  651.      st4 = tmp312 - st2;
  652.      st2 = tmp312 + st2;
  653.      st3 = tmp311 - st1;
  654.      c_re(input[8]) = st4;
  655.      c_re(input[0]) = st2;
  656.      c_im(input[4]) = st3;
  657.      st1 = st1 + tmp311;
  658.      st4 = c_im(input[3]);
  659.      c_im(input[3]) = tmp328;
  660.      c_im(input[12]) = st1;
  661.      st4 = st4 + c_im(input[11]);
  662.      c_im(input[11]) = tmp327;
  663.      st2 = c_im(input[7]);
  664.      c_im(input[7]) = tmp326;
  665.      st2 = st2 + c_im(input[15]);
  666.      c_im(input[15]) = tmp325;
  667.      st3 = st4 - st2;
  668.      st4 = st4 + st2;
  669.      st1 = st7 - st3;
  670.      st2 = st5 - st4;
  671.      tmp332 = st2;
  672.      st2 = st8 - st1;
  673.      st8 = st8 + st1;
  674.      st2 = st2 * K707106781[0];
  675.      st5 = st5 + st4;
  676.      st8 = st8 * K707106781[0];
  677.      st7 = st7 + st3;
  678.      st3 = tmp310 + st5;
  679.      st4 = st6 - st7;
  680.      st6 = st6 + st7;
  681.      c_im(input[0]) = st3;
  682.      st4 = st4 * K707106781[0];
  683.      st5 = tmp310 - st5;
  684.      st6 = st6 * K707106781[0];
  685.      st1 = tmp309 - st2;
  686.      c_im(input[8]) = st5;
  687.      c_re(input[14]) = st1;
  688.      st2 = tmp309 + st2;
  689.      st7 = tmp308 + tmp332;
  690.      st3 = tmp308 - tmp332;
  691.      c_re(input[6]) = st2;
  692.      c_re(input[4]) = st7;
  693.      c_re(input[12]) = st3;
  694.      st5 = tmp306 - st4;
  695.      st4 = tmp306 + st4;
  696.      st1 = tmp307 - st6;
  697.      c_im(input[10]) = st5;
  698.      c_im(input[2]) = st4;
  699.      c_re(input[10]) = st1;
  700.      st6 = tmp307 + st6;
  701.      st2 = tmp305 - st8;
  702.      st8 = tmp305 + st8;
  703.      c_re(input[2]) = st6;
  704.      c_im(input[6]) = st2;
  705.      c_im(input[14]) = st8;
  706. }
  708. void PFFTW(32) (fftw_complex * input) {
  709.      fftw_real tmp714;
  710.      fftw_real tmp713;
  711.      fftw_real tmp712;
  712.      fftw_real tmp711;
  713.      fftw_real tmp710;
  714.      fftw_real tmp709;
  715.      fftw_real tmp708;
  716.      fftw_real tmp707;
  717.      fftw_real tmp706;
  718.      fftw_real tmp705;
  719.      fftw_real tmp704;
  720.      fftw_real tmp703;
  721.      fftw_real tmp702;
  722.      fftw_real tmp701;
  723.      fftw_real tmp700;
  724.      fftw_real tmp699;
  725.      fftw_real tmp698;
  726.      fftw_real tmp697;
  727.      fftw_real tmp696;
  728.      fftw_real tmp695;
  729.      fftw_real tmp694;
  730.      fftw_real tmp693;
  731.      fftw_real tmp692;
  732.      fftw_real tmp691;
  733.      fftw_real tmp690;
  734.      fftw_real tmp689;
  735.      fftw_real tmp688;
  736.      fftw_real tmp687;
  737.      fftw_real tmp686;
  738.      fftw_real tmp685;
  739.      fftw_real tmp684;
  740.      fftw_real tmp683;
  741.      fftw_real tmp682;
  742.      fftw_real tmp681;
  743.      fftw_real tmp680;
  744.      fftw_real tmp679;
  745.      fftw_real tmp678;
  746.      fftw_real tmp677;
  747.      fftw_real tmp676;
  748.      fftw_real tmp675;
  749.      fftw_real tmp674;
  750.      fftw_real tmp673;
  751.      fftw_real tmp672;
  752.      fftw_real tmp671;
  753.      fftw_real tmp670;
  754.      fftw_real tmp669;
  755.      fftw_real tmp668;
  756.      fftw_real tmp667;
  757.      fftw_real tmp666;
  758.      fftw_real tmp665;
  759.      fftw_real tmp664;
  760.      fftw_real tmp663;
  761.      fftw_real tmp662;
  762.      fftw_real tmp661;
  763.      fftw_real tmp660;
  764.      fftw_real tmp659;
  765.      fftw_real tmp658;
  766.      fftw_real tmp657;
  767.      fftw_real tmp656;
  768.      fftw_real tmp655;
  769.      fftw_real tmp654;
  770.      fftw_real tmp653;
  771.      fftw_real tmp652;
  772.      fftw_real tmp651;
  773.      fftw_real tmp650;
  774.      fftw_real tmp649;
  775.      fftw_real tmp648;
  776.      fftw_real tmp647;
  777.      fftw_real tmp646;
  778.      fftw_real tmp645;
  779.      fftw_real tmp644;
  780.      fftw_real tmp643;
  781.      fftw_real tmp642;
  782.      fftw_real tmp641;
  783.      fftw_real tmp640;
  784.      fftw_real tmp639;
  785.      fftw_real tmp638;
  786.      fftw_real tmp637;
  787.      fftw_real tmp636;
  788.      fftw_real tmp635;
  789.      fftw_real tmp634;
  790.      fftw_real tmp633;
  791.      fftw_real tmp632;
  792.      fftw_real tmp631;
  793.      fftw_real tmp630;
  794.      fftw_real tmp629;
  795.      fftw_real tmp628;
  796.      fftw_real tmp627;
  797.      fftw_real tmp626;
  798.      fftw_real tmp625;
  799.      fftw_real tmp624;
  800.      fftw_real tmp623;
  801.      fftw_real tmp622;
  802.      fftw_real tmp621;
  803.      fftw_real st1;
  804.      fftw_real st2;
  805.      fftw_real st3;
  806.      fftw_real st4;
  807.      fftw_real st5;
  808.      fftw_real st6;
  809.      fftw_real st7;
  810.      fftw_real st8;
  811.      st8 = c_re(input[0]);
  812.      st8 = st8 - c_re(input[16]);
  813.      st7 = c_im(input[8]);
  814.      st7 = st7 - c_im(input[24]);
  815.      st6 = st8 - st7;
  816.      st8 = st8 + st7;
  817.      st5 = c_re(input[0]);
  818.      st5 = st5 + c_re(input[16]);
  819.      st4 = c_re(input[8]);
  820.      st4 = st4 + c_re(input[24]);
  821.      st3 = st5 + st4;
  822.      st5 = st5 - st4;
  823.      st2 = c_im(input[0]);
  824.      st2 = st2 + c_im(input[16]);
  825.      st1 = c_im(input[8]);
  826.      st1 = st1 + c_im(input[24]);
  827.      st7 = st2 + st1;
  828.      st2 = st2 - st1;
  829.      st4 = c_re(input[4]);
  830.      st4 = st4 + c_re(input[20]);
  831.      st1 = c_re(input[28]);
  832.      st1 = st1 + c_re(input[12]);
  833.      tmp621 = st6;
  834.      st6 = st4 + st1;
  835.      st1 = st1 - st4;
  836.      st4 = st3 + st6;
  837.      st3 = st3 - st6;
  838.      st6 = st2 - st1;
  839.      st1 = st1 + st2;
  840.      st2 = c_im(input[4]);
  841.      st2 = st2 + c_im(input[20]);
  842.      tmp622 = st1;
  843.      st1 = c_im(input[28]);
  844.      st1 = st1 + c_im(input[12]);
  845.      tmp623 = st6;
  846.      st6 = st2 + st1;
  847.      st2 = st2 - st1;
  848.      st1 = st7 + st6;
  849.      st7 = st7 - st6;
  850.      st6 = st5 - st2;
  851.      st5 = st5 + st2;
  852.      st2 = c_re(input[8]);
  853.      st2 = st2 - c_re(input[24]);
  854.      tmp624 = st5;
  855.      st5 = c_im(input[0]);
  856.      st5 = st5 - c_im(input[16]);
  857.      tmp625 = st6;
  858.      st6 = st2 + st5;
  859.      st5 = st5 - st2;
  860.      st2 = c_im(input[4]);
  861.      st2 = st2 - c_im(input[20]);
  862.      tmp626 = st3;
  863.      st3 = c_re(input[4]);
  864.      st3 = st3 - c_re(input[20]);
  865.      tmp627 = st1;
  866.      st1 = st2 - st3;
  867.      st3 = st3 + st2;
  868.      st2 = c_re(input[28]);
  869.      st2 = st2 - c_re(input[12]);
  870.      tmp628 = st7;
  871.      st7 = c_im(input[28]);
  872.      st7 = st7 - c_im(input[12]);
  873.      tmp629 = st4;
  874.      st4 = st2 + st7;
  875.      st2 = st2 - st7;
  876.      st7 = st1 - st4;
  877.      st7 = st7 * K707106781[0];
  878.      st1 = st1 + st4;
  879.      st1 = st1 * K707106781[0];
  880.      st4 = st2 - st3;
  881.      st4 = st4 * K707106781[0];
  882.      st3 = st3 + st2;
  883.      st3 = st3 * K707106781[0];
  884.      st2 = st8 - st3;
  885.      tmp630 = st2;
  886.      st2 = st5 - st1;
  887.      tmp631 = st2;
  888.      st2 = tmp621 - st7;
  889.      tmp632 = st2;
  890.      st2 = st6 - st4;
  891.      st7 = tmp621 + st7;
  892.      st6 = st6 + st4;
  893.      st8 = st8 + st3;
  894.      st5 = st5 + st1;
  895.      st1 = c_re(input[2]);
  896.      st1 = st1 + c_re(input[18]);
  897.      st4 = c_re(input[10]);
  898.      st4 = st4 + c_re(input[26]);
  899.      st3 = st1 + st4;
  900.      st1 = st1 - st4;
  901.      st4 = c_im(input[2]);
  902.      st4 = st4 + c_im(input[18]);
  903.      tmp633 = st5;
  904.      st5 = c_im(input[10]);
  905.      st5 = st5 + c_im(input[26]);
  906.      tmp634 = st8;
  907.      st8 = st4 + st5;
  908.      st4 = st4 - st5;
  909.      st5 = st1 + st4;
  910.      st4 = st4 - st1;
  911.      st1 = c_re(input[30]);
  912.      st1 = st1 + c_re(input[14]);
  913.      tmp635 = st6;
  914.      st6 = c_re(input[6]);
  915.      st6 = st6 + c_re(input[22]);
  916.      tmp636 = st7;
  917.      st7 = st1 + st6;
  918.      st1 = st1 - st6;
  919.      st6 = st3 + st7;
  920.      st7 = st7 - st3;
  921.      st3 = tmp629 + st6;
  922.      st6 = tmp629 - st6;
  923.      tmp637 = st6;
  924.      st6 = st7 + tmp628;
  925.      st7 = tmp628 - st7;
  926.      tmp638 = st7;
  927.      st7 = c_im(input[30]);
  928.      st7 = st7 + c_im(input[14]);
  929.      tmp639 = st6;
  930.      st6 = c_im(input[6]);
  931.      st6 = st6 + c_im(input[22]);
  932.      tmp640 = st3;
  933.      st3 = st7 + st6;
  934.      st7 = st7 - st6;
  935.      st6 = st8 + st3;
  936.      st8 = st8 - st3;
  937.      st3 = st1 - st7;
  938.      tmp641 = st2;
  939.      st2 = st3 - st5;
  940.      st2 = st2 * K707106781[0];
  941.      st5 = st5 + st3;
  942.      st5 = st5 * K707106781[0];
  943.      st1 = st1 + st7;
  944.      st7 = st4 - st1;
  945.      st7 = st7 * K707106781[0];
  946.      st4 = st4 + st1;
  947.      st4 = st4 * K707106781[0];
  948.      st3 = tmp627 - st6;
  949.      st6 = tmp627 + st6;
  950.      st1 = tmp626 - st8;
  951.      st8 = tmp626 + st8;
  952.      tmp642 = st8;
  953.      st8 = tmp625 + st7;
  954.      st7 = tmp625 - st7;
  955.      tmp643 = st7;
  956.      st7 = tmp623 - st2;
  957.      st2 = tmp623 + st2;
  958.      tmp644 = st2;
  959.      st2 = tmp624 + st5;
  960.      st5 = tmp624 - st5;
  961.      tmp645 = st5;
  962.      st5 = tmp622 - st4;
  963.      st4 = tmp622 + st4;
  964.      tmp646 = st4;
  965.      st4 = c_re(input[6]);
  966.      st4 = st4 - c_re(input[22]);
  967.      tmp647 = st5;
  968.      st5 = c_im(input[30]);
  969.      st5 = st5 - c_im(input[14]);
  970.      tmp648 = st2;
  971.      st2 = st4 + st5;
  972.      tmp649 = st7;
  973.      st7 = st2 * K382683432[0];
  974.      st5 = st5 - st4;
  975.      st2 = st2 * K923879532[0];
  976.      st4 = c_re(input[30]);
  977.      tmp650 = st8;
  978.      st8 = st5 * K923879532[0];
  979.      st4 = st4 - c_re(input[14]);
  980.      st5 = st5 * K382683432[0];
  981.      tmp651 = st1;
  982.      st1 = c_im(input[6]);
  983.      st1 = st1 - c_im(input[22]);
  984.      tmp652 = st6;
  985.      st6 = st4 - st1;
  986.      tmp653 = st3;
  987.      st3 = st6 * K923879532[0];
  988.      st4 = st4 + st1;
  989.      st6 = st6 * K382683432[0];
  990.      st7 = st7 + st3;
  991.      st1 = st4 * K382683432[0];
  992.      st6 = st6 - st2;
  993.      st4 = st4 * K923879532[0];
  994.      st8 = st8 + st1;
  995.      st4 = st4 - st5;
  996.      st2 = c_re(input[10]);
  997.      st2 = st2 - c_re(input[26]);
  998.      st5 = c_im(input[2]);
  999.      st5 = st5 - c_im(input[18]);
  1000.      st3 = st2 + st5;
  1001.      st1 = st3 * K382683432[0];
  1002.      st5 = st5 - st2;
  1003.      st3 = st3 * K923879532[0];
  1004.      st2 = c_re(input[2]);
  1005.      tmp654 = st6;
  1006.      st6 = st5 * K923879532[0];
  1007.      st2 = st2 - c_re(input[18]);
  1008.      st5 = st5 * K382683432[0];
  1009.      tmp655 = st7;
  1010.      st7 = c_im(input[10]);
  1011.      st7 = st7 - c_im(input[26]);
  1012.      tmp656 = st4;
  1013.      st4 = st2 - st7;
  1014.      tmp657 = st8;
  1015.      st8 = st4 * K923879532[0];
  1016.      st2 = st2 + st7;
  1017.      st4 = st4 * K382683432[0];
  1018.      st1 = st1 - st8;
  1019.      st7 = st2 * K382683432[0];
  1020.      st3 = st3 + st4;
  1021.      st2 = st2 * K923879532[0];
  1022.      st6 = st6 - st7;
  1023.      st5 = st5 + st2;
  1024.      st8 = st6 - tmp657;
  1025.      st4 = tmp630 + st8;
  1026.      st8 = tmp630 - st8;
  1027.      st7 = tmp656 - st5;
  1028.      st2 = tmp631 - st7;
  1029.      st7 = tmp631 + st7;
  1030.      tmp658 = st7;
  1031.      st7 = st1 - tmp655;
  1032.      tmp659 = st8;
  1033.      st8 = tmp632 + st7;
  1034.      st7 = tmp632 - st7;
  1035.      tmp660 = st7;
  1036.      st7 = tmp654 - st3;
  1037.      tmp661 = st8;
  1038.      st8 = tmp641 - st7;
  1039.      st7 = tmp641 + st7;
  1040.      st3 = st3 + tmp654;
  1041.      tmp662 = st7;
  1042.      st7 = tmp636 + st3;
  1043.      st3 = tmp636 - st3;
  1044.      st1 = st1 + tmp655;
  1045.      tmp663 = st3;
  1046.      st3 = tmp635 - st1;
  1047.      st1 = tmp635 + st1;
  1048.      st5 = st5 + tmp656;
  1049.      tmp664 = st1;
  1050.      st1 = tmp634 + st5;
  1051.      st5 = tmp634 - st5;
  1052.      st6 = st6 + tmp657;
  1053.      tmp665 = st5;
  1054.      st5 = tmp633 - st6;
  1055.      st6 = tmp633 + st6;
  1056.      tmp666 = st6;
  1057.      st6 = c_re(input[1]);
  1058.      st6 = st6 + c_re(input[17]);
  1059.      tmp667 = st5;
  1060.      st5 = c_re(input[9]);
  1061.      st5 = st5 + c_re(input[25]);
  1062.      tmp668 = st1;
  1063.      st1 = st6 + st5;
  1064.      st6 = st6 - st5;
  1065.      st5 = c_im(input[1]);
  1066.      st5 = st5 + c_im(input[17]);
  1067.      tmp669 = st3;
  1068.      st3 = c_im(input[9]);
  1069.      st3 = st3 + c_im(input[25]);
  1070.      tmp670 = st7;
  1071.      st7 = st5 - st3;
  1072.      st5 = st5 + st3;
  1073.      st3 = c_re(input[5]);
  1074.      st3 = st3 + c_re(input[21]);
  1075.      tmp671 = st8;
  1076.      st8 = c_re(input[29]);
  1077.      st8 = st8 + c_re(input[13]);
  1078.      tmp672 = st2;
  1079.      st2 = st3 + st8;
  1080.      st8 = st8 - st3;
  1081.      st3 = st1 + st2;
  1082.      st1 = st1 - st2;
  1083.      st2 = st8 + st7;
  1084.      tmp673 = st4;
  1085.      st4 = st2 * K382683432[0];
  1086.      st7 = st7 - st8;
  1087.      st2 = st2 * K923879532[0];
  1088.      st8 = c_im(input[5]);
  1089.      tmp674 = st3;
  1090.      st3 = st7 * K923879532[0];
  1091.      st8 = st8 + c_im(input[21]);
  1092.      st7 = st7 * K382683432[0];
  1093.      tmp675 = st4;
  1094.      st4 = c_im(input[29]);
  1095.      st4 = st4 + c_im(input[13]);
  1096.      tmp676 = st7;
  1097.      st7 = st8 - st4;
  1098.      st8 = st8 + st4;
  1099.      st4 = st5 + st8;
  1100.      st5 = st5 - st8;
  1101.      st8 = st1 + st5;
  1102.      st5 = st5 - st1;
  1103.      st1 = st6 + st7;
  1104.      tmp677 = st8;
  1105.      st8 = st1 * K923879532[0];
  1106.      st6 = st6 - st7;
  1107.      st1 = st1 * K382683432[0];
  1108.      st2 = st2 - st1;
  1109.      st7 = st6 * K382683432[0];
  1110.      st3 = st3 + st7;
  1111.      st6 = st6 * K923879532[0];
  1112.      st6 = tmp676 - st6;
  1113.      st8 = tmp675 + st8;
  1114.      st1 = c_re(input[31]);
  1115.      st1 = st1 + c_re(input[15]);
  1116.      st7 = c_re(input[7]);
  1117.      st7 = st7 + c_re(input[23]);
  1118.      tmp678 = st8;
  1119.      st8 = st1 + st7;
  1120.      st1 = st1 - st7;
  1121.      st7 = c_im(input[31]);
  1122.      st7 = st7 + c_im(input[15]);
  1123.      tmp679 = st2;
  1124.      st2 = c_im(input[7]);
  1125.      st2 = st2 + c_im(input[23]);
  1126.      tmp680 = st3;
  1127.      st3 = st7 - st2;
  1128.      st7 = st7 + st2;
  1129.      st2 = c_re(input[3]);
  1130.      st2 = st2 + c_re(input[19]);
  1131.      tmp681 = st6;
  1132.      st6 = c_re(input[27]);
  1133.      st6 = st6 + c_re(input[11]);
  1134.      tmp682 = st5;
  1135.      st5 = st2 + st6;
  1136.      st6 = st6 - st2;
  1137.      st2 = st8 + st5;
  1138.      st8 = st8 - st5;
  1139.      st5 = st6 + st3;
  1140.      tmp683 = st4;
  1141.      st4 = st5 * K382683432[0];
  1142.      st3 = st3 - st6;
  1143.      st5 = st5 * K923879532[0];
  1144.      st6 = tmp674 + st2;
  1145.      tmp684 = st4;
  1146.      st4 = st3 * K923879532[0];
  1147.      tmp685 = st5;
  1148.      st5 = tmp640 - st6;
  1149.      st3 = st3 * K382683432[0];
  1150.      st6 = tmp640 + st6;
  1151.      st2 = st2 - tmp674;
  1152.      c_re(input[16]) = st5;
  1153.      st5 = st2 + tmp653;
  1154.      st2 = tmp653 - st2;
  1155.      c_re(input[0]) = st6;
  1156.      st6 = c_im(input[3]);
  1157.      st6 = st6 + c_im(input[19]);
  1158.      c_im(input[8]) = st5;
  1159.      st5 = c_im(input[27]);
  1160.      st5 = st5 + c_im(input[11]);
  1161.      c_im(input[24]) = st2;
  1162.      st2 = st6 - st5;
  1163.      st6 = st6 + st5;
  1164.      st5 = st7 + st6;
  1165.      st7 = st7 - st6;
  1166.      st6 = st8 - st7;
  1167.      st8 = st8 + st7;
  1168.      st7 = st1 + st2;
  1169.      tmp686 = st4;
  1170.      st4 = st7 * K923879532[0];
  1171.      st1 = st1 - st2;
  1172.      st7 = st7 * K382683432[0];
  1173.      st2 = tmp683 + st5;
  1174.      tmp687 = st4;
  1175.      st4 = st1 * K382683432[0];
  1176.      tmp688 = st7;
  1177.      st7 = tmp652 - st2;
  1178.      st1 = st1 * K923879532[0];
  1179.      st2 = tmp652 + st2;
  1180.      st5 = tmp683 - st5;
  1181.      c_im(input[16]) = st7;
  1182.      st7 = tmp637 - st5;
  1183.      st5 = tmp637 + st5;
  1184.      c_im(input[0]) = st2;
  1185.      st2 = tmp682 + st8;
  1186.      st2 = st2 * K707106781[0];
  1187.      c_re(input[24]) = st7;
  1188.      st7 = tmp639 - st2;
  1189.      st2 = tmp639 + st2;
  1190.      st8 = tmp682 - st8;
  1191.      st8 = st8 * K707106781[0];
  1192.      c_re(input[8]) = st5;
  1193.      st5 = tmp651 - st8;
  1194.      st8 = tmp651 + st8;
  1195.      c_im(input[20]) = st7;
  1196.      st7 = tmp677 + st6;
  1197.      st7 = st7 * K707106781[0];
  1198.      c_im(input[4]) = st2;
  1199.      st2 = tmp642 - st7;
  1200.      st7 = tmp642 + st7;
  1201.      st6 = st6 - tmp677;
  1202.      st6 = st6 * K707106781[0];
  1203.      c_re(input[28]) = st5;
  1204.      st5 = tmp638 - st6;
  1205.      st6 = tmp638 + st6;
  1206.      st3 = st3 + st1;
  1207.      st1 = tmp681 - st3;
  1208.      st3 = tmp681 + st3;
  1209.      st4 = st4 - tmp686;
  1210.      c_re(input[12]) = st8;
  1211.      st8 = tmp680 + st4;
  1212.      st4 = st4 - tmp680;
  1213.      c_re(input[20]) = st2;
  1214.      st2 = tmp650 - st8;
  1215.      st8 = tmp650 + st8;
  1216.      c_re(input[4]) = st7;
  1217.      st7 = tmp649 - st4;
  1218.      st4 = tmp649 + st4;
  1219.      c_im(input[28]) = st5;
  1220.      st5 = tmp643 - st1;
  1221.      st1 = tmp643 + st1;
  1222.      c_im(input[12]) = st6;
  1223.      st6 = tmp644 - st3;
  1224.      st3 = tmp644 + st3;
  1225.      c_re(input[22]) = st2;
  1226.      st2 = tmp685 + tmp688;
  1227.      c_re(input[6]) = st8;
  1228.      st8 = tmp679 - st2;
  1229.      st2 = tmp679 + st2;
  1230.      c_im(input[30]) = st7;
  1231.      st7 = tmp687 - tmp684;
  1232.      c_im(input[14]) = st4;
  1233.      st4 = tmp678 + st7;
  1234.      st7 = st7 - tmp678;
  1235.      c_re(input[30]) = st5;
  1236.      st5 = tmp648 - st4;
  1237.      st4 = tmp648 + st4;
  1238.      c_re(input[14]) = st1;
  1239.      st1 = tmp647 - st7;
  1240.      st7 = tmp647 + st7;
  1241.      c_im(input[22]) = st6;
  1242.      st6 = tmp645 - st8;
  1243.      st8 = tmp645 + st8;
  1244.      c_im(input[6]) = st3;
  1245.      st3 = tmp646 - st2;
  1246.      st2 = tmp646 + st2;
  1247.      c_re(input[18]) = st5;
  1248.      st5 = c_re(input[31]);
  1249.      st5 = st5 - c_re(input[15]);
  1250.      c_re(input[2]) = st4;
  1251.      st4 = c_im(input[7]);
  1252.      st4 = st4 - c_im(input[23]);
  1253.      c_im(input[26]) = st1;
  1254.      st1 = st5 - st4;
  1255.      st5 = st5 + st4;
  1256.      c_im(input[10]) = st7;
  1257.      st7 = c_re(input[7]);
  1258.      st7 = st7 - c_re(input[23]);
  1259.      c_re(input[26]) = st6;
  1260.      st6 = c_im(input[31]);
  1261.      st6 = st6 - c_im(input[15]);
  1262.      c_re(input[10]) = st8;
  1263.      st8 = st7 + st6;
  1264.      st6 = st6 - st7;
  1265.      c_im(input[18]) = st3;
  1266.      st3 = c_im(input[3]);
  1267.      st3 = st3 - c_im(input[19]);
  1268.      c_im(input[2]) = st2;
  1269.      st2 = c_re(input[3]);
  1270.      st2 = st2 - c_re(input[19]);
  1271.      st4 = st3 - st2;
  1272.      st2 = st2 + st3;
  1273.      st7 = c_re(input[27]);
  1274.      st7 = st7 - c_re(input[11]);
  1275.      st3 = c_im(input[27]);
  1276.      st3 = st3 - c_im(input[11]);
  1277.      tmp689 = st5;
  1278.      st5 = st7 + st3;
  1279.      st7 = st7 - st3;
  1280.      st3 = st4 - st5;
  1281.      st3 = st3 * K707106781[0];
  1282.      st4 = st4 + st5;
  1283.      st4 = st4 * K707106781[0];
  1284.      st5 = st7 - st2;
  1285.      st5 = st5 * K707106781[0];
  1286.      st2 = st2 + st7;
  1287.      st2 = st2 * K707106781[0];
  1288.      st7 = st1 - st3;
  1289.      tmp690 = st2;
  1290.      st2 = st7 * K980785280[0];
  1291.      st1 = st1 + st3;
  1292.      st7 = st7 * K195090322[0];
  1293.      st3 = st6 - st4;
  1294.      tmp691 = st7;
  1295.      st7 = st3 * K555570233[0];
  1296.      st6 = st6 + st4;
  1297.      st3 = st3 * K831469612[0];
  1298.      st4 = st8 - st5;
  1299.      tmp692 = st6;
  1300.      st6 = st4 * K195090322[0];
  1301.      st8 = st8 + st5;
  1302.      st4 = st4 * K980785280[0];
  1303.      st5 = tmp689 - tmp690;
  1304.      tmp693 = st4;
  1305.      st4 = st5 * K831469612[0];
  1306.      tmp694 = st2;
  1307.      st2 = tmp689 + tmp690;
  1308.      st5 = st5 * K555570233[0];
  1309.      st4 = st4 + st7;
  1310.      st7 = st8 * K831469612[0];
  1311.      st5 = st5 - st3;
  1312.      st3 = st1 * K555570233[0];
  1313.      st6 = st6 + tmp694;
  1314.      st1 = st1 * K831469612[0];
  1315.      tmp695 = st6;
  1316.      st6 = tmp691 - tmp693;
  1317.      st8 = st8 * K555570233[0];
  1318.      st7 = st7 + st3;
  1319.      st3 = st2 * K195090322[0];
  1320.      st1 = st1 - st8;
  1321.      st8 = tmp692 * K980785280[0];
  1322.      st3 = st3 + st8;
  1323.      st2 = st2 * K980785280[0];
  1324.      st8 = c_re(input[9]);
  1325.      tmp696 = st3;
  1326.      st3 = tmp692 * K195090322[0];
  1327.      st8 = st8 - c_re(input[25]);
  1328.      st2 = st2 - st3;
  1329.      st3 = c_im(input[1]);
  1330.      st3 = st3 - c_im(input[17]);
  1331.      tmp697 = st2;
  1332.      st2 = st8 + st3;
  1333.      st3 = st3 - st8;
  1334.      st8 = c_re(input[1]);
  1335.      st8 = st8 - c_re(input[17]);
  1336.      tmp698 = st1;
  1337.      st1 = c_im(input[9]);
  1338.      st1 = st1 - c_im(input[25]);
  1339.      tmp699 = st7;
  1340.      st7 = st8 - st1;
  1341.      st8 = st8 + st1;
  1342.      st1 = c_re(input[29]);
  1343.      st1 = st1 - c_re(input[13]);
  1344.      tmp700 = st6;
  1345.      st6 = c_im(input[29]);
  1346.      st6 = st6 - c_im(input[13]);
  1347.      tmp701 = st5;
  1348.      st5 = st1 - st6;
  1349.      st1 = st1 + st6;
  1350.      st6 = c_re(input[5]);
  1351.      st6 = st6 - c_re(input[21]);
  1352.      tmp702 = st4;
  1353.      st4 = c_im(input[5]);
  1354.      st4 = st4 - c_im(input[21]);
  1355.      tmp703 = st8;
  1356.      st8 = st6 + st4;
  1357.      st4 = st4 - st6;
  1358.      st6 = st5 - st8;
  1359.      st6 = st6 * K707106781[0];
  1360.      st8 = st8 + st5;
  1361.      st8 = st8 * K707106781[0];
  1362.      st5 = st4 - st1;
  1363.      st5 = st5 * K707106781[0];
  1364.      st4 = st4 + st1;
  1365.      st4 = st4 * K707106781[0];
  1366.      st1 = st2 - st6;
  1367.      tmp704 = st8;
  1368.      st8 = st1 * K195090322[0];
  1369.      st2 = st2 + st6;
  1370.      st1 = st1 * K980785280[0];
  1371.      st6 = st3 - st4;
  1372.      tmp705 = st1;
  1373.      st1 = st6 * K555570233[0];
  1374.      st3 = st3 + st4;
  1375.      st6 = st6 * K831469612[0];
  1376.      st4 = st7 - st5;
  1377.      tmp706 = st8;
  1378.      st8 = st4 * K980785280[0];
  1379.      st7 = st7 + st5;
  1380.      st4 = st4 * K195090322[0];
  1381.      st5 = tmp703 - tmp704;
  1382.      tmp707 = st4;
  1383.      st4 = st5 * K831469612[0];
  1384.      tmp708 = st8;
  1385.      st8 = tmp703 + tmp704;
  1386.      st5 = st5 * K555570233[0];
  1387.      st1 = st1 - st4;
  1388.      st4 = st2 * K831469612[0];
  1389.      tmp709 = st4;
  1390.      st4 = st1 - tmp702;
  1391.      tmp710 = st4;
  1392.      st4 = st7 * K555570233[0];
  1393.      st1 = st1 + tmp702;
  1394.      st2 = st2 * K555570233[0];
  1395.      st5 = st5 + st6;
  1396.      st7 = st7 * K831469612[0];
  1397.      st6 = st5 + tmp701;
  1398.      tmp711 = st2;
  1399.      st2 = st3 * K980785280[0];
  1400.      st5 = tmp701 - st5;
  1401.      tmp712 = st2;
  1402.      st2 = st8 * K195090322[0];
  1403.      tmp713 = st2;
  1404.      st2 = tmp673 - st6;
  1405.      st8 = st8 * K980785280[0];
  1406.      st3 = st3 * K195090322[0];
  1407.      st6 = tmp673 + st6;
  1408.      c_re(input[21]) = st2;
  1409.      st2 = tmp672 - st5;
  1410.      st5 = tmp672 + st5;
  1411.      c_re(input[5]) = st6;
  1412.      st6 = tmp659 - tmp710;
  1413.      c_im(input[29]) = st2;
  1414.      st2 = tmp659 + tmp710;
  1415.      c_im(input[13]) = st5;
  1416.      st5 = tmp658 - st1;
  1417.      st1 = tmp658 + st1;
  1418.      c_re(input[29]) = st6;
  1419.      st6 = tmp706 - tmp708;
  1420.      c_re(input[13]) = st2;
  1421.      st2 = st6 - tmp695;
  1422.      st6 = st6 + tmp695;
  1423.      c_im(input[21]) = st5;
  1424.      st5 = tmp705 + tmp707;
  1425.      c_im(input[5]) = st1;
  1426.      st1 = st5 + tmp700;
  1427.      st5 = tmp700 - st5;
  1428.      tmp714 = st8;
  1429.      st8 = tmp661 - st1;
  1430.      st1 = tmp661 + st1;
  1431.      c_re(input[23]) = st8;
  1432.      st8 = tmp671 - st5;
  1433.      st5 = tmp671 + st5;
  1434.      c_re(input[7]) = st1;
  1435.      st1 = tmp660 - st2;
  1436.      st2 = tmp660 + st2;
  1437.      c_im(input[31]) = st8;
  1438.      st8 = tmp662 - st6;
  1439.      st6 = tmp662 + st6;
  1440.      st4 = tmp709 - st4;
  1441.      c_im(input[15]) = st5;
  1442.      st5 = st4 - tmp699;
  1443.      st4 = st4 + tmp699;
  1444.      st7 = tmp711 + st7;
  1445.      c_re(input[31]) = st1;
  1446.      st1 = st7 + tmp698;
  1447.      st7 = tmp698 - st7;
  1448.      c_re(input[15]) = st2;
  1449.      st2 = tmp670 - st1;
  1450.      st1 = tmp670 + st1;
  1451.      c_im(input[23]) = st8;
  1452.      st8 = tmp669 - st7;
  1453.      st7 = tmp669 + st7;
  1454.      c_im(input[7]) = st6;
  1455.      st6 = tmp663 - st5;
  1456.      st5 = tmp663 + st5;
  1457.      c_re(input[19]) = st2;
  1458.      st2 = tmp664 - st4;
  1459.      st4 = tmp664 + st4;
  1460.      c_re(input[3]) = st1;
  1461.      st1 = tmp712 - tmp713;
  1462.      c_im(input[27]) = st8;
  1463.      st8 = st1 - tmp696;
  1464.      st1 = st1 + tmp696;
  1465.      st3 = tmp714 + st3;
  1466.      c_im(input[11]) = st7;
  1467.      st7 = st3 + tmp697;
  1468.      st3 = tmp697 - st3;
  1469.      c_re(input[27]) = st6;
  1470.      st6 = tmp668 - st7;
  1471.      st7 = tmp668 + st7;
  1472.      c_re(input[11]) = st5;
  1473.      st5 = tmp667 - st3;
  1474.      st3 = tmp667 + st3;
  1475.      c_im(input[19]) = st2;
  1476.      st2 = tmp665 - st8;
  1477.      st8 = tmp665 + st8;
  1478.      c_im(input[3]) = st4;
  1479.      st4 = tmp666 - st1;
  1480.      st1 = tmp666 + st1;
  1481.      c_re(input[17]) = st6;
  1482.      c_re(input[1]) = st7;
  1483.      c_im(input[25]) = st5;
  1484.      c_im(input[9]) = st3;
  1485.      c_re(input[25]) = st2;
  1486.      c_re(input[9]) = st8;
  1487.      c_im(input[17]) = st4;
  1488.      c_im(input[1]) = st1;
  1489. }
  1491. void  PFFTW(64)(fftw_complex *input)
  1492. {
  1493.      PFFTW(twiddle_4)(input, PFFTW(W_64), 16);
  1494.      PFFTW(16)(input );
  1495.      PFFTW(16)(input + 16);
  1496.      PFFTW(16)(input + 32);
  1497.      PFFTW(16)(input + 48);
  1498. }
  1500. void PFFTW(128)(fftw_complex *input)
  1501. {
  1502.      PFFTW(twiddle_4)(input, PFFTW(W_128), 32);
  1503.      PFFTW(32)(input );
  1504.      PFFTW(32)(input + 32);
  1505.      PFFTW(32)(input + 64);
  1506.      PFFTW(32)(input + 96);
  1507. }
  1509. void PFFTW(512)(fftw_complex *input)
  1510. {
  1511.      PFFTW(twiddle_4)(input, PFFTW(W_512), 128);
  1512.      PFFTW(128)(input );
  1513.      PFFTW(128)(input + 128);
  1514.      PFFTW(128)(input + 256);
  1515.      PFFTW(128)(input + 384);
  1516. }
  1517. ///////////////////////////////////////////////////////////////
  1518. void PFFTWI(16) (fftw_complex * input) {
  1519.      fftw_real tmp333;
  1520.      fftw_real tmp332;
  1521.      fftw_real tmp331;
  1522.      fftw_real tmp330;
  1523.      fftw_real tmp329;
  1524.      fftw_real tmp328;
  1525.      fftw_real tmp327;
  1526.      fftw_real tmp326;
  1527.      fftw_real tmp325;
  1528.      fftw_real tmp324;
  1529.      fftw_real tmp323;
  1530.      fftw_real tmp322;
  1531.      fftw_real tmp321;
  1532.      fftw_real tmp320;
  1533.      fftw_real tmp319;
  1534.      fftw_real tmp318;
  1535.      fftw_real tmp317;
  1536.      fftw_real tmp316;
  1537.      fftw_real tmp315;
  1538.      fftw_real tmp314;
  1539.      fftw_real tmp313;
  1540.      fftw_real tmp312;
  1541.      fftw_real tmp311;
  1542.      fftw_real tmp310;
  1543.      fftw_real tmp309;
  1544.      fftw_real tmp308;
  1545.      fftw_real tmp307;
  1546.      fftw_real tmp306;
  1547.      fftw_real tmp305;
  1548.      fftw_real tmp304;
  1549.      fftw_real tmp303;
  1550.      fftw_real tmp302;
  1551.      fftw_real tmp301;
  1552.      fftw_real st1;
  1553.      fftw_real st2;
  1554.      fftw_real st3;
  1555.      fftw_real st4;
  1556.      fftw_real st5;
  1557.      fftw_real st6;
  1558.      fftw_real st7;
  1559.      fftw_real st8;
  1560.      st8 = c_re(input[4]);
  1561.      st8 = st8 - c_re(input[12]);
  1562.      st7 = c_im(input[0]);
  1563.      st7 = st7 - c_im(input[8]);
  1564.      st6 = c_re(input[0]);
  1565.      st6 = st6 - c_re(input[8]);
  1566.      st5 = st8 + st7;
  1567.      st7 = st7 - st8;
  1568.      st4 = c_im(input[4]);
  1569.      st4 = st4 - c_im(input[12]);
  1570.      st3 = c_re(input[0]);
  1571.      st3 = st3 + c_re(input[8]);
  1572.      st2 = st6 - st4;
  1573.      st6 = st6 + st4;
  1574.      st1 = c_re(input[4]);
  1575.      st1 = st1 + c_re(input[12]);
  1576.      st8 = c_im(input[0]);
  1577.      st8 = st8 + c_im(input[8]);
  1578.      st4 = st3 + st1;
  1579.      st3 = st3 - st1;
  1580.      st1 = c_im(input[4]);
  1581.      st1 = st1 + c_im(input[12]);
  1582.      tmp301 = st6;
  1583.      st6 = c_re(input[2]);
  1584.      st6 = st6 + c_re(input[10]);
  1585.      tmp302 = st7;
  1586.      st7 = st8 + st1;
  1587.      st8 = st8 - st1;
  1588.      st1 = c_re(input[6]);
  1589.      st1 = st1 + c_re(input[14]);
  1590.      tmp303 = st2;
  1591.      st2 = c_im(input[2]);
  1592.      st2 = st2 + c_im(input[10]);
  1593.      tmp304 = st5;
  1594.      st5 = st6 + st1;
  1595.      st6 = st6 - st1;
  1596.      st1 = c_im(input[6]);
  1597.      st1 = st1 + c_im(input[14]);
  1598.      tmp305 = st3;
  1599.      st3 = st4 + st5;
  1600.      st4 = st4 - st5;
  1601.      st5 = st2 + st1;
  1602.      st2 = st2 - st1;
  1603.      st1 = st6 + st8;
  1604.      tmp306 = st1;
  1605.      st1 = st7 - st5;
  1606.      st7 = st7 + st5;
  1607.      st5 = tmp305 - st2;
  1608.      st2 = tmp305 + st2;
  1609.      st8 = st8 - st6;
  1610.      st6 = c_re(input[6]);
  1611.      st6 = st6 - c_re(input[14]);
  1612.      tmp307 = st8;
  1613.      st8 = c_im(input[2]);
  1614.      st8 = st8 - c_im(input[10]);
  1615.      tmp308 = st2;
  1616.      st2 = c_re(input[2]);
  1617.      st2 = st2 - c_re(input[10]);
  1618.      tmp309 = st4;
  1619.      st4 = st6 + st8;
  1620.      st8 = st8 - st6;
  1621.      st6 = c_im(input[6]);
  1622.      st6 = st6 - c_im(input[14]);
  1623.      tmp310 = st5;
  1624.      st5 = c_re(input[1]);
  1625.      st5 = st5 - c_re(input[9]);
  1626.      tmp311 = st7;
  1627.      st7 = st2 - st6;
  1628.      st2 = st2 + st6;
  1629.      st6 = c_im(input[5]);
  1630.      tmp312 = st1;
  1631.      st1 = st4 + st7;
  1632.      st7 = st7 - st4;
  1633.      st4 = st2 - st8;
  1634.      st1 = st1 * K707106781[0];
  1635.      st8 = st8 + st2;
  1636.      st7 = st7 * K707106781[0];
  1637.      st6 = st6 - c_im(input[13]);
  1638.      st4 = st4 * K707106781[0];
  1639.      st2 = tmp304 - st1;
  1640.      st8 = st8 * K707106781[0];
  1641.      tmp313 = st3;
  1642.      st3 = st5 - st6;
  1643.      st5 = st5 + st6;
  1644.      st6 = tmp303 + st7;
  1645.      tmp314 = st6;
  1646.      st6 = st3 * K923879532[0];
  1647.      st7 = tmp303 - st7;
  1648.      st3 = st3 * K382683432[0];
  1649.      st1 = tmp304 + st1;
  1650.      tmp315 = st1;
  1651.      st1 = st5 * K382683432[0];
  1652.      tmp316 = st7;
  1653.      st7 = tmp302 - st4;
  1654.      st5 = st5 * K923879532[0];
  1655.      st4 = tmp302 + st4;
  1656.      tmp317 = st4;
  1657.      st4 = tmp301 - st8;
  1658.      st8 = tmp301 + st8;
  1659.      tmp318 = st8;
  1660.      st8 = c_re(input[5]);
  1661.      st8 = st8 - c_re(input[13]);
  1662.      tmp319 = st4;
  1663.      st4 = c_im(input[1]);
  1664.      st4 = st4 - c_im(input[9]);
  1665.      tmp320 = st7;
  1666.      st7 = c_re(input[3]);
  1667.      st7 = st7 - c_re(input[11]);
  1668.      tmp321 = st2;
  1669.      st2 = st8 + st4;
  1670.      st4 = st4 - st8;
  1671.      st8 = c_im(input[7]);
  1672.      tmp322 = st5;
  1673.      st5 = st2 * K382683432[0];
  1674.      st8 = st8 - c_im(input[15]);
  1675.      st2 = st2 * K923879532[0];
  1676.      st6 = st6 - st5;
  1677.      st5 = st4 * K923879532[0];
  1678.      st2 = st2 + st3;
  1679.      st4 = st4 * K382683432[0];
  1680.      st1 = st1 - st5;
  1681.      st3 = st7 - st8;
  1682.      st4 = st4 + tmp322;
  1683.      st7 = st7 + st8;
  1684.      st8 = st3 * K382683432[0];
  1685.      st5 = c_re(input[7]);
  1686.      st3 = st3 * K923879532[0];
  1687.      st5 = st5 - c_re(input[15]);
  1688.      tmp323 = st4;
  1689.      st4 = st7 * K923879532[0];
  1690.      tmp324 = st1;
  1691.      st1 = c_im(input[3]);
  1692.      st7 = st7 * K382683432[0];
  1693.      st1 = st1 - c_im(input[11]);
  1694.      tmp325 = st2;
  1695.      st2 = c_re(input[1]);
  1696.      st2 = st2 + c_re(input[9]);
  1697.      tmp326 = st6;
  1698.      st6 = st5 + st1;
  1699.      st1 = st1 - st5;
  1700.      st5 = c_re(input[5]);
  1701.      tmp327 = st7;
  1702.      st7 = st6 * K923879532[0];
  1703.      st5 = st5 + c_re(input[13]);
  1704.      st6 = st6 * K382683432[0];
  1705.      st8 = st8 - st7;
  1706.      st7 = st1 * K382683432[0];
  1707.      st6 = st6 + st3;
  1708.      st1 = st1 * K923879532[0];
  1709.      st7 = st7 - st4;
  1710.      st3 = st2 + st5;
  1711.      st1 = st1 + tmp327;
  1712.      st2 = st2 - st5;
  1713.      st4 = tmp326 - st8;
  1714.      st8 = tmp326 + st8;
  1715.      st5 = tmp325 - st6;
  1716.      tmp328 = st2;
  1717.      st2 = tmp321 - st4;
  1718.      st4 = tmp321 + st4;
  1719.      tmp329 = st3;
  1720.      st3 = tmp314 - st8;
  1721.      st8 = tmp314 + st8;
  1722.      tmp330 = st2;
  1723.      st2 = tmp316 - st5;
  1724.      st5 = tmp316 + st5;
  1725.      c_re(input[9]) = st3;
  1726.      c_re(input[1]) = st8;
  1727.      c_re(input[5]) = st2;
  1728.      c_re(input[13]) = st5;
  1729.      st6 = tmp325 + st6;
  1730.      st3 = tmp324 - st7;
  1731.      st8 = tmp323 - st1;
  1732.      st2 = tmp315 - st6;
  1733.      st6 = tmp315 + st6;
  1734.      st5 = tmp320 - st3;
  1735.      st3 = tmp320 + st3;
  1736.      tmp331 = st5;
  1737.      st5 = tmp317 - st8;
  1738.      st8 = tmp317 + st8;
  1739.      st7 = tmp324 + st7;
  1740.      st1 = tmp323 + st1;
  1741.      tmp332 = st3;
  1742.      st3 = c_im(input[1]);
  1743.      c_im(input[1]) = st6;
  1744.      st3 = st3 + c_im(input[9]);
  1745.      c_im(input[9]) = st2;
  1746.      st2 = tmp319 - st7;
  1747.      st7 = tmp319 + st7;
  1748.      st6 = tmp318 - st1;
  1749.      st1 = tmp318 + st1;
  1750.      tmp333 = st5;
  1751.      st5 = c_im(input[5]);
  1752.      c_im(input[5]) = st4;
  1753.      st5 = st5 + c_im(input[13]);
  1754.      c_im(input[13]) = tmp330;
  1755.      st4 = st3 - st5;
  1756.      st3 = st3 + st5;
  1757.      st5 = c_re(input[3]);
  1758.      c_re(input[3]) = st7;
  1759.      st5 = st5 + c_re(input[11]);
  1760.      c_re(input[11]) = st2;
  1761.      st2 = c_re(input[7]);
  1762.      c_re(input[7]) = st6;
  1763.      st2 = st2 + c_re(input[15]);
  1764.      c_re(input[15]) = st1;
  1765.      st7 = st5 + st2;
  1766.      st5 = st5 - st2;
  1767.      st6 = c_im(input[3]);
  1768.      c_im(input[3]) = st8;
  1769.      st6 = st6 + c_im(input[11]);
  1770.      c_im(input[11]) = tmp333;
  1771.      st8 = tmp329 + st7;
  1772.      st7 = tmp329 - st7;
  1773.      st1 = st5 + st4;
  1774.      st4 = st4 - st5;
  1775.      st2 = tmp313 - st8;
  1776.      st8 = tmp313 + st8;
  1777.      st5 = st7 + tmp312;
  1778.      st7 = tmp312 - st7;
  1779.      c_re(input[8]) = st2;
  1780.      c_re(input[0]) = st8;
  1781.      c_im(input[4]) = st5;
  1782.      c_im(input[12]) = st7;
  1783.      st2 = c_im(input[7]);
  1784.      c_im(input[7]) = tmp332;
  1785.      st2 = st2 + c_im(input[15]);
  1786.      c_im(input[15]) = tmp331;
  1787.      st8 = st6 - st2;
  1788.      st6 = st6 + st2;
  1789.      st5 = tmp328 - st8;
  1790.      st8 = tmp328 + st8;
  1791.      st7 = st3 - st6;
  1792.      st2 = st5 - st1;
  1793.      st1 = st1 + st5;
  1794.      st3 = st3 + st6;
  1795.      st2 = st2 * K707106781[0];
  1796.      st6 = st4 + st8;
  1797.      st1 = st1 * K707106781[0];
  1798.      st8 = st8 - st4;
  1799.      st6 = st6 * K707106781[0];
  1800.      st4 = tmp311 + st3;
  1801.      st8 = st8 * K707106781[0];
  1802.      st3 = tmp311 - st3;
  1803.      st5 = tmp310 - st2;
  1804.      c_im(input[0]) = st4;
  1805.      c_im(input[8]) = st3;
  1806.      c_re(input[10]) = st5;
  1807.      st2 = tmp310 + st2;
  1808.      st4 = tmp309 + st7;
  1809.      st7 = tmp309 - st7;
  1810.      st3 = tmp308 - st6;
  1811.      c_re(input[2]) = st2;
  1812.      c_re(input[12]) = st4;
  1813.      c_re(input[4]) = st7;
  1814.      c_re(input[6]) = st3;
  1815.      st6 = tmp308 + st6;
  1816.      st5 = tmp307 - st8;
  1817.      st8 = tmp307 + st8;
  1818.      st2 = tmp306 - st1;
  1819.      c_re(input[14]) = st6;
  1820.      c_im(input[14]) = st5;
  1821.      c_im(input[6]) = st8;
  1822.      c_im(input[10]) = st2;
  1823.      st1 = tmp306 + st1;
  1824.      c_im(input[2]) = st1;
  1825. }
  1827. void PFFTWI(32) (fftw_complex * input) {
  1828.      fftw_real tmp714;
  1829.      fftw_real tmp713;
  1830.      fftw_real tmp712;
  1831.      fftw_real tmp711;
  1832.      fftw_real tmp710;
  1833.      fftw_real tmp709;
  1834.      fftw_real tmp708;
  1835.      fftw_real tmp707;
  1836.      fftw_real tmp706;
  1837.      fftw_real tmp705;
  1838.      fftw_real tmp704;
  1839.      fftw_real tmp703;
  1840.      fftw_real tmp702;
  1841.      fftw_real tmp701;
  1842.      fftw_real tmp700;
  1843.      fftw_real tmp699;
  1844.      fftw_real tmp698;
  1845.      fftw_real tmp697;
  1846.      fftw_real tmp696;
  1847.      fftw_real tmp695;
  1848.      fftw_real tmp694;
  1849.      fftw_real tmp693;
  1850.      fftw_real tmp692;
  1851.      fftw_real tmp691;
  1852.      fftw_real tmp690;
  1853.      fftw_real tmp689;
  1854.      fftw_real tmp688;
  1855.      fftw_real tmp687;
  1856.      fftw_real tmp686;
  1857.      fftw_real tmp685;
  1858.      fftw_real tmp684;
  1859.      fftw_real tmp683;
  1860.      fftw_real tmp682;
  1861.      fftw_real tmp681;
  1862.      fftw_real tmp680;
  1863.      fftw_real tmp679;
  1864.      fftw_real tmp678;
  1865.      fftw_real tmp677;
  1866.      fftw_real tmp676;
  1867.      fftw_real tmp675;
  1868.      fftw_real tmp674;
  1869.      fftw_real tmp673;
  1870.      fftw_real tmp672;
  1871.      fftw_real tmp671;
  1872.      fftw_real tmp670;
  1873.      fftw_real tmp669;
  1874.      fftw_real tmp668;
  1875.      fftw_real tmp667;
  1876.      fftw_real tmp666;
  1877.      fftw_real tmp665;
  1878.      fftw_real tmp664;
  1879.      fftw_real tmp663;
  1880.      fftw_real tmp662;
  1881.      fftw_real tmp661;
  1882.      fftw_real tmp660;
  1883.      fftw_real tmp659;
  1884.      fftw_real tmp658;
  1885.      fftw_real tmp657;
  1886.      fftw_real tmp656;
  1887.      fftw_real tmp655;
  1888.      fftw_real tmp654;
  1889.      fftw_real tmp653;
  1890.      fftw_real tmp652;
  1891.      fftw_real tmp651;
  1892.      fftw_real tmp650;
  1893.      fftw_real tmp649;
  1894.      fftw_real tmp648;
  1895.      fftw_real tmp647;
  1896.      fftw_real tmp646;
  1897.      fftw_real tmp645;
  1898.      fftw_real tmp644;
  1899.      fftw_real tmp643;
  1900.      fftw_real tmp642;
  1901.      fftw_real tmp641;
  1902.      fftw_real tmp640;
  1903.      fftw_real tmp639;
  1904.      fftw_real tmp638;
  1905.      fftw_real tmp637;
  1906.      fftw_real tmp636;
  1907.      fftw_real tmp635;
  1908.      fftw_real tmp634;
  1909.      fftw_real tmp633;
  1910.      fftw_real tmp632;
  1911.      fftw_real tmp631;
  1912.      fftw_real tmp630;
  1913.      fftw_real tmp629;
  1914.      fftw_real tmp628;
  1915.      fftw_real tmp627;
  1916.      fftw_real tmp626;
  1917.      fftw_real tmp625;
  1918.      fftw_real tmp624;
  1919.      fftw_real tmp623;
  1920.      fftw_real tmp622;
  1921.      fftw_real tmp621;
  1922.      fftw_real st1;
  1923.      fftw_real st2;
  1924.      fftw_real st3;
  1925.      fftw_real st4;
  1926.      fftw_real st5;
  1927.      fftw_real st6;
  1928.      fftw_real st7;
  1929.      fftw_real st8;
  1930.      st8 = c_re(input[8]);
  1931.      st8 = st8 - c_re(input[24]);
  1932.      st7 = c_im(input[0]);
  1933.      st7 = st7 - c_im(input[16]);
  1934.      st6 = st8 + st7;
  1935.      st7 = st7 - st8;
  1936.      st5 = c_re(input[0]);
  1937.      st5 = st5 + c_re(input[16]);
  1938.      st4 = c_re(input[8]);
  1939.      st4 = st4 + c_re(input[24]);
  1940.      st3 = st5 + st4;
  1941.      st5 = st5 - st4;
  1942.      st2 = c_im(input[0]);
  1943.      st2 = st2 + c_im(input[16]);
  1944.      st1 = c_im(input[8]);
  1945.      st1 = st1 + c_im(input[24]);
  1946.      st8 = st2 + st1;
  1947.      st2 = st2 - st1;
  1948.      st4 = c_re(input[4]);
  1949.      st4 = st4 + c_re(input[20]);
  1950.      st1 = c_re(input[28]);
  1951.      st1 = st1 + c_re(input[12]);
  1952.      tmp621 = st6;
  1953.      st6 = st4 + st1;
  1954.      st4 = st4 - st1;
  1955.      st1 = st3 + st6;
  1956.      st3 = st3 - st6;
  1957.      st6 = st2 - st4;
  1958.      st4 = st4 + st2;
  1959.      st2 = c_im(input[4]);
  1960.      st2 = st2 + c_im(input[20]);
  1961.      tmp622 = st4;
  1962.      st4 = c_im(input[28]);
  1963.      st4 = st4 + c_im(input[12]);
  1964.      tmp623 = st6;
  1965.      st6 = st2 + st4;
  1966.      st4 = st4 - st2;
  1967.      st2 = st8 + st6;
  1968.      st8 = st8 - st6;
  1969.      st6 = st5 - st4;
  1970.      st5 = st5 + st4;
  1971.      st4 = c_re(input[0]);
  1972.      st4 = st4 - c_re(input[16]);
  1973.      tmp624 = st5;
  1974.      st5 = c_im(input[8]);
  1975.      st5 = st5 - c_im(input[24]);
  1976.      tmp625 = st6;
  1977.      st6 = st4 - st5;
  1978.      st4 = st4 + st5;
  1979.      st5 = c_re(input[4]);
  1980.      st5 = st5 - c_re(input[20]);
  1981.      tmp626 = st3;
  1982.      st3 = c_im(input[4]);
  1983.      st3 = st3 - c_im(input[20]);
  1984.      tmp627 = st2;
  1985.      st2 = st5 + st3;
  1986.      st5 = st5 - st3;
  1987.      st3 = c_im(input[28]);
  1988.      st3 = st3 - c_im(input[12]);
  1989.      tmp628 = st8;
  1990.      st8 = c_re(input[28]);
  1991.      st8 = st8 - c_re(input[12]);
  1992.      tmp629 = st1;
  1993.      st1 = st3 - st8;
  1994.      st8 = st8 + st3;
  1995.      st3 = st2 + st1;
  1996.      st3 = st3 * K707106781[0];
  1997.      st1 = st1 - st2;
  1998.      st1 = st1 * K707106781[0];
  1999.      st2 = st5 + st8;
  2000.      st2 = st2 * K707106781[0];
  2001.      st5 = st5 - st8;
  2002.      st5 = st5 * K707106781[0];
  2003.      st8 = st7 - st5;
  2004.      tmp630 = st8;
  2005.      st8 = st4 - st1;
  2006.      tmp631 = st8;
  2007.      st8 = tmp621 - st3;
  2008.      tmp632 = st8;
  2009.      st8 = st6 - st2;
  2010.      st3 = tmp621 + st3;
  2011.      st6 = st6 + st2;
  2012.      st7 = st7 + st5;
  2013.      st4 = st4 + st1;
  2014.      st1 = c_re(input[2]);
  2015.      st1 = st1 + c_re(input[18]);
  2016.      st2 = c_re(input[10]);
  2017.      st2 = st2 + c_re(input[26]);
  2018.      st5 = st1 + st2;
  2019.      st1 = st1 - st2;
  2020.      st2 = c_im(input[2]);
  2021.      st2 = st2 + c_im(input[18]);
  2022.      tmp633 = st4;
  2023.      st4 = c_im(input[10]);
  2024.      st4 = st4 + c_im(input[26]);
  2025.      tmp634 = st7;
  2026.      st7 = st2 + st4;
  2027.      st2 = st2 - st4;
  2028.      st4 = st1 + st2;
  2029.      st1 = st1 - st2;
  2030.      st2 = c_re(input[30]);
  2031.      st2 = st2 + c_re(input[14]);
  2032.      tmp635 = st6;
  2033.      st6 = c_re(input[6]);
  2034.      st6 = st6 + c_re(input[22]);
  2035.      tmp636 = st3;
  2036.      st3 = st2 + st6;
  2037.      st2 = st2 - st6;
  2038.      st6 = st5 + st3;
  2039.      st5 = st5 - st3;
  2040.      st3 = tmp629 + st6;
  2041.      st6 = tmp629 - st6;
  2042.      tmp637 = st6;
  2043.      st6 = tmp628 - st5;
  2044.      st5 = st5 + tmp628;
  2045.      tmp638 = st5;
  2046.      st5 = c_im(input[30]);
  2047.      st5 = st5 + c_im(input[14]);
  2048.      tmp639 = st6;
  2049.      st6 = c_im(input[6]);
  2050.      st6 = st6 + c_im(input[22]);
  2051.      tmp640 = st3;
  2052.      st3 = st5 + st6;
  2053.      st5 = st5 - st6;
  2054.      st6 = st7 + st3;
  2055.      st3 = st3 - st7;
  2056.      st7 = st5 - st2;
  2057.      tmp641 = st8;
  2058.      st8 = st7 - st4;
  2059.      st8 = st8 * K707106781[0];
  2060.      st4 = st4 + st7;
  2061.      st4 = st4 * K707106781[0];
  2062.      st2 = st2 + st5;
  2063.      st5 = st1 - st2;
  2064.      st5 = st5 * K707106781[0];
  2065.      st1 = st1 + st2;
  2066.      st1 = st1 * K707106781[0];
  2067.      st7 = tmp627 - st6;
  2068.      st6 = tmp627 + st6;
  2069.      st2 = tmp626 + st3;
  2070.      st3 = tmp626 - st3;
  2071.      tmp642 = st3;
  2072.      st3 = tmp623 + st5;
  2073.      st5 = tmp623 - st5;
  2074.      tmp643 = st5;
  2075.      st5 = tmp625 - st8;
  2076.      st8 = tmp625 + st8;
  2077.      tmp644 = st8;
  2078.      st8 = tmp622 + st4;
  2079.      st4 = tmp622 - st4;
  2080.      tmp645 = st4;
  2081.      st4 = tmp624 - st1;
  2082.      st1 = tmp624 + st1;
  2083.      tmp646 = st1;
  2084.      st1 = c_re(input[6]);
  2085.      st1 = st1 - c_re(input[22]);
  2086.      tmp647 = st4;
  2087.      st4 = c_im(input[30]);
  2088.      st4 = st4 - c_im(input[14]);
  2089.      tmp648 = st8;
  2090.      st8 = st1 + st4;
  2091.      tmp649 = st5;
  2092.      st5 = st8 * K382683432[0];
  2093.      st4 = st4 - st1;
  2094.      st8 = st8 * K923879532[0];
  2095.      st1 = c_re(input[30]);
  2096.      tmp650 = st3;
  2097.      st3 = st4 * K923879532[0];
  2098.      st1 = st1 - c_re(input[14]);
  2099.      st4 = st4 * K382683432[0];
  2100.      tmp651 = st2;
  2101.      st2 = c_im(input[6]);
  2102.      st2 = st2 - c_im(input[22]);
  2103.      tmp652 = st6;
  2104.      st6 = st1 - st2;
  2105.      tmp653 = st7;
  2106.      st7 = st6 * K923879532[0];
  2107.      st1 = st1 + st2;
  2108.      st6 = st6 * K382683432[0];
  2109.      st5 = st5 + st7;
  2110.      st2 = st1 * K382683432[0];
  2111.      st8 = st8 - st6;
  2112.      st1 = st1 * K923879532[0];
  2113.      st3 = st3 + st2;
  2114.      st4 = st4 - st1;
  2115.      st7 = c_re(input[2]);
  2116.      st7 = st7 - c_re(input[18]);
  2117.      st6 = c_im(input[10]);
  2118.      st6 = st6 - c_im(input[26]);
  2119.      st2 = st7 - st6;
  2120.      st1 = st2 * K923879532[0];
  2121.      st7 = st7 + st6;
  2122.      st2 = st2 * K382683432[0];
  2123.      st6 = c_re(input[10]);
  2124.      tmp654 = st8;
  2125.      st8 = st7 * K382683432[0];
  2126.      st6 = st6 - c_re(input[26]);
  2127.      st7 = st7 * K923879532[0];
  2128.      tmp655 = st5;
  2129.      st5 = c_im(input[2]);
  2130.      st5 = st5 - c_im(input[18]);
  2131.      tmp656 = st4;
  2132.      st4 = st6 + st5;
  2133.      tmp657 = st3;
  2134.      st3 = st4 * K382683432[0];
  2135.      st5 = st5 - st6;
  2136.      st4 = st4 * K923879532[0];
  2137.      st1 = st1 - st3;
  2138.      st6 = st5 * K923879532[0];
  2139.      st4 = st4 + st2;
  2140.      st5 = st5 * K382683432[0];
  2141.      st8 = st8 - st6;
  2142.      st5 = st5 + st7;
  2143.      st2 = st8 - tmp657;
  2144.      st7 = tmp630 + st2;
  2145.      st2 = tmp630 - st2;
  2146.      st3 = tmp656 - st5;
  2147.      st6 = tmp631 - st3;
  2148.      st3 = tmp631 + st3;
  2149.      tmp658 = st3;
  2150.      st3 = st1 - tmp655;
  2151.      tmp659 = st2;
  2152.      st2 = tmp632 - st3;
  2153.      st3 = tmp632 + st3;
  2154.      tmp660 = st3;
  2155.      st3 = tmp654 - st4;
  2156.      tmp661 = st2;
  2157.      st2 = tmp641 + st3;
  2158.      st3 = tmp641 - st3;
  2159.      st4 = st4 + tmp654;
  2160.      tmp662 = st3;
  2161.      st3 = tmp636 - st4;
  2162.      st4 = tmp636 + st4;
  2163.      st1 = st1 + tmp655;
  2164.      tmp663 = st4;
  2165.      st4 = tmp635 + st1;
  2166.      st1 = tmp635 - st1;
  2167.      st5 = st5 + tmp656;
  2168.      tmp664 = st1;
  2169.      st1 = tmp634 + st5;
  2170.      st5 = tmp634 - st5;
  2171.      st8 = st8 + tmp657;
  2172.      tmp665 = st5;
  2173.      st5 = tmp633 - st8;
  2174.      st8 = tmp633 + st8;
  2175.      tmp666 = st8;
  2176.      st8 = c_re(input[1]);
  2177.      st8 = st8 + c_re(input[17]);
  2178.      tmp667 = st5;
  2179.      st5 = c_re(input[9]);
  2180.      st5 = st5 + c_re(input[25]);
  2181.      tmp668 = st1;
  2182.      st1 = st8 + st5;
  2183.      st8 = st8 - st5;
  2184.      st5 = c_im(input[1]);
  2185.      st5 = st5 + c_im(input[17]);
  2186.      tmp669 = st4;
  2187.      st4 = c_im(input[9]);
  2188.      st4 = st4 + c_im(input[25]);
  2189.      tmp670 = st3;
  2190.      st3 = st5 + st4;
  2191.      st5 = st5 - st4;
  2192.      st4 = c_re(input[5]);
  2193.      st4 = st4 + c_re(input[21]);
  2194.      tmp671 = st2;
  2195.      st2 = c_re(input[29]);
  2196.      st2 = st2 + c_re(input[13]);
  2197.      tmp672 = st6;
  2198.      st6 = st4 + st2;
  2199.      st4 = st4 - st2;
  2200.      st2 = st1 + st6;
  2201.      st1 = st1 - st6;
  2202.      st6 = st4 + st5;
  2203.      tmp673 = st7;
  2204.      st7 = st6 * K923879532[0];
  2205.      st5 = st5 - st4;
  2206.      st6 = st6 * K382683432[0];
  2207.      st4 = c_im(input[5]);
  2208.      tmp674 = st2;
  2209.      st2 = st5 * K382683432[0];
  2210.      st4 = st4 + c_im(input[21]);
  2211.      st5 = st5 * K923879532[0];
  2212.      tmp675 = st7;
  2213.      st7 = c_im(input[29]);
  2214.      st7 = st7 + c_im(input[13]);
  2215.      tmp676 = st5;
  2216.      st5 = st4 + st7;
  2217.      st7 = st7 - st4;
  2218.      st4 = st3 + st5;
  2219.      st3 = st3 - st5;
  2220.      st5 = st1 - st3;
  2221.      st1 = st1 + st3;
  2222.      st3 = st8 + st7;
  2223.      tmp677 = st1;
  2224.      st1 = st3 * K382683432[0];
  2225.      st8 = st8 - st7;
  2226.      st3 = st3 * K923879532[0];
  2227.      st3 = st3 - st6;
  2228.      st6 = st8 * K923879532[0];
  2229.      st2 = st2 + st6;
  2230.      st8 = st8 * K382683432[0];
  2231.      st8 = st8 - tmp676;
  2232.      st1 = tmp675 + st1;
  2233.      st7 = c_re(input[31]);
  2234.      st7 = st7 + c_re(input[15]);
  2235.      st6 = c_re(input[7]);
  2236.      st6 = st6 + c_re(input[23]);
  2237.      tmp678 = st1;
  2238.      st1 = st7 + st6;
  2239.      st7 = st7 - st6;
  2240.      st6 = c_im(input[31]);
  2241.      st6 = st6 + c_im(input[15]);
  2242.      tmp679 = st3;
  2243.      st3 = c_im(input[7]);
  2244.      st3 = st3 + c_im(input[23]);
  2245.      tmp680 = st2;
  2246.      st2 = st6 + st3;
  2247.      st6 = st6 - st3;
  2248.      st3 = c_re(input[3]);
  2249.      st3 = st3 + c_re(input[19]);
  2250.      tmp681 = st8;
  2251.      st8 = c_re(input[27]);
  2252.      st8 = st8 + c_re(input[11]);
  2253.      tmp682 = st5;
  2254.      st5 = st3 + st8;
  2255.      st3 = st3 - st8;
  2256.      st8 = st1 + st5;
  2257.      st1 = st1 - st5;
  2258.      st5 = st3 + st6;
  2259.      tmp683 = st4;
  2260.      st4 = st5 * K923879532[0];
  2261.      st6 = st6 - st3;
  2262.      st5 = st5 * K382683432[0];
  2263.      st3 = tmp674 + st8;
  2264.      tmp684 = st4;
  2265.      st4 = st6 * K382683432[0];
  2266.      tmp685 = st5;
  2267.      st5 = tmp640 - st3;
  2268.      st6 = st6 * K923879532[0];
  2269.      st3 = tmp640 + st3;
  2270.      st8 = tmp674 - st8;
  2271.      c_re(input[16]) = st5;
  2272.      st5 = st8 + tmp653;
  2273.      st8 = tmp653 - st8;
  2274.      c_re(input[0]) = st3;
  2275.      st3 = c_im(input[3]);
  2276.      st3 = st3 + c_im(input[19]);
  2277.      c_im(input[8]) = st5;
  2278.      st5 = c_im(input[27]);
  2279.      st5 = st5 + c_im(input[11]);
  2280.      c_im(input[24]) = st8;
  2281.      st8 = st3 + st5;
  2282.      st5 = st5 - st3;
  2283.      st3 = st2 + st8;
  2284.      st2 = st2 - st8;
  2285.      st8 = st1 + st2;
  2286.      st2 = st2 - st1;
  2287.      st1 = st7 + st5;
  2288.      tmp686 = st4;
  2289.      st4 = st1 * K382683432[0];
  2290.      st7 = st7 - st5;
  2291.      st1 = st1 * K923879532[0];
  2292.      st5 = tmp683 + st3;
  2293.      tmp687 = st4;
  2294.      st4 = st7 * K923879532[0];
  2295.      tmp688 = st1;
  2296.      st1 = tmp652 - st5;
  2297.      st7 = st7 * K382683432[0];
  2298.      st5 = tmp652 + st5;
  2299.      st3 = st3 - tmp683;
  2300.      c_im(input[16]) = st1;
  2301.      st1 = tmp637 - st3;
  2302.      st3 = tmp637 + st3;
  2303.      c_im(input[0]) = st5;
  2304.      st5 = tmp682 - st8;
  2305.      st5 = st5 * K707106781[0];
  2306.      c_re(input[24]) = st1;
  2307.      st1 = tmp639 - st5;
  2308.      st5 = tmp639 + st5;
  2309.      st8 = tmp682 + st8;
  2310.      st8 = st8 * K707106781[0];
  2311.      c_re(input[8]) = st3;
  2312.      st3 = tmp651 - st8;
  2313.      st8 = tmp651 + st8;
  2314.      c_im(input[28]) = st1;
  2315.      st1 = st2 - tmp677;
  2316.      st1 = st1 * K707106781[0];
  2317.      c_im(input[12]) = st5;
  2318.      st5 = tmp642 - st1;
  2319.      st1 = tmp642 + st1;
  2320.      st2 = tmp677 + st2;
  2321.      st2 = st2 * K707106781[0];
  2322.      c_re(input[20]) = st3;
  2323.      st3 = tmp638 - st2;
  2324.      st2 = tmp638 + st2;
  2325.      st6 = st6 + st7;
  2326.      st7 = tmp681 - st6;
  2327.      st6 = tmp681 + st6;
  2328.      st4 = tmp686 - st4;
  2329.      c_re(input[4]) = st8;
  2330.      st8 = tmp680 + st4;
  2331.      st4 = st4 - tmp680;
  2332.      c_re(input[28]) = st5;
  2333.      st5 = tmp650 - st8;
  2334.      st8 = tmp650 + st8;
  2335.      c_re(input[12]) = st1;
  2336.      st1 = tmp649 - st4;
  2337.      st4 = tmp649 + st4;
  2338.      c_im(input[20]) = st3;
  2339.      st3 = tmp643 - st7;
  2340.      st7 = tmp643 + st7;
  2341.      c_im(input[4]) = st2;
  2342.      st2 = tmp644 - st6;
  2343.      st6 = tmp644 + st6;
  2344.      c_im(input[22]) = st5;
  2345.      st5 = tmp685 + tmp688;
  2346.      c_im(input[6]) = st8;
  2347.      st8 = tmp679 - st5;
  2348.      st5 = tmp679 + st5;
  2349.      c_re(input[30]) = st1;
  2350.      st1 = tmp684 - tmp687;
  2351.      c_re(input[14]) = st4;
  2352.      st4 = tmp678 + st1;
  2353.      st1 = st1 - tmp678;
  2354.      c_im(input[30]) = st3;
  2355.      st3 = tmp648 - st4;
  2356.      st4 = tmp648 + st4;
  2357.      c_im(input[14]) = st7;
  2358.      st7 = tmp647 - st1;
  2359.      st1 = tmp647 + st1;
  2360.      c_re(input[22]) = st2;
  2361.      st2 = tmp645 - st8;
  2362.      st8 = tmp645 + st8;
  2363.      c_re(input[6]) = st6;
  2364.      st6 = tmp646 - st5;
  2365.      st5 = tmp646 + st5;
  2366.      c_im(input[18]) = st3;
  2367.      st3 = c_re(input[31]);
  2368.      st3 = st3 - c_re(input[15]);
  2369.      c_im(input[2]) = st4;
  2370.      st4 = c_im(input[7]);
  2371.      st4 = st4 - c_im(input[23]);
  2372.      c_re(input[26]) = st7;
  2373.      st7 = st3 - st4;
  2374.      st3 = st3 + st4;
  2375.      c_re(input[10]) = st1;
  2376.      st1 = c_re(input[7]);
  2377.      st1 = st1 - c_re(input[23]);
  2378.      c_im(input[26]) = st2;
  2379.      st2 = c_im(input[31]);
  2380.      st2 = st2 - c_im(input[15]);
  2381.      c_im(input[10]) = st8;
  2382.      st8 = st1 + st2;
  2383.      st2 = st2 - st1;
  2384.      c_re(input[18]) = st6;
  2385.      st6 = c_re(input[3]);
  2386.      st6 = st6 - c_re(input[19]);
  2387.      c_re(input[2]) = st5;
  2388.      st5 = c_im(input[3]);
  2389.      st5 = st5 - c_im(input[19]);
  2390.      st4 = st6 - st5;
  2391.      st6 = st6 + st5;
  2392.      st1 = c_re(input[27]);
  2393.      st1 = st1 - c_re(input[11]);
  2394.      st5 = c_im(input[27]);
  2395.      st5 = st5 - c_im(input[11]);
  2396.      tmp689 = st2;
  2397.      st2 = st1 + st5;
  2398.      st5 = st5 - st1;
  2399.      st1 = st4 + st2;
  2400.      st1 = st1 * K707106781[0];
  2401.      st4 = st4 - st2;
  2402.      st4 = st4 * K707106781[0];
  2403.      st2 = st6 + st5;
  2404.      st2 = st2 * K707106781[0];
  2405.      st5 = st5 - st6;
  2406.      st5 = st5 * K707106781[0];
  2407.      st6 = st7 - st1;
  2408.      tmp690 = st4;
  2409.      st4 = st6 * K831469612[0];
  2410.      st7 = st7 + st1;
  2411.      st6 = st6 * K555570233[0];
  2412.      st1 = st3 - st5;
  2413.      tmp691 = st6;
  2414.      st6 = st1 * K195090322[0];
  2415.      st3 = st3 + st5;
  2416.      st1 = st1 * K980785280[0];
  2417.      st5 = st8 - st2;
  2418.      tmp692 = st3;
  2419.      st3 = st5 * K555570233[0];
  2420.      st8 = st8 + st2;
  2421.      st5 = st5 * K831469612[0];
  2422.      st2 = tmp689 - tmp690;
  2423.      tmp693 = st5;
  2424.      st5 = st2 * K980785280[0];
  2425.      tmp694 = st4;
  2426.      st4 = tmp689 + tmp690;
  2427.      st2 = st2 * K195090322[0];
  2428.      st5 = st5 + st6;
  2429.      st6 = st8 * K980785280[0];
  2430.      st2 = st2 - st1;
  2431.      st1 = st7 * K195090322[0];
  2432.      st3 = st3 - tmp694;
  2433.      st7 = st7 * K980785280[0];
  2434.      tmp695 = st3;
  2435.      st3 = tmp691 + tmp693;
  2436.      st8 = st8 * K195090322[0];
  2437.      st6 = st6 - st1;
  2438.      st1 = st4 * K555570233[0];
  2439.      st7 = st7 + st8;
  2440.      st8 = tmp692 * K831469612[0];
  2441.      st1 = st1 + st8;
  2442.      st4 = st4 * K831469612[0];
  2443.      st8 = c_re(input[1]);
  2444.      tmp696 = st1;
  2445.      st1 = tmp692 * K555570233[0];
  2446.      st8 = st8 - c_re(input[17]);
  2447.      st4 = st4 - st1;
  2448.      st1 = c_im(input[9]);
  2449.      st1 = st1 - c_im(input[25]);
  2450.      tmp697 = st4;
  2451.      st4 = st8 - st1;
  2452.      st8 = st8 + st1;
  2453.      st1 = c_re(input[9]);
  2454.      st1 = st1 - c_re(input[25]);
  2455.      tmp698 = st7;
  2456.      st7 = c_im(input[1]);
  2457.      st7 = st7 - c_im(input[17]);
  2458.      tmp699 = st6;
  2459.      st6 = st1 + st7;
  2460.      st7 = st7 - st1;
  2461.      st1 = c_re(input[5]);
  2462.      st1 = st1 - c_re(input[21]);
  2463.      tmp700 = st3;
  2464.      st3 = c_im(input[5]);
  2465.      st3 = st3 - c_im(input[21]);
  2466.      tmp701 = st2;
  2467.      st2 = st1 - st3;
  2468.      st1 = st1 + st3;
  2469.      st3 = c_re(input[29]);
  2470.      st3 = st3 - c_re(input[13]);
  2471.      tmp702 = st5;
  2472.      st5 = c_im(input[29]);
  2473.      st5 = st5 - c_im(input[13]);
  2474.      tmp703 = st7;
  2475.      st7 = st3 + st5;
  2476.      st5 = st5 - st3;
  2477.      st3 = st2 + st7;
  2478.      st3 = st3 * K707106781[0];
  2479.      st2 = st2 - st7;
  2480.      st2 = st2 * K707106781[0];
  2481.      st7 = st1 + st5;
  2482.      st7 = st7 * K707106781[0];
  2483.      st5 = st5 - st1;
  2484.      st5 = st5 * K707106781[0];
  2485.      st1 = st4 - st3;
  2486.      tmp704 = st2;
  2487.      st2 = st1 * K831469612[0];
  2488.      st4 = st4 + st3;
  2489.      st1 = st1 * K555570233[0];
  2490.      st3 = st8 - st5;
  2491.      tmp705 = st1;
  2492.      st1 = st3 * K195090322[0];
  2493.      st8 = st8 + st5;
  2494.      st3 = st3 * K980785280[0];
  2495.      st5 = st6 - st7;
  2496.      tmp706 = st2;
  2497.      st2 = st5 * K555570233[0];
  2498.      st6 = st6 + st7;
  2499.      st5 = st5 * K831469612[0];
  2500.      st7 = tmp703 - tmp704;
  2501.      tmp707 = st5;
  2502.      st5 = st7 * K980785280[0];
  2503.      tmp708 = st2;
  2504.      st2 = tmp703 + tmp704;
  2505.      st7 = st7 * K195090322[0];
  2506.      st1 = st1 - st5;
  2507.      st5 = st4 * K195090322[0];
  2508.      tmp709 = st5;
  2509.      st5 = st1 - tmp702;
  2510.      tmp710 = st5;
  2511.      st5 = st6 * K980785280[0];
  2512.      st1 = st1 + tmp702;
  2513.      st4 = st4 * K980785280[0];
  2514.      st7 = st7 + st3;
  2515.      st6 = st6 * K195090322[0];
  2516.      st3 = st7 + tmp701;
  2517.      tmp711 = st4;
  2518.      st4 = st8 * K831469612[0];
  2519.      st7 = tmp701 - st7;
  2520.      tmp712 = st4;
  2521.      st4 = st2 * K555570233[0];
  2522.      tmp713 = st4;
  2523.      st4 = tmp673 - st3;
  2524.      st2 = st2 * K831469612[0];
  2525.      st8 = st8 * K555570233[0];
  2526.      st3 = tmp673 + st3;
  2527.      c_im(input[23]) = st4;
  2528.      st4 = tmp672 - st7;
  2529.      st7 = tmp672 + st7;
  2530.      c_im(input[7]) = st3;
  2531.      st3 = tmp659 - tmp710;
  2532.      c_re(input[31]) = st4;
  2533.      st4 = tmp659 + tmp710;
  2534.      c_re(input[15]) = st7;
  2535.      st7 = tmp658 - st1;
  2536.      st1 = tmp658 + st1;
  2537.      c_im(input[31]) = st3;
  2538.      st3 = tmp706 + tmp708;
  2539.      c_im(input[15]) = st4;
  2540.      st4 = st3 + tmp695;
  2541.      st3 = tmp695 - st3;
  2542.      c_re(input[23]) = st7;
  2543.      st7 = tmp705 - tmp707;
  2544.      c_re(input[7]) = st1;
  2545.      st1 = st7 - tmp700;
  2546.      st7 = st7 + tmp700;
  2547.      tmp714 = st2;
  2548.      st2 = tmp661 - st1;
  2549.      st1 = tmp661 + st1;
  2550.      c_im(input[29]) = st2;
  2551.      st2 = tmp671 - st7;
  2552.      st7 = tmp671 + st7;
  2553.      c_im(input[13]) = st1;
  2554.      st1 = tmp660 - st4;
  2555.      st4 = tmp660 + st4;
  2556.      c_re(input[21]) = st2;
  2557.      st2 = tmp662 - st3;
  2558.      st3 = tmp662 + st3;
  2559.      st5 = tmp709 + st5;
  2560.      c_re(input[5]) = st7;
  2561.      st7 = st5 + tmp699;
  2562.      st5 = tmp699 - st5;
  2563.      st6 = tmp711 - st6;
  2564.      c_im(input[21]) = st1;
  2565.      st1 = st6 - tmp698;
  2566.      st6 = st6 + tmp698;
  2567.      c_im(input[5]) = st4;
  2568.      st4 = tmp670 - st1;
  2569.      st1 = tmp670 + st1;
  2570.      c_re(input[29]) = st2;
  2571.      st2 = tmp669 - st6;
  2572.      st6 = tmp669 + st6;
  2573.      c_re(input[13]) = st3;
  2574.      st3 = tmp663 - st7;
  2575.      st7 = tmp663 + st7;
  2576.      c_im(input[25]) = st4;
  2577.      st4 = tmp664 - st5;
  2578.      st5 = tmp664 + st5;
  2579.      c_im(input[9]) = st1;
  2580.      st1 = tmp712 - tmp713;
  2581.      c_re(input[17]) = st2;
  2582.      st2 = st1 - tmp696;
  2583.      st1 = st1 + tmp696;
  2584.      st8 = tmp714 + st8;
  2585.      c_re(input[1]) = st6;
  2586.      st6 = st8 + tmp697;
  2587.      st8 = tmp697 - st8;
  2588.      c_im(input[17]) = st3;
  2589.      st3 = tmp668 - st6;
  2590.      st6 = tmp668 + st6;
  2591.      c_im(input[1]) = st7;
  2592.      st7 = tmp667 - st8;
  2593.      st8 = tmp667 + st8;
  2594.      c_re(input[25]) = st4;
  2595.      st4 = tmp665 - st2;
  2596.      st2 = tmp665 + st2;
  2597.      c_re(input[9]) = st5;
  2598.      st5 = tmp666 - st1;
  2599.      st1 = tmp666 + st1;
  2600.      c_im(input[19]) = st3;
  2601.      c_im(input[3]) = st6;
  2602.      c_re(input[27]) = st7;
  2603.      c_re(input[11]) = st8;
  2604.      c_im(input[27]) = st4;
  2605.      c_im(input[11]) = st2;
  2606.      c_re(input[19]) = st5;
  2607.      c_re(input[3]) = st1;
  2608. }
  2610. void PFFTWI(64)(fftw_complex *input)
  2611. {
  2612.      PFFTWI(16)(input );
  2613.      PFFTWI(16)(input + 16);
  2614.      PFFTWI(16)(input + 32);
  2615.      PFFTWI(16)(input + 48);
  2616.      PFFTWI(twiddle_4)(input, PFFTW(W_64), 16);
  2617. }
  2619. void PFFTWI(128)(fftw_complex *input)
  2620. {
  2621.      PFFTWI(32)(input );
  2622.      PFFTWI(32)(input + 32);
  2623.      PFFTWI(32)(input + 64);
  2624.      PFFTWI(32)(input + 96);
  2625.      PFFTWI(twiddle_4)(input, PFFTW(W_128), 32);
  2626. }
  2628. void PFFTWI(512)(fftw_complex *input)
  2629. {
  2630.      PFFTWI(128)(input );
  2631.      PFFTWI(128)(input + 128);
  2632.      PFFTWI(128)(input + 256);
  2633.      PFFTWI(128)(input + 384);
  2634.      PFFTWI(twiddle_4)(input, PFFTW(W_512), 128);
  2635. }
  2637. ///////////////////////////////////////////////////////////////////////////
  2638. void  PFFTW(twiddle_4) (fftw_complex * A, const fftw_complex * W, int iostride) {
  2639.      int i;
  2640.      fftw_complex *inout;
  2641.      inout = A;
  2642.      {
  2643.       fftw_real st1;
  2644.       fftw_real st2;
  2645.       fftw_real st3;
  2646.       fftw_real st4;
  2647.       fftw_real st5;
  2648.       fftw_real st6;
  2649.       fftw_real st7;
  2650.       fftw_real st8;
  2651.       st8 = c_re(inout[0]);
  2652.       st8 = st8 + c_re(inout[2 * iostride]);
  2653.       st7 = c_re(inout[iostride]);
  2654.       st7 = st7 + c_re(inout[3 * iostride]);
  2655.       st6 = st8 - st7;
  2656.       st8 = st8 + st7;
  2657.       st5 = c_im(inout[0]);
  2658.       st5 = st5 + c_im(inout[2 * iostride]);
  2659.       st4 = c_im(inout[iostride]);
  2660.       st4 = st4 + c_im(inout[3 * iostride]);
  2661.       st3 = st5 - st4;
  2662.       st5 = st5 + st4;
  2663.       st2 = c_im(inout[0]);
  2664.       st2 = st2 - c_im(inout[2 * iostride]);
  2665.       st1 = c_re(inout[iostride]);
  2666.       st1 = st1 - c_re(inout[3 * iostride]);
  2667.       st7 = st2 - st1;
  2668.       st1 = st1 + st2;
  2669.       st4 = c_re(inout[0]);
  2670.       st4 = st4 - c_re(inout[2 * iostride]);
  2671.       c_re(inout[2 * iostride]) = st6;
  2672.       st6 = c_im(inout[iostride]);
  2673.       st6 = st6 - c_im(inout[3 * iostride]);
  2674.       c_re(inout[0]) = st8;
  2675.       st8 = st4 - st6;
  2676.       st4 = st4 + st6;
  2677.       c_im(inout[0]) = st5;
  2678.       c_im(inout[2 * iostride]) = st3;
  2679.       c_im(inout[iostride]) = st7;
  2680.       c_im(inout[3 * iostride]) = st1;
  2681.       c_re(inout[3 * iostride]) = st8;
  2682.       c_re(inout[iostride]) = st4;
  2683.      }
  2684.      inout = inout + 1;
  2685.      i = iostride - 1;
  2686.      do {
  2687.       {
  2688.            fftw_real st1;
  2689.            fftw_real st2;
  2690.            fftw_real st3;
  2691.            fftw_real st4;
  2692.            fftw_real st5;
  2693.            fftw_real st6;
  2694.            fftw_real st7;
  2695.            fftw_real st8;
  2696.            st8 = c_re(inout[0]);
  2697.            st8 = st8 + c_re(inout[2 * iostride]);
  2698.            st7 = c_re(inout[iostride]);
  2699.            st7 = st7 + c_re(inout[3 * iostride]);
  2700.            st6 = st8 - st7;
  2701.            st5 = st6 * c_im(W[1]);
  2702.            st8 = st8 + st7;
  2703.            st6 = st6 * c_re(W[1]);
  2704.            st4 = c_im(inout[0]);
  2705.            st4 = st4 + c_im(inout[2 * iostride]);
  2706.            st3 = c_im(inout[iostride]);
  2707.            st3 = st3 + c_im(inout[3 * iostride]);
  2708.            st2 = st4 - st3;
  2709.            st1 = st2 * c_im(W[1]);
  2710.            st4 = st4 + st3;
  2711.            st2 = st2 * c_re(W[1]);
  2712.            st2 = st2 - st5;
  2713.            st6 = st6 + st1;
  2714.            st7 = c_re(inout[0]);
  2715.            st7 = st7 - c_re(inout[2 * iostride]);
  2716.            st5 = c_im(inout[iostride]);
  2717.            st5 = st5 - c_im(inout[3 * iostride]);
  2718.            c_re(inout[0]) = st8;
  2719.            st8 = st7 - st5;
  2720.            st3 = st8 * c_re(W[0]);
  2721.            st7 = st7 + st5;
  2722.            st8 = st8 * c_im(W[0]);
  2723.            st1 = c_re(inout[iostride]);
  2724.            c_re(inout[2 * iostride]) = st6;
  2725.            st6 = st7 * c_im(W[0]);
  2726.            st1 = st1 - c_re(inout[3 * iostride]);
  2727.            st7 = st7 * c_re(W[0]);
  2728.            st5 = c_im(inout[0]);
  2729.            st5 = st5 - c_im(inout[2 * iostride]);
  2730.            c_im(inout[0]) = st4;
  2731.            st4 = st1 + st5;
  2732.            c_im(inout[2 * iostride]) = st2;
  2733.            st2 = st4 * c_im(W[0]);
  2734.            st5 = st5 - st1;
  2735.            st4 = st4 * c_re(W[0]);
  2736.            st3 = st3 - st2;
  2737.            st1 = st5 * c_re(W[0]);
  2738.            st5 = st5 * c_im(W[0]);
  2739.            st4 = st4 + st8;
  2740.            st5 = st5 + st7;
  2741.            st1 = st1 - st6;
  2742.            c_re(inout[3 * iostride]) = st3;
  2743.            c_im(inout[3 * iostride]) = st4;
  2744.            c_re(inout[iostride]) = st5;
  2745.            c_im(inout[iostride]) = st1;
  2746.       }
  2747.       i = i - 1, inout = inout + 1, W = W + 2;
  2748.      } while (i > 0);
  2749. }
  2751. void PFFTWI(twiddle_4) (fftw_complex * A, const fftw_complex * W, int iostride) {
  2752.      int i;
  2753.      fftw_complex *inout;
  2754.      inout = A;
  2755.      {
  2756.       fftw_real st1;
  2757.       fftw_real st2;
  2758.       fftw_real st3;
  2759.       fftw_real st4;
  2760.       fftw_real st5;
  2761.       fftw_real st6;
  2762.       fftw_real st7;
  2763.       fftw_real st8;
  2764.       st8 = c_re(inout[0]);
  2765.       st8 = st8 + c_re(inout[2 * iostride]);
  2766.       st7 = c_re(inout[iostride]);
  2767.       st7 = st7 + c_re(inout[3 * iostride]);
  2768.       st6 = st8 - st7;
  2769.       st8 = st8 + st7;
  2770.       st5 = c_im(inout[0]);
  2771.       st5 = st5 + c_im(inout[2 * iostride]);
  2772.       st4 = c_im(inout[iostride]);
  2773.       st4 = st4 + c_im(inout[3 * iostride]);
  2774.       st3 = st5 - st4;
  2775.       st5 = st5 + st4;
  2776.       st2 = c_re(inout[iostride]);
  2777.       st2 = st2 - c_re(inout[3 * iostride]);
  2778.       st1 = c_im(inout[0]);
  2779.       st1 = st1 - c_im(inout[2 * iostride]);
  2780.       st7 = st2 + st1;
  2781.       st1 = st1 - st2;
  2782.       st4 = c_re(inout[0]);
  2783.       st4 = st4 - c_re(inout[2 * iostride]);
  2784.       c_re(inout[2 * iostride]) = st6;
  2785.       st6 = c_im(inout[iostride]);
  2786.       st6 = st6 - c_im(inout[3 * iostride]);
  2787.       c_re(inout[0]) = st8;
  2788.       st8 = st4 - st6;
  2789.       st4 = st4 + st6;
  2790.       c_im(inout[0]) = st5;
  2791.       c_im(inout[2 * iostride]) = st3;
  2792.       c_im(inout[iostride]) = st7;
  2793.       c_im(inout[3 * iostride]) = st1;
  2794.       c_re(inout[iostride]) = st8;
  2795.       c_re(inout[3 * iostride]) = st4;
  2796.      }
  2797.      inout = inout + 1;
  2798.      i = iostride - 1;
  2799.      do {
  2800.       {
  2801.            fftw_real st1;
  2802.            fftw_real st2;
  2803.            fftw_real st3;
  2804.            fftw_real st4;
  2805.            fftw_real st5;
  2806.            fftw_real st6;
  2807.            fftw_real st7;
  2808.            fftw_real st8;
  2809.            st8 = c_re(inout[2 * iostride]);
  2810.            st8 = st8 * c_re(W[1]);
  2811.            st7 = c_im(inout[2 * iostride]);
  2812.            st7 = st7 * c_im(W[1]);
  2813.            st8 = st8 - st7;
  2814.            st6 = st8 + c_re(inout[0]);
  2815.            st8 = c_re(inout[0]) - st8;
  2816.            st5 = c_re(inout[2 * iostride]);
  2817.            st5 = st5 * c_im(W[1]);
  2818.            st4 = c_im(inout[2 * iostride]);
  2819.            st4 = st4 * c_re(W[1]);
  2820.            st5 = st5 + st4;
  2821.            st3 = st5 + c_im(inout[0]);
  2822.            st5 = c_im(inout[0]) - st5;
  2823.            st2 = c_re(inout[iostride]);
  2824.            st2 = st2 * c_re(W[0]);
  2825.            st1 = c_im(inout[iostride]);
  2826.            st1 = st1 * c_im(W[0]);
  2827.            st2 = st2 - st1;
  2828.            st7 = c_re(inout[3 * iostride]);
  2829.            st7 = st7 * c_re(W[0]);
  2830.            st4 = c_im(inout[3 * iostride]);
  2831.            st4 = st4 * c_im(W[0]);
  2832.            st7 = st7 + st4;
  2833.            st1 = st2 + st7;
  2834.            st2 = st2 - st7;
  2835.            st4 = st6 - st1;
  2836.            st6 = st6 + st1;
  2837.            st7 = st2 + st5;
  2838.            st5 = st5 - st2;
  2839.            st1 = c_re(inout[iostride]);
  2840.            st1 = st1 * c_im(W[0]);
  2841.            st2 = c_im(inout[iostride]);
  2842.            st2 = st2 * c_re(W[0]);
  2843.            st1 = st1 + st2;
  2844.            c_re(inout[2 * iostride]) = st4;
  2845.            st4 = c_im(inout[3 * iostride]);
  2846.            st4 = st4 * c_re(W[0]);
  2847.            c_re(inout[0]) = st6;
  2848.            st6 = c_re(inout[3 * iostride]);
  2849.            st6 = st6 * c_im(W[0]);
  2850.            st4 = st4 - st6;
  2851.            c_im(inout[iostride]) = st7;
  2852.            st7 = st1 - st4;
  2853.            st1 = st1 + st4;
  2854.            c_im(inout[3 * iostride]) = st5;
  2855.            st5 = st8 - st7;
  2856.            st8 = st8 + st7;
  2857.            st2 = st1 + st3;
  2858.            st3 = st3 - st1;
  2859.            c_re(inout[iostride]) = st5;
  2860.            c_re(inout[3 * iostride]) = st8;
  2861.            c_im(inout[0]) = st2;
  2862.            c_im(inout[2 * iostride]) = st3;
  2863.       }
  2864.       i = i - 1, inout = inout + 1, W = W + 2;
  2865.      } while (i > 0);
  2866. }
  2867. //////////////////////////////////////////////////////////////////////
  2868. int PFFTW(permutation_64)(int i)
  2869. {
  2870.     int i1 = i % 4;
  2871.     int i2 = i / 4;
  2872.     if (i1 <= 2)
  2873.        return (i1 * 16 + i2);
  2874.     else
  2875.        return (i1 * 16 + ((i2 + 1) % 16));
  2876. }
  2878. int PFFTW(permutation_128)(int i)
  2879. {
  2880.     int i1 = i % 4;
  2881.     int i2 = i / 4;
  2882.     if (i1 <= 2)
  2883.        return (i1 * 32 + i2);
  2884.     else
  2885.        return (i1 * 32 + ((i2 + 1) % 32));
  2886. }
  2888. int PFFTW(permutation_512)(int i)
  2889. {
  2890.     int i1 = i % 4;
  2891.     int i2 = i / 4;
  2892.     if (i1 <= 2)
  2893.        return (i1 * 128 + PFFTW(permutation_128)(i2));
  2894.     else
  2895.        return (i1 * 128 + PFFTW(permutation_128)((i2 + 1) % 128));
  2896. }
  2897. /////////////////////////////////////////////////////////////////
  2898. void MakeFFTOrder(faacDecHandle hDecoder)
  2899. {
  2900.     int i;
  2902.     for (i=0 ; i < 512 ; i++)
  2903.     {
  2904.         if (i < 64)
  2905.             hDecoder->unscambled64[i] = PFFTW(permutation_64)(i);
  2906.         hDecoder->unscambled512[i] = PFFTW(permutation_512)(i);
  2907.     }
  2908. }