21
21
/* This file was automatically generated --- DO NOT EDIT */
22
/* Generated on Sat Jul 5 22:12:01 EDT 2003 */
22
/* Generated on Tue Mar 7 11:55:08 EST 2006 */
24
24
#include "codelet-rdft.h"
26
/* Generated by: /homee/stevenj/cvs/fftw3.0.1/genfft/gen_hc2r -compact -variables 4 -sign 1 -n 8 -name hc2rIII_8 -dft-III -include hc2rIII.h */
28
/* Generated by: ../../../genfft/gen_hc2r -fma -reorder-insns -schedule-for-pipeline -compact -variables 4 -pipeline-latency 4 -sign 1 -n 8 -name hc2rIII_8 -dft-III -include hc2rIII.h */
31
* This function contains 22 FP additions, 12 FP multiplications,
32
* (or, 18 additions, 8 multiplications, 4 fused multiply/add),
33
* 23 stack variables, and 16 memory accesses
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_hc2r.ml,v 1.19 2006-02-12 23:34:12 athena Exp $
44
static void hc2rIII_8(const R *ri, const R *ii, R *O, stride ris, stride iis, stride os, INT v, INT ivs, INT ovs)
46
DK(KP414213562, +0.414213562373095048801688724209698078569671875);
47
DK(KP1_847759065, +1.847759065022573512256366378793576573644833252);
48
DK(KP1_414213562, +1.414213562373095048801688724209698078569671875);
49
DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
51
for (i = v; i > 0; i = i - 1, ri = ri + ivs, ii = ii + ivs, O = O + ovs, MAKE_VOLATILE_STRIDE(ris), MAKE_VOLATILE_STRIDE(iis), MAKE_VOLATILE_STRIDE(os)) {
52
E T4, T7, T3, Tl, Tf, T5, T8, T9, T6, Tc;
79
O[0] = KP2_000000000 * (T3 + T6);
81
O[WS(os, 4)] = KP2_000000000 * (Tl - Tk);
84
O[WS(os, 6)] = KP1_414213562 * (Tm - Tj);
85
O[WS(os, 2)] = KP1_414213562 * (Tj + Tm);
86
O[WS(os, 7)] = -(KP1_847759065 * (FNMS(KP414213562, Th, Ti)));
87
O[WS(os, 3)] = KP1_847759065 * (FMA(KP414213562, Ti, Th));
88
O[WS(os, 5)] = -(KP1_847759065 * (FMA(KP414213562, Tb, Tg)));
89
O[WS(os, 1)] = KP1_847759065 * (FNMS(KP414213562, Tg, Tb));
95
static const khc2r_desc desc = { 8, "hc2rIII_8", {18, 8, 4, 0}, &GENUS, 0, 0, 0, 0, 0 };
97
void X(codelet_hc2rIII_8) (planner *p) {
98
X(khc2rIII_register) (p, hc2rIII_8, &desc);
103
/* Generated by: ../../../genfft/gen_hc2r -compact -variables 4 -pipeline-latency 4 -sign 1 -n 8 -name hc2rIII_8 -dft-III -include hc2rIII.h */
29
106
* This function contains 22 FP additions, 12 FP multiplications,
34
111
* Generator Id's :
35
* $Id: algsimp.ml,v 1.7 2003/03/15 20:29:42 stevenj Exp $
36
* $Id: fft.ml,v 1.2 2003/03/15 20:29:42 stevenj Exp $
37
* $Id: gen_hc2r.ml,v 1.14 2003/04/17 19:25:50 athena Exp $
112
* $Id: algsimp.ml,v 1.9 2006-02-12 23:34:12 athena Exp $
113
* $Id: fft.ml,v 1.4 2006-01-05 03:04:27 stevenj Exp $
114
* $Id: gen_hc2r.ml,v 1.19 2006-02-12 23:34:12 athena Exp $
40
117
#include "hc2rIII.h"
42
static void hc2rIII_8(const R *ri, const R *ii, R *O, stride ris, stride iis, stride os, int v, int ivs, int ovs)
119
static void hc2rIII_8(const R *ri, const R *ii, R *O, stride ris, stride iis, stride os, INT v, INT ivs, INT ovs)
44
121
DK(KP1_414213562, +1.414213562373095048801688724209698078569671875);
45
122
DK(KP765366864, +0.765366864730179543456919968060797733522689125);
46
123
DK(KP1_847759065, +1.847759065022573512256366378793576573644833252);
47
124
DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
49
for (i = v; i > 0; i = i - 1, ri = ri + ivs, ii = ii + ivs, O = O + ovs) {
126
for (i = v; i > 0; i = i - 1, ri = ri + ivs, ii = ii + ivs, O = O + ovs, MAKE_VOLATILE_STRIDE(ris), MAKE_VOLATILE_STRIDE(iis), MAKE_VOLATILE_STRIDE(os)) {
50
127
E T3, T7, Tf, Tl, T6, Tc, Ta, Tk, Tb, Tg;