annotate src/fftw-3.3.8/rdft/simd/common/hc2cbdftv_32.c @ 169:223a55898ab9 tip default

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