annotate src/fftw-3.3.5/dft/simd/common/t3bv_32.c @ 165:7e6e71a29886

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