annotate src/fftw-3.3.8/dft/simd/common/n2sv_32.c @ 168:ceec0dd9ec9c

Replace these with versions built using an older toolset (so as to avoid ABI compatibilities when linking on Ubuntu 14.04 for packaging purposes)
author Chris Cannam <cannam@all-day-breakfast.com>
date Fri, 07 Feb 2020 11:51:13 +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:05:20 EDT 2018 */
cannam@167 23
cannam@167 24 #include "dft/codelet-dft.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_notw.native -fma -simd -compact -variables 4 -pipeline-latency 8 -n 32 -name n2sv_32 -with-ostride 1 -include dft/simd/n2s.h -store-multiple 4 */
cannam@167 29
cannam@167 30 /*
cannam@167 31 * This function contains 372 FP additions, 136 FP multiplications,
cannam@167 32 * (or, 236 additions, 0 multiplications, 136 fused multiply/add),
cannam@167 33 * 138 stack variables, 7 constants, and 144 memory accesses
cannam@167 34 */
cannam@167 35 #include "dft/simd/n2s.h"
cannam@167 36
cannam@167 37 static void n2sv_32(const R *ri, const R *ii, R *ro, R *io, stride is, stride os, INT v, INT ivs, INT ovs)
cannam@167 38 {
cannam@167 39 DVK(KP980785280, +0.980785280403230449126182236134239036973933731);
cannam@167 40 DVK(KP198912367, +0.198912367379658006911597622644676228597850501);
cannam@167 41 DVK(KP831469612, +0.831469612302545237078788377617905756738560812);
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 i;
cannam@167 48 for (i = v; i > 0; i = i - (2 * VL), ri = ri + ((2 * VL) * ivs), ii = ii + ((2 * VL) * ivs), ro = ro + ((2 * VL) * ovs), io = io + ((2 * VL) * ovs), MAKE_VOLATILE_STRIDE(128, is), MAKE_VOLATILE_STRIDE(128, os)) {
cannam@167 49 V T7, T4r, T4Z, T18, T1z, T3t, T3T, T2T, Te, T1f, T50, T4s, T2W, T3u, T1G;
cannam@167 50 V T3U, Tm, T1n, T1O, T2Z, T3y, T3X, T4w, T53, Tt, T1u, T1V, T2Y, T3B, T3W;
cannam@167 51 V T4z, T52, T2t, T3L, T3O, T2K, TR, TY, T5F, T5G, T5H, T5I, T4R, T5k, T2E;
cannam@167 52 V T3M, T4W, T5j, T2N, T3P, T22, T3E, T3H, T2j, TC, TJ, T5A, T5B, T5C, T5D;
cannam@167 53 V T4G, T5h, T2d, T3F, T4L, T5g, T2m, T3I;
cannam@167 54 {
cannam@167 55 V T3, T1x, T14, T2R, T6, T2S, T17, T1y;
cannam@167 56 {
cannam@167 57 V T1, T2, T12, T13;
cannam@167 58 T1 = LD(&(ri[0]), ivs, &(ri[0]));
cannam@167 59 T2 = LD(&(ri[WS(is, 16)]), ivs, &(ri[0]));
cannam@167 60 T3 = VADD(T1, T2);
cannam@167 61 T1x = VSUB(T1, T2);
cannam@167 62 T12 = LD(&(ii[0]), ivs, &(ii[0]));
cannam@167 63 T13 = LD(&(ii[WS(is, 16)]), ivs, &(ii[0]));
cannam@167 64 T14 = VADD(T12, T13);
cannam@167 65 T2R = VSUB(T12, T13);
cannam@167 66 }
cannam@167 67 {
cannam@167 68 V T4, T5, T15, T16;
cannam@167 69 T4 = LD(&(ri[WS(is, 8)]), ivs, &(ri[0]));
cannam@167 70 T5 = LD(&(ri[WS(is, 24)]), ivs, &(ri[0]));
cannam@167 71 T6 = VADD(T4, T5);
cannam@167 72 T2S = VSUB(T4, T5);
cannam@167 73 T15 = LD(&(ii[WS(is, 8)]), ivs, &(ii[0]));
cannam@167 74 T16 = LD(&(ii[WS(is, 24)]), ivs, &(ii[0]));
cannam@167 75 T17 = VADD(T15, T16);
cannam@167 76 T1y = VSUB(T15, T16);
cannam@167 77 }
cannam@167 78 T7 = VADD(T3, T6);
cannam@167 79 T4r = VSUB(T3, T6);
cannam@167 80 T4Z = VSUB(T14, T17);
cannam@167 81 T18 = VADD(T14, T17);
cannam@167 82 T1z = VADD(T1x, T1y);
cannam@167 83 T3t = VSUB(T1x, T1y);
cannam@167 84 T3T = VADD(T2S, T2R);
cannam@167 85 T2T = VSUB(T2R, T2S);
cannam@167 86 }
cannam@167 87 {
cannam@167 88 V Ta, T1A, T1b, T1B, Td, T1D, T1e, T1E;
cannam@167 89 {
cannam@167 90 V T8, T9, T19, T1a;
cannam@167 91 T8 = LD(&(ri[WS(is, 4)]), ivs, &(ri[0]));
cannam@167 92 T9 = LD(&(ri[WS(is, 20)]), ivs, &(ri[0]));
cannam@167 93 Ta = VADD(T8, T9);
cannam@167 94 T1A = VSUB(T8, T9);
cannam@167 95 T19 = LD(&(ii[WS(is, 4)]), ivs, &(ii[0]));
cannam@167 96 T1a = LD(&(ii[WS(is, 20)]), ivs, &(ii[0]));
cannam@167 97 T1b = VADD(T19, T1a);
cannam@167 98 T1B = VSUB(T19, T1a);
cannam@167 99 }
cannam@167 100 {
cannam@167 101 V Tb, Tc, T1c, T1d;
cannam@167 102 Tb = LD(&(ri[WS(is, 28)]), ivs, &(ri[0]));
cannam@167 103 Tc = LD(&(ri[WS(is, 12)]), ivs, &(ri[0]));
cannam@167 104 Td = VADD(Tb, Tc);
cannam@167 105 T1D = VSUB(Tb, Tc);
cannam@167 106 T1c = LD(&(ii[WS(is, 28)]), ivs, &(ii[0]));
cannam@167 107 T1d = LD(&(ii[WS(is, 12)]), ivs, &(ii[0]));
cannam@167 108 T1e = VADD(T1c, T1d);
cannam@167 109 T1E = VSUB(T1c, T1d);
cannam@167 110 }
cannam@167 111 Te = VADD(Ta, Td);
cannam@167 112 T1f = VADD(T1b, T1e);
cannam@167 113 T50 = VSUB(Td, Ta);
cannam@167 114 T4s = VSUB(T1b, T1e);
cannam@167 115 {
cannam@167 116 V T2U, T2V, T1C, T1F;
cannam@167 117 T2U = VSUB(T1B, T1A);
cannam@167 118 T2V = VADD(T1D, T1E);
cannam@167 119 T2W = VADD(T2U, T2V);
cannam@167 120 T3u = VSUB(T2U, T2V);
cannam@167 121 T1C = VADD(T1A, T1B);
cannam@167 122 T1F = VSUB(T1D, T1E);
cannam@167 123 T1G = VADD(T1C, T1F);
cannam@167 124 T3U = VSUB(T1F, T1C);
cannam@167 125 }
cannam@167 126 }
cannam@167 127 {
cannam@167 128 V Ti, T1L, T1j, T1I, Tl, T1J, T1m, T1M, T1K, T1N;
cannam@167 129 {
cannam@167 130 V Tg, Th, T1h, T1i;
cannam@167 131 Tg = LD(&(ri[WS(is, 2)]), ivs, &(ri[0]));
cannam@167 132 Th = LD(&(ri[WS(is, 18)]), ivs, &(ri[0]));
cannam@167 133 Ti = VADD(Tg, Th);
cannam@167 134 T1L = VSUB(Tg, Th);
cannam@167 135 T1h = LD(&(ii[WS(is, 2)]), ivs, &(ii[0]));
cannam@167 136 T1i = LD(&(ii[WS(is, 18)]), ivs, &(ii[0]));
cannam@167 137 T1j = VADD(T1h, T1i);
cannam@167 138 T1I = VSUB(T1h, T1i);
cannam@167 139 }
cannam@167 140 {
cannam@167 141 V Tj, Tk, T1k, T1l;
cannam@167 142 Tj = LD(&(ri[WS(is, 10)]), ivs, &(ri[0]));
cannam@167 143 Tk = LD(&(ri[WS(is, 26)]), ivs, &(ri[0]));
cannam@167 144 Tl = VADD(Tj, Tk);
cannam@167 145 T1J = VSUB(Tj, Tk);
cannam@167 146 T1k = LD(&(ii[WS(is, 10)]), ivs, &(ii[0]));
cannam@167 147 T1l = LD(&(ii[WS(is, 26)]), ivs, &(ii[0]));
cannam@167 148 T1m = VADD(T1k, T1l);
cannam@167 149 T1M = VSUB(T1k, T1l);
cannam@167 150 }
cannam@167 151 Tm = VADD(Ti, Tl);
cannam@167 152 T1n = VADD(T1j, T1m);
cannam@167 153 T1K = VSUB(T1I, T1J);
cannam@167 154 T1N = VADD(T1L, T1M);
cannam@167 155 T1O = VFNMS(LDK(KP414213562), T1N, T1K);
cannam@167 156 T2Z = VFMA(LDK(KP414213562), T1K, T1N);
cannam@167 157 {
cannam@167 158 V T3w, T3x, T4u, T4v;
cannam@167 159 T3w = VADD(T1J, T1I);
cannam@167 160 T3x = VSUB(T1L, T1M);
cannam@167 161 T3y = VFMA(LDK(KP414213562), T3x, T3w);
cannam@167 162 T3X = VFNMS(LDK(KP414213562), T3w, T3x);
cannam@167 163 T4u = VSUB(T1j, T1m);
cannam@167 164 T4v = VSUB(Ti, Tl);
cannam@167 165 T4w = VSUB(T4u, T4v);
cannam@167 166 T53 = VADD(T4v, T4u);
cannam@167 167 }
cannam@167 168 }
cannam@167 169 {
cannam@167 170 V Tp, T1S, T1q, T1P, Ts, T1Q, T1t, T1T, T1R, T1U;
cannam@167 171 {
cannam@167 172 V Tn, To, T1o, T1p;
cannam@167 173 Tn = LD(&(ri[WS(is, 30)]), ivs, &(ri[0]));
cannam@167 174 To = LD(&(ri[WS(is, 14)]), ivs, &(ri[0]));
cannam@167 175 Tp = VADD(Tn, To);
cannam@167 176 T1S = VSUB(Tn, To);
cannam@167 177 T1o = LD(&(ii[WS(is, 30)]), ivs, &(ii[0]));
cannam@167 178 T1p = LD(&(ii[WS(is, 14)]), ivs, &(ii[0]));
cannam@167 179 T1q = VADD(T1o, T1p);
cannam@167 180 T1P = VSUB(T1o, T1p);
cannam@167 181 }
cannam@167 182 {
cannam@167 183 V Tq, Tr, T1r, T1s;
cannam@167 184 Tq = LD(&(ri[WS(is, 6)]), ivs, &(ri[0]));
cannam@167 185 Tr = LD(&(ri[WS(is, 22)]), ivs, &(ri[0]));
cannam@167 186 Ts = VADD(Tq, Tr);
cannam@167 187 T1Q = VSUB(Tq, Tr);
cannam@167 188 T1r = LD(&(ii[WS(is, 6)]), ivs, &(ii[0]));
cannam@167 189 T1s = LD(&(ii[WS(is, 22)]), ivs, &(ii[0]));
cannam@167 190 T1t = VADD(T1r, T1s);
cannam@167 191 T1T = VSUB(T1r, T1s);
cannam@167 192 }
cannam@167 193 Tt = VADD(Tp, Ts);
cannam@167 194 T1u = VADD(T1q, T1t);
cannam@167 195 T1R = VSUB(T1P, T1Q);
cannam@167 196 T1U = VADD(T1S, T1T);
cannam@167 197 T1V = VFMA(LDK(KP414213562), T1U, T1R);
cannam@167 198 T2Y = VFNMS(LDK(KP414213562), T1R, T1U);
cannam@167 199 {
cannam@167 200 V T3z, T3A, T4x, T4y;
cannam@167 201 T3z = VADD(T1Q, T1P);
cannam@167 202 T3A = VSUB(T1S, T1T);
cannam@167 203 T3B = VFNMS(LDK(KP414213562), T3A, T3z);
cannam@167 204 T3W = VFMA(LDK(KP414213562), T3z, T3A);
cannam@167 205 T4x = VSUB(Tp, Ts);
cannam@167 206 T4y = VSUB(T1q, T1t);
cannam@167 207 T4z = VADD(T4x, T4y);
cannam@167 208 T52 = VSUB(T4x, T4y);
cannam@167 209 }
cannam@167 210 }
cannam@167 211 {
cannam@167 212 V TN, T2G, T2r, T4N, TQ, T2s, T2J, T4O, TU, T2x, T2w, T4T, TX, T2z, T2C;
cannam@167 213 V T4U;
cannam@167 214 {
cannam@167 215 V TL, TM, T2p, T2q;
cannam@167 216 TL = LD(&(ri[WS(is, 31)]), ivs, &(ri[WS(is, 1)]));
cannam@167 217 TM = LD(&(ri[WS(is, 15)]), ivs, &(ri[WS(is, 1)]));
cannam@167 218 TN = VADD(TL, TM);
cannam@167 219 T2G = VSUB(TL, TM);
cannam@167 220 T2p = LD(&(ii[WS(is, 31)]), ivs, &(ii[WS(is, 1)]));
cannam@167 221 T2q = LD(&(ii[WS(is, 15)]), ivs, &(ii[WS(is, 1)]));
cannam@167 222 T2r = VSUB(T2p, T2q);
cannam@167 223 T4N = VADD(T2p, T2q);
cannam@167 224 }
cannam@167 225 {
cannam@167 226 V TO, TP, T2H, T2I;
cannam@167 227 TO = LD(&(ri[WS(is, 7)]), ivs, &(ri[WS(is, 1)]));
cannam@167 228 TP = LD(&(ri[WS(is, 23)]), ivs, &(ri[WS(is, 1)]));
cannam@167 229 TQ = VADD(TO, TP);
cannam@167 230 T2s = VSUB(TO, TP);
cannam@167 231 T2H = LD(&(ii[WS(is, 7)]), ivs, &(ii[WS(is, 1)]));
cannam@167 232 T2I = LD(&(ii[WS(is, 23)]), ivs, &(ii[WS(is, 1)]));
cannam@167 233 T2J = VSUB(T2H, T2I);
cannam@167 234 T4O = VADD(T2H, T2I);
cannam@167 235 }
cannam@167 236 {
cannam@167 237 V TS, TT, T2u, T2v;
cannam@167 238 TS = LD(&(ri[WS(is, 3)]), ivs, &(ri[WS(is, 1)]));
cannam@167 239 TT = LD(&(ri[WS(is, 19)]), ivs, &(ri[WS(is, 1)]));
cannam@167 240 TU = VADD(TS, TT);
cannam@167 241 T2x = VSUB(TS, TT);
cannam@167 242 T2u = LD(&(ii[WS(is, 3)]), ivs, &(ii[WS(is, 1)]));
cannam@167 243 T2v = LD(&(ii[WS(is, 19)]), ivs, &(ii[WS(is, 1)]));
cannam@167 244 T2w = VSUB(T2u, T2v);
cannam@167 245 T4T = VADD(T2u, T2v);
cannam@167 246 }
cannam@167 247 {
cannam@167 248 V TV, TW, T2A, T2B;
cannam@167 249 TV = LD(&(ri[WS(is, 27)]), ivs, &(ri[WS(is, 1)]));
cannam@167 250 TW = LD(&(ri[WS(is, 11)]), ivs, &(ri[WS(is, 1)]));
cannam@167 251 TX = VADD(TV, TW);
cannam@167 252 T2z = VSUB(TV, TW);
cannam@167 253 T2A = LD(&(ii[WS(is, 27)]), ivs, &(ii[WS(is, 1)]));
cannam@167 254 T2B = LD(&(ii[WS(is, 11)]), ivs, &(ii[WS(is, 1)]));
cannam@167 255 T2C = VSUB(T2A, T2B);
cannam@167 256 T4U = VADD(T2A, T2B);
cannam@167 257 }
cannam@167 258 T2t = VSUB(T2r, T2s);
cannam@167 259 T3L = VSUB(T2G, T2J);
cannam@167 260 T3O = VADD(T2s, T2r);
cannam@167 261 T2K = VADD(T2G, T2J);
cannam@167 262 TR = VADD(TN, TQ);
cannam@167 263 TY = VADD(TU, TX);
cannam@167 264 T5F = VSUB(TR, TY);
cannam@167 265 {
cannam@167 266 V T4P, T4Q, T2y, T2D;
cannam@167 267 T5G = VADD(T4N, T4O);
cannam@167 268 T5H = VADD(T4T, T4U);
cannam@167 269 T5I = VSUB(T5G, T5H);
cannam@167 270 T4P = VSUB(T4N, T4O);
cannam@167 271 T4Q = VSUB(TX, TU);
cannam@167 272 T4R = VSUB(T4P, T4Q);
cannam@167 273 T5k = VADD(T4Q, T4P);
cannam@167 274 T2y = VSUB(T2w, T2x);
cannam@167 275 T2D = VADD(T2z, T2C);
cannam@167 276 T2E = VADD(T2y, T2D);
cannam@167 277 T3M = VSUB(T2D, T2y);
cannam@167 278 {
cannam@167 279 V T4S, T4V, T2L, T2M;
cannam@167 280 T4S = VSUB(TN, TQ);
cannam@167 281 T4V = VSUB(T4T, T4U);
cannam@167 282 T4W = VSUB(T4S, T4V);
cannam@167 283 T5j = VADD(T4S, T4V);
cannam@167 284 T2L = VADD(T2x, T2w);
cannam@167 285 T2M = VSUB(T2z, T2C);
cannam@167 286 T2N = VADD(T2L, T2M);
cannam@167 287 T3P = VSUB(T2L, T2M);
cannam@167 288 }
cannam@167 289 }
cannam@167 290 }
cannam@167 291 {
cannam@167 292 V Ty, T2f, T20, T4C, TB, T21, T2i, T4D, TF, T26, T25, T4I, TI, T28, T2b;
cannam@167 293 V T4J;
cannam@167 294 {
cannam@167 295 V Tw, Tx, T1Y, T1Z;
cannam@167 296 Tw = LD(&(ri[WS(is, 1)]), ivs, &(ri[WS(is, 1)]));
cannam@167 297 Tx = LD(&(ri[WS(is, 17)]), ivs, &(ri[WS(is, 1)]));
cannam@167 298 Ty = VADD(Tw, Tx);
cannam@167 299 T2f = VSUB(Tw, Tx);
cannam@167 300 T1Y = LD(&(ii[WS(is, 1)]), ivs, &(ii[WS(is, 1)]));
cannam@167 301 T1Z = LD(&(ii[WS(is, 17)]), ivs, &(ii[WS(is, 1)]));
cannam@167 302 T20 = VSUB(T1Y, T1Z);
cannam@167 303 T4C = VADD(T1Y, T1Z);
cannam@167 304 }
cannam@167 305 {
cannam@167 306 V Tz, TA, T2g, T2h;
cannam@167 307 Tz = LD(&(ri[WS(is, 9)]), ivs, &(ri[WS(is, 1)]));
cannam@167 308 TA = LD(&(ri[WS(is, 25)]), ivs, &(ri[WS(is, 1)]));
cannam@167 309 TB = VADD(Tz, TA);
cannam@167 310 T21 = VSUB(Tz, TA);
cannam@167 311 T2g = LD(&(ii[WS(is, 9)]), ivs, &(ii[WS(is, 1)]));
cannam@167 312 T2h = LD(&(ii[WS(is, 25)]), ivs, &(ii[WS(is, 1)]));
cannam@167 313 T2i = VSUB(T2g, T2h);
cannam@167 314 T4D = VADD(T2g, T2h);
cannam@167 315 }
cannam@167 316 {
cannam@167 317 V TD, TE, T23, T24;
cannam@167 318 TD = LD(&(ri[WS(is, 5)]), ivs, &(ri[WS(is, 1)]));
cannam@167 319 TE = LD(&(ri[WS(is, 21)]), ivs, &(ri[WS(is, 1)]));
cannam@167 320 TF = VADD(TD, TE);
cannam@167 321 T26 = VSUB(TD, TE);
cannam@167 322 T23 = LD(&(ii[WS(is, 5)]), ivs, &(ii[WS(is, 1)]));
cannam@167 323 T24 = LD(&(ii[WS(is, 21)]), ivs, &(ii[WS(is, 1)]));
cannam@167 324 T25 = VSUB(T23, T24);
cannam@167 325 T4I = VADD(T23, T24);
cannam@167 326 }
cannam@167 327 {
cannam@167 328 V TG, TH, T29, T2a;
cannam@167 329 TG = LD(&(ri[WS(is, 29)]), ivs, &(ri[WS(is, 1)]));
cannam@167 330 TH = LD(&(ri[WS(is, 13)]), ivs, &(ri[WS(is, 1)]));
cannam@167 331 TI = VADD(TG, TH);
cannam@167 332 T28 = VSUB(TG, TH);
cannam@167 333 T29 = LD(&(ii[WS(is, 29)]), ivs, &(ii[WS(is, 1)]));
cannam@167 334 T2a = LD(&(ii[WS(is, 13)]), ivs, &(ii[WS(is, 1)]));
cannam@167 335 T2b = VSUB(T29, T2a);
cannam@167 336 T4J = VADD(T29, T2a);
cannam@167 337 }
cannam@167 338 T22 = VSUB(T20, T21);
cannam@167 339 T3E = VSUB(T2f, T2i);
cannam@167 340 T3H = VADD(T21, T20);
cannam@167 341 T2j = VADD(T2f, T2i);
cannam@167 342 TC = VADD(Ty, TB);
cannam@167 343 TJ = VADD(TF, TI);
cannam@167 344 T5A = VSUB(TC, TJ);
cannam@167 345 {
cannam@167 346 V T4E, T4F, T27, T2c;
cannam@167 347 T5B = VADD(T4C, T4D);
cannam@167 348 T5C = VADD(T4I, T4J);
cannam@167 349 T5D = VSUB(T5B, T5C);
cannam@167 350 T4E = VSUB(T4C, T4D);
cannam@167 351 T4F = VSUB(TI, TF);
cannam@167 352 T4G = VSUB(T4E, T4F);
cannam@167 353 T5h = VADD(T4F, T4E);
cannam@167 354 T27 = VSUB(T25, T26);
cannam@167 355 T2c = VADD(T28, T2b);
cannam@167 356 T2d = VADD(T27, T2c);
cannam@167 357 T3F = VSUB(T2c, T27);
cannam@167 358 {
cannam@167 359 V T4H, T4K, T2k, T2l;
cannam@167 360 T4H = VSUB(Ty, TB);
cannam@167 361 T4K = VSUB(T4I, T4J);
cannam@167 362 T4L = VSUB(T4H, T4K);
cannam@167 363 T5g = VADD(T4H, T4K);
cannam@167 364 T2k = VADD(T26, T25);
cannam@167 365 T2l = VSUB(T28, T2b);
cannam@167 366 T2m = VADD(T2k, T2l);
cannam@167 367 T3I = VSUB(T2k, T2l);
cannam@167 368 }
cannam@167 369 }
cannam@167 370 }
cannam@167 371 {
cannam@167 372 V T61, T62, T63, T64, T65, T66, T67, T68, T69, T6a, T6b, T6c, T6d, T6e, T6f;
cannam@167 373 V T6g, T6h, T6i, T6j, T6k, T6l, T6m, T6n, T6o, T6p, T6q, T6r, T6s, T6t, T6u;
cannam@167 374 V T6v, T6w;
cannam@167 375 {
cannam@167 376 V T4B, T5b, T5a, T5c, T4Y, T56, T55, T57;
cannam@167 377 {
cannam@167 378 V T4t, T4A, T58, T59;
cannam@167 379 T4t = VSUB(T4r, T4s);
cannam@167 380 T4A = VSUB(T4w, T4z);
cannam@167 381 T4B = VFMA(LDK(KP707106781), T4A, T4t);
cannam@167 382 T5b = VFNMS(LDK(KP707106781), T4A, T4t);
cannam@167 383 T58 = VFMA(LDK(KP414213562), T4R, T4W);
cannam@167 384 T59 = VFNMS(LDK(KP414213562), T4G, T4L);
cannam@167 385 T5a = VSUB(T58, T59);
cannam@167 386 T5c = VADD(T59, T58);
cannam@167 387 }
cannam@167 388 {
cannam@167 389 V T4M, T4X, T51, T54;
cannam@167 390 T4M = VFMA(LDK(KP414213562), T4L, T4G);
cannam@167 391 T4X = VFNMS(LDK(KP414213562), T4W, T4R);
cannam@167 392 T4Y = VSUB(T4M, T4X);
cannam@167 393 T56 = VADD(T4M, T4X);
cannam@167 394 T51 = VSUB(T4Z, T50);
cannam@167 395 T54 = VSUB(T52, T53);
cannam@167 396 T55 = VFNMS(LDK(KP707106781), T54, T51);
cannam@167 397 T57 = VFMA(LDK(KP707106781), T54, T51);
cannam@167 398 }
cannam@167 399 T61 = VFNMS(LDK(KP923879532), T4Y, T4B);
cannam@167 400 STM4(&(ro[22]), T61, ovs, &(ro[0]));
cannam@167 401 T62 = VFNMS(LDK(KP923879532), T5a, T57);
cannam@167 402 STM4(&(io[22]), T62, ovs, &(io[0]));
cannam@167 403 T63 = VFMA(LDK(KP923879532), T4Y, T4B);
cannam@167 404 STM4(&(ro[6]), T63, ovs, &(ro[0]));
cannam@167 405 T64 = VFMA(LDK(KP923879532), T5a, T57);
cannam@167 406 STM4(&(io[6]), T64, ovs, &(io[0]));
cannam@167 407 T65 = VFNMS(LDK(KP923879532), T56, T55);
cannam@167 408 STM4(&(io[14]), T65, ovs, &(io[0]));
cannam@167 409 T66 = VFNMS(LDK(KP923879532), T5c, T5b);
cannam@167 410 STM4(&(ro[14]), T66, ovs, &(ro[0]));
cannam@167 411 T67 = VFMA(LDK(KP923879532), T56, T55);
cannam@167 412 STM4(&(io[30]), T67, ovs, &(io[0]));
cannam@167 413 T68 = VFMA(LDK(KP923879532), T5c, T5b);
cannam@167 414 STM4(&(ro[30]), T68, ovs, &(ro[0]));
cannam@167 415 }
cannam@167 416 {
cannam@167 417 V T5f, T5r, T5u, T5w, T5m, T5q, T5p, T5v;
cannam@167 418 {
cannam@167 419 V T5d, T5e, T5s, T5t;
cannam@167 420 T5d = VADD(T4r, T4s);
cannam@167 421 T5e = VADD(T53, T52);
cannam@167 422 T5f = VFMA(LDK(KP707106781), T5e, T5d);
cannam@167 423 T5r = VFNMS(LDK(KP707106781), T5e, T5d);
cannam@167 424 T5s = VFNMS(LDK(KP414213562), T5g, T5h);
cannam@167 425 T5t = VFMA(LDK(KP414213562), T5j, T5k);
cannam@167 426 T5u = VSUB(T5s, T5t);
cannam@167 427 T5w = VADD(T5s, T5t);
cannam@167 428 }
cannam@167 429 {
cannam@167 430 V T5i, T5l, T5n, T5o;
cannam@167 431 T5i = VFMA(LDK(KP414213562), T5h, T5g);
cannam@167 432 T5l = VFNMS(LDK(KP414213562), T5k, T5j);
cannam@167 433 T5m = VADD(T5i, T5l);
cannam@167 434 T5q = VSUB(T5l, T5i);
cannam@167 435 T5n = VADD(T50, T4Z);
cannam@167 436 T5o = VADD(T4w, T4z);
cannam@167 437 T5p = VFNMS(LDK(KP707106781), T5o, T5n);
cannam@167 438 T5v = VFMA(LDK(KP707106781), T5o, T5n);
cannam@167 439 }
cannam@167 440 T69 = VFNMS(LDK(KP923879532), T5m, T5f);
cannam@167 441 STM4(&(ro[18]), T69, ovs, &(ro[0]));
cannam@167 442 T6a = VFNMS(LDK(KP923879532), T5w, T5v);
cannam@167 443 STM4(&(io[18]), T6a, ovs, &(io[0]));
cannam@167 444 T6b = VFMA(LDK(KP923879532), T5m, T5f);
cannam@167 445 STM4(&(ro[2]), T6b, ovs, &(ro[0]));
cannam@167 446 T6c = VFMA(LDK(KP923879532), T5w, T5v);
cannam@167 447 STM4(&(io[2]), T6c, ovs, &(io[0]));
cannam@167 448 T6d = VFNMS(LDK(KP923879532), T5q, T5p);
cannam@167 449 STM4(&(io[26]), T6d, ovs, &(io[0]));
cannam@167 450 T6e = VFNMS(LDK(KP923879532), T5u, T5r);
cannam@167 451 STM4(&(ro[26]), T6e, ovs, &(ro[0]));
cannam@167 452 T6f = VFMA(LDK(KP923879532), T5q, T5p);
cannam@167 453 STM4(&(io[10]), T6f, ovs, &(io[0]));
cannam@167 454 T6g = VFMA(LDK(KP923879532), T5u, T5r);
cannam@167 455 STM4(&(ro[10]), T6g, ovs, &(ro[0]));
cannam@167 456 }
cannam@167 457 {
cannam@167 458 V T5z, T5P, T5S, T5U, T5K, T5O, T5N, T5T;
cannam@167 459 {
cannam@167 460 V T5x, T5y, T5Q, T5R;
cannam@167 461 T5x = VSUB(T7, Te);
cannam@167 462 T5y = VSUB(T1n, T1u);
cannam@167 463 T5z = VADD(T5x, T5y);
cannam@167 464 T5P = VSUB(T5x, T5y);
cannam@167 465 T5Q = VSUB(T5D, T5A);
cannam@167 466 T5R = VADD(T5F, T5I);
cannam@167 467 T5S = VSUB(T5Q, T5R);
cannam@167 468 T5U = VADD(T5Q, T5R);
cannam@167 469 }
cannam@167 470 {
cannam@167 471 V T5E, T5J, T5L, T5M;
cannam@167 472 T5E = VADD(T5A, T5D);
cannam@167 473 T5J = VSUB(T5F, T5I);
cannam@167 474 T5K = VADD(T5E, T5J);
cannam@167 475 T5O = VSUB(T5J, T5E);
cannam@167 476 T5L = VSUB(T18, T1f);
cannam@167 477 T5M = VSUB(Tt, Tm);
cannam@167 478 T5N = VSUB(T5L, T5M);
cannam@167 479 T5T = VADD(T5M, T5L);
cannam@167 480 }
cannam@167 481 T6h = VFNMS(LDK(KP707106781), T5K, T5z);
cannam@167 482 STM4(&(ro[20]), T6h, ovs, &(ro[0]));
cannam@167 483 T6i = VFNMS(LDK(KP707106781), T5U, T5T);
cannam@167 484 STM4(&(io[20]), T6i, ovs, &(io[0]));
cannam@167 485 T6j = VFMA(LDK(KP707106781), T5K, T5z);
cannam@167 486 STM4(&(ro[4]), T6j, ovs, &(ro[0]));
cannam@167 487 T6k = VFMA(LDK(KP707106781), T5U, T5T);
cannam@167 488 STM4(&(io[4]), T6k, ovs, &(io[0]));
cannam@167 489 T6l = VFNMS(LDK(KP707106781), T5O, T5N);
cannam@167 490 STM4(&(io[28]), T6l, ovs, &(io[0]));
cannam@167 491 T6m = VFNMS(LDK(KP707106781), T5S, T5P);
cannam@167 492 STM4(&(ro[28]), T6m, ovs, &(ro[0]));
cannam@167 493 T6n = VFMA(LDK(KP707106781), T5O, T5N);
cannam@167 494 STM4(&(io[12]), T6n, ovs, &(io[0]));
cannam@167 495 T6o = VFMA(LDK(KP707106781), T5S, T5P);
cannam@167 496 STM4(&(ro[12]), T6o, ovs, &(ro[0]));
cannam@167 497 }
cannam@167 498 {
cannam@167 499 V Tv, T5V, T5Y, T60, T10, T11, T1w, T5Z;
cannam@167 500 {
cannam@167 501 V Tf, Tu, T5W, T5X;
cannam@167 502 Tf = VADD(T7, Te);
cannam@167 503 Tu = VADD(Tm, Tt);
cannam@167 504 Tv = VADD(Tf, Tu);
cannam@167 505 T5V = VSUB(Tf, Tu);
cannam@167 506 T5W = VADD(T5B, T5C);
cannam@167 507 T5X = VADD(T5G, T5H);
cannam@167 508 T5Y = VSUB(T5W, T5X);
cannam@167 509 T60 = VADD(T5W, T5X);
cannam@167 510 }
cannam@167 511 {
cannam@167 512 V TK, TZ, T1g, T1v;
cannam@167 513 TK = VADD(TC, TJ);
cannam@167 514 TZ = VADD(TR, TY);
cannam@167 515 T10 = VADD(TK, TZ);
cannam@167 516 T11 = VSUB(TZ, TK);
cannam@167 517 T1g = VADD(T18, T1f);
cannam@167 518 T1v = VADD(T1n, T1u);
cannam@167 519 T1w = VSUB(T1g, T1v);
cannam@167 520 T5Z = VADD(T1g, T1v);
cannam@167 521 }
cannam@167 522 T6p = VSUB(Tv, T10);
cannam@167 523 STM4(&(ro[16]), T6p, ovs, &(ro[0]));
cannam@167 524 T6q = VSUB(T5Z, T60);
cannam@167 525 STM4(&(io[16]), T6q, ovs, &(io[0]));
cannam@167 526 T6r = VADD(Tv, T10);
cannam@167 527 STM4(&(ro[0]), T6r, ovs, &(ro[0]));
cannam@167 528 T6s = VADD(T5Z, T60);
cannam@167 529 STM4(&(io[0]), T6s, ovs, &(io[0]));
cannam@167 530 T6t = VADD(T11, T1w);
cannam@167 531 STM4(&(io[8]), T6t, ovs, &(io[0]));
cannam@167 532 T6u = VADD(T5V, T5Y);
cannam@167 533 STM4(&(ro[8]), T6u, ovs, &(ro[0]));
cannam@167 534 T6v = VSUB(T1w, T11);
cannam@167 535 STM4(&(io[24]), T6v, ovs, &(io[0]));
cannam@167 536 T6w = VSUB(T5V, T5Y);
cannam@167 537 STM4(&(ro[24]), T6w, ovs, &(ro[0]));
cannam@167 538 }
cannam@167 539 {
cannam@167 540 V T6x, T6y, T6z, T6A, T6B, T6C, T6D, T6E, T6F, T6G, T6H, T6I, T6J, T6K, T6L;
cannam@167 541 V T6M;
cannam@167 542 {
cannam@167 543 V T1X, T37, T31, T33, T2o, T35, T2P, T34;
cannam@167 544 {
cannam@167 545 V T1H, T1W, T2X, T30;
cannam@167 546 T1H = VFNMS(LDK(KP707106781), T1G, T1z);
cannam@167 547 T1W = VSUB(T1O, T1V);
cannam@167 548 T1X = VFMA(LDK(KP923879532), T1W, T1H);
cannam@167 549 T37 = VFNMS(LDK(KP923879532), T1W, T1H);
cannam@167 550 T2X = VFNMS(LDK(KP707106781), T2W, T2T);
cannam@167 551 T30 = VSUB(T2Y, T2Z);
cannam@167 552 T31 = VFNMS(LDK(KP923879532), T30, T2X);
cannam@167 553 T33 = VFMA(LDK(KP923879532), T30, T2X);
cannam@167 554 }
cannam@167 555 {
cannam@167 556 V T2e, T2n, T2F, T2O;
cannam@167 557 T2e = VFNMS(LDK(KP707106781), T2d, T22);
cannam@167 558 T2n = VFNMS(LDK(KP707106781), T2m, T2j);
cannam@167 559 T2o = VFMA(LDK(KP668178637), T2n, T2e);
cannam@167 560 T35 = VFNMS(LDK(KP668178637), T2e, T2n);
cannam@167 561 T2F = VFNMS(LDK(KP707106781), T2E, T2t);
cannam@167 562 T2O = VFNMS(LDK(KP707106781), T2N, T2K);
cannam@167 563 T2P = VFNMS(LDK(KP668178637), T2O, T2F);
cannam@167 564 T34 = VFMA(LDK(KP668178637), T2F, T2O);
cannam@167 565 }
cannam@167 566 {
cannam@167 567 V T2Q, T36, T32, T38;
cannam@167 568 T2Q = VSUB(T2o, T2P);
cannam@167 569 T6x = VFNMS(LDK(KP831469612), T2Q, T1X);
cannam@167 570 STM4(&(ro[21]), T6x, ovs, &(ro[1]));
cannam@167 571 T6y = VFMA(LDK(KP831469612), T2Q, T1X);
cannam@167 572 STM4(&(ro[5]), T6y, ovs, &(ro[1]));
cannam@167 573 T36 = VSUB(T34, T35);
cannam@167 574 T6z = VFNMS(LDK(KP831469612), T36, T33);
cannam@167 575 STM4(&(io[21]), T6z, ovs, &(io[1]));
cannam@167 576 T6A = VFMA(LDK(KP831469612), T36, T33);
cannam@167 577 STM4(&(io[5]), T6A, ovs, &(io[1]));
cannam@167 578 T32 = VADD(T2o, T2P);
cannam@167 579 T6B = VFNMS(LDK(KP831469612), T32, T31);
cannam@167 580 STM4(&(io[13]), T6B, ovs, &(io[1]));
cannam@167 581 T6C = VFMA(LDK(KP831469612), T32, T31);
cannam@167 582 STM4(&(io[29]), T6C, ovs, &(io[1]));
cannam@167 583 T38 = VADD(T35, T34);
cannam@167 584 T6D = VFNMS(LDK(KP831469612), T38, T37);
cannam@167 585 STM4(&(ro[13]), T6D, ovs, &(ro[1]));
cannam@167 586 T6E = VFMA(LDK(KP831469612), T38, T37);
cannam@167 587 STM4(&(ro[29]), T6E, ovs, &(ro[1]));
cannam@167 588 }
cannam@167 589 }
cannam@167 590 {
cannam@167 591 V T3D, T41, T3Z, T45, T3K, T42, T3R, T43;
cannam@167 592 {
cannam@167 593 V T3v, T3C, T3V, T3Y;
cannam@167 594 T3v = VFMA(LDK(KP707106781), T3u, T3t);
cannam@167 595 T3C = VSUB(T3y, T3B);
cannam@167 596 T3D = VFMA(LDK(KP923879532), T3C, T3v);
cannam@167 597 T41 = VFNMS(LDK(KP923879532), T3C, T3v);
cannam@167 598 T3V = VFMA(LDK(KP707106781), T3U, T3T);
cannam@167 599 T3Y = VSUB(T3W, T3X);
cannam@167 600 T3Z = VFNMS(LDK(KP923879532), T3Y, T3V);
cannam@167 601 T45 = VFMA(LDK(KP923879532), T3Y, T3V);
cannam@167 602 }
cannam@167 603 {
cannam@167 604 V T3G, T3J, T3N, T3Q;
cannam@167 605 T3G = VFNMS(LDK(KP707106781), T3F, T3E);
cannam@167 606 T3J = VFNMS(LDK(KP707106781), T3I, T3H);
cannam@167 607 T3K = VFMA(LDK(KP668178637), T3J, T3G);
cannam@167 608 T42 = VFNMS(LDK(KP668178637), T3G, T3J);
cannam@167 609 T3N = VFNMS(LDK(KP707106781), T3M, T3L);
cannam@167 610 T3Q = VFNMS(LDK(KP707106781), T3P, T3O);
cannam@167 611 T3R = VFNMS(LDK(KP668178637), T3Q, T3N);
cannam@167 612 T43 = VFMA(LDK(KP668178637), T3N, T3Q);
cannam@167 613 }
cannam@167 614 {
cannam@167 615 V T3S, T46, T40, T44;
cannam@167 616 T3S = VADD(T3K, T3R);
cannam@167 617 T6F = VFNMS(LDK(KP831469612), T3S, T3D);
cannam@167 618 STM4(&(ro[19]), T6F, ovs, &(ro[1]));
cannam@167 619 T6G = VFMA(LDK(KP831469612), T3S, T3D);
cannam@167 620 STM4(&(ro[3]), T6G, ovs, &(ro[1]));
cannam@167 621 T46 = VADD(T42, T43);
cannam@167 622 T6H = VFNMS(LDK(KP831469612), T46, T45);
cannam@167 623 STM4(&(io[19]), T6H, ovs, &(io[1]));
cannam@167 624 T6I = VFMA(LDK(KP831469612), T46, T45);
cannam@167 625 STM4(&(io[3]), T6I, ovs, &(io[1]));
cannam@167 626 T40 = VSUB(T3R, T3K);
cannam@167 627 T6J = VFNMS(LDK(KP831469612), T40, T3Z);
cannam@167 628 STM4(&(io[27]), T6J, ovs, &(io[1]));
cannam@167 629 T6K = VFMA(LDK(KP831469612), T40, T3Z);
cannam@167 630 STM4(&(io[11]), T6K, ovs, &(io[1]));
cannam@167 631 T44 = VSUB(T42, T43);
cannam@167 632 T6L = VFNMS(LDK(KP831469612), T44, T41);
cannam@167 633 STM4(&(ro[27]), T6L, ovs, &(ro[1]));
cannam@167 634 T6M = VFMA(LDK(KP831469612), T44, T41);
cannam@167 635 STM4(&(ro[11]), T6M, ovs, &(ro[1]));
cannam@167 636 }
cannam@167 637 }
cannam@167 638 {
cannam@167 639 V T49, T4p, T4j, T4l, T4c, T4n, T4f, T4m;
cannam@167 640 {
cannam@167 641 V T47, T48, T4h, T4i;
cannam@167 642 T47 = VFNMS(LDK(KP707106781), T3u, T3t);
cannam@167 643 T48 = VADD(T3X, T3W);
cannam@167 644 T49 = VFNMS(LDK(KP923879532), T48, T47);
cannam@167 645 T4p = VFMA(LDK(KP923879532), T48, T47);
cannam@167 646 T4h = VFNMS(LDK(KP707106781), T3U, T3T);
cannam@167 647 T4i = VADD(T3y, T3B);
cannam@167 648 T4j = VFMA(LDK(KP923879532), T4i, T4h);
cannam@167 649 T4l = VFNMS(LDK(KP923879532), T4i, T4h);
cannam@167 650 }
cannam@167 651 {
cannam@167 652 V T4a, T4b, T4d, T4e;
cannam@167 653 T4a = VFMA(LDK(KP707106781), T3I, T3H);
cannam@167 654 T4b = VFMA(LDK(KP707106781), T3F, T3E);
cannam@167 655 T4c = VFMA(LDK(KP198912367), T4b, T4a);
cannam@167 656 T4n = VFNMS(LDK(KP198912367), T4a, T4b);
cannam@167 657 T4d = VFMA(LDK(KP707106781), T3P, T3O);
cannam@167 658 T4e = VFMA(LDK(KP707106781), T3M, T3L);
cannam@167 659 T4f = VFNMS(LDK(KP198912367), T4e, T4d);
cannam@167 660 T4m = VFMA(LDK(KP198912367), T4d, T4e);
cannam@167 661 }
cannam@167 662 {
cannam@167 663 V T4g, T6N, T6O, T4o, T6P, T6Q;
cannam@167 664 T4g = VSUB(T4c, T4f);
cannam@167 665 T6N = VFNMS(LDK(KP980785280), T4g, T49);
cannam@167 666 STM4(&(ro[23]), T6N, ovs, &(ro[1]));
cannam@167 667 STN4(&(ro[20]), T6h, T6x, T61, T6N, ovs);
cannam@167 668 T6O = VFMA(LDK(KP980785280), T4g, T49);
cannam@167 669 STM4(&(ro[7]), T6O, ovs, &(ro[1]));
cannam@167 670 STN4(&(ro[4]), T6j, T6y, T63, T6O, ovs);
cannam@167 671 T4o = VSUB(T4m, T4n);
cannam@167 672 T6P = VFNMS(LDK(KP980785280), T4o, T4l);
cannam@167 673 STM4(&(io[23]), T6P, ovs, &(io[1]));
cannam@167 674 STN4(&(io[20]), T6i, T6z, T62, T6P, ovs);
cannam@167 675 T6Q = VFMA(LDK(KP980785280), T4o, T4l);
cannam@167 676 STM4(&(io[7]), T6Q, ovs, &(io[1]));
cannam@167 677 STN4(&(io[4]), T6k, T6A, T64, T6Q, ovs);
cannam@167 678 }
cannam@167 679 {
cannam@167 680 V T4k, T6R, T6S, T4q, T6T, T6U;
cannam@167 681 T4k = VADD(T4c, T4f);
cannam@167 682 T6R = VFNMS(LDK(KP980785280), T4k, T4j);
cannam@167 683 STM4(&(io[15]), T6R, ovs, &(io[1]));
cannam@167 684 STN4(&(io[12]), T6n, T6B, T65, T6R, ovs);
cannam@167 685 T6S = VFMA(LDK(KP980785280), T4k, T4j);
cannam@167 686 STM4(&(io[31]), T6S, ovs, &(io[1]));
cannam@167 687 STN4(&(io[28]), T6l, T6C, T67, T6S, ovs);
cannam@167 688 T4q = VADD(T4n, T4m);
cannam@167 689 T6T = VFNMS(LDK(KP980785280), T4q, T4p);
cannam@167 690 STM4(&(ro[15]), T6T, ovs, &(ro[1]));
cannam@167 691 STN4(&(ro[12]), T6o, T6D, T66, T6T, ovs);
cannam@167 692 T6U = VFMA(LDK(KP980785280), T4q, T4p);
cannam@167 693 STM4(&(ro[31]), T6U, ovs, &(ro[1]));
cannam@167 694 STN4(&(ro[28]), T6m, T6E, T68, T6U, ovs);
cannam@167 695 }
cannam@167 696 }
cannam@167 697 {
cannam@167 698 V T3b, T3n, T3l, T3r, T3e, T3o, T3h, T3p;
cannam@167 699 {
cannam@167 700 V T39, T3a, T3j, T3k;
cannam@167 701 T39 = VFMA(LDK(KP707106781), T1G, T1z);
cannam@167 702 T3a = VADD(T2Z, T2Y);
cannam@167 703 T3b = VFMA(LDK(KP923879532), T3a, T39);
cannam@167 704 T3n = VFNMS(LDK(KP923879532), T3a, T39);
cannam@167 705 T3j = VFMA(LDK(KP707106781), T2W, T2T);
cannam@167 706 T3k = VADD(T1O, T1V);
cannam@167 707 T3l = VFNMS(LDK(KP923879532), T3k, T3j);
cannam@167 708 T3r = VFMA(LDK(KP923879532), T3k, T3j);
cannam@167 709 }
cannam@167 710 {
cannam@167 711 V T3c, T3d, T3f, T3g;
cannam@167 712 T3c = VFMA(LDK(KP707106781), T2m, T2j);
cannam@167 713 T3d = VFMA(LDK(KP707106781), T2d, T22);
cannam@167 714 T3e = VFMA(LDK(KP198912367), T3d, T3c);
cannam@167 715 T3o = VFNMS(LDK(KP198912367), T3c, T3d);
cannam@167 716 T3f = VFMA(LDK(KP707106781), T2N, T2K);
cannam@167 717 T3g = VFMA(LDK(KP707106781), T2E, T2t);
cannam@167 718 T3h = VFNMS(LDK(KP198912367), T3g, T3f);
cannam@167 719 T3p = VFMA(LDK(KP198912367), T3f, T3g);
cannam@167 720 }
cannam@167 721 {
cannam@167 722 V T3i, T6V, T6W, T3s, T6X, T6Y;
cannam@167 723 T3i = VADD(T3e, T3h);
cannam@167 724 T6V = VFNMS(LDK(KP980785280), T3i, T3b);
cannam@167 725 STM4(&(ro[17]), T6V, ovs, &(ro[1]));
cannam@167 726 STN4(&(ro[16]), T6p, T6V, T69, T6F, ovs);
cannam@167 727 T6W = VFMA(LDK(KP980785280), T3i, T3b);
cannam@167 728 STM4(&(ro[1]), T6W, ovs, &(ro[1]));
cannam@167 729 STN4(&(ro[0]), T6r, T6W, T6b, T6G, ovs);
cannam@167 730 T3s = VADD(T3o, T3p);
cannam@167 731 T6X = VFNMS(LDK(KP980785280), T3s, T3r);
cannam@167 732 STM4(&(io[17]), T6X, ovs, &(io[1]));
cannam@167 733 STN4(&(io[16]), T6q, T6X, T6a, T6H, ovs);
cannam@167 734 T6Y = VFMA(LDK(KP980785280), T3s, T3r);
cannam@167 735 STM4(&(io[1]), T6Y, ovs, &(io[1]));
cannam@167 736 STN4(&(io[0]), T6s, T6Y, T6c, T6I, ovs);
cannam@167 737 }
cannam@167 738 {
cannam@167 739 V T3m, T6Z, T70, T3q, T71, T72;
cannam@167 740 T3m = VSUB(T3h, T3e);
cannam@167 741 T6Z = VFNMS(LDK(KP980785280), T3m, T3l);
cannam@167 742 STM4(&(io[25]), T6Z, ovs, &(io[1]));
cannam@167 743 STN4(&(io[24]), T6v, T6Z, T6d, T6J, ovs);
cannam@167 744 T70 = VFMA(LDK(KP980785280), T3m, T3l);
cannam@167 745 STM4(&(io[9]), T70, ovs, &(io[1]));
cannam@167 746 STN4(&(io[8]), T6t, T70, T6f, T6K, ovs);
cannam@167 747 T3q = VSUB(T3o, T3p);
cannam@167 748 T71 = VFNMS(LDK(KP980785280), T3q, T3n);
cannam@167 749 STM4(&(ro[25]), T71, ovs, &(ro[1]));
cannam@167 750 STN4(&(ro[24]), T6w, T71, T6e, T6L, ovs);
cannam@167 751 T72 = VFMA(LDK(KP980785280), T3q, T3n);
cannam@167 752 STM4(&(ro[9]), T72, ovs, &(ro[1]));
cannam@167 753 STN4(&(ro[8]), T6u, T72, T6g, T6M, ovs);
cannam@167 754 }
cannam@167 755 }
cannam@167 756 }
cannam@167 757 }
cannam@167 758 }
cannam@167 759 }
cannam@167 760 VLEAVE();
cannam@167 761 }
cannam@167 762
cannam@167 763 static const kdft_desc desc = { 32, XSIMD_STRING("n2sv_32"), {236, 0, 136, 0}, &GENUS, 0, 1, 0, 0 };
cannam@167 764
cannam@167 765 void XSIMD(codelet_n2sv_32) (planner *p) {
cannam@167 766 X(kdft_register) (p, n2sv_32, &desc);
cannam@167 767 }
cannam@167 768
cannam@167 769 #else
cannam@167 770
cannam@167 771 /* Generated by: ../../../genfft/gen_notw.native -simd -compact -variables 4 -pipeline-latency 8 -n 32 -name n2sv_32 -with-ostride 1 -include dft/simd/n2s.h -store-multiple 4 */
cannam@167 772
cannam@167 773 /*
cannam@167 774 * This function contains 372 FP additions, 84 FP multiplications,
cannam@167 775 * (or, 340 additions, 52 multiplications, 32 fused multiply/add),
cannam@167 776 * 130 stack variables, 7 constants, and 144 memory accesses
cannam@167 777 */
cannam@167 778 #include "dft/simd/n2s.h"
cannam@167 779
cannam@167 780 static void n2sv_32(const R *ri, const R *ii, R *ro, R *io, stride is, stride os, INT v, INT ivs, INT ovs)
cannam@167 781 {
cannam@167 782 DVK(KP831469612, +0.831469612302545237078788377617905756738560812);
cannam@167 783 DVK(KP555570233, +0.555570233019602224742830813948532874374937191);
cannam@167 784 DVK(KP195090322, +0.195090322016128267848284868477022240927691618);
cannam@167 785 DVK(KP980785280, +0.980785280403230449126182236134239036973933731);
cannam@167 786 DVK(KP923879532, +0.923879532511286756128183189396788286822416626);
cannam@167 787 DVK(KP382683432, +0.382683432365089771728459984030398866761344562);
cannam@167 788 DVK(KP707106781, +0.707106781186547524400844362104849039284835938);
cannam@167 789 {
cannam@167 790 INT i;
cannam@167 791 for (i = v; i > 0; i = i - (2 * VL), ri = ri + ((2 * VL) * ivs), ii = ii + ((2 * VL) * ivs), ro = ro + ((2 * VL) * ovs), io = io + ((2 * VL) * ovs), MAKE_VOLATILE_STRIDE(128, is), MAKE_VOLATILE_STRIDE(128, os)) {
cannam@167 792 V T7, T4r, T4Z, T18, T1z, T3t, T3T, T2T, Te, T1f, T50, T4s, T2W, T3u, T1G;
cannam@167 793 V T3U, Tm, T1n, T1O, T2Z, T3y, T3X, T4w, T53, Tt, T1u, T1V, T2Y, T3B, T3W;
cannam@167 794 V T4z, T52, T2t, T3L, T3O, T2K, TR, TY, T5F, T5G, T5H, T5I, T4R, T5j, T2E;
cannam@167 795 V T3P, T4W, T5k, T2N, T3M, T22, T3E, T3H, T2j, TC, TJ, T5A, T5B, T5C, T5D;
cannam@167 796 V T4G, T5g, T2d, T3F, T4L, T5h, T2m, T3I;
cannam@167 797 {
cannam@167 798 V T3, T1x, T14, T2S, T6, T2R, T17, T1y;
cannam@167 799 {
cannam@167 800 V T1, T2, T12, T13;
cannam@167 801 T1 = LD(&(ri[0]), ivs, &(ri[0]));
cannam@167 802 T2 = LD(&(ri[WS(is, 16)]), ivs, &(ri[0]));
cannam@167 803 T3 = VADD(T1, T2);
cannam@167 804 T1x = VSUB(T1, T2);
cannam@167 805 T12 = LD(&(ii[0]), ivs, &(ii[0]));
cannam@167 806 T13 = LD(&(ii[WS(is, 16)]), ivs, &(ii[0]));
cannam@167 807 T14 = VADD(T12, T13);
cannam@167 808 T2S = VSUB(T12, T13);
cannam@167 809 }
cannam@167 810 {
cannam@167 811 V T4, T5, T15, T16;
cannam@167 812 T4 = LD(&(ri[WS(is, 8)]), ivs, &(ri[0]));
cannam@167 813 T5 = LD(&(ri[WS(is, 24)]), ivs, &(ri[0]));
cannam@167 814 T6 = VADD(T4, T5);
cannam@167 815 T2R = VSUB(T4, T5);
cannam@167 816 T15 = LD(&(ii[WS(is, 8)]), ivs, &(ii[0]));
cannam@167 817 T16 = LD(&(ii[WS(is, 24)]), ivs, &(ii[0]));
cannam@167 818 T17 = VADD(T15, T16);
cannam@167 819 T1y = VSUB(T15, T16);
cannam@167 820 }
cannam@167 821 T7 = VADD(T3, T6);
cannam@167 822 T4r = VSUB(T3, T6);
cannam@167 823 T4Z = VSUB(T14, T17);
cannam@167 824 T18 = VADD(T14, T17);
cannam@167 825 T1z = VSUB(T1x, T1y);
cannam@167 826 T3t = VADD(T1x, T1y);
cannam@167 827 T3T = VSUB(T2S, T2R);
cannam@167 828 T2T = VADD(T2R, T2S);
cannam@167 829 }
cannam@167 830 {
cannam@167 831 V Ta, T1B, T1b, T1A, Td, T1D, T1e, T1E;
cannam@167 832 {
cannam@167 833 V T8, T9, T19, T1a;
cannam@167 834 T8 = LD(&(ri[WS(is, 4)]), ivs, &(ri[0]));
cannam@167 835 T9 = LD(&(ri[WS(is, 20)]), ivs, &(ri[0]));
cannam@167 836 Ta = VADD(T8, T9);
cannam@167 837 T1B = VSUB(T8, T9);
cannam@167 838 T19 = LD(&(ii[WS(is, 4)]), ivs, &(ii[0]));
cannam@167 839 T1a = LD(&(ii[WS(is, 20)]), ivs, &(ii[0]));
cannam@167 840 T1b = VADD(T19, T1a);
cannam@167 841 T1A = VSUB(T19, T1a);
cannam@167 842 }
cannam@167 843 {
cannam@167 844 V Tb, Tc, T1c, T1d;
cannam@167 845 Tb = LD(&(ri[WS(is, 28)]), ivs, &(ri[0]));
cannam@167 846 Tc = LD(&(ri[WS(is, 12)]), ivs, &(ri[0]));
cannam@167 847 Td = VADD(Tb, Tc);
cannam@167 848 T1D = VSUB(Tb, Tc);
cannam@167 849 T1c = LD(&(ii[WS(is, 28)]), ivs, &(ii[0]));
cannam@167 850 T1d = LD(&(ii[WS(is, 12)]), ivs, &(ii[0]));
cannam@167 851 T1e = VADD(T1c, T1d);
cannam@167 852 T1E = VSUB(T1c, T1d);
cannam@167 853 }
cannam@167 854 Te = VADD(Ta, Td);
cannam@167 855 T1f = VADD(T1b, T1e);
cannam@167 856 T50 = VSUB(Td, Ta);
cannam@167 857 T4s = VSUB(T1b, T1e);
cannam@167 858 {
cannam@167 859 V T2U, T2V, T1C, T1F;
cannam@167 860 T2U = VSUB(T1D, T1E);
cannam@167 861 T2V = VADD(T1B, T1A);
cannam@167 862 T2W = VMUL(LDK(KP707106781), VSUB(T2U, T2V));
cannam@167 863 T3u = VMUL(LDK(KP707106781), VADD(T2V, T2U));
cannam@167 864 T1C = VSUB(T1A, T1B);
cannam@167 865 T1F = VADD(T1D, T1E);
cannam@167 866 T1G = VMUL(LDK(KP707106781), VSUB(T1C, T1F));
cannam@167 867 T3U = VMUL(LDK(KP707106781), VADD(T1C, T1F));
cannam@167 868 }
cannam@167 869 }
cannam@167 870 {
cannam@167 871 V Ti, T1L, T1j, T1J, Tl, T1I, T1m, T1M, T1K, T1N;
cannam@167 872 {
cannam@167 873 V Tg, Th, T1h, T1i;
cannam@167 874 Tg = LD(&(ri[WS(is, 2)]), ivs, &(ri[0]));
cannam@167 875 Th = LD(&(ri[WS(is, 18)]), ivs, &(ri[0]));
cannam@167 876 Ti = VADD(Tg, Th);
cannam@167 877 T1L = VSUB(Tg, Th);
cannam@167 878 T1h = LD(&(ii[WS(is, 2)]), ivs, &(ii[0]));
cannam@167 879 T1i = LD(&(ii[WS(is, 18)]), ivs, &(ii[0]));
cannam@167 880 T1j = VADD(T1h, T1i);
cannam@167 881 T1J = VSUB(T1h, T1i);
cannam@167 882 }
cannam@167 883 {
cannam@167 884 V Tj, Tk, T1k, T1l;
cannam@167 885 Tj = LD(&(ri[WS(is, 10)]), ivs, &(ri[0]));
cannam@167 886 Tk = LD(&(ri[WS(is, 26)]), ivs, &(ri[0]));
cannam@167 887 Tl = VADD(Tj, Tk);
cannam@167 888 T1I = VSUB(Tj, Tk);
cannam@167 889 T1k = LD(&(ii[WS(is, 10)]), ivs, &(ii[0]));
cannam@167 890 T1l = LD(&(ii[WS(is, 26)]), ivs, &(ii[0]));
cannam@167 891 T1m = VADD(T1k, T1l);
cannam@167 892 T1M = VSUB(T1k, T1l);
cannam@167 893 }
cannam@167 894 Tm = VADD(Ti, Tl);
cannam@167 895 T1n = VADD(T1j, T1m);
cannam@167 896 T1K = VADD(T1I, T1J);
cannam@167 897 T1N = VSUB(T1L, T1M);
cannam@167 898 T1O = VFNMS(LDK(KP923879532), T1N, VMUL(LDK(KP382683432), T1K));
cannam@167 899 T2Z = VFMA(LDK(KP923879532), T1K, VMUL(LDK(KP382683432), T1N));
cannam@167 900 {
cannam@167 901 V T3w, T3x, T4u, T4v;
cannam@167 902 T3w = VSUB(T1J, T1I);
cannam@167 903 T3x = VADD(T1L, T1M);
cannam@167 904 T3y = VFNMS(LDK(KP382683432), T3x, VMUL(LDK(KP923879532), T3w));
cannam@167 905 T3X = VFMA(LDK(KP382683432), T3w, VMUL(LDK(KP923879532), T3x));
cannam@167 906 T4u = VSUB(T1j, T1m);
cannam@167 907 T4v = VSUB(Ti, Tl);
cannam@167 908 T4w = VSUB(T4u, T4v);
cannam@167 909 T53 = VADD(T4v, T4u);
cannam@167 910 }
cannam@167 911 }
cannam@167 912 {
cannam@167 913 V Tp, T1S, T1q, T1Q, Ts, T1P, T1t, T1T, T1R, T1U;
cannam@167 914 {
cannam@167 915 V Tn, To, T1o, T1p;
cannam@167 916 Tn = LD(&(ri[WS(is, 30)]), ivs, &(ri[0]));
cannam@167 917 To = LD(&(ri[WS(is, 14)]), ivs, &(ri[0]));
cannam@167 918 Tp = VADD(Tn, To);
cannam@167 919 T1S = VSUB(Tn, To);
cannam@167 920 T1o = LD(&(ii[WS(is, 30)]), ivs, &(ii[0]));
cannam@167 921 T1p = LD(&(ii[WS(is, 14)]), ivs, &(ii[0]));
cannam@167 922 T1q = VADD(T1o, T1p);
cannam@167 923 T1Q = VSUB(T1o, T1p);
cannam@167 924 }
cannam@167 925 {
cannam@167 926 V Tq, Tr, T1r, T1s;
cannam@167 927 Tq = LD(&(ri[WS(is, 6)]), ivs, &(ri[0]));
cannam@167 928 Tr = LD(&(ri[WS(is, 22)]), ivs, &(ri[0]));
cannam@167 929 Ts = VADD(Tq, Tr);
cannam@167 930 T1P = VSUB(Tq, Tr);
cannam@167 931 T1r = LD(&(ii[WS(is, 6)]), ivs, &(ii[0]));
cannam@167 932 T1s = LD(&(ii[WS(is, 22)]), ivs, &(ii[0]));
cannam@167 933 T1t = VADD(T1r, T1s);
cannam@167 934 T1T = VSUB(T1r, T1s);
cannam@167 935 }
cannam@167 936 Tt = VADD(Tp, Ts);
cannam@167 937 T1u = VADD(T1q, T1t);
cannam@167 938 T1R = VADD(T1P, T1Q);
cannam@167 939 T1U = VSUB(T1S, T1T);
cannam@167 940 T1V = VFMA(LDK(KP382683432), T1R, VMUL(LDK(KP923879532), T1U));
cannam@167 941 T2Y = VFNMS(LDK(KP923879532), T1R, VMUL(LDK(KP382683432), T1U));
cannam@167 942 {
cannam@167 943 V T3z, T3A, T4x, T4y;
cannam@167 944 T3z = VSUB(T1Q, T1P);
cannam@167 945 T3A = VADD(T1S, T1T);
cannam@167 946 T3B = VFMA(LDK(KP923879532), T3z, VMUL(LDK(KP382683432), T3A));
cannam@167 947 T3W = VFNMS(LDK(KP382683432), T3z, VMUL(LDK(KP923879532), T3A));
cannam@167 948 T4x = VSUB(Tp, Ts);
cannam@167 949 T4y = VSUB(T1q, T1t);
cannam@167 950 T4z = VADD(T4x, T4y);
cannam@167 951 T52 = VSUB(T4x, T4y);
cannam@167 952 }
cannam@167 953 }
cannam@167 954 {
cannam@167 955 V TN, T2p, T2J, T4S, TQ, T2G, T2s, T4T, TU, T2x, T2w, T4O, TX, T2z, T2C;
cannam@167 956 V T4P;
cannam@167 957 {
cannam@167 958 V TL, TM, T2H, T2I;
cannam@167 959 TL = LD(&(ri[WS(is, 31)]), ivs, &(ri[WS(is, 1)]));
cannam@167 960 TM = LD(&(ri[WS(is, 15)]), ivs, &(ri[WS(is, 1)]));
cannam@167 961 TN = VADD(TL, TM);
cannam@167 962 T2p = VSUB(TL, TM);
cannam@167 963 T2H = LD(&(ii[WS(is, 31)]), ivs, &(ii[WS(is, 1)]));
cannam@167 964 T2I = LD(&(ii[WS(is, 15)]), ivs, &(ii[WS(is, 1)]));
cannam@167 965 T2J = VSUB(T2H, T2I);
cannam@167 966 T4S = VADD(T2H, T2I);
cannam@167 967 }
cannam@167 968 {
cannam@167 969 V TO, TP, T2q, T2r;
cannam@167 970 TO = LD(&(ri[WS(is, 7)]), ivs, &(ri[WS(is, 1)]));
cannam@167 971 TP = LD(&(ri[WS(is, 23)]), ivs, &(ri[WS(is, 1)]));
cannam@167 972 TQ = VADD(TO, TP);
cannam@167 973 T2G = VSUB(TO, TP);
cannam@167 974 T2q = LD(&(ii[WS(is, 7)]), ivs, &(ii[WS(is, 1)]));
cannam@167 975 T2r = LD(&(ii[WS(is, 23)]), ivs, &(ii[WS(is, 1)]));
cannam@167 976 T2s = VSUB(T2q, T2r);
cannam@167 977 T4T = VADD(T2q, T2r);
cannam@167 978 }
cannam@167 979 {
cannam@167 980 V TS, TT, T2u, T2v;
cannam@167 981 TS = LD(&(ri[WS(is, 3)]), ivs, &(ri[WS(is, 1)]));
cannam@167 982 TT = LD(&(ri[WS(is, 19)]), ivs, &(ri[WS(is, 1)]));
cannam@167 983 TU = VADD(TS, TT);
cannam@167 984 T2x = VSUB(TS, TT);
cannam@167 985 T2u = LD(&(ii[WS(is, 3)]), ivs, &(ii[WS(is, 1)]));
cannam@167 986 T2v = LD(&(ii[WS(is, 19)]), ivs, &(ii[WS(is, 1)]));
cannam@167 987 T2w = VSUB(T2u, T2v);
cannam@167 988 T4O = VADD(T2u, T2v);
cannam@167 989 }
cannam@167 990 {
cannam@167 991 V TV, TW, T2A, T2B;
cannam@167 992 TV = LD(&(ri[WS(is, 27)]), ivs, &(ri[WS(is, 1)]));
cannam@167 993 TW = LD(&(ri[WS(is, 11)]), ivs, &(ri[WS(is, 1)]));
cannam@167 994 TX = VADD(TV, TW);
cannam@167 995 T2z = VSUB(TV, TW);
cannam@167 996 T2A = LD(&(ii[WS(is, 27)]), ivs, &(ii[WS(is, 1)]));
cannam@167 997 T2B = LD(&(ii[WS(is, 11)]), ivs, &(ii[WS(is, 1)]));
cannam@167 998 T2C = VSUB(T2A, T2B);
cannam@167 999 T4P = VADD(T2A, T2B);
cannam@167 1000 }
cannam@167 1001 T2t = VSUB(T2p, T2s);
cannam@167 1002 T3L = VADD(T2p, T2s);
cannam@167 1003 T3O = VSUB(T2J, T2G);
cannam@167 1004 T2K = VADD(T2G, T2J);
cannam@167 1005 TR = VADD(TN, TQ);
cannam@167 1006 TY = VADD(TU, TX);
cannam@167 1007 T5F = VSUB(TR, TY);
cannam@167 1008 {
cannam@167 1009 V T4N, T4Q, T2y, T2D;
cannam@167 1010 T5G = VADD(T4S, T4T);
cannam@167 1011 T5H = VADD(T4O, T4P);
cannam@167 1012 T5I = VSUB(T5G, T5H);
cannam@167 1013 T4N = VSUB(TN, TQ);
cannam@167 1014 T4Q = VSUB(T4O, T4P);
cannam@167 1015 T4R = VSUB(T4N, T4Q);
cannam@167 1016 T5j = VADD(T4N, T4Q);
cannam@167 1017 T2y = VSUB(T2w, T2x);
cannam@167 1018 T2D = VADD(T2z, T2C);
cannam@167 1019 T2E = VMUL(LDK(KP707106781), VSUB(T2y, T2D));
cannam@167 1020 T3P = VMUL(LDK(KP707106781), VADD(T2y, T2D));
cannam@167 1021 {
cannam@167 1022 V T4U, T4V, T2L, T2M;
cannam@167 1023 T4U = VSUB(T4S, T4T);
cannam@167 1024 T4V = VSUB(TX, TU);
cannam@167 1025 T4W = VSUB(T4U, T4V);
cannam@167 1026 T5k = VADD(T4V, T4U);
cannam@167 1027 T2L = VSUB(T2z, T2C);
cannam@167 1028 T2M = VADD(T2x, T2w);
cannam@167 1029 T2N = VMUL(LDK(KP707106781), VSUB(T2L, T2M));
cannam@167 1030 T3M = VMUL(LDK(KP707106781), VADD(T2M, T2L));
cannam@167 1031 }
cannam@167 1032 }
cannam@167 1033 }
cannam@167 1034 {
cannam@167 1035 V Ty, T2f, T21, T4C, TB, T1Y, T2i, T4D, TF, T28, T2b, T4I, TI, T23, T26;
cannam@167 1036 V T4J;
cannam@167 1037 {
cannam@167 1038 V Tw, Tx, T1Z, T20;
cannam@167 1039 Tw = LD(&(ri[WS(is, 1)]), ivs, &(ri[WS(is, 1)]));
cannam@167 1040 Tx = LD(&(ri[WS(is, 17)]), ivs, &(ri[WS(is, 1)]));
cannam@167 1041 Ty = VADD(Tw, Tx);
cannam@167 1042 T2f = VSUB(Tw, Tx);
cannam@167 1043 T1Z = LD(&(ii[WS(is, 1)]), ivs, &(ii[WS(is, 1)]));
cannam@167 1044 T20 = LD(&(ii[WS(is, 17)]), ivs, &(ii[WS(is, 1)]));
cannam@167 1045 T21 = VSUB(T1Z, T20);
cannam@167 1046 T4C = VADD(T1Z, T20);
cannam@167 1047 }
cannam@167 1048 {
cannam@167 1049 V Tz, TA, T2g, T2h;
cannam@167 1050 Tz = LD(&(ri[WS(is, 9)]), ivs, &(ri[WS(is, 1)]));
cannam@167 1051 TA = LD(&(ri[WS(is, 25)]), ivs, &(ri[WS(is, 1)]));
cannam@167 1052 TB = VADD(Tz, TA);
cannam@167 1053 T1Y = VSUB(Tz, TA);
cannam@167 1054 T2g = LD(&(ii[WS(is, 9)]), ivs, &(ii[WS(is, 1)]));
cannam@167 1055 T2h = LD(&(ii[WS(is, 25)]), ivs, &(ii[WS(is, 1)]));
cannam@167 1056 T2i = VSUB(T2g, T2h);
cannam@167 1057 T4D = VADD(T2g, T2h);
cannam@167 1058 }
cannam@167 1059 {
cannam@167 1060 V TD, TE, T29, T2a;
cannam@167 1061 TD = LD(&(ri[WS(is, 5)]), ivs, &(ri[WS(is, 1)]));
cannam@167 1062 TE = LD(&(ri[WS(is, 21)]), ivs, &(ri[WS(is, 1)]));
cannam@167 1063 TF = VADD(TD, TE);
cannam@167 1064 T28 = VSUB(TD, TE);
cannam@167 1065 T29 = LD(&(ii[WS(is, 5)]), ivs, &(ii[WS(is, 1)]));
cannam@167 1066 T2a = LD(&(ii[WS(is, 21)]), ivs, &(ii[WS(is, 1)]));
cannam@167 1067 T2b = VSUB(T29, T2a);
cannam@167 1068 T4I = VADD(T29, T2a);
cannam@167 1069 }
cannam@167 1070 {
cannam@167 1071 V TG, TH, T24, T25;
cannam@167 1072 TG = LD(&(ri[WS(is, 29)]), ivs, &(ri[WS(is, 1)]));
cannam@167 1073 TH = LD(&(ri[WS(is, 13)]), ivs, &(ri[WS(is, 1)]));
cannam@167 1074 TI = VADD(TG, TH);
cannam@167 1075 T23 = VSUB(TG, TH);
cannam@167 1076 T24 = LD(&(ii[WS(is, 29)]), ivs, &(ii[WS(is, 1)]));
cannam@167 1077 T25 = LD(&(ii[WS(is, 13)]), ivs, &(ii[WS(is, 1)]));
cannam@167 1078 T26 = VSUB(T24, T25);
cannam@167 1079 T4J = VADD(T24, T25);
cannam@167 1080 }
cannam@167 1081 T22 = VADD(T1Y, T21);
cannam@167 1082 T3E = VADD(T2f, T2i);
cannam@167 1083 T3H = VSUB(T21, T1Y);
cannam@167 1084 T2j = VSUB(T2f, T2i);
cannam@167 1085 TC = VADD(Ty, TB);
cannam@167 1086 TJ = VADD(TF, TI);
cannam@167 1087 T5A = VSUB(TC, TJ);
cannam@167 1088 {
cannam@167 1089 V T4E, T4F, T27, T2c;
cannam@167 1090 T5B = VADD(T4C, T4D);
cannam@167 1091 T5C = VADD(T4I, T4J);
cannam@167 1092 T5D = VSUB(T5B, T5C);
cannam@167 1093 T4E = VSUB(T4C, T4D);
cannam@167 1094 T4F = VSUB(TI, TF);
cannam@167 1095 T4G = VSUB(T4E, T4F);
cannam@167 1096 T5g = VADD(T4F, T4E);
cannam@167 1097 T27 = VSUB(T23, T26);
cannam@167 1098 T2c = VADD(T28, T2b);
cannam@167 1099 T2d = VMUL(LDK(KP707106781), VSUB(T27, T2c));
cannam@167 1100 T3F = VMUL(LDK(KP707106781), VADD(T2c, T27));
cannam@167 1101 {
cannam@167 1102 V T4H, T4K, T2k, T2l;
cannam@167 1103 T4H = VSUB(Ty, TB);
cannam@167 1104 T4K = VSUB(T4I, T4J);
cannam@167 1105 T4L = VSUB(T4H, T4K);
cannam@167 1106 T5h = VADD(T4H, T4K);
cannam@167 1107 T2k = VSUB(T2b, T28);
cannam@167 1108 T2l = VADD(T23, T26);
cannam@167 1109 T2m = VMUL(LDK(KP707106781), VSUB(T2k, T2l));
cannam@167 1110 T3I = VMUL(LDK(KP707106781), VADD(T2k, T2l));
cannam@167 1111 }
cannam@167 1112 }
cannam@167 1113 }
cannam@167 1114 {
cannam@167 1115 V T61, T62, T63, T64, T65, T66, T67, T68, T69, T6a, T6b, T6c, T6d, T6e, T6f;
cannam@167 1116 V T6g, T6h, T6i, T6j, T6k, T6l, T6m, T6n, T6o, T6p, T6q, T6r, T6s, T6t, T6u;
cannam@167 1117 V T6v, T6w;
cannam@167 1118 {
cannam@167 1119 V T4B, T57, T5a, T5c, T4Y, T56, T55, T5b;
cannam@167 1120 {
cannam@167 1121 V T4t, T4A, T58, T59;
cannam@167 1122 T4t = VSUB(T4r, T4s);
cannam@167 1123 T4A = VMUL(LDK(KP707106781), VSUB(T4w, T4z));
cannam@167 1124 T4B = VADD(T4t, T4A);
cannam@167 1125 T57 = VSUB(T4t, T4A);
cannam@167 1126 T58 = VFNMS(LDK(KP923879532), T4L, VMUL(LDK(KP382683432), T4G));
cannam@167 1127 T59 = VFMA(LDK(KP382683432), T4W, VMUL(LDK(KP923879532), T4R));
cannam@167 1128 T5a = VSUB(T58, T59);
cannam@167 1129 T5c = VADD(T58, T59);
cannam@167 1130 }
cannam@167 1131 {
cannam@167 1132 V T4M, T4X, T51, T54;
cannam@167 1133 T4M = VFMA(LDK(KP923879532), T4G, VMUL(LDK(KP382683432), T4L));
cannam@167 1134 T4X = VFNMS(LDK(KP923879532), T4W, VMUL(LDK(KP382683432), T4R));
cannam@167 1135 T4Y = VADD(T4M, T4X);
cannam@167 1136 T56 = VSUB(T4X, T4M);
cannam@167 1137 T51 = VSUB(T4Z, T50);
cannam@167 1138 T54 = VMUL(LDK(KP707106781), VSUB(T52, T53));
cannam@167 1139 T55 = VSUB(T51, T54);
cannam@167 1140 T5b = VADD(T51, T54);
cannam@167 1141 }
cannam@167 1142 T61 = VSUB(T4B, T4Y);
cannam@167 1143 STM4(&(ro[22]), T61, ovs, &(ro[0]));
cannam@167 1144 T62 = VSUB(T5b, T5c);
cannam@167 1145 STM4(&(io[22]), T62, ovs, &(io[0]));
cannam@167 1146 T63 = VADD(T4B, T4Y);
cannam@167 1147 STM4(&(ro[6]), T63, ovs, &(ro[0]));
cannam@167 1148 T64 = VADD(T5b, T5c);
cannam@167 1149 STM4(&(io[6]), T64, ovs, &(io[0]));
cannam@167 1150 T65 = VSUB(T55, T56);
cannam@167 1151 STM4(&(io[30]), T65, ovs, &(io[0]));
cannam@167 1152 T66 = VSUB(T57, T5a);
cannam@167 1153 STM4(&(ro[30]), T66, ovs, &(ro[0]));
cannam@167 1154 T67 = VADD(T55, T56);
cannam@167 1155 STM4(&(io[14]), T67, ovs, &(io[0]));
cannam@167 1156 T68 = VADD(T57, T5a);
cannam@167 1157 STM4(&(ro[14]), T68, ovs, &(ro[0]));
cannam@167 1158 }
cannam@167 1159 {
cannam@167 1160 V T5f, T5r, T5u, T5w, T5m, T5q, T5p, T5v;
cannam@167 1161 {
cannam@167 1162 V T5d, T5e, T5s, T5t;
cannam@167 1163 T5d = VADD(T4r, T4s);
cannam@167 1164 T5e = VMUL(LDK(KP707106781), VADD(T53, T52));
cannam@167 1165 T5f = VADD(T5d, T5e);
cannam@167 1166 T5r = VSUB(T5d, T5e);
cannam@167 1167 T5s = VFNMS(LDK(KP382683432), T5h, VMUL(LDK(KP923879532), T5g));
cannam@167 1168 T5t = VFMA(LDK(KP923879532), T5k, VMUL(LDK(KP382683432), T5j));
cannam@167 1169 T5u = VSUB(T5s, T5t);
cannam@167 1170 T5w = VADD(T5s, T5t);
cannam@167 1171 }
cannam@167 1172 {
cannam@167 1173 V T5i, T5l, T5n, T5o;
cannam@167 1174 T5i = VFMA(LDK(KP382683432), T5g, VMUL(LDK(KP923879532), T5h));
cannam@167 1175 T5l = VFNMS(LDK(KP382683432), T5k, VMUL(LDK(KP923879532), T5j));
cannam@167 1176 T5m = VADD(T5i, T5l);
cannam@167 1177 T5q = VSUB(T5l, T5i);
cannam@167 1178 T5n = VADD(T50, T4Z);
cannam@167 1179 T5o = VMUL(LDK(KP707106781), VADD(T4w, T4z));
cannam@167 1180 T5p = VSUB(T5n, T5o);
cannam@167 1181 T5v = VADD(T5n, T5o);
cannam@167 1182 }
cannam@167 1183 T69 = VSUB(T5f, T5m);
cannam@167 1184 STM4(&(ro[18]), T69, ovs, &(ro[0]));
cannam@167 1185 T6a = VSUB(T5v, T5w);
cannam@167 1186 STM4(&(io[18]), T6a, ovs, &(io[0]));
cannam@167 1187 T6b = VADD(T5f, T5m);
cannam@167 1188 STM4(&(ro[2]), T6b, ovs, &(ro[0]));
cannam@167 1189 T6c = VADD(T5v, T5w);
cannam@167 1190 STM4(&(io[2]), T6c, ovs, &(io[0]));
cannam@167 1191 T6d = VSUB(T5p, T5q);
cannam@167 1192 STM4(&(io[26]), T6d, ovs, &(io[0]));
cannam@167 1193 T6e = VSUB(T5r, T5u);
cannam@167 1194 STM4(&(ro[26]), T6e, ovs, &(ro[0]));
cannam@167 1195 T6f = VADD(T5p, T5q);
cannam@167 1196 STM4(&(io[10]), T6f, ovs, &(io[0]));
cannam@167 1197 T6g = VADD(T5r, T5u);
cannam@167 1198 STM4(&(ro[10]), T6g, ovs, &(ro[0]));
cannam@167 1199 }
cannam@167 1200 {
cannam@167 1201 V T5z, T5P, T5S, T5U, T5K, T5O, T5N, T5T;
cannam@167 1202 {
cannam@167 1203 V T5x, T5y, T5Q, T5R;
cannam@167 1204 T5x = VSUB(T7, Te);
cannam@167 1205 T5y = VSUB(T1n, T1u);
cannam@167 1206 T5z = VADD(T5x, T5y);
cannam@167 1207 T5P = VSUB(T5x, T5y);
cannam@167 1208 T5Q = VSUB(T5D, T5A);
cannam@167 1209 T5R = VADD(T5F, T5I);
cannam@167 1210 T5S = VMUL(LDK(KP707106781), VSUB(T5Q, T5R));
cannam@167 1211 T5U = VMUL(LDK(KP707106781), VADD(T5Q, T5R));
cannam@167 1212 }
cannam@167 1213 {
cannam@167 1214 V T5E, T5J, T5L, T5M;
cannam@167 1215 T5E = VADD(T5A, T5D);
cannam@167 1216 T5J = VSUB(T5F, T5I);
cannam@167 1217 T5K = VMUL(LDK(KP707106781), VADD(T5E, T5J));
cannam@167 1218 T5O = VMUL(LDK(KP707106781), VSUB(T5J, T5E));
cannam@167 1219 T5L = VSUB(T18, T1f);
cannam@167 1220 T5M = VSUB(Tt, Tm);
cannam@167 1221 T5N = VSUB(T5L, T5M);
cannam@167 1222 T5T = VADD(T5M, T5L);
cannam@167 1223 }
cannam@167 1224 T6h = VSUB(T5z, T5K);
cannam@167 1225 STM4(&(ro[20]), T6h, ovs, &(ro[0]));
cannam@167 1226 T6i = VSUB(T5T, T5U);
cannam@167 1227 STM4(&(io[20]), T6i, ovs, &(io[0]));
cannam@167 1228 T6j = VADD(T5z, T5K);
cannam@167 1229 STM4(&(ro[4]), T6j, ovs, &(ro[0]));
cannam@167 1230 T6k = VADD(T5T, T5U);
cannam@167 1231 STM4(&(io[4]), T6k, ovs, &(io[0]));
cannam@167 1232 T6l = VSUB(T5N, T5O);
cannam@167 1233 STM4(&(io[28]), T6l, ovs, &(io[0]));
cannam@167 1234 T6m = VSUB(T5P, T5S);
cannam@167 1235 STM4(&(ro[28]), T6m, ovs, &(ro[0]));
cannam@167 1236 T6n = VADD(T5N, T5O);
cannam@167 1237 STM4(&(io[12]), T6n, ovs, &(io[0]));
cannam@167 1238 T6o = VADD(T5P, T5S);
cannam@167 1239 STM4(&(ro[12]), T6o, ovs, &(ro[0]));
cannam@167 1240 }
cannam@167 1241 {
cannam@167 1242 V Tv, T5V, T5Y, T60, T10, T11, T1w, T5Z;
cannam@167 1243 {
cannam@167 1244 V Tf, Tu, T5W, T5X;
cannam@167 1245 Tf = VADD(T7, Te);
cannam@167 1246 Tu = VADD(Tm, Tt);
cannam@167 1247 Tv = VADD(Tf, Tu);
cannam@167 1248 T5V = VSUB(Tf, Tu);
cannam@167 1249 T5W = VADD(T5B, T5C);
cannam@167 1250 T5X = VADD(T5G, T5H);
cannam@167 1251 T5Y = VSUB(T5W, T5X);
cannam@167 1252 T60 = VADD(T5W, T5X);
cannam@167 1253 }
cannam@167 1254 {
cannam@167 1255 V TK, TZ, T1g, T1v;
cannam@167 1256 TK = VADD(TC, TJ);
cannam@167 1257 TZ = VADD(TR, TY);
cannam@167 1258 T10 = VADD(TK, TZ);
cannam@167 1259 T11 = VSUB(TZ, TK);
cannam@167 1260 T1g = VADD(T18, T1f);
cannam@167 1261 T1v = VADD(T1n, T1u);
cannam@167 1262 T1w = VSUB(T1g, T1v);
cannam@167 1263 T5Z = VADD(T1g, T1v);
cannam@167 1264 }
cannam@167 1265 T6p = VSUB(Tv, T10);
cannam@167 1266 STM4(&(ro[16]), T6p, ovs, &(ro[0]));
cannam@167 1267 T6q = VSUB(T5Z, T60);
cannam@167 1268 STM4(&(io[16]), T6q, ovs, &(io[0]));
cannam@167 1269 T6r = VADD(Tv, T10);
cannam@167 1270 STM4(&(ro[0]), T6r, ovs, &(ro[0]));
cannam@167 1271 T6s = VADD(T5Z, T60);
cannam@167 1272 STM4(&(io[0]), T6s, ovs, &(io[0]));
cannam@167 1273 T6t = VADD(T11, T1w);
cannam@167 1274 STM4(&(io[8]), T6t, ovs, &(io[0]));
cannam@167 1275 T6u = VADD(T5V, T5Y);
cannam@167 1276 STM4(&(ro[8]), T6u, ovs, &(ro[0]));
cannam@167 1277 T6v = VSUB(T1w, T11);
cannam@167 1278 STM4(&(io[24]), T6v, ovs, &(io[0]));
cannam@167 1279 T6w = VSUB(T5V, T5Y);
cannam@167 1280 STM4(&(ro[24]), T6w, ovs, &(ro[0]));
cannam@167 1281 }
cannam@167 1282 {
cannam@167 1283 V T6x, T6y, T6z, T6A, T6B, T6C, T6D, T6E;
cannam@167 1284 {
cannam@167 1285 V T1X, T33, T31, T37, T2o, T34, T2P, T35;
cannam@167 1286 {
cannam@167 1287 V T1H, T1W, T2X, T30;
cannam@167 1288 T1H = VSUB(T1z, T1G);
cannam@167 1289 T1W = VSUB(T1O, T1V);
cannam@167 1290 T1X = VADD(T1H, T1W);
cannam@167 1291 T33 = VSUB(T1H, T1W);
cannam@167 1292 T2X = VSUB(T2T, T2W);
cannam@167 1293 T30 = VSUB(T2Y, T2Z);
cannam@167 1294 T31 = VSUB(T2X, T30);
cannam@167 1295 T37 = VADD(T2X, T30);
cannam@167 1296 }
cannam@167 1297 {
cannam@167 1298 V T2e, T2n, T2F, T2O;
cannam@167 1299 T2e = VSUB(T22, T2d);
cannam@167 1300 T2n = VSUB(T2j, T2m);
cannam@167 1301 T2o = VFMA(LDK(KP980785280), T2e, VMUL(LDK(KP195090322), T2n));
cannam@167 1302 T34 = VFNMS(LDK(KP980785280), T2n, VMUL(LDK(KP195090322), T2e));
cannam@167 1303 T2F = VSUB(T2t, T2E);
cannam@167 1304 T2O = VSUB(T2K, T2N);
cannam@167 1305 T2P = VFNMS(LDK(KP980785280), T2O, VMUL(LDK(KP195090322), T2F));
cannam@167 1306 T35 = VFMA(LDK(KP195090322), T2O, VMUL(LDK(KP980785280), T2F));
cannam@167 1307 }
cannam@167 1308 {
cannam@167 1309 V T2Q, T38, T32, T36;
cannam@167 1310 T2Q = VADD(T2o, T2P);
cannam@167 1311 T6x = VSUB(T1X, T2Q);
cannam@167 1312 STM4(&(ro[23]), T6x, ovs, &(ro[1]));
cannam@167 1313 T6y = VADD(T1X, T2Q);
cannam@167 1314 STM4(&(ro[7]), T6y, ovs, &(ro[1]));
cannam@167 1315 T38 = VADD(T34, T35);
cannam@167 1316 T6z = VSUB(T37, T38);
cannam@167 1317 STM4(&(io[23]), T6z, ovs, &(io[1]));
cannam@167 1318 T6A = VADD(T37, T38);
cannam@167 1319 STM4(&(io[7]), T6A, ovs, &(io[1]));
cannam@167 1320 T32 = VSUB(T2P, T2o);
cannam@167 1321 T6B = VSUB(T31, T32);
cannam@167 1322 STM4(&(io[31]), T6B, ovs, &(io[1]));
cannam@167 1323 T6C = VADD(T31, T32);
cannam@167 1324 STM4(&(io[15]), T6C, ovs, &(io[1]));
cannam@167 1325 T36 = VSUB(T34, T35);
cannam@167 1326 T6D = VSUB(T33, T36);
cannam@167 1327 STM4(&(ro[31]), T6D, ovs, &(ro[1]));
cannam@167 1328 T6E = VADD(T33, T36);
cannam@167 1329 STM4(&(ro[15]), T6E, ovs, &(ro[1]));
cannam@167 1330 }
cannam@167 1331 }
cannam@167 1332 {
cannam@167 1333 V T3D, T41, T3Z, T45, T3K, T42, T3R, T43;
cannam@167 1334 {
cannam@167 1335 V T3v, T3C, T3V, T3Y;
cannam@167 1336 T3v = VSUB(T3t, T3u);
cannam@167 1337 T3C = VSUB(T3y, T3B);
cannam@167 1338 T3D = VADD(T3v, T3C);
cannam@167 1339 T41 = VSUB(T3v, T3C);
cannam@167 1340 T3V = VSUB(T3T, T3U);
cannam@167 1341 T3Y = VSUB(T3W, T3X);
cannam@167 1342 T3Z = VSUB(T3V, T3Y);
cannam@167 1343 T45 = VADD(T3V, T3Y);
cannam@167 1344 }
cannam@167 1345 {
cannam@167 1346 V T3G, T3J, T3N, T3Q;
cannam@167 1347 T3G = VSUB(T3E, T3F);
cannam@167 1348 T3J = VSUB(T3H, T3I);
cannam@167 1349 T3K = VFMA(LDK(KP555570233), T3G, VMUL(LDK(KP831469612), T3J));
cannam@167 1350 T42 = VFNMS(LDK(KP831469612), T3G, VMUL(LDK(KP555570233), T3J));
cannam@167 1351 T3N = VSUB(T3L, T3M);
cannam@167 1352 T3Q = VSUB(T3O, T3P);
cannam@167 1353 T3R = VFNMS(LDK(KP831469612), T3Q, VMUL(LDK(KP555570233), T3N));
cannam@167 1354 T43 = VFMA(LDK(KP831469612), T3N, VMUL(LDK(KP555570233), T3Q));
cannam@167 1355 }
cannam@167 1356 {
cannam@167 1357 V T3S, T6F, T6G, T46, T6H, T6I;
cannam@167 1358 T3S = VADD(T3K, T3R);
cannam@167 1359 T6F = VSUB(T3D, T3S);
cannam@167 1360 STM4(&(ro[21]), T6F, ovs, &(ro[1]));
cannam@167 1361 STN4(&(ro[20]), T6h, T6F, T61, T6x, ovs);
cannam@167 1362 T6G = VADD(T3D, T3S);
cannam@167 1363 STM4(&(ro[5]), T6G, ovs, &(ro[1]));
cannam@167 1364 STN4(&(ro[4]), T6j, T6G, T63, T6y, ovs);
cannam@167 1365 T46 = VADD(T42, T43);
cannam@167 1366 T6H = VSUB(T45, T46);
cannam@167 1367 STM4(&(io[21]), T6H, ovs, &(io[1]));
cannam@167 1368 STN4(&(io[20]), T6i, T6H, T62, T6z, ovs);
cannam@167 1369 T6I = VADD(T45, T46);
cannam@167 1370 STM4(&(io[5]), T6I, ovs, &(io[1]));
cannam@167 1371 STN4(&(io[4]), T6k, T6I, T64, T6A, ovs);
cannam@167 1372 }
cannam@167 1373 {
cannam@167 1374 V T40, T6J, T6K, T44, T6L, T6M;
cannam@167 1375 T40 = VSUB(T3R, T3K);
cannam@167 1376 T6J = VSUB(T3Z, T40);
cannam@167 1377 STM4(&(io[29]), T6J, ovs, &(io[1]));
cannam@167 1378 STN4(&(io[28]), T6l, T6J, T65, T6B, ovs);
cannam@167 1379 T6K = VADD(T3Z, T40);
cannam@167 1380 STM4(&(io[13]), T6K, ovs, &(io[1]));
cannam@167 1381 STN4(&(io[12]), T6n, T6K, T67, T6C, ovs);
cannam@167 1382 T44 = VSUB(T42, T43);
cannam@167 1383 T6L = VSUB(T41, T44);
cannam@167 1384 STM4(&(ro[29]), T6L, ovs, &(ro[1]));
cannam@167 1385 STN4(&(ro[28]), T6m, T6L, T66, T6D, ovs);
cannam@167 1386 T6M = VADD(T41, T44);
cannam@167 1387 STM4(&(ro[13]), T6M, ovs, &(ro[1]));
cannam@167 1388 STN4(&(ro[12]), T6o, T6M, T68, T6E, ovs);
cannam@167 1389 }
cannam@167 1390 }
cannam@167 1391 }
cannam@167 1392 {
cannam@167 1393 V T6N, T6O, T6P, T6Q, T6R, T6S, T6T, T6U;
cannam@167 1394 {
cannam@167 1395 V T49, T4l, T4j, T4p, T4c, T4m, T4f, T4n;
cannam@167 1396 {
cannam@167 1397 V T47, T48, T4h, T4i;
cannam@167 1398 T47 = VADD(T3t, T3u);
cannam@167 1399 T48 = VADD(T3X, T3W);
cannam@167 1400 T49 = VADD(T47, T48);
cannam@167 1401 T4l = VSUB(T47, T48);
cannam@167 1402 T4h = VADD(T3T, T3U);
cannam@167 1403 T4i = VADD(T3y, T3B);
cannam@167 1404 T4j = VSUB(T4h, T4i);
cannam@167 1405 T4p = VADD(T4h, T4i);
cannam@167 1406 }
cannam@167 1407 {
cannam@167 1408 V T4a, T4b, T4d, T4e;
cannam@167 1409 T4a = VADD(T3E, T3F);
cannam@167 1410 T4b = VADD(T3H, T3I);
cannam@167 1411 T4c = VFMA(LDK(KP980785280), T4a, VMUL(LDK(KP195090322), T4b));
cannam@167 1412 T4m = VFNMS(LDK(KP195090322), T4a, VMUL(LDK(KP980785280), T4b));
cannam@167 1413 T4d = VADD(T3L, T3M);
cannam@167 1414 T4e = VADD(T3O, T3P);
cannam@167 1415 T4f = VFNMS(LDK(KP195090322), T4e, VMUL(LDK(KP980785280), T4d));
cannam@167 1416 T4n = VFMA(LDK(KP195090322), T4d, VMUL(LDK(KP980785280), T4e));
cannam@167 1417 }
cannam@167 1418 {
cannam@167 1419 V T4g, T4q, T4k, T4o;
cannam@167 1420 T4g = VADD(T4c, T4f);
cannam@167 1421 T6N = VSUB(T49, T4g);
cannam@167 1422 STM4(&(ro[17]), T6N, ovs, &(ro[1]));
cannam@167 1423 T6O = VADD(T49, T4g);
cannam@167 1424 STM4(&(ro[1]), T6O, ovs, &(ro[1]));
cannam@167 1425 T4q = VADD(T4m, T4n);
cannam@167 1426 T6P = VSUB(T4p, T4q);
cannam@167 1427 STM4(&(io[17]), T6P, ovs, &(io[1]));
cannam@167 1428 T6Q = VADD(T4p, T4q);
cannam@167 1429 STM4(&(io[1]), T6Q, ovs, &(io[1]));
cannam@167 1430 T4k = VSUB(T4f, T4c);
cannam@167 1431 T6R = VSUB(T4j, T4k);
cannam@167 1432 STM4(&(io[25]), T6R, ovs, &(io[1]));
cannam@167 1433 T6S = VADD(T4j, T4k);
cannam@167 1434 STM4(&(io[9]), T6S, ovs, &(io[1]));
cannam@167 1435 T4o = VSUB(T4m, T4n);
cannam@167 1436 T6T = VSUB(T4l, T4o);
cannam@167 1437 STM4(&(ro[25]), T6T, ovs, &(ro[1]));
cannam@167 1438 T6U = VADD(T4l, T4o);
cannam@167 1439 STM4(&(ro[9]), T6U, ovs, &(ro[1]));
cannam@167 1440 }
cannam@167 1441 }
cannam@167 1442 {
cannam@167 1443 V T3b, T3n, T3l, T3r, T3e, T3o, T3h, T3p;
cannam@167 1444 {
cannam@167 1445 V T39, T3a, T3j, T3k;
cannam@167 1446 T39 = VADD(T1z, T1G);
cannam@167 1447 T3a = VADD(T2Z, T2Y);
cannam@167 1448 T3b = VADD(T39, T3a);
cannam@167 1449 T3n = VSUB(T39, T3a);
cannam@167 1450 T3j = VADD(T2T, T2W);
cannam@167 1451 T3k = VADD(T1O, T1V);
cannam@167 1452 T3l = VSUB(T3j, T3k);
cannam@167 1453 T3r = VADD(T3j, T3k);
cannam@167 1454 }
cannam@167 1455 {
cannam@167 1456 V T3c, T3d, T3f, T3g;
cannam@167 1457 T3c = VADD(T22, T2d);
cannam@167 1458 T3d = VADD(T2j, T2m);
cannam@167 1459 T3e = VFMA(LDK(KP555570233), T3c, VMUL(LDK(KP831469612), T3d));
cannam@167 1460 T3o = VFNMS(LDK(KP555570233), T3d, VMUL(LDK(KP831469612), T3c));
cannam@167 1461 T3f = VADD(T2t, T2E);
cannam@167 1462 T3g = VADD(T2K, T2N);
cannam@167 1463 T3h = VFNMS(LDK(KP555570233), T3g, VMUL(LDK(KP831469612), T3f));
cannam@167 1464 T3p = VFMA(LDK(KP831469612), T3g, VMUL(LDK(KP555570233), T3f));
cannam@167 1465 }
cannam@167 1466 {
cannam@167 1467 V T3i, T6V, T6W, T3s, T6X, T6Y;
cannam@167 1468 T3i = VADD(T3e, T3h);
cannam@167 1469 T6V = VSUB(T3b, T3i);
cannam@167 1470 STM4(&(ro[19]), T6V, ovs, &(ro[1]));
cannam@167 1471 STN4(&(ro[16]), T6p, T6N, T69, T6V, ovs);
cannam@167 1472 T6W = VADD(T3b, T3i);
cannam@167 1473 STM4(&(ro[3]), T6W, ovs, &(ro[1]));
cannam@167 1474 STN4(&(ro[0]), T6r, T6O, T6b, T6W, ovs);
cannam@167 1475 T3s = VADD(T3o, T3p);
cannam@167 1476 T6X = VSUB(T3r, T3s);
cannam@167 1477 STM4(&(io[19]), T6X, ovs, &(io[1]));
cannam@167 1478 STN4(&(io[16]), T6q, T6P, T6a, T6X, ovs);
cannam@167 1479 T6Y = VADD(T3r, T3s);
cannam@167 1480 STM4(&(io[3]), T6Y, ovs, &(io[1]));
cannam@167 1481 STN4(&(io[0]), T6s, T6Q, T6c, T6Y, ovs);
cannam@167 1482 }
cannam@167 1483 {
cannam@167 1484 V T3m, T6Z, T70, T3q, T71, T72;
cannam@167 1485 T3m = VSUB(T3h, T3e);
cannam@167 1486 T6Z = VSUB(T3l, T3m);
cannam@167 1487 STM4(&(io[27]), T6Z, ovs, &(io[1]));
cannam@167 1488 STN4(&(io[24]), T6v, T6R, T6d, T6Z, ovs);
cannam@167 1489 T70 = VADD(T3l, T3m);
cannam@167 1490 STM4(&(io[11]), T70, ovs, &(io[1]));
cannam@167 1491 STN4(&(io[8]), T6t, T6S, T6f, T70, ovs);
cannam@167 1492 T3q = VSUB(T3o, T3p);
cannam@167 1493 T71 = VSUB(T3n, T3q);
cannam@167 1494 STM4(&(ro[27]), T71, ovs, &(ro[1]));
cannam@167 1495 STN4(&(ro[24]), T6w, T6T, T6e, T71, ovs);
cannam@167 1496 T72 = VADD(T3n, T3q);
cannam@167 1497 STM4(&(ro[11]), T72, ovs, &(ro[1]));
cannam@167 1498 STN4(&(ro[8]), T6u, T6U, T6g, T72, ovs);
cannam@167 1499 }
cannam@167 1500 }
cannam@167 1501 }
cannam@167 1502 }
cannam@167 1503 }
cannam@167 1504 }
cannam@167 1505 VLEAVE();
cannam@167 1506 }
cannam@167 1507
cannam@167 1508 static const kdft_desc desc = { 32, XSIMD_STRING("n2sv_32"), {340, 52, 32, 0}, &GENUS, 0, 1, 0, 0 };
cannam@167 1509
cannam@167 1510 void XSIMD(codelet_n2sv_32) (planner *p) {
cannam@167 1511 X(kdft_register) (p, n2sv_32, &desc);
cannam@167 1512 }
cannam@167 1513
cannam@167 1514 #endif