← Back to branch summary

~ubuntu-branches/ubuntu/maverick/blender/maverick

« back to all changes in this revision

Viewing changes to extern/fftw/rdft/codelets/hc2r/hb_16.c

  • Committer: Bazaar Package Importer
  • Author(s): Khashayar Naderehvandi, Khashayar Naderehvandi, Alessio Treglia
  • Date: 2009-01-22 16:53:59 UTC
  • mfrom: (14.1.1 experimental)
  • Revision ID: james.westby@ubuntu.com-20090122165359-v0996tn7fbit64ni
Tags: 2.48a+dfsg-1ubuntu1
https://launchpad.net/bugs/320045
[ Khashayar Naderehvandi ]
* Merge from debian experimental (LP: #320045), Ubuntu remaining changes:
  - Add patch correcting header file locations.
  - Add libvorbis-dev and libgsm1-dev to Build-Depends.
  - Use avcodec_decode_audio2() in source/blender/src/hddaudio.c

[ Alessio Treglia ]
* Add missing previous changelog entries.

Show diffs side-by-side

added added

removed removed

Lines of Context:
extern/fftw/rdft/codelets/hc2r/hb_16.c
 
1
/*
 
2
 * Copyright (c) 2003, 2006 Matteo Frigo
 
3
 * Copyright (c) 2003, 2006 Massachusetts Institute of Technology
 
4
 *
 
5
 * This program is free software; you can redistribute it and/or modify
 
6
 * it under the terms of the GNU General Public License as published by
 
7
 * the Free Software Foundation; either version 2 of the License, or
 
8
 * (at your option) any later version.
 
9
 *
 
10
 * This program 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
 
13
 * GNU General Public License for more details.
 
14
 *
 
15
 * You should have received a copy of the GNU General Public License
 
16
 * along with this program; if not, write to the Free Software
 
17
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 
18
 *
 
19
 */
 
20
 
 
21
/* This file was automatically generated --- DO NOT EDIT */
 
22
/* Generated on Sun Jul  2 16:31:31 EDT 2006 */
 
23
 
 
24
#include "codelet-rdft.h"
 
25
 
 
26
#ifdef HAVE_FMA
 
27
 
 
28
/* Generated by: ../../../genfft/gen_hc2hc -fma -reorder-insns -schedule-for-pipeline -compact -variables 4 -pipeline-latency 4 -sign 1 -n 16 -dif -name hb_16 -include hb.h */
 
29
 
 
30
/*
 
31
 * This function contains 174 FP additions, 100 FP multiplications,
 
32
 * (or, 104 additions, 30 multiplications, 70 fused multiply/add),
 
33
 * 83 stack variables, and 64 memory accesses
 
34
 */
 
35
/*
 
36
 * Generator Id's : 
 
37
 * $Id: algsimp.ml,v 1.9 2006-02-12 23:34:12 athena Exp $
 
38
 * $Id: fft.ml,v 1.4 2006-01-05 03:04:27 stevenj Exp $
 
39
 * $Id: gen_hc2hc.ml,v 1.16 2006-02-12 23:34:12 athena Exp $
 
40
 */
 
41
 
 
42
#include "hb.h"
 
43
 
 
44
static const R *hb_16(R *rio, R *iio, const R *W, stride ios, INT m, INT dist)
 
45
{
 
46
     DK(KP923879532, +0.923879532511286756128183189396788286822416626);
 
47
     DK(KP707106781, +0.707106781186547524400844362104849039284835938);
 
48
     DK(KP414213562, +0.414213562373095048801688724209698078569671875);
 
49
     INT i;
 
50
     for (i = m - 2; i > 0; i = i - 2, rio = rio + dist, iio = iio - dist, W = W + 30, MAKE_VOLATILE_STRIDE(ios)) {
 
51
          E T3v, T3s, T3u, T3w, T3t;
 
52
          {
 
53
               E T26, T3j, T2z, T36, T11, T1K, T18, T1L, T1C, Tf, T37, T2d, T1m, TE, T3k;
 
54
               E T2C, T1J, Tu, T20, T1F, T1a, TN, T3n, T3e, T1b, TW, T2k, T2h, T2F, T2s;
 
55
               E T3m, T3b;
 
56
               {
 
57
                    E TD, Tw, T2B, T2A;
 
58
                    {
 
59
                         E T24, T3, T2y, T14, T2x, T6, T25, T17, Tb, T2b, Ta, T2a, Tz, Tc, TA;
 
60
                         E TB;
 
61
                         {
 
62
                              E T4, T5, T15, T16;
 
63
                              {
 
64
                                   E T1, T2, T12, T13;
 
65
                                   T1 = rio[0];
 
66
                                   T2 = iio[-WS(ios, 8)];
 
67
                                   T12 = iio[0];
 
68
                                   T13 = rio[WS(ios, 8)];
 
69
                                   T4 = rio[WS(ios, 4)];
 
70
                                   T24 = T1 - T2;
 
71
                                   T3 = T1 + T2;
 
72
                                   T2y = T12 + T13;
 
73
                                   T14 = T12 - T13;
 
74
                                   T5 = iio[-WS(ios, 12)];
 
75
                                   T15 = iio[-WS(ios, 4)];
 
76
                                   T16 = rio[WS(ios, 12)];
 
77
                              }
 
78
                              {
 
79
                                   E T8, T9, Tx, Ty;
 
80
                                   T8 = rio[WS(ios, 2)];
 
81
                                   T2x = T4 - T5;
 
82
                                   T6 = T4 + T5;
 
83
                                   T25 = T15 + T16;
 
84
                                   T17 = T15 - T16;
 
85
                                   T9 = iio[-WS(ios, 10)];
 
86
                                   Tx = iio[-WS(ios, 2)];
 
87
                                   Ty = rio[WS(ios, 10)];
 
88
                                   Tb = iio[-WS(ios, 14)];
 
89
                                   T2b = T8 - T9;
 
90
                                   Ta = T8 + T9;
 
91
                                   T2a = Tx + Ty;
 
92
                                   Tz = Tx - Ty;
 
93
                                   Tc = rio[WS(ios, 6)];
 
94
                                   TA = iio[-WS(ios, 6)];
 
95
                                   TB = rio[WS(ios, 14)];
 
96
                              }
 
97
                         }
 
98
                         {
 
99
                              E T27, T28, TC, Te, Td, T7, T29, T2c;
 
100
                              T26 = T24 - T25;
 
101
                              T3j = T24 + T25;
 
102
                              T27 = Tb - Tc;
 
103
                              Td = Tb + Tc;
 
104
                              T28 = TB + TA;
 
105
                              TC = TA - TB;
 
106
                              T2z = T2x + T2y;
 
107
                              T36 = T2y - T2x;
 
108
                              Te = Ta + Td;
 
109
                              T11 = Td - Ta;
 
110
                              T1K = T14 + T17;
 
111
                              T18 = T14 - T17;
 
112
                              TD = Tz - TC;
 
113
                              T1L = Tz + TC;
 
114
                              Tw = T3 - T6;
 
115
                              T7 = T3 + T6;
 
116
                              T2B = T27 + T28;
 
117
                              T29 = T27 - T28;
 
118
                              T2c = T2a - T2b;
 
119
                              T2A = T2b + T2a;
 
120
                              T1C = T7 - Te;
 
121
                              Tf = T7 + Te;
 
122
                              T37 = T2c + T29;
 
123
                              T2d = T29 - T2c;
 
124
                         }
 
125
                    }
 
126
                    {
 
127
                         E T2f, Ti, T2j, TI, T2i, Tl, T2g, TL, Tq, T2m, Tp, T2q, TR, Tr, TS;
 
128
                         E TT;
 
129
                         {
 
130
                              E Tj, Tk, TJ, TK;
 
131
                              {
 
132
                                   E Tg, Th, TG, TH;
 
133
                                   Tg = rio[WS(ios, 1)];
 
134
                                   T1m = Tw - TD;
 
135
                                   TE = Tw + TD;
 
136
                                   T3k = T2A + T2B;
 
137
                                   T2C = T2A - T2B;
 
138
                                   Th = iio[-WS(ios, 9)];
 
139
                                   TG = iio[-WS(ios, 1)];
 
140
                                   TH = rio[WS(ios, 9)];
 
141
                                   Tj = rio[WS(ios, 5)];
 
142
                                   T2f = Tg - Th;
 
143
                                   Ti = Tg + Th;
 
144
                                   T2j = TG + TH;
 
145
                                   TI = TG - TH;
 
146
                                   Tk = iio[-WS(ios, 13)];
 
147
                                   TJ = iio[-WS(ios, 5)];
 
148
                                   TK = rio[WS(ios, 13)];
 
149
                              }
 
150
                              {
 
151
                                   E Tn, To, TP, TQ;
 
152
                                   Tn = iio[-WS(ios, 15)];
 
153
                                   T2i = Tj - Tk;
 
154
                                   Tl = Tj + Tk;
 
155
                                   T2g = TJ + TK;
 
156
                                   TL = TJ - TK;
 
157
                                   To = rio[WS(ios, 7)];
 
158
                                   TP = iio[-WS(ios, 7)];
 
159
                                   TQ = rio[WS(ios, 15)];
 
160
                                   Tq = rio[WS(ios, 3)];
 
161
                                   T2m = Tn - To;
 
162
                                   Tp = Tn + To;
 
163
                                   T2q = TQ + TP;
 
164
                                   TR = TP - TQ;
 
165
                                   Tr = iio[-WS(ios, 11)];
 
166
                                   TS = iio[-WS(ios, 3)];
 
167
                                   TT = rio[WS(ios, 11)];
 
168
                              }
 
169
                         }
 
170
                         {
 
171
                              E TO, TV, T3c, T2r, T3d, T2o, T39, T3a;
 
172
                              {
 
173
                                   E TF, Tm, T2p, T2n, TU, T1D, TM, Tt, Ts, T1E;
 
174
                                   TF = Ti - Tl;
 
175
                                   Tm = Ti + Tl;
 
176
                                   T2p = Tq - Tr;
 
177
                                   Ts = Tq + Tr;
 
178
                                   T2n = TS + TT;
 
179
                                   TU = TS - TT;
 
180
                                   T1D = TI + TL;
 
181
                                   TM = TI - TL;
 
182
                                   TO = Tp - Ts;
 
183
                                   Tt = Tp + Ts;
 
184
                                   TV = TR - TU;
 
185
                                   T1E = TR + TU;
 
186
                                   T1J = Tt - Tm;
 
187
                                   Tu = Tm + Tt;
 
188
                                   T3c = T2p + T2q;
 
189
                                   T2r = T2p - T2q;
 
190
                                   T20 = T1D + T1E;
 
191
                                   T1F = T1D - T1E;
 
192
                                   T1a = TM - TF;
 
193
                                   TN = TF + TM;
 
194
                                   T3d = T2m + T2n;
 
195
                                   T2o = T2m - T2n;
 
196
                              }
 
197
                              T3n = FMA(KP414213562, T3c, T3d);
 
198
                              T3e = FNMS(KP414213562, T3d, T3c);
 
199
                              T1b = TO + TV;
 
200
                              TW = TO - TV;
 
201
                              T2k = T2i + T2j;
 
202
                              T39 = T2j - T2i;
 
203
                              T3a = T2f + T2g;
 
204
                              T2h = T2f - T2g;
 
205
                              T2F = FNMS(KP414213562, T2o, T2r);
 
206
                              T2s = FMA(KP414213562, T2r, T2o);
 
207
                              T3m = FMA(KP414213562, T39, T3a);
 
208
                              T3b = FNMS(KP414213562, T3a, T39);
 
209
                         }
 
210
                    }
 
211
               }
 
212
               {
 
213
                    E T2E, T2l, T1c, T19, TX, T1z, T1v, T1y, T1x, T1A;
 
214
                    {
 
215
                         E T1M, T1W, T21, T1V, T1Y, T1Z;
 
216
                         rio[0] = Tf + Tu;
 
217
                         T1Z = T1L + T1K;
 
218
                         T1M = T1K - T1L;
 
219
                         T2E = FMA(KP414213562, T2h, T2k);
 
220
                         T2l = FNMS(KP414213562, T2k, T2h);
 
221
                         T1W = Tf - Tu;
 
222
                         T21 = T1Z - T20;
 
223
                         T1V = W[14];
 
224
                         T1Y = W[15];
 
225
                         iio[-WS(ios, 15)] = T20 + T1Z;
 
226
                         {
 
227
                              E T1G, T1T, T1N, T1P, T1B, T1U, T1I, T1H, T1O;
 
228
                              {
 
229
                                   E T1S, T1R, T1X, T22, T1Q;
 
230
                                   T1X = T1V * T1W;
 
231
                                   T22 = T1Y * T1W;
 
232
                                   T1G = T1C + T1F;
 
233
                                   T1Q = T1C - T1F;
 
234
                                   rio[WS(ios, 8)] = FNMS(T1Y, T21, T1X);
 
235
                                   iio[-WS(ios, 7)] = FMA(T1V, T21, T22);
 
236
                                   T1T = T1M - T1J;
 
237
                                   T1N = T1J + T1M;
 
238
                                   T1P = W[6];
 
239
                                   T1S = W[7];
 
240
                                   T1B = W[22];
 
241
                                   T1R = T1P * T1Q;
 
242
                                   T1U = T1S * T1Q;
 
243
                                   T1I = W[23];
 
244
                                   T1H = T1B * T1G;
 
245
                                   rio[WS(ios, 4)] = FNMS(T1S, T1T, T1R);
 
246
                              }
 
247
                              iio[-WS(ios, 11)] = FMA(T1P, T1T, T1U);
 
248
                              T1O = T1I * T1G;
 
249
                              rio[WS(ios, 12)] = FNMS(T1I, T1N, T1H);
 
250
                              {
 
251
                                   E T1r, T1s, T1w, T1o, T1n;
 
252
                                   T1n = T1b - T1a;
 
253
                                   T1c = T1a + T1b;
 
254
                                   T19 = T11 + T18;
 
255
                                   T1r = T18 - T11;
 
256
                                   iio[-WS(ios, 3)] = FMA(T1B, T1N, T1O);
 
257
                                   TX = TN + TW;
 
258
                                   T1s = TN - TW;
 
259
                                   T1w = FNMS(KP707106781, T1n, T1m);
 
260
                                   T1o = FMA(KP707106781, T1n, T1m);
 
261
                                   {
 
262
                                        E T1l, T1t, T1q, T1p, T1u;
 
263
                                        T1l = W[2];
 
264
                                        T1t = FMA(KP707106781, T1s, T1r);
 
265
                                        T1z = FNMS(KP707106781, T1s, T1r);
 
266
                                        T1q = W[3];
 
267
                                        T1p = T1l * T1o;
 
268
                                        T1v = W[18];
 
269
                                        T1y = W[19];
 
270
                                        T1u = T1q * T1o;
 
271
                                        rio[WS(ios, 2)] = FNMS(T1q, T1t, T1p);
 
272
                                        T1x = T1v * T1w;
 
273
                                        T1A = T1y * T1w;
 
274
                                        iio[-WS(ios, 13)] = FMA(T1l, T1t, T1u);
 
275
                                   }
 
276
                              }
 
277
                         }
 
278
                    }
 
279
                    {
 
280
                         E T2V, T2R, T2Q, T2W, T2N, T2M, T2L;
 
281
                         {
 
282
                              E T1g, T1f, T1j, T1h, T1i, TY;
 
283
                              rio[WS(ios, 10)] = FNMS(T1y, T1z, T1x);
 
284
                              iio[-WS(ios, 5)] = FMA(T1v, T1z, T1A);
 
285
                              T1g = FNMS(KP707106781, TX, TE);
 
286
                              TY = FMA(KP707106781, TX, TE);
 
287
                              {
 
288
                                   E Tv, T10, T1d, TZ, T1e;
 
289
                                   Tv = W[26];
 
290
                                   T10 = W[27];
 
291
                                   T1f = W[10];
 
292
                                   T1j = FNMS(KP707106781, T1c, T19);
 
293
                                   T1d = FMA(KP707106781, T1c, T19);
 
294
                                   TZ = Tv * TY;
 
295
                                   T1e = T10 * TY;
 
296
                                   T1h = T1f * T1g;
 
297
                                   T1i = W[11];
 
298
                                   rio[WS(ios, 14)] = FNMS(T10, T1d, TZ);
 
299
                                   iio[-WS(ios, 1)] = FMA(Tv, T1d, T1e);
 
300
                              }
 
301
                              {
 
302
                                   E T2u, T2K, T2H, T23, T2w;
 
303
                                   {
 
304
                                        E T2e, T1k, T2t, T2D, T2G;
 
305
                                        T2e = FMA(KP707106781, T2d, T26);
 
306
                                        T2V = FNMS(KP707106781, T2d, T26);
 
307
                                        rio[WS(ios, 6)] = FNMS(T1i, T1j, T1h);
 
308
                                        T1k = T1i * T1g;
 
309
                                        T2R = T2s - T2l;
 
310
                                        T2t = T2l + T2s;
 
311
                                        T2D = FMA(KP707106781, T2C, T2z);
 
312
                                        T2Q = FNMS(KP707106781, T2C, T2z);
 
313
                                        T2W = T2E - T2F;
 
314
                                        T2G = T2E + T2F;
 
315
                                        iio[-WS(ios, 9)] = FMA(T1f, T1j, T1k);
 
316
                                        T2u = FMA(KP923879532, T2t, T2e);
 
317
                                        T2N = FNMS(KP923879532, T2t, T2e);
 
318
                                        T2K = FNMS(KP923879532, T2G, T2D);
 
319
                                        T2H = FMA(KP923879532, T2G, T2D);
 
320
                                   }
 
321
                                   T23 = W[0];
 
322
                                   T2w = W[1];
 
323
                                   {
 
324
                                        E T2J, T2I, T2v, T2O;
 
325
                                        T2J = W[16];
 
326
                                        T2M = W[17];
 
327
                                        T2I = T23 * T2H;
 
328
                                        T2v = T23 * T2u;
 
329
                                        T2O = T2J * T2N;
 
330
                                        T2L = T2J * T2K;
 
331
                                        iio[-WS(ios, 14)] = FMA(T2w, T2u, T2I);
 
332
                                        rio[WS(ios, 1)] = FNMS(T2w, T2H, T2v);
 
333
                                        rio[WS(ios, 9)] = FNMS(T2M, T2K, T2O);
 
334
                                   }
 
335
                              }
 
336
                         }
 
337
                         iio[-WS(ios, 6)] = FMA(T2M, T2N, T2L);
 
338
                         {
 
339
                              E T33, T30, T32, T34, T31;
 
340
                              {
 
341
                                   E T2P, T2S, T2X, T2U, T2T, T2Z, T2Y;
 
342
                                   T2P = W[24];
 
343
                                   T33 = FNMS(KP923879532, T2R, T2Q);
 
344
                                   T2S = FMA(KP923879532, T2R, T2Q);
 
345
                                   T30 = FNMS(KP923879532, T2W, T2V);
 
346
                                   T2X = FMA(KP923879532, T2W, T2V);
 
347
                                   T2U = W[25];
 
348
                                   T2T = T2P * T2S;
 
349
                                   T2Z = W[8];
 
350
                                   T2Y = T2P * T2X;
 
351
                                   T32 = W[9];
 
352
                                   iio[-WS(ios, 2)] = FMA(T2U, T2X, T2T);
 
353
                                   T34 = T2Z * T33;
 
354
                                   T31 = T2Z * T30;
 
355
                                   rio[WS(ios, 13)] = FNMS(T2U, T2S, T2Y);
 
356
                              }
 
357
                              {
 
358
                                   E T3l, T3f, T38, T3o, T3L, T3I, T3K, T3M, T3J;
 
359
                                   {
 
360
                                        E T3y, T3z, T3D, T3E;
 
361
                                        T3l = FMA(KP707106781, T3k, T3j);
 
362
                                        T3y = FNMS(KP707106781, T3k, T3j);
 
363
                                        iio[-WS(ios, 10)] = FMA(T32, T30, T34);
 
364
                                        rio[WS(ios, 5)] = FNMS(T32, T33, T31);
 
365
                                        T3z = T3b + T3e;
 
366
                                        T3f = T3b - T3e;
 
367
                                        T38 = FMA(KP707106781, T37, T36);
 
368
                                        T3D = FNMS(KP707106781, T37, T36);
 
369
                                        T3E = T3m - T3n;
 
370
                                        T3o = T3m + T3n;
 
371
                                        {
 
372
                                             E T3x, T3A, T3F, T3C, T3B, T3H, T3G;
 
373
                                             T3x = W[4];
 
374
                                             T3L = FMA(KP923879532, T3z, T3y);
 
375
                                             T3A = FNMS(KP923879532, T3z, T3y);
 
376
                                             T3I = FNMS(KP923879532, T3E, T3D);
 
377
                                             T3F = FMA(KP923879532, T3E, T3D);
 
378
                                             T3C = W[5];
 
379
                                             T3B = T3x * T3A;
 
380
                                             T3H = W[20];
 
381
                                             T3G = T3x * T3F;
 
382
                                             T3K = W[21];
 
383
                                             rio[WS(ios, 3)] = FNMS(T3C, T3F, T3B);
 
384
                                             T3M = T3H * T3L;
 
385
                                             T3J = T3H * T3I;
 
386
                                             iio[-WS(ios, 12)] = FMA(T3C, T3A, T3G);
 
387
                                        }
 
388
                                   }
 
389
                                   rio[WS(ios, 11)] = FNMS(T3K, T3I, T3M);
 
390
                                   iio[-WS(ios, 4)] = FMA(T3K, T3L, T3J);
 
391
                                   {
 
392
                                        E T35, T3g, T3p, T3i, T3h, T3r, T3q;
 
393
                                        T35 = W[28];
 
394
                                        T3v = FNMS(KP923879532, T3f, T38);
 
395
                                        T3g = FMA(KP923879532, T3f, T38);
 
396
                                        T3s = FNMS(KP923879532, T3o, T3l);
 
397
                                        T3p = FMA(KP923879532, T3o, T3l);
 
398
                                        T3i = W[29];
 
399
                                        T3h = T35 * T3g;
 
400
                                        T3r = W[12];
 
401
                                        T3q = T35 * T3p;
 
402
                                        T3u = W[13];
 
403
                                        iio[0] = FMA(T3i, T3p, T3h);
 
404
                                        T3w = T3r * T3v;
 
405
                                        T3t = T3r * T3s;
 
406
                                        rio[WS(ios, 15)] = FNMS(T3i, T3g, T3q);
 
407
                                   }
 
408
                              }
 
409
                         }
 
410
                    }
 
411
               }
 
412
          }
 
413
          iio[-WS(ios, 8)] = FMA(T3u, T3s, T3w);
 
414
          rio[WS(ios, 7)] = FNMS(T3u, T3v, T3t);
 
415
     }
 
416
     return W;
 
417
}
 
418
 
 
419
static const tw_instr twinstr[] = {
 
420
     {TW_FULL, 0, 16},
 
421
     {TW_NEXT, 1, 0}
 
422
};
 
423
 
 
424
static const hc2hc_desc desc = { 16, "hb_16", twinstr, &GENUS, {104, 30, 70, 0}, 0, 0, 0 };
 
425
 
 
426
void X(codelet_hb_16) (planner *p) {
 
427
     X(khc2hc_register) (p, hb_16, &desc);
 
428
}
 
429
#else                           /* HAVE_FMA */
 
430
 
 
431
/* Generated by: ../../../genfft/gen_hc2hc -compact -variables 4 -pipeline-latency 4 -sign 1 -n 16 -dif -name hb_16 -include hb.h */
 
432
 
 
433
/*
 
434
 * This function contains 174 FP additions, 84 FP multiplications,
 
435
 * (or, 136 additions, 46 multiplications, 38 fused multiply/add),
 
436
 * 50 stack variables, and 64 memory accesses
 
437
 */
 
438
/*
 
439
 * Generator Id's : 
 
440
 * $Id: algsimp.ml,v 1.9 2006-02-12 23:34:12 athena Exp $
 
441
 * $Id: fft.ml,v 1.4 2006-01-05 03:04:27 stevenj Exp $
 
442
 * $Id: gen_hc2hc.ml,v 1.16 2006-02-12 23:34:12 athena Exp $
 
443
 */
 
444
 
 
445
#include "hb.h"
 
446
 
 
447
static const R *hb_16(R *rio, R *iio, const R *W, stride ios, INT m, INT dist)
 
448
{
 
449
     DK(KP382683432, +0.382683432365089771728459984030398866761344562);
 
450
     DK(KP923879532, +0.923879532511286756128183189396788286822416626);
 
451
     DK(KP707106781, +0.707106781186547524400844362104849039284835938);
 
452
     INT i;
 
453
     for (i = m - 2; i > 0; i = i - 2, rio = rio + dist, iio = iio - dist, W = W + 30, MAKE_VOLATILE_STRIDE(ios)) {
 
454
          E T7, T2K, T30, Tw, T1a, T2e, T2k, T1B, Te, TD, T1C, T13, T2n, T2Z, T2b;
 
455
          E T2L, Tm, T1v, TN, T10, T1W, T2p, T2P, T2W, Tt, T1w, TW, T11, T23, T2q;
 
456
          E T2S, T2X;
 
457
          {
 
458
               E T3, T2c, T16, T2j, T6, T2i, T19, T2d;
 
459
               {
 
460
                    E T1, T2, T14, T15;
 
461
                    T1 = rio[0];
 
462
                    T2 = iio[-WS(ios, 8)];
 
463
                    T3 = T1 + T2;
 
464
                    T2c = T1 - T2;
 
465
                    T14 = iio[0];
 
466
                    T15 = rio[WS(ios, 8)];
 
467
                    T16 = T14 - T15;
 
468
                    T2j = T14 + T15;
 
469
               }
 
470
               {
 
471
                    E T4, T5, T17, T18;
 
472
                    T4 = rio[WS(ios, 4)];
 
473
                    T5 = iio[-WS(ios, 12)];
 
474
                    T6 = T4 + T5;
 
475
                    T2i = T4 - T5;
 
476
                    T17 = iio[-WS(ios, 4)];
 
477
                    T18 = rio[WS(ios, 12)];
 
478
                    T19 = T17 - T18;
 
479
                    T2d = T17 + T18;
 
480
               }
 
481
               T7 = T3 + T6;
 
482
               T2K = T2j - T2i;
 
483
               T30 = T2c + T2d;
 
484
               Tw = T3 - T6;
 
485
               T1a = T16 - T19;
 
486
               T2e = T2c - T2d;
 
487
               T2k = T2i + T2j;
 
488
               T1B = T16 + T19;
 
489
          }
 
490
          {
 
491
               E Ta, T29, Tz, T28, Td, T25, TC, T26;
 
492
               {
 
493
                    E T8, T9, Tx, Ty;
 
494
                    T8 = rio[WS(ios, 2)];
 
495
                    T9 = iio[-WS(ios, 10)];
 
496
                    Ta = T8 + T9;
 
497
                    T29 = T8 - T9;
 
498
                    Tx = iio[-WS(ios, 2)];
 
499
                    Ty = rio[WS(ios, 10)];
 
500
                    Tz = Tx - Ty;
 
501
                    T28 = Tx + Ty;
 
502
               }
 
503
               {
 
504
                    E Tb, Tc, TA, TB;
 
505
                    Tb = iio[-WS(ios, 14)];
 
506
                    Tc = rio[WS(ios, 6)];
 
507
                    Td = Tb + Tc;
 
508
                    T25 = Tb - Tc;
 
509
                    TA = iio[-WS(ios, 6)];
 
510
                    TB = rio[WS(ios, 14)];
 
511
                    TC = TA - TB;
 
512
                    T26 = TB + TA;
 
513
               }
 
514
               Te = Ta + Td;
 
515
               TD = Tz - TC;
 
516
               T1C = Tz + TC;
 
517
               T13 = Td - Ta;
 
518
               {
 
519
                    E T2l, T2m, T27, T2a;
 
520
                    T2l = T29 + T28;
 
521
                    T2m = T25 + T26;
 
522
                    T2n = KP707106781 * (T2l - T2m);
 
523
                    T2Z = KP707106781 * (T2l + T2m);
 
524
                    T27 = T25 - T26;
 
525
                    T2a = T28 - T29;
 
526
                    T2b = KP707106781 * (T27 - T2a);
 
527
                    T2L = KP707106781 * (T2a + T27);
 
528
               }
 
529
          }
 
530
          {
 
531
               E Ti, T1Q, TI, T1U, Tl, T1T, TL, T1R, TF, TM;
 
532
               {
 
533
                    E Tg, Th, TG, TH;
 
534
                    Tg = rio[WS(ios, 1)];
 
535
                    Th = iio[-WS(ios, 9)];
 
536
                    Ti = Tg + Th;
 
537
                    T1Q = Tg - Th;
 
538
                    TG = iio[-WS(ios, 1)];
 
539
                    TH = rio[WS(ios, 9)];
 
540
                    TI = TG - TH;
 
541
                    T1U = TG + TH;
 
542
               }
 
543
               {
 
544
                    E Tj, Tk, TJ, TK;
 
545
                    Tj = rio[WS(ios, 5)];
 
546
                    Tk = iio[-WS(ios, 13)];
 
547
                    Tl = Tj + Tk;
 
548
                    T1T = Tj - Tk;
 
549
                    TJ = iio[-WS(ios, 5)];
 
550
                    TK = rio[WS(ios, 13)];
 
551
                    TL = TJ - TK;
 
552
                    T1R = TJ + TK;
 
553
               }
 
554
               Tm = Ti + Tl;
 
555
               T1v = TI + TL;
 
556
               TF = Ti - Tl;
 
557
               TM = TI - TL;
 
558
               TN = TF + TM;
 
559
               T10 = TM - TF;
 
560
               {
 
561
                    E T1S, T1V, T2N, T2O;
 
562
                    T1S = T1Q - T1R;
 
563
                    T1V = T1T + T1U;
 
564
                    T1W = FNMS(KP382683432, T1V, KP923879532 * T1S);
 
565
                    T2p = FMA(KP923879532, T1V, KP382683432 * T1S);
 
566
                    T2N = T1U - T1T;
 
567
                    T2O = T1Q + T1R;
 
568
                    T2P = FNMS(KP382683432, T2O, KP923879532 * T2N);
 
569
                    T2W = FMA(KP382683432, T2N, KP923879532 * T2O);
 
570
               }
 
571
          }
 
572
          {
 
573
               E Tp, T1X, TR, T21, Ts, T20, TU, T1Y, TO, TV;
 
574
               {
 
575
                    E Tn, To, TP, TQ;
 
576
                    Tn = iio[-WS(ios, 15)];
 
577
                    To = rio[WS(ios, 7)];
 
578
                    Tp = Tn + To;
 
579
                    T1X = Tn - To;
 
580
                    TP = iio[-WS(ios, 7)];
 
581
                    TQ = rio[WS(ios, 15)];
 
582
                    TR = TP - TQ;
 
583
                    T21 = TQ + TP;
 
584
               }
 
585
               {
 
586
                    E Tq, Tr, TS, TT;
 
587
                    Tq = rio[WS(ios, 3)];
 
588
                    Tr = iio[-WS(ios, 11)];
 
589
                    Ts = Tq + Tr;
 
590
                    T20 = Tq - Tr;
 
591
                    TS = iio[-WS(ios, 3)];
 
592
                    TT = rio[WS(ios, 11)];
 
593
                    TU = TS - TT;
 
594
                    T1Y = TS + TT;
 
595
               }
 
596
               Tt = Tp + Ts;
 
597
               T1w = TU + TR;
 
598
               TO = Tp - Ts;
 
599
               TV = TR - TU;
 
600
               TW = TO - TV;
 
601
               T11 = TO + TV;
 
602
               {
 
603
                    E T1Z, T22, T2Q, T2R;
 
604
                    T1Z = T1X - T1Y;
 
605
                    T22 = T20 - T21;
 
606
                    T23 = FMA(KP923879532, T1Z, KP382683432 * T22);
 
607
                    T2q = FNMS(KP382683432, T1Z, KP923879532 * T22);
 
608
                    T2Q = T1X + T1Y;
 
609
                    T2R = T20 + T21;
 
610
                    T2S = FNMS(KP923879532, T2R, KP382683432 * T2Q);
 
611
                    T2X = FMA(KP923879532, T2Q, KP382683432 * T2R);
 
612
               }
 
613
          }
 
614
          {
 
615
               E Tf, Tu, T1K, T1M, T1N, T1O, T1J, T1L;
 
616
               Tf = T7 + Te;
 
617
               Tu = Tm + Tt;
 
618
               T1K = Tf - Tu;
 
619
               T1M = T1C + T1B;
 
620
               T1N = T1v + T1w;
 
621
               T1O = T1M - T1N;
 
622
               rio[0] = Tf + Tu;
 
623
               iio[-WS(ios, 15)] = T1N + T1M;
 
624
               T1J = W[14];
 
625
               T1L = W[15];
 
626
               rio[WS(ios, 8)] = FNMS(T1L, T1O, T1J * T1K);
 
627
               iio[-WS(ios, 7)] = FMA(T1L, T1K, T1J * T1O);
 
628
          }
 
629
          {
 
630
               E T2U, T36, T32, T34;
 
631
               {
 
632
                    E T2M, T2T, T2Y, T31;
 
633
                    T2M = T2K + T2L;
 
634
                    T2T = T2P + T2S;
 
635
                    T2U = T2M + T2T;
 
636
                    T36 = T2M - T2T;
 
637
                    T2Y = T2W + T2X;
 
638
                    T31 = T2Z + T30;
 
639
                    T32 = T2Y + T31;
 
640
                    T34 = T31 - T2Y;
 
641
               }
 
642
               {
 
643
                    E T2J, T2V, T33, T35;
 
644
                    T2J = W[28];
 
645
                    T2V = W[29];
 
646
                    iio[0] = FMA(T2J, T2U, T2V * T32);
 
647
                    rio[WS(ios, 15)] = FNMS(T2V, T2U, T2J * T32);
 
648
                    T33 = W[12];
 
649
                    T35 = W[13];
 
650
                    rio[WS(ios, 7)] = FNMS(T35, T36, T33 * T34);
 
651
                    iio[-WS(ios, 8)] = FMA(T33, T36, T35 * T34);
 
652
               }
 
653
          }
 
654
          {
 
655
               E TY, T1e, T1c, T1g;
 
656
               {
 
657
                    E TE, TX, T12, T1b;
 
658
                    TE = Tw + TD;
 
659
                    TX = KP707106781 * (TN + TW);
 
660
                    TY = TE + TX;
 
661
                    T1e = TE - TX;
 
662
                    T12 = KP707106781 * (T10 + T11);
 
663
                    T1b = T13 + T1a;
 
664
                    T1c = T12 + T1b;
 
665
                    T1g = T1b - T12;
 
666
               }
 
667
               {
 
668
                    E Tv, TZ, T1d, T1f;
 
669
                    Tv = W[26];
 
670
                    TZ = W[27];
 
671
                    rio[WS(ios, 14)] = FNMS(TZ, T1c, Tv * TY);
 
672
                    iio[-WS(ios, 1)] = FMA(TZ, TY, Tv * T1c);
 
673
                    T1d = W[10];
 
674
                    T1f = W[11];
 
675
                    rio[WS(ios, 6)] = FNMS(T1f, T1g, T1d * T1e);
 
676
                    iio[-WS(ios, 9)] = FMA(T1f, T1e, T1d * T1g);
 
677
               }
 
678
          }
 
679
          {
 
680
               E T2g, T2w, T2s, T2u;
 
681
               {
 
682
                    E T24, T2f, T2o, T2r;
 
683
                    T24 = T1W + T23;
 
684
                    T2f = T2b + T2e;
 
685
                    T2g = T24 + T2f;
 
686
                    T2w = T2f - T24;
 
687
                    T2o = T2k + T2n;
 
688
                    T2r = T2p + T2q;
 
689
                    T2s = T2o + T2r;
 
690
                    T2u = T2o - T2r;
 
691
               }
 
692
               {
 
693
                    E T1P, T2h, T2t, T2v;
 
694
                    T1P = W[0];
 
695
                    T2h = W[1];
 
696
                    rio[WS(ios, 1)] = FNMS(T2h, T2s, T1P * T2g);
 
697
                    iio[-WS(ios, 14)] = FMA(T1P, T2s, T2h * T2g);
 
698
                    T2t = W[16];
 
699
                    T2v = W[17];
 
700
                    iio[-WS(ios, 6)] = FMA(T2t, T2u, T2v * T2w);
 
701
                    rio[WS(ios, 9)] = FNMS(T2v, T2u, T2t * T2w);
 
702
               }
 
703
          }
 
704
          {
 
705
               E T1k, T1q, T1o, T1s;
 
706
               {
 
707
                    E T1i, T1j, T1m, T1n;
 
708
                    T1i = Tw - TD;
 
709
                    T1j = KP707106781 * (T11 - T10);
 
710
                    T1k = T1i + T1j;
 
711
                    T1q = T1i - T1j;
 
712
                    T1m = KP707106781 * (TN - TW);
 
713
                    T1n = T1a - T13;
 
714
                    T1o = T1m + T1n;
 
715
                    T1s = T1n - T1m;
 
716
               }
 
717
               {
 
718
                    E T1h, T1l, T1p, T1r;
 
719
                    T1h = W[2];
 
720
                    T1l = W[3];
 
721
                    rio[WS(ios, 2)] = FNMS(T1l, T1o, T1h * T1k);
 
722
                    iio[-WS(ios, 13)] = FMA(T1l, T1k, T1h * T1o);
 
723
                    T1p = W[18];
 
724
                    T1r = W[19];
 
725
                    rio[WS(ios, 10)] = FNMS(T1r, T1s, T1p * T1q);
 
726
                    iio[-WS(ios, 5)] = FMA(T1r, T1q, T1p * T1s);
 
727
               }
 
728
          }
 
729
          {
 
730
               E T2A, T2I, T2E, T2G;
 
731
               {
 
732
                    E T2y, T2z, T2C, T2D;
 
733
                    T2y = T2k - T2n;
 
734
                    T2z = T23 - T1W;
 
735
                    T2A = T2y + T2z;
 
736
                    T2I = T2y - T2z;
 
737
                    T2C = T2p - T2q;
 
738
                    T2D = T2e - T2b;
 
739
                    T2E = T2C + T2D;
 
740
                    T2G = T2D - T2C;
 
741
               }
 
742
               {
 
743
                    E T2x, T2B, T2F, T2H;
 
744
                    T2x = W[24];
 
745
                    T2B = W[25];
 
746
                    iio[-WS(ios, 2)] = FMA(T2x, T2A, T2B * T2E);
 
747
                    rio[WS(ios, 13)] = FNMS(T2B, T2A, T2x * T2E);
 
748
                    T2F = W[8];
 
749
                    T2H = W[9];
 
750
                    rio[WS(ios, 5)] = FNMS(T2H, T2I, T2F * T2G);
 
751
                    iio[-WS(ios, 10)] = FMA(T2F, T2I, T2H * T2G);
 
752
               }
 
753
          }
 
754
          {
 
755
               E T1y, T1G, T1E, T1I;
 
756
               {
 
757
                    E T1u, T1x, T1A, T1D;
 
758
                    T1u = T7 - Te;
 
759
                    T1x = T1v - T1w;
 
760
                    T1y = T1u + T1x;
 
761
                    T1G = T1u - T1x;
 
762
                    T1A = Tt - Tm;
 
763
                    T1D = T1B - T1C;
 
764
                    T1E = T1A + T1D;
 
765
                    T1I = T1D - T1A;
 
766
               }
 
767
               {
 
768
                    E T1t, T1z, T1F, T1H;
 
769
                    T1t = W[22];
 
770
                    T1z = W[23];
 
771
                    rio[WS(ios, 12)] = FNMS(T1z, T1E, T1t * T1y);
 
772
                    iio[-WS(ios, 3)] = FMA(T1z, T1y, T1t * T1E);
 
773
                    T1F = W[6];
 
774
                    T1H = W[7];
 
775
                    rio[WS(ios, 4)] = FNMS(T1H, T1I, T1F * T1G);
 
776
                    iio[-WS(ios, 11)] = FMA(T1H, T1G, T1F * T1I);
 
777
               }
 
778
          }
 
779
          {
 
780
               E T3a, T3i, T3e, T3g;
 
781
               {
 
782
                    E T38, T39, T3c, T3d;
 
783
                    T38 = T2S - T2P;
 
784
                    T39 = T30 - T2Z;
 
785
                    T3a = T38 + T39;
 
786
                    T3i = T39 - T38;
 
787
                    T3c = T2K - T2L;
 
788
                    T3d = T2W - T2X;
 
789
                    T3e = T3c + T3d;
 
790
                    T3g = T3c - T3d;
 
791
               }
 
792
               {
 
793
                    E T37, T3b, T3f, T3h;
 
794
                    T37 = W[4];
 
795
                    T3b = W[5];
 
796
                    rio[WS(ios, 3)] = FNMS(T3b, T3e, T37 * T3a);
 
797
                    iio[-WS(ios, 12)] = FMA(T37, T3e, T3b * T3a);
 
798
                    T3f = W[20];
 
799
                    T3h = W[21];
 
800
                    iio[-WS(ios, 4)] = FMA(T3f, T3g, T3h * T3i);
 
801
                    rio[WS(ios, 11)] = FNMS(T3h, T3g, T3f * T3i);
 
802
               }
 
803
          }
 
804
     }
 
805
     return W;
 
806
}
 
807
 
 
808
static const tw_instr twinstr[] = {
 
809
     {TW_FULL, 0, 16},
 
810
     {TW_NEXT, 1, 0}
 
811
};
 
812
 
 
813
static const hc2hc_desc desc = { 16, "hb_16", twinstr, &GENUS, {136, 46, 38, 0}, 0, 0, 0 };
 
814
 
 
815
void X(codelet_hb_16) (planner *p) {
 
816
     X(khc2hc_register) (p, hb_16, &desc);
 
817
}
 
818
#endif                          /* HAVE_FMA */