sv-dependency-builds: src/fftw-3.3.8/rdft/scalar/r2cf/r2cf

annotate src/fftw-3.3.8/rdft/scalar/r2cf/r2cf_128.c @ 83:ae30d91d2ffe

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
date	Fri, 07 Feb 2020 11:51:13 +0000
parents	d0c2a83c1364
children

rev	line source
Chris@82	1 /*
Chris@82	2 * Copyright (c) 2003, 2007-14 Matteo Frigo
Chris@82	3 * Copyright (c) 2003, 2007-14 Massachusetts Institute of Technology
Chris@82	4 *
Chris@82	5 * This program is free software; you can redistribute it and/or modify
Chris@82	6 * it under the terms of the GNU General Public License as published by
Chris@82	7 * the Free Software Foundation; either version 2 of the License, or
Chris@82	8 * (at your option) any later version.
Chris@82	9 *
Chris@82	10 * This program is distributed in the hope that it will be useful,
Chris@82	11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
Chris@82	12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
Chris@82	13 * GNU General Public License for more details.
Chris@82	14 *
Chris@82	15 * You should have received a copy of the GNU General Public License
Chris@82	16 * along with this program; if not, write to the Free Software
Chris@82	17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
Chris@82	18 *
Chris@82	19 */
Chris@82	20
Chris@82	21 /* This file was automatically generated --- DO NOT EDIT */
Chris@82	22 /* Generated on Thu May 24 08:06:27 EDT 2018 */
Chris@82	23
Chris@82	24 #include "rdft/codelet-rdft.h"
Chris@82	25
Chris@82	26 #if defined(ARCH_PREFERS_FMA) \|\| defined(ISA_EXTENSION_PREFERS_FMA)
Chris@82	27
Chris@82	28 /* Generated by: ../../../genfft/gen_r2cf.native -fma -compact -variables 4 -pipeline-latency 4 -n 128 -name r2cf_128 -include rdft/scalar/r2cf.h */
Chris@82	29
Chris@82	30 /*
Chris@82	31 * This function contains 956 FP additions, 516 FP multiplications,
Chris@82	32 * (or, 440 additions, 0 multiplications, 516 fused multiply/add),
Chris@82	33 * 186 stack variables, 31 constants, and 256 memory accesses
Chris@82	34 */
Chris@82	35 #include "rdft/scalar/r2cf.h"
Chris@82	36
Chris@82	37 static void r2cf_128(R R0, R R1, R Cr, R Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
Chris@82	38 {
Chris@82	39 DK(KP803207531, +0.803207531480644909806676512963141923879569427);
Chris@82	40 DK(KP989176509, +0.989176509964780973451673738016243063983689533);
Chris@82	41 DK(KP741650546, +0.741650546272035369581266691172079863842265220);
Chris@82	42 DK(KP148335987, +0.148335987538347428753676511486911367000625355);
Chris@82	43 DK(KP740951125, +0.740951125354959091175616897495162729728955309);
Chris@82	44 DK(KP998795456, +0.998795456205172392714771604759100694443203615);
Chris@82	45 DK(KP906347169, +0.906347169019147157946142717268914412664134293);
Chris@82	46 DK(KP049126849, +0.049126849769467254105343321271313617079695752);
Chris@82	47 DK(KP970031253, +0.970031253194543992603984207286100251456865962);
Chris@82	48 DK(KP857728610, +0.857728610000272069902269984284770137042490799);
Chris@82	49 DK(KP250486960, +0.250486960191305461595702160124721208578685568);
Chris@82	50 DK(KP599376933, +0.599376933681923766271389869014404232837890546);
Chris@82	51 DK(KP941544065, +0.941544065183020778412509402599502357185589796);
Chris@82	52 DK(KP903989293, +0.903989293123443331586200297230537048710132025);
Chris@82	53 DK(KP357805721, +0.357805721314524104672487743774474392487532769);
Chris@82	54 DK(KP472964775, +0.472964775891319928124438237972992463904131113);
Chris@82	55 DK(KP773010453, +0.773010453362736960810906609758469800971041293);
Chris@82	56 DK(KP995184726, +0.995184726672196886244836953109479921575474869);
Chris@82	57 DK(KP098491403, +0.098491403357164253077197521291327432293052451);
Chris@82	58 DK(KP820678790, +0.820678790828660330972281985331011598767386482);
Chris@82	59 DK(KP956940335, +0.956940335732208864935797886980269969482849206);
Chris@82	60 DK(KP881921264, +0.881921264348355029712756863660388349508442621);
Chris@82	61 DK(KP534511135, +0.534511135950791641089685961295362908582039528);
Chris@82	62 DK(KP303346683, +0.303346683607342391675883946941299872384187453);
Chris@82	63 DK(KP831469612, +0.831469612302545237078788377617905756738560812);
Chris@82	64 DK(KP980785280, +0.980785280403230449126182236134239036973933731);
Chris@82	65 DK(KP923879532, +0.923879532511286756128183189396788286822416626);
Chris@82	66 DK(KP668178637, +0.668178637919298919997757686523080761552472251);
Chris@82	67 DK(KP198912367, +0.198912367379658006911597622644676228597850501);
Chris@82	68 DK(KP414213562, +0.414213562373095048801688724209698078569671875);
Chris@82	69 DK(KP707106781, +0.707106781186547524400844362104849039284835938);
Chris@82	70 {
Chris@82	71 INT i;
Chris@82	72 for (i = v; i > 0; i = i - 1, R0 = R0 + ivs, R1 = R1 + ivs, Cr = Cr + ovs, Ci = Ci + ovs, MAKE_VOLATILE_STRIDE(512, rs), MAKE_VOLATILE_STRIDE(512, csr), MAKE_VOLATILE_STRIDE(512, csi)) {
Chris@82	73 E TcD, TdR, T27, T7r, T5P, T8v, Tf, Ta5, Tu, Tbn, TcG, TdS, T2e, T8w, T5S;
Chris@82	74 E T7s, TK, Ta6, TcK, TdU, T2o, T5U, T7w, T8y, TZ, Ta7, TcN, TdV, T2x, T5V;
Chris@82	75 E T7z, T8z, T1g, Taa, TcU, Tex, TcX, Tew, T1v, Tab, T2M, T6A, T7E, T9b, T7H;
Chris@82	76 E T9a, T2T, T6z, T4X, T6L, Tdz, TeL, TdK, TeO, T5G, T6O, T8d, T9p, TaV, Tc3;
Chris@82	77 E Tbi, Tc4, T8o, T9s, T3I, T6E, Tde, TeE, Tdp, TeH, T4r, T6H, T7U, T9i, Tao;
Chris@82	78 E TbW, TaL, TbX, T85, T9l, T1L, Tad, Td3, TeA, Td6, Tez, T20, Tae, T37, T6x;
Chris@82	79 E T7L, T9e, T7O, T9d, T3e, T6w, TbZ, Tc0, T3Z, T4s, Tds, TeF, T4g, T4t, T80;
Chris@82	80 E T87, Tdl, TeI, T7X, T86, TaD, TaM, Tc6, Tc7, T5e, T5H, TdN, TeM, T5v, T5I;
Chris@82	81 E T8j, T8q, TdG, TeP, T8g, T8p, Tba, Tbj;
Chris@82	82 {
Chris@82	83 E T3, T23, Td, T25, T6, T5N, Ta, T24;
Chris@82	84 {
Chris@82	85 E T1, T2, Tb, Tc;
Chris@82	86 T1 = R0[0];
Chris@82	87 T2 = R0[WS(rs, 32)];
Chris@82	88 T3 = T1 + T2;
Chris@82	89 T23 = T1 - T2;
Chris@82	90 Tb = R0[WS(rs, 56)];
Chris@82	91 Tc = R0[WS(rs, 24)];
Chris@82	92 Td = Tb + Tc;
Chris@82	93 T25 = Tb - Tc;
Chris@82	94 }
Chris@82	95 {
Chris@82	96 E T4, T5, T8, T9;
Chris@82	97 T4 = R0[WS(rs, 16)];
Chris@82	98 T5 = R0[WS(rs, 48)];
Chris@82	99 T6 = T4 + T5;
Chris@82	100 T5N = T4 - T5;
Chris@82	101 T8 = R0[WS(rs, 8)];
Chris@82	102 T9 = R0[WS(rs, 40)];
Chris@82	103 Ta = T8 + T9;
Chris@82	104 T24 = T8 - T9;
Chris@82	105 }
Chris@82	106 TcD = T3 - T6;
Chris@82	107 TdR = Td - Ta;
Chris@82	108 {
Chris@82	109 E T26, T5O, T7, Te;
Chris@82	110 T26 = T24 + T25;
Chris@82	111 T27 = FMA(KP707106781, T26, T23);
Chris@82	112 T7r = FNMS(KP707106781, T26, T23);
Chris@82	113 T5O = T25 - T24;
Chris@82	114 T5P = FNMS(KP707106781, T5O, T5N);
Chris@82	115 T8v = FMA(KP707106781, T5O, T5N);
Chris@82	116 T7 = T3 + T6;
Chris@82	117 Te = Ta + Td;
Chris@82	118 Tf = T7 + Te;
Chris@82	119 Ta5 = T7 - Te;
Chris@82	120 }
Chris@82	121 }
Chris@82	122 {
Chris@82	123 E Ti, T28, Ts, T2c, Tl, T29, Tp, T2b;
Chris@82	124 {
Chris@82	125 E Tg, Th, Tq, Tr;
Chris@82	126 Tg = R0[WS(rs, 4)];
Chris@82	127 Th = R0[WS(rs, 36)];
Chris@82	128 Ti = Tg + Th;
Chris@82	129 T28 = Tg - Th;
Chris@82	130 Tq = R0[WS(rs, 12)];
Chris@82	131 Tr = R0[WS(rs, 44)];
Chris@82	132 Ts = Tq + Tr;
Chris@82	133 T2c = Tq - Tr;
Chris@82	134 }
Chris@82	135 {
Chris@82	136 E Tj, Tk, Tn, To;
Chris@82	137 Tj = R0[WS(rs, 20)];
Chris@82	138 Tk = R0[WS(rs, 52)];
Chris@82	139 Tl = Tj + Tk;
Chris@82	140 T29 = Tj - Tk;
Chris@82	141 Tn = R0[WS(rs, 60)];
Chris@82	142 To = R0[WS(rs, 28)];
Chris@82	143 Tp = Tn + To;
Chris@82	144 T2b = Tn - To;
Chris@82	145 }
Chris@82	146 {
Chris@82	147 E Tm, Tt, TcE, TcF;
Chris@82	148 Tm = Ti + Tl;
Chris@82	149 Tt = Tp + Ts;
Chris@82	150 Tu = Tm + Tt;
Chris@82	151 Tbn = Tt - Tm;
Chris@82	152 TcE = Ti - Tl;
Chris@82	153 TcF = Tp - Ts;
Chris@82	154 TcG = TcE + TcF;
Chris@82	155 TdS = TcF - TcE;
Chris@82	156 }
Chris@82	157 {
Chris@82	158 E T2a, T2d, T5Q, T5R;
Chris@82	159 T2a = FNMS(KP414213562, T29, T28);
Chris@82	160 T2d = FMA(KP414213562, T2c, T2b);
Chris@82	161 T2e = T2a + T2d;
Chris@82	162 T8w = T2d - T2a;
Chris@82	163 T5Q = FMA(KP414213562, T28, T29);
Chris@82	164 T5R = FNMS(KP414213562, T2b, T2c);
Chris@82	165 T5S = T5Q + T5R;
Chris@82	166 T7s = T5Q - T5R;
Chris@82	167 }
Chris@82	168 }
Chris@82	169 {
Chris@82	170 E Ty, T2g, TB, T2l, TF, T2m, TI, T2j;
Chris@82	171 {
Chris@82	172 E Tw, Tx, Tz, TA;
Chris@82	173 Tw = R0[WS(rs, 2)];
Chris@82	174 Tx = R0[WS(rs, 34)];
Chris@82	175 Ty = Tw + Tx;
Chris@82	176 T2g = Tw - Tx;
Chris@82	177 Tz = R0[WS(rs, 18)];
Chris@82	178 TA = R0[WS(rs, 50)];
Chris@82	179 TB = Tz + TA;
Chris@82	180 T2l = Tz - TA;
Chris@82	181 {
Chris@82	182 E TD, TE, T2h, TG, TH, T2i;
Chris@82	183 TD = R0[WS(rs, 10)];
Chris@82	184 TE = R0[WS(rs, 42)];
Chris@82	185 T2h = TD - TE;
Chris@82	186 TG = R0[WS(rs, 58)];
Chris@82	187 TH = R0[WS(rs, 26)];
Chris@82	188 T2i = TG - TH;
Chris@82	189 TF = TD + TE;
Chris@82	190 T2m = T2h - T2i;
Chris@82	191 TI = TG + TH;
Chris@82	192 T2j = T2h + T2i;
Chris@82	193 }
Chris@82	194 }
Chris@82	195 {
Chris@82	196 E TC, TJ, TcI, TcJ;
Chris@82	197 TC = Ty + TB;
Chris@82	198 TJ = TF + TI;
Chris@82	199 TK = TC + TJ;
Chris@82	200 Ta6 = TC - TJ;
Chris@82	201 TcI = Ty - TB;
Chris@82	202 TcJ = TI - TF;
Chris@82	203 TcK = FMA(KP414213562, TcJ, TcI);
Chris@82	204 TdU = FNMS(KP414213562, TcI, TcJ);
Chris@82	205 }
Chris@82	206 {
Chris@82	207 E T2k, T2n, T7u, T7v;
Chris@82	208 T2k = FMA(KP707106781, T2j, T2g);
Chris@82	209 T2n = FMA(KP707106781, T2m, T2l);
Chris@82	210 T2o = FNMS(KP198912367, T2n, T2k);
Chris@82	211 T5U = FMA(KP198912367, T2k, T2n);
Chris@82	212 T7u = FNMS(KP707106781, T2j, T2g);
Chris@82	213 T7v = FNMS(KP707106781, T2m, T2l);
Chris@82	214 T7w = FMA(KP668178637, T7v, T7u);
Chris@82	215 T8y = FNMS(KP668178637, T7u, T7v);
Chris@82	216 }
Chris@82	217 }
Chris@82	218 {
Chris@82	219 E TN, T2p, TQ, T2u, TU, T2v, TX, T2s;
Chris@82	220 {
Chris@82	221 E TL, TM, TO, TP;
Chris@82	222 TL = R0[WS(rs, 62)];
Chris@82	223 TM = R0[WS(rs, 30)];
Chris@82	224 TN = TL + TM;
Chris@82	225 T2p = TL - TM;
Chris@82	226 TO = R0[WS(rs, 14)];
Chris@82	227 TP = R0[WS(rs, 46)];
Chris@82	228 TQ = TO + TP;
Chris@82	229 T2u = TO - TP;
Chris@82	230 {
Chris@82	231 E TS, TT, T2q, TV, TW, T2r;
Chris@82	232 TS = R0[WS(rs, 6)];
Chris@82	233 TT = R0[WS(rs, 38)];
Chris@82	234 T2q = TS - TT;
Chris@82	235 TV = R0[WS(rs, 54)];
Chris@82	236 TW = R0[WS(rs, 22)];
Chris@82	237 T2r = TV - TW;
Chris@82	238 TU = TS + TT;
Chris@82	239 T2v = T2q - T2r;
Chris@82	240 TX = TV + TW;
Chris@82	241 T2s = T2q + T2r;
Chris@82	242 }
Chris@82	243 }
Chris@82	244 {
Chris@82	245 E TR, TY, TcL, TcM;
Chris@82	246 TR = TN + TQ;
Chris@82	247 TY = TU + TX;
Chris@82	248 TZ = TR + TY;
Chris@82	249 Ta7 = TR - TY;
Chris@82	250 TcL = TN - TQ;
Chris@82	251 TcM = TX - TU;
Chris@82	252 TcN = FNMS(KP414213562, TcM, TcL);
Chris@82	253 TdV = FMA(KP414213562, TcL, TcM);
Chris@82	254 }
Chris@82	255 {
Chris@82	256 E T2t, T2w, T7x, T7y;
Chris@82	257 T2t = FMA(KP707106781, T2s, T2p);
Chris@82	258 T2w = FMA(KP707106781, T2v, T2u);
Chris@82	259 T2x = FMA(KP198912367, T2w, T2t);
Chris@82	260 T5V = FNMS(KP198912367, T2t, T2w);
Chris@82	261 T7x = FNMS(KP707106781, T2s, T2p);
Chris@82	262 T7y = FNMS(KP707106781, T2v, T2u);
Chris@82	263 T7z = FNMS(KP668178637, T7y, T7x);
Chris@82	264 T8z = FMA(KP668178637, T7x, T7y);
Chris@82	265 }
Chris@82	266 }
Chris@82	267 {
Chris@82	268 E T14, T2A, T17, T2N, T1b, T2O, T1e, T2D, T1j, T1m, T2H, TcR, T2Q, T1q, T1t;
Chris@82	269 E T2K, TcS, T2R;
Chris@82	270 {
Chris@82	271 E T12, T13, T15, T16;
Chris@82	272 T12 = R0[WS(rs, 1)];
Chris@82	273 T13 = R0[WS(rs, 33)];
Chris@82	274 T14 = T12 + T13;
Chris@82	275 T2A = T12 - T13;
Chris@82	276 T15 = R0[WS(rs, 17)];
Chris@82	277 T16 = R0[WS(rs, 49)];
Chris@82	278 T17 = T15 + T16;
Chris@82	279 T2N = T15 - T16;
Chris@82	280 }
Chris@82	281 {
Chris@82	282 E T19, T1a, T2B, T1c, T1d, T2C;
Chris@82	283 T19 = R0[WS(rs, 9)];
Chris@82	284 T1a = R0[WS(rs, 41)];
Chris@82	285 T2B = T19 - T1a;
Chris@82	286 T1c = R0[WS(rs, 57)];
Chris@82	287 T1d = R0[WS(rs, 25)];
Chris@82	288 T2C = T1c - T1d;
Chris@82	289 T1b = T19 + T1a;
Chris@82	290 T2O = T2B - T2C;
Chris@82	291 T1e = T1c + T1d;
Chris@82	292 T2D = T2B + T2C;
Chris@82	293 }
Chris@82	294 {
Chris@82	295 E T2F, T2G, T2I, T2J;
Chris@82	296 {
Chris@82	297 E T1h, T1i, T1k, T1l;
Chris@82	298 T1h = R0[WS(rs, 5)];
Chris@82	299 T1i = R0[WS(rs, 37)];
Chris@82	300 T1j = T1h + T1i;
Chris@82	301 T2F = T1h - T1i;
Chris@82	302 T1k = R0[WS(rs, 21)];
Chris@82	303 T1l = R0[WS(rs, 53)];
Chris@82	304 T1m = T1k + T1l;
Chris@82	305 T2G = T1k - T1l;
Chris@82	306 }
Chris@82	307 T2H = FNMS(KP414213562, T2G, T2F);
Chris@82	308 TcR = T1j - T1m;
Chris@82	309 T2Q = FMA(KP414213562, T2F, T2G);
Chris@82	310 {
Chris@82	311 E T1o, T1p, T1r, T1s;
Chris@82	312 T1o = R0[WS(rs, 61)];
Chris@82	313 T1p = R0[WS(rs, 29)];
Chris@82	314 T1q = T1o + T1p;
Chris@82	315 T2I = T1o - T1p;
Chris@82	316 T1r = R0[WS(rs, 13)];
Chris@82	317 T1s = R0[WS(rs, 45)];
Chris@82	318 T1t = T1r + T1s;
Chris@82	319 T2J = T1r - T1s;
Chris@82	320 }
Chris@82	321 T2K = FMA(KP414213562, T2J, T2I);
Chris@82	322 TcS = T1q - T1t;
Chris@82	323 T2R = FNMS(KP414213562, T2I, T2J);
Chris@82	324 }
Chris@82	325 {
Chris@82	326 E T18, T1f, TcQ, TcT;
Chris@82	327 T18 = T14 + T17;
Chris@82	328 T1f = T1b + T1e;
Chris@82	329 T1g = T18 + T1f;
Chris@82	330 Taa = T18 - T1f;
Chris@82	331 TcQ = T14 - T17;
Chris@82	332 TcT = TcR + TcS;
Chris@82	333 TcU = FMA(KP707106781, TcT, TcQ);
Chris@82	334 Tex = FNMS(KP707106781, TcT, TcQ);
Chris@82	335 }
Chris@82	336 {
Chris@82	337 E TcV, TcW, T1n, T1u;
Chris@82	338 TcV = T1e - T1b;
Chris@82	339 TcW = TcS - TcR;
Chris@82	340 TcX = FMA(KP707106781, TcW, TcV);
Chris@82	341 Tew = FNMS(KP707106781, TcW, TcV);
Chris@82	342 T1n = T1j + T1m;
Chris@82	343 T1u = T1q + T1t;
Chris@82	344 T1v = T1n + T1u;
Chris@82	345 Tab = T1u - T1n;
Chris@82	346 }
Chris@82	347 {
Chris@82	348 E T2E, T2L, T7C, T7D;
Chris@82	349 T2E = FMA(KP707106781, T2D, T2A);
Chris@82	350 T2L = T2H + T2K;
Chris@82	351 T2M = FMA(KP923879532, T2L, T2E);
Chris@82	352 T6A = FNMS(KP923879532, T2L, T2E);
Chris@82	353 T7C = FNMS(KP707106781, T2D, T2A);
Chris@82	354 T7D = T2Q - T2R;
Chris@82	355 T7E = FMA(KP923879532, T7D, T7C);
Chris@82	356 T9b = FNMS(KP923879532, T7D, T7C);
Chris@82	357 }
Chris@82	358 {
Chris@82	359 E T7F, T7G, T2P, T2S;
Chris@82	360 T7F = FNMS(KP707106781, T2O, T2N);
Chris@82	361 T7G = T2K - T2H;
Chris@82	362 T7H = FMA(KP923879532, T7G, T7F);
Chris@82	363 T9a = FNMS(KP923879532, T7G, T7F);
Chris@82	364 T2P = FMA(KP707106781, T2O, T2N);
Chris@82	365 T2S = T2Q + T2R;
Chris@82	366 T2T = FMA(KP923879532, T2S, T2P);
Chris@82	367 T6z = FNMS(KP923879532, T2S, T2P);
Chris@82	368 }
Chris@82	369 }
Chris@82	370 {
Chris@82	371 E T4z, TaP, T5A, TaQ, T4G, TaT, T5B, TaS, Tbf, Tbg, T4O, Tdw, T5D, Tbc, Tbd;
Chris@82	372 E T4V, Tdx, T5E;
Chris@82	373 {
Chris@82	374 E T4x, T4y, T5y, T5z;
Chris@82	375 T4x = R1[WS(rs, 63)];
Chris@82	376 T4y = R1[WS(rs, 31)];
Chris@82	377 T4z = T4x - T4y;
Chris@82	378 TaP = T4x + T4y;
Chris@82	379 T5y = R1[WS(rs, 47)];
Chris@82	380 T5z = R1[WS(rs, 15)];
Chris@82	381 T5A = T5y - T5z;
Chris@82	382 TaQ = T5z + T5y;
Chris@82	383 }
Chris@82	384 {
Chris@82	385 E T4A, T4B, T4C, T4D, T4E, T4F;
Chris@82	386 T4A = R1[WS(rs, 7)];
Chris@82	387 T4B = R1[WS(rs, 39)];
Chris@82	388 T4C = T4A - T4B;
Chris@82	389 T4D = R1[WS(rs, 55)];
Chris@82	390 T4E = R1[WS(rs, 23)];
Chris@82	391 T4F = T4D - T4E;
Chris@82	392 T4G = T4C + T4F;
Chris@82	393 TaT = T4D + T4E;
Chris@82	394 T5B = T4F - T4C;
Chris@82	395 TaS = T4A + T4B;
Chris@82	396 }
Chris@82	397 {
Chris@82	398 E T4K, T4N, T4R, T4U;
Chris@82	399 {
Chris@82	400 E T4I, T4J, T4L, T4M;
Chris@82	401 T4I = R1[WS(rs, 3)];
Chris@82	402 T4J = R1[WS(rs, 35)];
Chris@82	403 T4K = T4I - T4J;
Chris@82	404 Tbf = T4I + T4J;
Chris@82	405 T4L = R1[WS(rs, 51)];
Chris@82	406 T4M = R1[WS(rs, 19)];
Chris@82	407 T4N = T4L - T4M;
Chris@82	408 Tbg = T4M + T4L;
Chris@82	409 }
Chris@82	410 T4O = FMA(KP414213562, T4N, T4K);
Chris@82	411 Tdw = Tbf - Tbg;
Chris@82	412 T5D = FNMS(KP414213562, T4K, T4N);
Chris@82	413 {
Chris@82	414 E T4P, T4Q, T4S, T4T;
Chris@82	415 T4P = R1[WS(rs, 59)];
Chris@82	416 T4Q = R1[WS(rs, 27)];
Chris@82	417 T4R = T4P - T4Q;
Chris@82	418 Tbc = T4P + T4Q;
Chris@82	419 T4S = R1[WS(rs, 43)];
Chris@82	420 T4T = R1[WS(rs, 11)];
Chris@82	421 T4U = T4S - T4T;
Chris@82	422 Tbd = T4T + T4S;
Chris@82	423 }
Chris@82	424 T4V = FNMS(KP414213562, T4U, T4R);
Chris@82	425 Tdx = Tbc - Tbd;
Chris@82	426 T5E = FMA(KP414213562, T4R, T4U);
Chris@82	427 }
Chris@82	428 {
Chris@82	429 E T4H, T4W, Tdv, Tdy;
Chris@82	430 T4H = FMA(KP707106781, T4G, T4z);
Chris@82	431 T4W = T4O + T4V;
Chris@82	432 T4X = FMA(KP923879532, T4W, T4H);
Chris@82	433 T6L = FNMS(KP923879532, T4W, T4H);
Chris@82	434 Tdv = TaP - TaQ;
Chris@82	435 Tdy = Tdw + Tdx;
Chris@82	436 Tdz = FMA(KP707106781, Tdy, Tdv);
Chris@82	437 TeL = FNMS(KP707106781, Tdy, Tdv);
Chris@82	438 }
Chris@82	439 {
Chris@82	440 E TdI, TdJ, T5C, T5F;
Chris@82	441 TdI = TaT - TaS;
Chris@82	442 TdJ = Tdx - Tdw;
Chris@82	443 TdK = FMA(KP707106781, TdJ, TdI);
Chris@82	444 TeO = FNMS(KP707106781, TdJ, TdI);
Chris@82	445 T5C = FMA(KP707106781, T5B, T5A);
Chris@82	446 T5F = T5D + T5E;
Chris@82	447 T5G = FMA(KP923879532, T5F, T5C);
Chris@82	448 T6O = FNMS(KP923879532, T5F, T5C);
Chris@82	449 }
Chris@82	450 {
Chris@82	451 E T8b, T8c, TaR, TaU;
Chris@82	452 T8b = FNMS(KP707106781, T4G, T4z);
Chris@82	453 T8c = T5E - T5D;
Chris@82	454 T8d = FMA(KP923879532, T8c, T8b);
Chris@82	455 T9p = FNMS(KP923879532, T8c, T8b);
Chris@82	456 TaR = TaP + TaQ;
Chris@82	457 TaU = TaS + TaT;
Chris@82	458 TaV = TaR - TaU;
Chris@82	459 Tc3 = TaR + TaU;
Chris@82	460 }
Chris@82	461 {
Chris@82	462 E Tbe, Tbh, T8m, T8n;
Chris@82	463 Tbe = Tbc + Tbd;
Chris@82	464 Tbh = Tbf + Tbg;
Chris@82	465 Tbi = Tbe - Tbh;
Chris@82	466 Tc4 = Tbh + Tbe;
Chris@82	467 T8m = FNMS(KP707106781, T5B, T5A);
Chris@82	468 T8n = T4V - T4O;
Chris@82	469 T8o = FNMS(KP923879532, T8n, T8m);
Chris@82	470 T9s = FMA(KP923879532, T8n, T8m);
Chris@82	471 }
Chris@82	472 }
Chris@82	473 {
Chris@82	474 E T3k, Tai, T4l, Taj, T3r, Tam, T4m, Tal, TaI, TaJ, T3z, Tdb, T4o, TaF, TaG;
Chris@82	475 E T3G, Tdc, T4p;
Chris@82	476 {
Chris@82	477 E T3i, T3j, T4j, T4k;
Chris@82	478 T3i = R1[0];
Chris@82	479 T3j = R1[WS(rs, 32)];
Chris@82	480 T3k = T3i - T3j;
Chris@82	481 Tai = T3i + T3j;
Chris@82	482 T4j = R1[WS(rs, 16)];
Chris@82	483 T4k = R1[WS(rs, 48)];
Chris@82	484 T4l = T4j - T4k;
Chris@82	485 Taj = T4j + T4k;
Chris@82	486 }
Chris@82	487 {
Chris@82	488 E T3l, T3m, T3n, T3o, T3p, T3q;
Chris@82	489 T3l = R1[WS(rs, 8)];
Chris@82	490 T3m = R1[WS(rs, 40)];
Chris@82	491 T3n = T3l - T3m;
Chris@82	492 T3o = R1[WS(rs, 56)];
Chris@82	493 T3p = R1[WS(rs, 24)];
Chris@82	494 T3q = T3o - T3p;
Chris@82	495 T3r = T3n + T3q;
Chris@82	496 Tam = T3o + T3p;
Chris@82	497 T4m = T3n - T3q;
Chris@82	498 Tal = T3l + T3m;
Chris@82	499 }
Chris@82	500 {
Chris@82	501 E T3v, T3y, T3C, T3F;
Chris@82	502 {
Chris@82	503 E T3t, T3u, T3w, T3x;
Chris@82	504 T3t = R1[WS(rs, 4)];
Chris@82	505 T3u = R1[WS(rs, 36)];
Chris@82	506 T3v = T3t - T3u;
Chris@82	507 TaI = T3t + T3u;
Chris@82	508 T3w = R1[WS(rs, 20)];
Chris@82	509 T3x = R1[WS(rs, 52)];
Chris@82	510 T3y = T3w - T3x;
Chris@82	511 TaJ = T3w + T3x;
Chris@82	512 }
Chris@82	513 T3z = FNMS(KP414213562, T3y, T3v);
Chris@82	514 Tdb = TaI - TaJ;
Chris@82	515 T4o = FMA(KP414213562, T3v, T3y);
Chris@82	516 {
Chris@82	517 E T3A, T3B, T3D, T3E;
Chris@82	518 T3A = R1[WS(rs, 60)];
Chris@82	519 T3B = R1[WS(rs, 28)];
Chris@82	520 T3C = T3A - T3B;
Chris@82	521 TaF = T3A + T3B;
Chris@82	522 T3D = R1[WS(rs, 12)];
Chris@82	523 T3E = R1[WS(rs, 44)];
Chris@82	524 T3F = T3D - T3E;
Chris@82	525 TaG = T3D + T3E;
Chris@82	526 }
Chris@82	527 T3G = FMA(KP414213562, T3F, T3C);
Chris@82	528 Tdc = TaF - TaG;
Chris@82	529 T4p = FNMS(KP414213562, T3C, T3F);
Chris@82	530 }
Chris@82	531 {
Chris@82	532 E T3s, T3H, Tda, Tdd;
Chris@82	533 T3s = FMA(KP707106781, T3r, T3k);
Chris@82	534 T3H = T3z + T3G;
Chris@82	535 T3I = FMA(KP923879532, T3H, T3s);
Chris@82	536 T6E = FNMS(KP923879532, T3H, T3s);
Chris@82	537 Tda = Tai - Taj;
Chris@82	538 Tdd = Tdb + Tdc;
Chris@82	539 Tde = FMA(KP707106781, Tdd, Tda);
Chris@82	540 TeE = FNMS(KP707106781, Tdd, Tda);
Chris@82	541 }
Chris@82	542 {
Chris@82	543 E Tdn, Tdo, T4n, T4q;
Chris@82	544 Tdn = Tam - Tal;
Chris@82	545 Tdo = Tdc - Tdb;
Chris@82	546 Tdp = FMA(KP707106781, Tdo, Tdn);
Chris@82	547 TeH = FNMS(KP707106781, Tdo, Tdn);
Chris@82	548 T4n = FMA(KP707106781, T4m, T4l);
Chris@82	549 T4q = T4o + T4p;
Chris@82	550 T4r = FMA(KP923879532, T4q, T4n);
Chris@82	551 T6H = FNMS(KP923879532, T4q, T4n);
Chris@82	552 }
Chris@82	553 {
Chris@82	554 E T7S, T7T, Tak, Tan;
Chris@82	555 T7S = FNMS(KP707106781, T3r, T3k);
Chris@82	556 T7T = T4o - T4p;
Chris@82	557 T7U = FMA(KP923879532, T7T, T7S);
Chris@82	558 T9i = FNMS(KP923879532, T7T, T7S);
Chris@82	559 Tak = Tai + Taj;
Chris@82	560 Tan = Tal + Tam;
Chris@82	561 Tao = Tak - Tan;
Chris@82	562 TbW = Tak + Tan;
Chris@82	563 }
Chris@82	564 {
Chris@82	565 E TaH, TaK, T83, T84;
Chris@82	566 TaH = TaF + TaG;
Chris@82	567 TaK = TaI + TaJ;
Chris@82	568 TaL = TaH - TaK;
Chris@82	569 TbX = TaK + TaH;
Chris@82	570 T83 = FNMS(KP707106781, T4m, T4l);
Chris@82	571 T84 = T3G - T3z;
Chris@82	572 T85 = FMA(KP923879532, T84, T83);
Chris@82	573 T9l = FNMS(KP923879532, T84, T83);
Chris@82	574 }
Chris@82	575 }
Chris@82	576 {
Chris@82	577 E T1z, T2V, T1C, T38, T1G, T39, T1J, T2Y, T1O, T1R, T32, Td0, T3b, T1V, T1Y;
Chris@82	578 E T35, Td1, T3c;
Chris@82	579 {
Chris@82	580 E T1x, T1y, T1A, T1B;
Chris@82	581 T1x = R0[WS(rs, 63)];
Chris@82	582 T1y = R0[WS(rs, 31)];
Chris@82	583 T1z = T1x + T1y;
Chris@82	584 T2V = T1x - T1y;
Chris@82	585 T1A = R0[WS(rs, 15)];
Chris@82	586 T1B = R0[WS(rs, 47)];
Chris@82	587 T1C = T1A + T1B;
Chris@82	588 T38 = T1A - T1B;
Chris@82	589 }
Chris@82	590 {
Chris@82	591 E T1E, T1F, T2W, T1H, T1I, T2X;
Chris@82	592 T1E = R0[WS(rs, 7)];
Chris@82	593 T1F = R0[WS(rs, 39)];
Chris@82	594 T2W = T1E - T1F;
Chris@82	595 T1H = R0[WS(rs, 55)];
Chris@82	596 T1I = R0[WS(rs, 23)];
Chris@82	597 T2X = T1H - T1I;
Chris@82	598 T1G = T1E + T1F;
Chris@82	599 T39 = T2W - T2X;
Chris@82	600 T1J = T1H + T1I;
Chris@82	601 T2Y = T2W + T2X;
Chris@82	602 }
Chris@82	603 {
Chris@82	604 E T30, T31, T33, T34;
Chris@82	605 {
Chris@82	606 E T1M, T1N, T1P, T1Q;
Chris@82	607 T1M = R0[WS(rs, 3)];
Chris@82	608 T1N = R0[WS(rs, 35)];
Chris@82	609 T1O = T1M + T1N;
Chris@82	610 T30 = T1M - T1N;
Chris@82	611 T1P = R0[WS(rs, 19)];
Chris@82	612 T1Q = R0[WS(rs, 51)];
Chris@82	613 T1R = T1P + T1Q;
Chris@82	614 T31 = T1P - T1Q;
Chris@82	615 }
Chris@82	616 T32 = FNMS(KP414213562, T31, T30);
Chris@82	617 Td0 = T1O - T1R;
Chris@82	618 T3b = FMA(KP414213562, T30, T31);
Chris@82	619 {
Chris@82	620 E T1T, T1U, T1W, T1X;
Chris@82	621 T1T = R0[WS(rs, 59)];
Chris@82	622 T1U = R0[WS(rs, 27)];
Chris@82	623 T1V = T1T + T1U;
Chris@82	624 T33 = T1T - T1U;
Chris@82	625 T1W = R0[WS(rs, 11)];
Chris@82	626 T1X = R0[WS(rs, 43)];
Chris@82	627 T1Y = T1W + T1X;
Chris@82	628 T34 = T1W - T1X;
Chris@82	629 }
Chris@82	630 T35 = FMA(KP414213562, T34, T33);
Chris@82	631 Td1 = T1V - T1Y;
Chris@82	632 T3c = FNMS(KP414213562, T33, T34);
Chris@82	633 }
Chris@82	634 {
Chris@82	635 E T1D, T1K, TcZ, Td2;
Chris@82	636 T1D = T1z + T1C;
Chris@82	637 T1K = T1G + T1J;
Chris@82	638 T1L = T1D + T1K;
Chris@82	639 Tad = T1D - T1K;
Chris@82	640 TcZ = T1z - T1C;
Chris@82	641 Td2 = Td0 + Td1;
Chris@82	642 Td3 = FMA(KP707106781, Td2, TcZ);
Chris@82	643 TeA = FNMS(KP707106781, Td2, TcZ);
Chris@82	644 }
Chris@82	645 {
Chris@82	646 E Td4, Td5, T1S, T1Z;
Chris@82	647 Td4 = T1J - T1G;
Chris@82	648 Td5 = Td1 - Td0;
Chris@82	649 Td6 = FMA(KP707106781, Td5, Td4);
Chris@82	650 Tez = FNMS(KP707106781, Td5, Td4);
Chris@82	651 T1S = T1O + T1R;
Chris@82	652 T1Z = T1V + T1Y;
Chris@82	653 T20 = T1S + T1Z;
Chris@82	654 Tae = T1Z - T1S;
Chris@82	655 }
Chris@82	656 {
Chris@82	657 E T2Z, T36, T7J, T7K;
Chris@82	658 T2Z = FMA(KP707106781, T2Y, T2V);
Chris@82	659 T36 = T32 + T35;
Chris@82	660 T37 = FMA(KP923879532, T36, T2Z);
Chris@82	661 T6x = FNMS(KP923879532, T36, T2Z);
Chris@82	662 T7J = FNMS(KP707106781, T2Y, T2V);
Chris@82	663 T7K = T3b - T3c;
Chris@82	664 T7L = FMA(KP923879532, T7K, T7J);
Chris@82	665 T9e = FNMS(KP923879532, T7K, T7J);
Chris@82	666 }
Chris@82	667 {
Chris@82	668 E T7M, T7N, T3a, T3d;
Chris@82	669 T7M = FNMS(KP707106781, T39, T38);
Chris@82	670 T7N = T35 - T32;
Chris@82	671 T7O = FMA(KP923879532, T7N, T7M);
Chris@82	672 T9d = FNMS(KP923879532, T7N, T7M);
Chris@82	673 T3a = FMA(KP707106781, T39, T38);
Chris@82	674 T3d = T3b + T3c;
Chris@82	675 T3e = FMA(KP923879532, T3d, T3a);
Chris@82	676 T6w = FNMS(KP923879532, T3d, T3a);
Chris@82	677 }
Chris@82	678 }
Chris@82	679 {
Chris@82	680 E T3L, Tdf, T3W, Tar, T42, Tdi, T4d, Tay, T3S, Tdg, T3X, Tau, T49, Tdj, T4e;
Chris@82	681 E TaB, Tdh, Tdk;
Chris@82	682 {
Chris@82	683 E T3J, T3K, Tap, T3U, T3V, Taq;
Chris@82	684 T3J = R1[WS(rs, 2)];
Chris@82	685 T3K = R1[WS(rs, 34)];
Chris@82	686 Tap = T3J + T3K;
Chris@82	687 T3U = R1[WS(rs, 18)];
Chris@82	688 T3V = R1[WS(rs, 50)];
Chris@82	689 Taq = T3U + T3V;
Chris@82	690 T3L = T3J - T3K;
Chris@82	691 Tdf = Tap - Taq;
Chris@82	692 T3W = T3U - T3V;
Chris@82	693 Tar = Tap + Taq;
Chris@82	694 }
Chris@82	695 {
Chris@82	696 E T40, T41, Taw, T4b, T4c, Tax;
Chris@82	697 T40 = R1[WS(rs, 62)];
Chris@82	698 T41 = R1[WS(rs, 30)];
Chris@82	699 Taw = T40 + T41;
Chris@82	700 T4b = R1[WS(rs, 14)];
Chris@82	701 T4c = R1[WS(rs, 46)];
Chris@82	702 Tax = T4b + T4c;
Chris@82	703 T42 = T40 - T41;
Chris@82	704 Tdi = Taw - Tax;
Chris@82	705 T4d = T4b - T4c;
Chris@82	706 Tay = Taw + Tax;
Chris@82	707 }
Chris@82	708 {
Chris@82	709 E T3O, Tas, T3R, Tat;
Chris@82	710 {
Chris@82	711 E T3M, T3N, T3P, T3Q;
Chris@82	712 T3M = R1[WS(rs, 10)];
Chris@82	713 T3N = R1[WS(rs, 42)];
Chris@82	714 T3O = T3M - T3N;
Chris@82	715 Tas = T3M + T3N;
Chris@82	716 T3P = R1[WS(rs, 58)];
Chris@82	717 T3Q = R1[WS(rs, 26)];
Chris@82	718 T3R = T3P - T3Q;
Chris@82	719 Tat = T3P + T3Q;
Chris@82	720 }
Chris@82	721 T3S = T3O + T3R;
Chris@82	722 Tdg = Tat - Tas;
Chris@82	723 T3X = T3O - T3R;
Chris@82	724 Tau = Tas + Tat;
Chris@82	725 }
Chris@82	726 {
Chris@82	727 E T45, Taz, T48, TaA;
Chris@82	728 {
Chris@82	729 E T43, T44, T46, T47;
Chris@82	730 T43 = R1[WS(rs, 6)];
Chris@82	731 T44 = R1[WS(rs, 38)];
Chris@82	732 T45 = T43 - T44;
Chris@82	733 Taz = T43 + T44;
Chris@82	734 T46 = R1[WS(rs, 54)];
Chris@82	735 T47 = R1[WS(rs, 22)];
Chris@82	736 T48 = T46 - T47;
Chris@82	737 TaA = T46 + T47;
Chris@82	738 }
Chris@82	739 T49 = T45 + T48;
Chris@82	740 Tdj = TaA - Taz;
Chris@82	741 T4e = T45 - T48;
Chris@82	742 TaB = Taz + TaA;
Chris@82	743 }
Chris@82	744 TbZ = Tar + Tau;
Chris@82	745 Tc0 = Tay + TaB;
Chris@82	746 {
Chris@82	747 E T3T, T3Y, Tdq, Tdr;
Chris@82	748 T3T = FMA(KP707106781, T3S, T3L);
Chris@82	749 T3Y = FMA(KP707106781, T3X, T3W);
Chris@82	750 T3Z = FNMS(KP198912367, T3Y, T3T);
Chris@82	751 T4s = FMA(KP198912367, T3T, T3Y);
Chris@82	752 Tdq = FNMS(KP414213562, Tdf, Tdg);
Chris@82	753 Tdr = FMA(KP414213562, Tdi, Tdj);
Chris@82	754 Tds = Tdq + Tdr;
Chris@82	755 TeF = Tdr - Tdq;
Chris@82	756 }
Chris@82	757 {
Chris@82	758 E T4a, T4f, T7Y, T7Z;
Chris@82	759 T4a = FMA(KP707106781, T49, T42);
Chris@82	760 T4f = FMA(KP707106781, T4e, T4d);
Chris@82	761 T4g = FMA(KP198912367, T4f, T4a);
Chris@82	762 T4t = FNMS(KP198912367, T4a, T4f);
Chris@82	763 T7Y = FNMS(KP707106781, T49, T42);
Chris@82	764 T7Z = FNMS(KP707106781, T4e, T4d);
Chris@82	765 T80 = FNMS(KP668178637, T7Z, T7Y);
Chris@82	766 T87 = FMA(KP668178637, T7Y, T7Z);
Chris@82	767 }
Chris@82	768 Tdh = FMA(KP414213562, Tdg, Tdf);
Chris@82	769 Tdk = FNMS(KP414213562, Tdj, Tdi);
Chris@82	770 Tdl = Tdh + Tdk;
Chris@82	771 TeI = Tdh - Tdk;
Chris@82	772 {
Chris@82	773 E T7V, T7W, Tav, TaC;
Chris@82	774 T7V = FNMS(KP707106781, T3S, T3L);
Chris@82	775 T7W = FNMS(KP707106781, T3X, T3W);
Chris@82	776 T7X = FMA(KP668178637, T7W, T7V);
Chris@82	777 T86 = FNMS(KP668178637, T7V, T7W);
Chris@82	778 Tav = Tar - Tau;
Chris@82	779 TaC = Tay - TaB;
Chris@82	780 TaD = Tav + TaC;
Chris@82	781 TaM = TaC - Tav;
Chris@82	782 }
Chris@82	783 }
Chris@82	784 {
Chris@82	785 E T50, TdA, T5b, TaY, T5h, TdD, T5s, Tb5, T57, TdB, T5c, Tb1, T5o, TdE, T5t;
Chris@82	786 E Tb8, TdC, TdF;
Chris@82	787 {
Chris@82	788 E T4Y, T4Z, TaW, T59, T5a, TaX;
Chris@82	789 T4Y = R1[WS(rs, 1)];
Chris@82	790 T4Z = R1[WS(rs, 33)];
Chris@82	791 TaW = T4Y + T4Z;
Chris@82	792 T59 = R1[WS(rs, 49)];
Chris@82	793 T5a = R1[WS(rs, 17)];
Chris@82	794 TaX = T5a + T59;
Chris@82	795 T50 = T4Y - T4Z;
Chris@82	796 TdA = TaW - TaX;
Chris@82	797 T5b = T59 - T5a;
Chris@82	798 TaY = TaW + TaX;
Chris@82	799 }
Chris@82	800 {
Chris@82	801 E T5f, T5g, Tb3, T5q, T5r, Tb4;
Chris@82	802 T5f = R1[WS(rs, 61)];
Chris@82	803 T5g = R1[WS(rs, 29)];
Chris@82	804 Tb3 = T5f + T5g;
Chris@82	805 T5q = R1[WS(rs, 45)];
Chris@82	806 T5r = R1[WS(rs, 13)];
Chris@82	807 Tb4 = T5r + T5q;
Chris@82	808 T5h = T5f - T5g;
Chris@82	809 TdD = Tb3 - Tb4;
Chris@82	810 T5s = T5q - T5r;
Chris@82	811 Tb5 = Tb3 + Tb4;
Chris@82	812 }
Chris@82	813 {
Chris@82	814 E T53, TaZ, T56, Tb0;
Chris@82	815 {
Chris@82	816 E T51, T52, T54, T55;
Chris@82	817 T51 = R1[WS(rs, 9)];
Chris@82	818 T52 = R1[WS(rs, 41)];
Chris@82	819 T53 = T51 - T52;
Chris@82	820 TaZ = T51 + T52;
Chris@82	821 T54 = R1[WS(rs, 57)];
Chris@82	822 T55 = R1[WS(rs, 25)];
Chris@82	823 T56 = T54 - T55;
Chris@82	824 Tb0 = T54 + T55;
Chris@82	825 }
Chris@82	826 T57 = T53 + T56;
Chris@82	827 TdB = Tb0 - TaZ;
Chris@82	828 T5c = T56 - T53;
Chris@82	829 Tb1 = TaZ + Tb0;
Chris@82	830 }
Chris@82	831 {
Chris@82	832 E T5k, Tb6, T5n, Tb7;
Chris@82	833 {
Chris@82	834 E T5i, T5j, T5l, T5m;
Chris@82	835 T5i = R1[WS(rs, 5)];
Chris@82	836 T5j = R1[WS(rs, 37)];
Chris@82	837 T5k = T5i - T5j;
Chris@82	838 Tb6 = T5i + T5j;
Chris@82	839 T5l = R1[WS(rs, 53)];
Chris@82	840 T5m = R1[WS(rs, 21)];
Chris@82	841 T5n = T5l - T5m;
Chris@82	842 Tb7 = T5l + T5m;
Chris@82	843 }
Chris@82	844 T5o = T5k + T5n;
Chris@82	845 TdE = Tb7 - Tb6;
Chris@82	846 T5t = T5n - T5k;
Chris@82	847 Tb8 = Tb6 + Tb7;
Chris@82	848 }
Chris@82	849 Tc6 = TaY + Tb1;
Chris@82	850 Tc7 = Tb5 + Tb8;
Chris@82	851 {
Chris@82	852 E T58, T5d, TdL, TdM;
Chris@82	853 T58 = FMA(KP707106781, T57, T50);
Chris@82	854 T5d = FMA(KP707106781, T5c, T5b);
Chris@82	855 T5e = FMA(KP198912367, T5d, T58);
Chris@82	856 T5H = FNMS(KP198912367, T58, T5d);
Chris@82	857 TdL = FNMS(KP414213562, TdA, TdB);
Chris@82	858 TdM = FMA(KP414213562, TdD, TdE);
Chris@82	859 TdN = TdL + TdM;
Chris@82	860 TeM = TdM - TdL;
Chris@82	861 }
Chris@82	862 {
Chris@82	863 E T5p, T5u, T8h, T8i;
Chris@82	864 T5p = FMA(KP707106781, T5o, T5h);
Chris@82	865 T5u = FMA(KP707106781, T5t, T5s);
Chris@82	866 T5v = FNMS(KP198912367, T5u, T5p);
Chris@82	867 T5I = FMA(KP198912367, T5p, T5u);
Chris@82	868 T8h = FNMS(KP707106781, T5o, T5h);
Chris@82	869 T8i = FNMS(KP707106781, T5t, T5s);
Chris@82	870 T8j = FMA(KP668178637, T8i, T8h);
Chris@82	871 T8q = FNMS(KP668178637, T8h, T8i);
Chris@82	872 }
Chris@82	873 TdC = FMA(KP414213562, TdB, TdA);
Chris@82	874 TdF = FNMS(KP414213562, TdE, TdD);
Chris@82	875 TdG = TdC + TdF;
Chris@82	876 TeP = TdF - TdC;
Chris@82	877 {
Chris@82	878 E T8e, T8f, Tb2, Tb9;
Chris@82	879 T8e = FNMS(KP707106781, T57, T50);
Chris@82	880 T8f = FNMS(KP707106781, T5c, T5b);
Chris@82	881 T8g = FNMS(KP668178637, T8f, T8e);
Chris@82	882 T8p = FMA(KP668178637, T8e, T8f);
Chris@82	883 Tb2 = TaY - Tb1;
Chris@82	884 Tb9 = Tb5 - Tb8;
Chris@82	885 Tba = Tb2 + Tb9;
Chris@82	886 Tbj = Tb9 - Tb2;
Chris@82	887 }
Chris@82	888 }
Chris@82	889 {
Chris@82	890 E T11, TbV, Tc9, Tcf, T22, Tcb, Tc2, Tce;
Chris@82	891 {
Chris@82	892 E Tv, T10, Tc5, Tc8;
Chris@82	893 Tv = Tf + Tu;
Chris@82	894 T10 = TK + TZ;
Chris@82	895 T11 = Tv + T10;
Chris@82	896 TbV = Tv - T10;
Chris@82	897 Tc5 = Tc3 + Tc4;
Chris@82	898 Tc8 = Tc6 + Tc7;
Chris@82	899 Tc9 = Tc5 - Tc8;
Chris@82	900 Tcf = Tc5 + Tc8;
Chris@82	901 }
Chris@82	902 {
Chris@82	903 E T1w, T21, TbY, Tc1;
Chris@82	904 T1w = T1g + T1v;
Chris@82	905 T21 = T1L + T20;
Chris@82	906 T22 = T1w + T21;
Chris@82	907 Tcb = T21 - T1w;
Chris@82	908 TbY = TbW + TbX;
Chris@82	909 Tc1 = TbZ + Tc0;
Chris@82	910 Tc2 = TbY - Tc1;
Chris@82	911 Tce = TbY + Tc1;
Chris@82	912 }
Chris@82	913 Cr[WS(csr, 32)] = T11 - T22;
Chris@82	914 Ci[WS(csi, 32)] = Tcf - Tce;
Chris@82	915 {
Chris@82	916 E Tca, Tcc, Tcd, Tcg;
Chris@82	917 Tca = Tc2 + Tc9;
Chris@82	918 Cr[WS(csr, 48)] = FNMS(KP707106781, Tca, TbV);
Chris@82	919 Cr[WS(csr, 16)] = FMA(KP707106781, Tca, TbV);
Chris@82	920 Tcc = Tc9 - Tc2;
Chris@82	921 Ci[WS(csi, 16)] = FMA(KP707106781, Tcc, Tcb);
Chris@82	922 Ci[WS(csi, 48)] = FMS(KP707106781, Tcc, Tcb);
Chris@82	923 Tcd = T11 + T22;
Chris@82	924 Tcg = Tce + Tcf;
Chris@82	925 Cr[WS(csr, 64)] = Tcd - Tcg;
Chris@82	926 Cr[0] = Tcd + Tcg;
Chris@82	927 }
Chris@82	928 }
Chris@82	929 {
Chris@82	930 E Tch, Tct, Tck, Tcu, Tco, Tcy, Tcr, Tcz, Tci, Tcj;
Chris@82	931 Tch = Tf - Tu;
Chris@82	932 Tct = TZ - TK;
Chris@82	933 Tci = T1g - T1v;
Chris@82	934 Tcj = T1L - T20;
Chris@82	935 Tck = Tci + Tcj;
Chris@82	936 Tcu = Tcj - Tci;
Chris@82	937 {
Chris@82	938 E Tcm, Tcn, Tcp, Tcq;
Chris@82	939 Tcm = TbW - TbX;
Chris@82	940 Tcn = Tc0 - TbZ;
Chris@82	941 Tco = FMA(KP414213562, Tcn, Tcm);
Chris@82	942 Tcy = FNMS(KP414213562, Tcm, Tcn);
Chris@82	943 Tcp = Tc3 - Tc4;
Chris@82	944 Tcq = Tc7 - Tc6;
Chris@82	945 Tcr = FNMS(KP414213562, Tcq, Tcp);
Chris@82	946 Tcz = FMA(KP414213562, Tcp, Tcq);
Chris@82	947 }
Chris@82	948 {
Chris@82	949 E Tcl, Tcs, Tcx, TcA;
Chris@82	950 Tcl = FMA(KP707106781, Tck, Tch);
Chris@82	951 Tcs = Tco + Tcr;
Chris@82	952 Cr[WS(csr, 56)] = FNMS(KP923879532, Tcs, Tcl);
Chris@82	953 Cr[WS(csr, 8)] = FMA(KP923879532, Tcs, Tcl);
Chris@82	954 Tcx = FMA(KP707106781, Tcu, Tct);
Chris@82	955 TcA = Tcy + Tcz;
Chris@82	956 Ci[WS(csi, 8)] = FMA(KP923879532, TcA, Tcx);
Chris@82	957 Ci[WS(csi, 56)] = FMS(KP923879532, TcA, Tcx);
Chris@82	958 }
Chris@82	959 {
Chris@82	960 E Tcv, Tcw, TcB, TcC;
Chris@82	961 Tcv = FNMS(KP707106781, Tcu, Tct);
Chris@82	962 Tcw = Tcr - Tco;
Chris@82	963 Ci[WS(csi, 24)] = FMS(KP923879532, Tcw, Tcv);
Chris@82	964 Ci[WS(csi, 40)] = FMA(KP923879532, Tcw, Tcv);
Chris@82	965 TcB = FNMS(KP707106781, Tck, Tch);
Chris@82	966 TcC = Tcz - Tcy;
Chris@82	967 Cr[WS(csr, 40)] = FNMS(KP923879532, TcC, TcB);
Chris@82	968 Cr[WS(csr, 24)] = FMA(KP923879532, TcC, TcB);
Chris@82	969 }
Chris@82	970 }
Chris@82	971 {
Chris@82	972 E Ta9, TbB, Tbp, TbL, Tag, TbM, TbJ, TbR, TaO, Tbw, Tbs, TbC, TbG, TbQ, Tbl;
Chris@82	973 E Tbx, Ta8, Tbo;
Chris@82	974 Ta8 = Ta6 + Ta7;
Chris@82	975 Ta9 = FMA(KP707106781, Ta8, Ta5);
Chris@82	976 TbB = FNMS(KP707106781, Ta8, Ta5);
Chris@82	977 Tbo = Ta7 - Ta6;
Chris@82	978 Tbp = FMA(KP707106781, Tbo, Tbn);
Chris@82	979 TbL = FNMS(KP707106781, Tbo, Tbn);
Chris@82	980 {
Chris@82	981 E Tac, Taf, TbH, TbI;
Chris@82	982 Tac = FMA(KP414213562, Tab, Taa);
Chris@82	983 Taf = FNMS(KP414213562, Tae, Tad);
Chris@82	984 Tag = Tac + Taf;
Chris@82	985 TbM = Taf - Tac;
Chris@82	986 TbH = FNMS(KP707106781, Tba, TaV);
Chris@82	987 TbI = FNMS(KP707106781, Tbj, Tbi);
Chris@82	988 TbJ = FMA(KP668178637, TbI, TbH);
Chris@82	989 TbR = FNMS(KP668178637, TbH, TbI);
Chris@82	990 }
Chris@82	991 {
Chris@82	992 E TaE, TaN, Tbq, Tbr;
Chris@82	993 TaE = FMA(KP707106781, TaD, Tao);
Chris@82	994 TaN = FMA(KP707106781, TaM, TaL);
Chris@82	995 TaO = FMA(KP198912367, TaN, TaE);
Chris@82	996 Tbw = FNMS(KP198912367, TaE, TaN);
Chris@82	997 Tbq = FNMS(KP414213562, Taa, Tab);
Chris@82	998 Tbr = FMA(KP414213562, Tad, Tae);
Chris@82	999 Tbs = Tbq + Tbr;
Chris@82	1000 TbC = Tbr - Tbq;
Chris@82	1001 }
Chris@82	1002 {
Chris@82	1003 E TbE, TbF, Tbb, Tbk;
Chris@82	1004 TbE = FNMS(KP707106781, TaD, Tao);
Chris@82	1005 TbF = FNMS(KP707106781, TaM, TaL);
Chris@82	1006 TbG = FNMS(KP668178637, TbF, TbE);
Chris@82	1007 TbQ = FMA(KP668178637, TbE, TbF);
Chris@82	1008 Tbb = FMA(KP707106781, Tba, TaV);
Chris@82	1009 Tbk = FMA(KP707106781, Tbj, Tbi);
Chris@82	1010 Tbl = FNMS(KP198912367, Tbk, Tbb);
Chris@82	1011 Tbx = FMA(KP198912367, Tbb, Tbk);
Chris@82	1012 }
Chris@82	1013 {
Chris@82	1014 E Tah, Tbm, Tbv, Tby;
Chris@82	1015 Tah = FMA(KP923879532, Tag, Ta9);
Chris@82	1016 Tbm = TaO + Tbl;
Chris@82	1017 Cr[WS(csr, 60)] = FNMS(KP980785280, Tbm, Tah);
Chris@82	1018 Cr[WS(csr, 4)] = FMA(KP980785280, Tbm, Tah);
Chris@82	1019 Tbv = FMA(KP923879532, Tbs, Tbp);
Chris@82	1020 Tby = Tbw + Tbx;
Chris@82	1021 Ci[WS(csi, 4)] = FMA(KP980785280, Tby, Tbv);
Chris@82	1022 Ci[WS(csi, 60)] = FMS(KP980785280, Tby, Tbv);
Chris@82	1023 }
Chris@82	1024 {
Chris@82	1025 E Tbt, Tbu, Tbz, TbA;
Chris@82	1026 Tbt = FNMS(KP923879532, Tbs, Tbp);
Chris@82	1027 Tbu = Tbl - TaO;
Chris@82	1028 Ci[WS(csi, 28)] = FMS(KP980785280, Tbu, Tbt);
Chris@82	1029 Ci[WS(csi, 36)] = FMA(KP980785280, Tbu, Tbt);
Chris@82	1030 Tbz = FNMS(KP923879532, Tag, Ta9);
Chris@82	1031 TbA = Tbx - Tbw;
Chris@82	1032 Cr[WS(csr, 36)] = FNMS(KP980785280, TbA, Tbz);
Chris@82	1033 Cr[WS(csr, 28)] = FMA(KP980785280, TbA, Tbz);
Chris@82	1034 }
Chris@82	1035 {
Chris@82	1036 E TbD, TbK, TbP, TbS;
Chris@82	1037 TbD = FMA(KP923879532, TbC, TbB);
Chris@82	1038 TbK = TbG + TbJ;
Chris@82	1039 Cr[WS(csr, 52)] = FNMS(KP831469612, TbK, TbD);
Chris@82	1040 Cr[WS(csr, 12)] = FMA(KP831469612, TbK, TbD);
Chris@82	1041 TbP = FNMS(KP923879532, TbM, TbL);
Chris@82	1042 TbS = TbQ + TbR;
Chris@82	1043 Ci[WS(csi, 12)] = -(FMA(KP831469612, TbS, TbP));
Chris@82	1044 Ci[WS(csi, 52)] = FNMS(KP831469612, TbS, TbP);
Chris@82	1045 }
Chris@82	1046 {
Chris@82	1047 E TbN, TbO, TbT, TbU;
Chris@82	1048 TbN = FMA(KP923879532, TbM, TbL);
Chris@82	1049 TbO = TbJ - TbG;
Chris@82	1050 Ci[WS(csi, 20)] = FMA(KP831469612, TbO, TbN);
Chris@82	1051 Ci[WS(csi, 44)] = FMS(KP831469612, TbO, TbN);
Chris@82	1052 TbT = FNMS(KP923879532, TbC, TbB);
Chris@82	1053 TbU = TbQ - TbR;
Chris@82	1054 Cr[WS(csr, 44)] = FNMS(KP831469612, TbU, TbT);
Chris@82	1055 Cr[WS(csr, 20)] = FMA(KP831469612, TbU, TbT);
Chris@82	1056 }
Chris@82	1057 }
Chris@82	1058 {
Chris@82	1059 E Tev, Tf7, Tfc, Tfm, Tff, Tfn, TeC, Tfi, TeK, Tf2, TeV, Tfh, TeY, Tf8, TeR;
Chris@82	1060 E Tf3;
Chris@82	1061 {
Chris@82	1062 E Tet, Teu, Tfa, Tfb;
Chris@82	1063 Tet = FNMS(KP707106781, TcG, TcD);
Chris@82	1064 Teu = TdV - TdU;
Chris@82	1065 Tev = FNMS(KP923879532, Teu, Tet);
Chris@82	1066 Tf7 = FMA(KP923879532, Teu, Tet);
Chris@82	1067 Tfa = FMA(KP923879532, TeF, TeE);
Chris@82	1068 Tfb = FMA(KP923879532, TeI, TeH);
Chris@82	1069 Tfc = FNMS(KP303346683, Tfb, Tfa);
Chris@82	1070 Tfm = FMA(KP303346683, Tfa, Tfb);
Chris@82	1071 }
Chris@82	1072 {
Chris@82	1073 E Tfd, Tfe, Tey, TeB;
Chris@82	1074 Tfd = FMA(KP923879532, TeM, TeL);
Chris@82	1075 Tfe = FNMS(KP923879532, TeP, TeO);
Chris@82	1076 Tff = FMA(KP303346683, Tfe, Tfd);
Chris@82	1077 Tfn = FNMS(KP303346683, Tfd, Tfe);
Chris@82	1078 Tey = FMA(KP668178637, Tex, Tew);
Chris@82	1079 TeB = FNMS(KP668178637, TeA, Tez);
Chris@82	1080 TeC = Tey - TeB;
Chris@82	1081 Tfi = Tey + TeB;
Chris@82	1082 }
Chris@82	1083 {
Chris@82	1084 E TeG, TeJ, TeT, TeU;
Chris@82	1085 TeG = FNMS(KP923879532, TeF, TeE);
Chris@82	1086 TeJ = FNMS(KP923879532, TeI, TeH);
Chris@82	1087 TeK = FMA(KP534511135, TeJ, TeG);
Chris@82	1088 Tf2 = FNMS(KP534511135, TeG, TeJ);
Chris@82	1089 TeT = FNMS(KP707106781, TdS, TdR);
Chris@82	1090 TeU = TcN - TcK;
Chris@82	1091 TeV = FMA(KP923879532, TeU, TeT);
Chris@82	1092 Tfh = FNMS(KP923879532, TeU, TeT);
Chris@82	1093 }
Chris@82	1094 {
Chris@82	1095 E TeW, TeX, TeN, TeQ;
Chris@82	1096 TeW = FMA(KP668178637, Tez, TeA);
Chris@82	1097 TeX = FNMS(KP668178637, Tew, Tex);
Chris@82	1098 TeY = TeW - TeX;
Chris@82	1099 Tf8 = TeX + TeW;
Chris@82	1100 TeN = FNMS(KP923879532, TeM, TeL);
Chris@82	1101 TeQ = FMA(KP923879532, TeP, TeO);
Chris@82	1102 TeR = FNMS(KP534511135, TeQ, TeN);
Chris@82	1103 Tf3 = FMA(KP534511135, TeN, TeQ);
Chris@82	1104 }
Chris@82	1105 {
Chris@82	1106 E TeD, TeS, Tf1, Tf4;
Chris@82	1107 TeD = FMA(KP831469612, TeC, Tev);
Chris@82	1108 TeS = TeK + TeR;
Chris@82	1109 Cr[WS(csr, 54)] = FNMS(KP881921264, TeS, TeD);
Chris@82	1110 Cr[WS(csr, 10)] = FMA(KP881921264, TeS, TeD);
Chris@82	1111 Tf1 = FMA(KP831469612, TeY, TeV);
Chris@82	1112 Tf4 = Tf2 + Tf3;
Chris@82	1113 Ci[WS(csi, 10)] = FMA(KP881921264, Tf4, Tf1);
Chris@82	1114 Ci[WS(csi, 54)] = FMS(KP881921264, Tf4, Tf1);
Chris@82	1115 }
Chris@82	1116 {
Chris@82	1117 E TeZ, Tf0, Tf5, Tf6;
Chris@82	1118 TeZ = FNMS(KP831469612, TeY, TeV);
Chris@82	1119 Tf0 = TeR - TeK;
Chris@82	1120 Ci[WS(csi, 22)] = FMS(KP881921264, Tf0, TeZ);
Chris@82	1121 Ci[WS(csi, 42)] = FMA(KP881921264, Tf0, TeZ);
Chris@82	1122 Tf5 = FNMS(KP831469612, TeC, Tev);
Chris@82	1123 Tf6 = Tf3 - Tf2;
Chris@82	1124 Cr[WS(csr, 42)] = FNMS(KP881921264, Tf6, Tf5);
Chris@82	1125 Cr[WS(csr, 22)] = FMA(KP881921264, Tf6, Tf5);
Chris@82	1126 }
Chris@82	1127 {
Chris@82	1128 E Tf9, Tfg, Tfl, Tfo;
Chris@82	1129 Tf9 = FMA(KP831469612, Tf8, Tf7);
Chris@82	1130 Tfg = Tfc + Tff;
Chris@82	1131 Cr[WS(csr, 58)] = FNMS(KP956940335, Tfg, Tf9);
Chris@82	1132 Cr[WS(csr, 6)] = FMA(KP956940335, Tfg, Tf9);
Chris@82	1133 Tfl = FMA(KP831469612, Tfi, Tfh);
Chris@82	1134 Tfo = Tfm + Tfn;
Chris@82	1135 Ci[WS(csi, 6)] = -(FMA(KP956940335, Tfo, Tfl));
Chris@82	1136 Ci[WS(csi, 58)] = FNMS(KP956940335, Tfo, Tfl);
Chris@82	1137 }
Chris@82	1138 {
Chris@82	1139 E Tfj, Tfk, Tfp, Tfq;
Chris@82	1140 Tfj = FNMS(KP831469612, Tfi, Tfh);
Chris@82	1141 Tfk = Tff - Tfc;
Chris@82	1142 Ci[WS(csi, 26)] = FMA(KP956940335, Tfk, Tfj);
Chris@82	1143 Ci[WS(csi, 38)] = FMS(KP956940335, Tfk, Tfj);
Chris@82	1144 Tfp = FNMS(KP831469612, Tf8, Tf7);
Chris@82	1145 Tfq = Tfm - Tfn;
Chris@82	1146 Cr[WS(csr, 38)] = FNMS(KP956940335, Tfq, Tfp);
Chris@82	1147 Cr[WS(csr, 26)] = FMA(KP956940335, Tfq, Tfp);
Chris@82	1148 }
Chris@82	1149 }
Chris@82	1150 {
Chris@82	1151 E TcP, Te9, Tee, Teo, Teh, Tep, Td8, Tek, Tdu, Te4, TdX, Tej, Te0, Tea, TdP;
Chris@82	1152 E Te5;
Chris@82	1153 {
Chris@82	1154 E TcH, TcO, Tec, Ted;
Chris@82	1155 TcH = FMA(KP707106781, TcG, TcD);
Chris@82	1156 TcO = TcK + TcN;
Chris@82	1157 TcP = FMA(KP923879532, TcO, TcH);
Chris@82	1158 Te9 = FNMS(KP923879532, TcO, TcH);
Chris@82	1159 Tec = FNMS(KP923879532, Tdl, Tde);
Chris@82	1160 Ted = FNMS(KP923879532, Tds, Tdp);
Chris@82	1161 Tee = FNMS(KP820678790, Ted, Tec);
Chris@82	1162 Teo = FMA(KP820678790, Tec, Ted);
Chris@82	1163 }
Chris@82	1164 {
Chris@82	1165 E Tef, Teg, TcY, Td7;
Chris@82	1166 Tef = FNMS(KP923879532, TdG, Tdz);
Chris@82	1167 Teg = FNMS(KP923879532, TdN, TdK);
Chris@82	1168 Teh = FMA(KP820678790, Teg, Tef);
Chris@82	1169 Tep = FNMS(KP820678790, Tef, Teg);
Chris@82	1170 TcY = FMA(KP198912367, TcX, TcU);
Chris@82	1171 Td7 = FNMS(KP198912367, Td6, Td3);
Chris@82	1172 Td8 = TcY + Td7;
Chris@82	1173 Tek = Td7 - TcY;
Chris@82	1174 }
Chris@82	1175 {
Chris@82	1176 E Tdm, Tdt, TdT, TdW;
Chris@82	1177 Tdm = FMA(KP923879532, Tdl, Tde);
Chris@82	1178 Tdt = FMA(KP923879532, Tds, Tdp);
Chris@82	1179 Tdu = FMA(KP098491403, Tdt, Tdm);
Chris@82	1180 Te4 = FNMS(KP098491403, Tdm, Tdt);
Chris@82	1181 TdT = FMA(KP707106781, TdS, TdR);
Chris@82	1182 TdW = TdU + TdV;
Chris@82	1183 TdX = FMA(KP923879532, TdW, TdT);
Chris@82	1184 Tej = FNMS(KP923879532, TdW, TdT);
Chris@82	1185 }
Chris@82	1186 {
Chris@82	1187 E TdY, TdZ, TdH, TdO;
Chris@82	1188 TdY = FNMS(KP198912367, TcU, TcX);
Chris@82	1189 TdZ = FMA(KP198912367, Td3, Td6);
Chris@82	1190 Te0 = TdY + TdZ;
Chris@82	1191 Tea = TdZ - TdY;
Chris@82	1192 TdH = FMA(KP923879532, TdG, Tdz);
Chris@82	1193 TdO = FMA(KP923879532, TdN, TdK);
Chris@82	1194 TdP = FNMS(KP098491403, TdO, TdH);
Chris@82	1195 Te5 = FMA(KP098491403, TdH, TdO);
Chris@82	1196 }
Chris@82	1197 {
Chris@82	1198 E Td9, TdQ, Te3, Te6;
Chris@82	1199 Td9 = FMA(KP980785280, Td8, TcP);
Chris@82	1200 TdQ = Tdu + TdP;
Chris@82	1201 Cr[WS(csr, 62)] = FNMS(KP995184726, TdQ, Td9);
Chris@82	1202 Cr[WS(csr, 2)] = FMA(KP995184726, TdQ, Td9);
Chris@82	1203 Te3 = FMA(KP980785280, Te0, TdX);
Chris@82	1204 Te6 = Te4 + Te5;
Chris@82	1205 Ci[WS(csi, 2)] = FMA(KP995184726, Te6, Te3);
Chris@82	1206 Ci[WS(csi, 62)] = FMS(KP995184726, Te6, Te3);
Chris@82	1207 }
Chris@82	1208 {
Chris@82	1209 E Te1, Te2, Te7, Te8;
Chris@82	1210 Te1 = FNMS(KP980785280, Te0, TdX);
Chris@82	1211 Te2 = TdP - Tdu;
Chris@82	1212 Ci[WS(csi, 30)] = FMS(KP995184726, Te2, Te1);
Chris@82	1213 Ci[WS(csi, 34)] = FMA(KP995184726, Te2, Te1);
Chris@82	1214 Te7 = FNMS(KP980785280, Td8, TcP);
Chris@82	1215 Te8 = Te5 - Te4;
Chris@82	1216 Cr[WS(csr, 34)] = FNMS(KP995184726, Te8, Te7);
Chris@82	1217 Cr[WS(csr, 30)] = FMA(KP995184726, Te8, Te7);
Chris@82	1218 }
Chris@82	1219 {
Chris@82	1220 E Teb, Tei, Ten, Teq;
Chris@82	1221 Teb = FMA(KP980785280, Tea, Te9);
Chris@82	1222 Tei = Tee + Teh;
Chris@82	1223 Cr[WS(csr, 50)] = FNMS(KP773010453, Tei, Teb);
Chris@82	1224 Cr[WS(csr, 14)] = FMA(KP773010453, Tei, Teb);
Chris@82	1225 Ten = FNMS(KP980785280, Tek, Tej);
Chris@82	1226 Teq = Teo + Tep;
Chris@82	1227 Ci[WS(csi, 14)] = -(FMA(KP773010453, Teq, Ten));
Chris@82	1228 Ci[WS(csi, 50)] = FNMS(KP773010453, Teq, Ten);
Chris@82	1229 }
Chris@82	1230 {
Chris@82	1231 E Tel, Tem, Ter, Tes;
Chris@82	1232 Tel = FMA(KP980785280, Tek, Tej);
Chris@82	1233 Tem = Teh - Tee;
Chris@82	1234 Ci[WS(csi, 18)] = FMA(KP773010453, Tem, Tel);
Chris@82	1235 Ci[WS(csi, 46)] = FMS(KP773010453, Tem, Tel);
Chris@82	1236 Ter = FNMS(KP980785280, Tea, Te9);
Chris@82	1237 Tes = Teo - Tep;
Chris@82	1238 Cr[WS(csr, 46)] = FNMS(KP773010453, Tes, Ter);
Chris@82	1239 Cr[WS(csr, 18)] = FMA(KP773010453, Tes, Ter);
Chris@82	1240 }
Chris@82	1241 }
Chris@82	1242 {
Chris@82	1243 E T6v, T77, T6C, T7i, T6Y, T78, T6V, T7h, T6R, T7n, T72, T7f, T6K, T7m, T73;
Chris@82	1244 E T7c;
Chris@82	1245 {
Chris@82	1246 E T6t, T6u, T6T, T6U;
Chris@82	1247 T6t = FNMS(KP923879532, T2e, T27);
Chris@82	1248 T6u = T5U - T5V;
Chris@82	1249 T6v = FNMS(KP980785280, T6u, T6t);
Chris@82	1250 T77 = FMA(KP980785280, T6u, T6t);
Chris@82	1251 {
Chris@82	1252 E T6y, T6B, T6W, T6X;
Chris@82	1253 T6y = FMA(KP820678790, T6x, T6w);
Chris@82	1254 T6B = FNMS(KP820678790, T6A, T6z);
Chris@82	1255 T6C = T6y - T6B;
Chris@82	1256 T7i = T6B + T6y;
Chris@82	1257 T6W = FNMS(KP820678790, T6w, T6x);
Chris@82	1258 T6X = FMA(KP820678790, T6z, T6A);
Chris@82	1259 T6Y = T6W - T6X;
Chris@82	1260 T78 = T6X + T6W;
Chris@82	1261 }
Chris@82	1262 T6T = FNMS(KP923879532, T5S, T5P);
Chris@82	1263 T6U = T2x - T2o;
Chris@82	1264 T6V = FNMS(KP980785280, T6U, T6T);
Chris@82	1265 T7h = FMA(KP980785280, T6U, T6T);
Chris@82	1266 {
Chris@82	1267 E T6N, T7d, T6Q, T7e, T6M, T6P;
Chris@82	1268 T6M = T5I - T5H;
Chris@82	1269 T6N = FNMS(KP980785280, T6M, T6L);
Chris@82	1270 T7d = FMA(KP980785280, T6M, T6L);
Chris@82	1271 T6P = T5v - T5e;
Chris@82	1272 T6Q = FMA(KP980785280, T6P, T6O);
Chris@82	1273 T7e = FNMS(KP980785280, T6P, T6O);
Chris@82	1274 T6R = FNMS(KP472964775, T6Q, T6N);
Chris@82	1275 T7n = FNMS(KP357805721, T7d, T7e);
Chris@82	1276 T72 = FMA(KP472964775, T6N, T6Q);
Chris@82	1277 T7f = FMA(KP357805721, T7e, T7d);
Chris@82	1278 }
Chris@82	1279 {
Chris@82	1280 E T6G, T7a, T6J, T7b, T6F, T6I;
Chris@82	1281 T6F = T4s - T4t;
Chris@82	1282 T6G = FNMS(KP980785280, T6F, T6E);
Chris@82	1283 T7a = FMA(KP980785280, T6F, T6E);
Chris@82	1284 T6I = T4g - T3Z;
Chris@82	1285 T6J = FNMS(KP980785280, T6I, T6H);
Chris@82	1286 T7b = FMA(KP980785280, T6I, T6H);
Chris@82	1287 T6K = FNMS(KP472964775, T6J, T6G);
Chris@82	1288 T7m = FNMS(KP357805721, T7a, T7b);
Chris@82	1289 T73 = FMA(KP472964775, T6G, T6J);
Chris@82	1290 T7c = FMA(KP357805721, T7b, T7a);
Chris@82	1291 }
Chris@82	1292 }
Chris@82	1293 {
Chris@82	1294 E T6D, T6S, T71, T74;
Chris@82	1295 T6D = FMA(KP773010453, T6C, T6v);
Chris@82	1296 T6S = T6K + T6R;
Chris@82	1297 Cr[WS(csr, 55)] = FNMS(KP903989293, T6S, T6D);
Chris@82	1298 Cr[WS(csr, 9)] = FMA(KP903989293, T6S, T6D);
Chris@82	1299 T71 = FNMS(KP773010453, T6Y, T6V);
Chris@82	1300 T74 = T72 - T73;
Chris@82	1301 Ci[WS(csi, 9)] = FMS(KP903989293, T74, T71);
Chris@82	1302 Ci[WS(csi, 55)] = FMA(KP903989293, T74, T71);
Chris@82	1303 }
Chris@82	1304 {
Chris@82	1305 E T6Z, T70, T75, T76;
Chris@82	1306 T6Z = FMA(KP773010453, T6Y, T6V);
Chris@82	1307 T70 = T6R - T6K;
Chris@82	1308 Ci[WS(csi, 23)] = FMA(KP903989293, T70, T6Z);
Chris@82	1309 Ci[WS(csi, 41)] = FMS(KP903989293, T70, T6Z);
Chris@82	1310 T75 = FNMS(KP773010453, T6C, T6v);
Chris@82	1311 T76 = T73 + T72;
Chris@82	1312 Cr[WS(csr, 41)] = FNMS(KP903989293, T76, T75);
Chris@82	1313 Cr[WS(csr, 23)] = FMA(KP903989293, T76, T75);
Chris@82	1314 }
Chris@82	1315 {
Chris@82	1316 E T79, T7g, T7l, T7o;
Chris@82	1317 T79 = FMA(KP773010453, T78, T77);
Chris@82	1318 T7g = T7c + T7f;
Chris@82	1319 Cr[WS(csr, 57)] = FNMS(KP941544065, T7g, T79);
Chris@82	1320 Cr[WS(csr, 7)] = FMA(KP941544065, T7g, T79);
Chris@82	1321 T7l = FMA(KP773010453, T7i, T7h);
Chris@82	1322 T7o = T7m - T7n;
Chris@82	1323 Ci[WS(csi, 7)] = FMA(KP941544065, T7o, T7l);
Chris@82	1324 Ci[WS(csi, 57)] = FMS(KP941544065, T7o, T7l);
Chris@82	1325 }
Chris@82	1326 {
Chris@82	1327 E T7j, T7k, T7p, T7q;
Chris@82	1328 T7j = FNMS(KP773010453, T7i, T7h);
Chris@82	1329 T7k = T7f - T7c;
Chris@82	1330 Ci[WS(csi, 25)] = FMS(KP941544065, T7k, T7j);
Chris@82	1331 Ci[WS(csi, 39)] = FMA(KP941544065, T7k, T7j);
Chris@82	1332 T7p = FNMS(KP773010453, T78, T77);
Chris@82	1333 T7q = T7m + T7n;
Chris@82	1334 Cr[WS(csr, 39)] = FMA(KP941544065, T7q, T7p);
Chris@82	1335 Cr[WS(csr, 25)] = FNMS(KP941544065, T7q, T7p);
Chris@82	1336 }
Chris@82	1337 }
Chris@82	1338 {
Chris@82	1339 E T99, T9L, T9g, T9W, T9C, T9M, T9z, T9V, T9v, Ta0, T9H, T9T, T9o, Ta1, T9G;
Chris@82	1340 E T9Q;
Chris@82	1341 {
Chris@82	1342 E T97, T98, T9x, T9y;
Chris@82	1343 T97 = FNMS(KP923879532, T7s, T7r);
Chris@82	1344 T98 = T8z - T8y;
Chris@82	1345 T99 = FNMS(KP831469612, T98, T97);
Chris@82	1346 T9L = FMA(KP831469612, T98, T97);
Chris@82	1347 {
Chris@82	1348 E T9c, T9f, T9A, T9B;
Chris@82	1349 T9c = FMA(KP534511135, T9b, T9a);
Chris@82	1350 T9f = FNMS(KP534511135, T9e, T9d);
Chris@82	1351 T9g = T9c - T9f;
Chris@82	1352 T9W = T9c + T9f;
Chris@82	1353 T9A = FMA(KP534511135, T9d, T9e);
Chris@82	1354 T9B = FNMS(KP534511135, T9a, T9b);
Chris@82	1355 T9C = T9A - T9B;
Chris@82	1356 T9M = T9B + T9A;
Chris@82	1357 }
Chris@82	1358 T9x = FNMS(KP923879532, T8w, T8v);
Chris@82	1359 T9y = T7z - T7w;
Chris@82	1360 T9z = FMA(KP831469612, T9y, T9x);
Chris@82	1361 T9V = FNMS(KP831469612, T9y, T9x);
Chris@82	1362 {
Chris@82	1363 E T9r, T9R, T9u, T9S, T9q, T9t;
Chris@82	1364 T9q = T8p - T8q;
Chris@82	1365 T9r = FNMS(KP831469612, T9q, T9p);
Chris@82	1366 T9R = FMA(KP831469612, T9q, T9p);
Chris@82	1367 T9t = T8j - T8g;
Chris@82	1368 T9u = FNMS(KP831469612, T9t, T9s);
Chris@82	1369 T9S = FMA(KP831469612, T9t, T9s);
Chris@82	1370 T9v = FMA(KP599376933, T9u, T9r);
Chris@82	1371 Ta0 = FMA(KP250486960, T9R, T9S);
Chris@82	1372 T9H = FNMS(KP599376933, T9r, T9u);
Chris@82	1373 T9T = FNMS(KP250486960, T9S, T9R);
Chris@82	1374 }
Chris@82	1375 {
Chris@82	1376 E T9k, T9O, T9n, T9P, T9j, T9m;
Chris@82	1377 T9j = T87 - T86;
Chris@82	1378 T9k = FNMS(KP831469612, T9j, T9i);
Chris@82	1379 T9O = FMA(KP831469612, T9j, T9i);
Chris@82	1380 T9m = T7X - T80;
Chris@82	1381 T9n = FNMS(KP831469612, T9m, T9l);
Chris@82	1382 T9P = FMA(KP831469612, T9m, T9l);
Chris@82	1383 T9o = FMA(KP599376933, T9n, T9k);
Chris@82	1384 Ta1 = FMA(KP250486960, T9O, T9P);
Chris@82	1385 T9G = FNMS(KP599376933, T9k, T9n);
Chris@82	1386 T9Q = FNMS(KP250486960, T9P, T9O);
Chris@82	1387 }
Chris@82	1388 }
Chris@82	1389 {
Chris@82	1390 E T9h, T9w, T9F, T9I;
Chris@82	1391 T9h = FMA(KP881921264, T9g, T99);
Chris@82	1392 T9w = T9o + T9v;
Chris@82	1393 Cr[WS(csr, 53)] = FNMS(KP857728610, T9w, T9h);
Chris@82	1394 Cr[WS(csr, 11)] = FMA(KP857728610, T9w, T9h);
Chris@82	1395 T9F = FMA(KP881921264, T9C, T9z);
Chris@82	1396 T9I = T9G - T9H;
Chris@82	1397 Ci[WS(csi, 11)] = FMA(KP857728610, T9I, T9F);
Chris@82	1398 Ci[WS(csi, 53)] = FMS(KP857728610, T9I, T9F);
Chris@82	1399 }
Chris@82	1400 {
Chris@82	1401 E T9D, T9E, T9J, T9K;
Chris@82	1402 T9D = FNMS(KP881921264, T9C, T9z);
Chris@82	1403 T9E = T9v - T9o;
Chris@82	1404 Ci[WS(csi, 21)] = FMS(KP857728610, T9E, T9D);
Chris@82	1405 Ci[WS(csi, 43)] = FMA(KP857728610, T9E, T9D);
Chris@82	1406 T9J = FNMS(KP881921264, T9g, T99);
Chris@82	1407 T9K = T9G + T9H;
Chris@82	1408 Cr[WS(csr, 43)] = FMA(KP857728610, T9K, T9J);
Chris@82	1409 Cr[WS(csr, 21)] = FNMS(KP857728610, T9K, T9J);
Chris@82	1410 }
Chris@82	1411 {
Chris@82	1412 E T9N, T9U, T9Z, Ta2;
Chris@82	1413 T9N = FMA(KP881921264, T9M, T9L);
Chris@82	1414 T9U = T9Q + T9T;
Chris@82	1415 Cr[WS(csr, 59)] = FNMS(KP970031253, T9U, T9N);
Chris@82	1416 Cr[WS(csr, 5)] = FMA(KP970031253, T9U, T9N);
Chris@82	1417 T9Z = FMA(KP881921264, T9W, T9V);
Chris@82	1418 Ta2 = Ta0 - Ta1;
Chris@82	1419 Ci[WS(csi, 5)] = FMS(KP970031253, Ta2, T9Z);
Chris@82	1420 Ci[WS(csi, 59)] = FMA(KP970031253, Ta2, T9Z);
Chris@82	1421 }
Chris@82	1422 {
Chris@82	1423 E T9X, T9Y, Ta3, Ta4;
Chris@82	1424 T9X = FNMS(KP881921264, T9W, T9V);
Chris@82	1425 T9Y = T9T - T9Q;
Chris@82	1426 Ci[WS(csi, 27)] = FMA(KP970031253, T9Y, T9X);
Chris@82	1427 Ci[WS(csi, 37)] = FMS(KP970031253, T9Y, T9X);
Chris@82	1428 Ta3 = FNMS(KP881921264, T9M, T9L);
Chris@82	1429 Ta4 = Ta1 + Ta0;
Chris@82	1430 Cr[WS(csr, 37)] = FNMS(KP970031253, Ta4, Ta3);
Chris@82	1431 Cr[WS(csr, 27)] = FMA(KP970031253, Ta4, Ta3);
Chris@82	1432 }
Chris@82	1433 }
Chris@82	1434 {
Chris@82	1435 E T2z, T69, T3g, T6k, T60, T6a, T5X, T6j, T5L, T6p, T64, T6h, T4w, T6o, T65;
Chris@82	1436 E T6e;
Chris@82	1437 {
Chris@82	1438 E T2f, T2y, T5T, T5W;
Chris@82	1439 T2f = FMA(KP923879532, T2e, T27);
Chris@82	1440 T2y = T2o + T2x;
Chris@82	1441 T2z = FMA(KP980785280, T2y, T2f);
Chris@82	1442 T69 = FNMS(KP980785280, T2y, T2f);
Chris@82	1443 {
Chris@82	1444 E T2U, T3f, T5Y, T5Z;
Chris@82	1445 T2U = FNMS(KP098491403, T2T, T2M);
Chris@82	1446 T3f = FMA(KP098491403, T3e, T37);
Chris@82	1447 T3g = T2U + T3f;
Chris@82	1448 T6k = T3f - T2U;
Chris@82	1449 T5Y = FMA(KP098491403, T2M, T2T);
Chris@82	1450 T5Z = FNMS(KP098491403, T37, T3e);
Chris@82	1451 T60 = T5Y + T5Z;
Chris@82	1452 T6a = T5Y - T5Z;
Chris@82	1453 }
Chris@82	1454 T5T = FMA(KP923879532, T5S, T5P);
Chris@82	1455 T5W = T5U + T5V;
Chris@82	1456 T5X = FMA(KP980785280, T5W, T5T);
Chris@82	1457 T6j = FNMS(KP980785280, T5W, T5T);
Chris@82	1458 {
Chris@82	1459 E T5x, T6f, T5K, T6g, T5w, T5J;
Chris@82	1460 T5w = T5e + T5v;
Chris@82	1461 T5x = FMA(KP980785280, T5w, T4X);
Chris@82	1462 T6f = FNMS(KP980785280, T5w, T4X);
Chris@82	1463 T5J = T5H + T5I;
Chris@82	1464 T5K = FMA(KP980785280, T5J, T5G);
Chris@82	1465 T6g = FNMS(KP980785280, T5J, T5G);
Chris@82	1466 T5L = FNMS(KP049126849, T5K, T5x);
Chris@82	1467 T6p = FNMS(KP906347169, T6f, T6g);
Chris@82	1468 T64 = FMA(KP049126849, T5x, T5K);
Chris@82	1469 T6h = FMA(KP906347169, T6g, T6f);
Chris@82	1470 }
Chris@82	1471 {
Chris@82	1472 E T4i, T6c, T4v, T6d, T4h, T4u;
Chris@82	1473 T4h = T3Z + T4g;
Chris@82	1474 T4i = FMA(KP980785280, T4h, T3I);
Chris@82	1475 T6c = FNMS(KP980785280, T4h, T3I);
Chris@82	1476 T4u = T4s + T4t;
Chris@82	1477 T4v = FMA(KP980785280, T4u, T4r);
Chris@82	1478 T6d = FNMS(KP980785280, T4u, T4r);
Chris@82	1479 T4w = FNMS(KP049126849, T4v, T4i);
Chris@82	1480 T6o = FNMS(KP906347169, T6c, T6d);
Chris@82	1481 T65 = FMA(KP049126849, T4i, T4v);
Chris@82	1482 T6e = FMA(KP906347169, T6d, T6c);
Chris@82	1483 }
Chris@82	1484 }
Chris@82	1485 {
Chris@82	1486 E T3h, T5M, T63, T66;
Chris@82	1487 T3h = FMA(KP995184726, T3g, T2z);
Chris@82	1488 T5M = T4w + T5L;
Chris@82	1489 Cr[WS(csr, 63)] = FNMS(KP998795456, T5M, T3h);
Chris@82	1490 Cr[WS(csr, 1)] = FMA(KP998795456, T5M, T3h);
Chris@82	1491 T63 = FMA(KP995184726, T60, T5X);
Chris@82	1492 T66 = T64 - T65;
Chris@82	1493 Ci[WS(csi, 1)] = FMS(KP998795456, T66, T63);
Chris@82	1494 Ci[WS(csi, 63)] = FMA(KP998795456, T66, T63);
Chris@82	1495 }
Chris@82	1496 {
Chris@82	1497 E T61, T62, T67, T68;
Chris@82	1498 T61 = FNMS(KP995184726, T60, T5X);
Chris@82	1499 T62 = T5L - T4w;
Chris@82	1500 Ci[WS(csi, 31)] = FMA(KP998795456, T62, T61);
Chris@82	1501 Ci[WS(csi, 33)] = FMS(KP998795456, T62, T61);
Chris@82	1502 T67 = FNMS(KP995184726, T3g, T2z);
Chris@82	1503 T68 = T65 + T64;
Chris@82	1504 Cr[WS(csr, 33)] = FNMS(KP998795456, T68, T67);
Chris@82	1505 Cr[WS(csr, 31)] = FMA(KP998795456, T68, T67);
Chris@82	1506 }
Chris@82	1507 {
Chris@82	1508 E T6b, T6i, T6n, T6q;
Chris@82	1509 T6b = FMA(KP995184726, T6a, T69);
Chris@82	1510 T6i = T6e + T6h;
Chris@82	1511 Cr[WS(csr, 49)] = FNMS(KP740951125, T6i, T6b);
Chris@82	1512 Cr[WS(csr, 15)] = FMA(KP740951125, T6i, T6b);
Chris@82	1513 T6n = FMA(KP995184726, T6k, T6j);
Chris@82	1514 T6q = T6o - T6p;
Chris@82	1515 Ci[WS(csi, 15)] = FMA(KP740951125, T6q, T6n);
Chris@82	1516 Ci[WS(csi, 49)] = FMS(KP740951125, T6q, T6n);
Chris@82	1517 }
Chris@82	1518 {
Chris@82	1519 E T6l, T6m, T6r, T6s;
Chris@82	1520 T6l = FNMS(KP995184726, T6k, T6j);
Chris@82	1521 T6m = T6h - T6e;
Chris@82	1522 Ci[WS(csi, 17)] = FMS(KP740951125, T6m, T6l);
Chris@82	1523 Ci[WS(csi, 47)] = FMA(KP740951125, T6m, T6l);
Chris@82	1524 T6r = FNMS(KP995184726, T6a, T69);
Chris@82	1525 T6s = T6o + T6p;
Chris@82	1526 Cr[WS(csr, 47)] = FMA(KP740951125, T6s, T6r);
Chris@82	1527 Cr[WS(csr, 17)] = FNMS(KP740951125, T6s, T6r);
Chris@82	1528 }
Chris@82	1529 }
Chris@82	1530 {
Chris@82	1531 E T7B, T8N, T7Q, T8Y, T8E, T8O, T8B, T8X, T8t, T92, T8J, T8V, T8a, T93, T8I;
Chris@82	1532 E T8S;
Chris@82	1533 {
Chris@82	1534 E T7t, T7A, T8x, T8A;
Chris@82	1535 T7t = FMA(KP923879532, T7s, T7r);
Chris@82	1536 T7A = T7w + T7z;
Chris@82	1537 T7B = FMA(KP831469612, T7A, T7t);
Chris@82	1538 T8N = FNMS(KP831469612, T7A, T7t);
Chris@82	1539 {
Chris@82	1540 E T7I, T7P, T8C, T8D;
Chris@82	1541 T7I = FMA(KP303346683, T7H, T7E);
Chris@82	1542 T7P = FNMS(KP303346683, T7O, T7L);
Chris@82	1543 T7Q = T7I + T7P;
Chris@82	1544 T8Y = T7P - T7I;
Chris@82	1545 T8C = FNMS(KP303346683, T7E, T7H);
Chris@82	1546 T8D = FMA(KP303346683, T7L, T7O);
Chris@82	1547 T8E = T8C + T8D;
Chris@82	1548 T8O = T8D - T8C;
Chris@82	1549 }
Chris@82	1550 T8x = FMA(KP923879532, T8w, T8v);
Chris@82	1551 T8A = T8y + T8z;
Chris@82	1552 T8B = FMA(KP831469612, T8A, T8x);
Chris@82	1553 T8X = FNMS(KP831469612, T8A, T8x);
Chris@82	1554 {
Chris@82	1555 E T8l, T8T, T8s, T8U, T8k, T8r;
Chris@82	1556 T8k = T8g + T8j;
Chris@82	1557 T8l = FMA(KP831469612, T8k, T8d);
Chris@82	1558 T8T = FNMS(KP831469612, T8k, T8d);
Chris@82	1559 T8r = T8p + T8q;
Chris@82	1560 T8s = FMA(KP831469612, T8r, T8o);
Chris@82	1561 T8U = FNMS(KP831469612, T8r, T8o);
Chris@82	1562 T8t = FMA(KP148335987, T8s, T8l);
Chris@82	1563 T92 = FMA(KP741650546, T8T, T8U);
Chris@82	1564 T8J = FNMS(KP148335987, T8l, T8s);
Chris@82	1565 T8V = FNMS(KP741650546, T8U, T8T);
Chris@82	1566 }
Chris@82	1567 {
Chris@82	1568 E T82, T8Q, T89, T8R, T81, T88;
Chris@82	1569 T81 = T7X + T80;
Chris@82	1570 T82 = FMA(KP831469612, T81, T7U);
Chris@82	1571 T8Q = FNMS(KP831469612, T81, T7U);
Chris@82	1572 T88 = T86 + T87;
Chris@82	1573 T89 = FMA(KP831469612, T88, T85);
Chris@82	1574 T8R = FNMS(KP831469612, T88, T85);
Chris@82	1575 T8a = FMA(KP148335987, T89, T82);
Chris@82	1576 T93 = FMA(KP741650546, T8Q, T8R);
Chris@82	1577 T8I = FNMS(KP148335987, T82, T89);
Chris@82	1578 T8S = FNMS(KP741650546, T8R, T8Q);
Chris@82	1579 }
Chris@82	1580 }
Chris@82	1581 {
Chris@82	1582 E T7R, T8u, T8H, T8K;
Chris@82	1583 T7R = FMA(KP956940335, T7Q, T7B);
Chris@82	1584 T8u = T8a + T8t;
Chris@82	1585 Cr[WS(csr, 61)] = FNMS(KP989176509, T8u, T7R);
Chris@82	1586 Cr[WS(csr, 3)] = FMA(KP989176509, T8u, T7R);
Chris@82	1587 T8H = FMA(KP956940335, T8E, T8B);
Chris@82	1588 T8K = T8I - T8J;
Chris@82	1589 Ci[WS(csi, 3)] = FMA(KP989176509, T8K, T8H);
Chris@82	1590 Ci[WS(csi, 61)] = FMS(KP989176509, T8K, T8H);
Chris@82	1591 }
Chris@82	1592 {
Chris@82	1593 E T8F, T8G, T8L, T8M;
Chris@82	1594 T8F = FNMS(KP956940335, T8E, T8B);
Chris@82	1595 T8G = T8t - T8a;
Chris@82	1596 Ci[WS(csi, 29)] = FMS(KP989176509, T8G, T8F);
Chris@82	1597 Ci[WS(csi, 35)] = FMA(KP989176509, T8G, T8F);
Chris@82	1598 T8L = FNMS(KP956940335, T7Q, T7B);
Chris@82	1599 T8M = T8I + T8J;
Chris@82	1600 Cr[WS(csr, 35)] = FMA(KP989176509, T8M, T8L);
Chris@82	1601 Cr[WS(csr, 29)] = FNMS(KP989176509, T8M, T8L);
Chris@82	1602 }
Chris@82	1603 {
Chris@82	1604 E T8P, T8W, T91, T94;
Chris@82	1605 T8P = FMA(KP956940335, T8O, T8N);
Chris@82	1606 T8W = T8S + T8V;
Chris@82	1607 Cr[WS(csr, 51)] = FNMS(KP803207531, T8W, T8P);
Chris@82	1608 Cr[WS(csr, 13)] = FMA(KP803207531, T8W, T8P);
Chris@82	1609 T91 = FNMS(KP956940335, T8Y, T8X);
Chris@82	1610 T94 = T92 - T93;
Chris@82	1611 Ci[WS(csi, 13)] = FMS(KP803207531, T94, T91);
Chris@82	1612 Ci[WS(csi, 51)] = FMA(KP803207531, T94, T91);
Chris@82	1613 }
Chris@82	1614 {
Chris@82	1615 E T8Z, T90, T95, T96;
Chris@82	1616 T8Z = FMA(KP956940335, T8Y, T8X);
Chris@82	1617 T90 = T8V - T8S;
Chris@82	1618 Ci[WS(csi, 19)] = FMA(KP803207531, T90, T8Z);
Chris@82	1619 Ci[WS(csi, 45)] = FMS(KP803207531, T90, T8Z);
Chris@82	1620 T95 = FNMS(KP956940335, T8O, T8N);
Chris@82	1621 T96 = T93 + T92;
Chris@82	1622 Cr[WS(csr, 45)] = FNMS(KP803207531, T96, T95);
Chris@82	1623 Cr[WS(csr, 19)] = FMA(KP803207531, T96, T95);
Chris@82	1624 }
Chris@82	1625 }
Chris@82	1626 }
Chris@82	1627 }
Chris@82	1628 }
Chris@82	1629
Chris@82	1630 static const kr2c_desc desc = { 128, "r2cf_128", {440, 0, 516, 0}, &GENUS };
Chris@82	1631
Chris@82	1632 void X(codelet_r2cf_128) (planner *p) {
Chris@82	1633 X(kr2c_register) (p, r2cf_128, &desc);
Chris@82	1634 }
Chris@82	1635
Chris@82	1636 #else
Chris@82	1637
Chris@82	1638 /* Generated by: ../../../genfft/gen_r2cf.native -compact -variables 4 -pipeline-latency 4 -n 128 -name r2cf_128 -include rdft/scalar/r2cf.h */
Chris@82	1639
Chris@82	1640 /*
Chris@82	1641 * This function contains 956 FP additions, 330 FP multiplications,
Chris@82	1642 * (or, 812 additions, 186 multiplications, 144 fused multiply/add),
Chris@82	1643 * 186 stack variables, 31 constants, and 256 memory accesses
Chris@82	1644 */
Chris@82	1645 #include "rdft/scalar/r2cf.h"
Chris@82	1646
Chris@82	1647 static void r2cf_128(R R0, R R1, R Cr, R Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
Chris@82	1648 {
Chris@82	1649 DK(KP803207531, +0.803207531480644909806676512963141923879569427);
Chris@82	1650 DK(KP595699304, +0.595699304492433343467036528829969889511926338);
Chris@82	1651 DK(KP146730474, +0.146730474455361751658850129646717819706215317);
Chris@82	1652 DK(KP989176509, +0.989176509964780973451673738016243063983689533);
Chris@82	1653 DK(KP740951125, +0.740951125354959091175616897495162729728955309);
Chris@82	1654 DK(KP671558954, +0.671558954847018400625376850427421803228750632);
Chris@82	1655 DK(KP049067674, +0.049067674327418014254954976942682658314745363);
Chris@82	1656 DK(KP998795456, +0.998795456205172392714771604759100694443203615);
Chris@82	1657 DK(KP242980179, +0.242980179903263889948274162077471118320990783);
Chris@82	1658 DK(KP970031253, +0.970031253194543992603984207286100251456865962);
Chris@82	1659 DK(KP514102744, +0.514102744193221726593693838968815772608049120);
Chris@82	1660 DK(KP857728610, +0.857728610000272069902269984284770137042490799);
Chris@82	1661 DK(KP336889853, +0.336889853392220050689253212619147570477766780);
Chris@82	1662 DK(KP941544065, +0.941544065183020778412509402599502357185589796);
Chris@82	1663 DK(KP427555093, +0.427555093430282094320966856888798534304578629);
Chris@82	1664 DK(KP903989293, +0.903989293123443331586200297230537048710132025);
Chris@82	1665 DK(KP098017140, +0.098017140329560601994195563888641845861136673);
Chris@82	1666 DK(KP995184726, +0.995184726672196886244836953109479921575474869);
Chris@82	1667 DK(KP634393284, +0.634393284163645498215171613225493370675687095);
Chris@82	1668 DK(KP773010453, +0.773010453362736960810906609758469800971041293);
Chris@82	1669 DK(KP881921264, +0.881921264348355029712756863660388349508442621);
Chris@82	1670 DK(KP471396736, +0.471396736825997648556387625905254377657460319);
Chris@82	1671 DK(KP956940335, +0.956940335732208864935797886980269969482849206);
Chris@82	1672 DK(KP290284677, +0.290284677254462367636192375817395274691476278);
Chris@82	1673 DK(KP555570233, +0.555570233019602224742830813948532874374937191);
Chris@82	1674 DK(KP831469612, +0.831469612302545237078788377617905756738560812);
Chris@82	1675 DK(KP195090322, +0.195090322016128267848284868477022240927691618);
Chris@82	1676 DK(KP980785280, +0.980785280403230449126182236134239036973933731);
Chris@82	1677 DK(KP382683432, +0.382683432365089771728459984030398866761344562);
Chris@82	1678 DK(KP923879532, +0.923879532511286756128183189396788286822416626);
Chris@82	1679 DK(KP707106781, +0.707106781186547524400844362104849039284835938);
Chris@82	1680 {
Chris@82	1681 INT i;
Chris@82	1682 for (i = v; i > 0; i = i - 1, R0 = R0 + ivs, R1 = R1 + ivs, Cr = Cr + ovs, Ci = Ci + ovs, MAKE_VOLATILE_STRIDE(512, rs), MAKE_VOLATILE_STRIDE(512, csr), MAKE_VOLATILE_STRIDE(512, csi)) {
Chris@82	1683 E TcD, TdU, T27, T7r, T5S, T8y, Tf, Ta5, Tu, Tbq, TcG, TdV, T2e, T8z, T5V;
Chris@82	1684 E T7s, TK, Ta6, TcK, TdX, T2o, T5X, T7w, T8B, TZ, Ta7, TcN, TdY, T2x, T5Y;
Chris@82	1685 E T7z, T8C, T1g, Taa, TcU, TeA, TcX, Tez, T1v, Tab, T2M, T6z, T7E, T9e, T7H;
Chris@82	1686 E T9d, T2T, T6A, T4X, T6L, Tdz, TeL, TdK, TeP, T5G, T6P, T8d, T9p, TaV, Tc3;
Chris@82	1687 E Tbi, Tc4, T8o, T9t, T3I, T6H, Tde, TeH, Tdp, TeF, T4r, T6F, T7U, T9l, Tao;
Chris@82	1688 E TbW, TaL, TbX, T85, T9j, T1L, Tad, Td3, Tew, Td6, Tex, T20, Tae, T37, T6x;
Chris@82	1689 E T7L, T9a, T7O, T9b, T3e, T6w, TbZ, Tc0, T3Z, T4s, Tds, TeI, T4g, T4t, T80;
Chris@82	1690 E T87, Tdl, TeE, T7X, T86, TaD, TaM, Tc6, Tc7, T5e, T5H, TdN, TeM, T5v, T5I;
Chris@82	1691 E T8j, T8q, TdG, TeO, T8g, T8p, Tba, Tbj;
Chris@82	1692 {
Chris@82	1693 E T3, T23, Td, T25, T6, T5R, Ta, T24;
Chris@82	1694 {
Chris@82	1695 E T1, T2, Tb, Tc;
Chris@82	1696 T1 = R0[0];
Chris@82	1697 T2 = R0[WS(rs, 32)];
Chris@82	1698 T3 = T1 + T2;
Chris@82	1699 T23 = T1 - T2;
Chris@82	1700 Tb = R0[WS(rs, 56)];
Chris@82	1701 Tc = R0[WS(rs, 24)];
Chris@82	1702 Td = Tb + Tc;
Chris@82	1703 T25 = Tb - Tc;
Chris@82	1704 }
Chris@82	1705 {
Chris@82	1706 E T4, T5, T8, T9;
Chris@82	1707 T4 = R0[WS(rs, 16)];
Chris@82	1708 T5 = R0[WS(rs, 48)];
Chris@82	1709 T6 = T4 + T5;
Chris@82	1710 T5R = T4 - T5;
Chris@82	1711 T8 = R0[WS(rs, 8)];
Chris@82	1712 T9 = R0[WS(rs, 40)];
Chris@82	1713 Ta = T8 + T9;
Chris@82	1714 T24 = T8 - T9;
Chris@82	1715 }
Chris@82	1716 TcD = T3 - T6;
Chris@82	1717 TdU = Td - Ta;
Chris@82	1718 {
Chris@82	1719 E T26, T5Q, T7, Te;
Chris@82	1720 T26 = KP707106781 * (T24 + T25);
Chris@82	1721 T27 = T23 + T26;
Chris@82	1722 T7r = T23 - T26;
Chris@82	1723 T5Q = KP707106781 * (T25 - T24);
Chris@82	1724 T5S = T5Q - T5R;
Chris@82	1725 T8y = T5R + T5Q;
Chris@82	1726 T7 = T3 + T6;
Chris@82	1727 Te = Ta + Td;
Chris@82	1728 Tf = T7 + Te;
Chris@82	1729 Ta5 = T7 - Te;
Chris@82	1730 }
Chris@82	1731 }
Chris@82	1732 {
Chris@82	1733 E Ti, T28, Ts, T2c, Tl, T29, Tp, T2b;
Chris@82	1734 {
Chris@82	1735 E Tg, Th, Tq, Tr;
Chris@82	1736 Tg = R0[WS(rs, 4)];
Chris@82	1737 Th = R0[WS(rs, 36)];
Chris@82	1738 Ti = Tg + Th;
Chris@82	1739 T28 = Tg - Th;
Chris@82	1740 Tq = R0[WS(rs, 12)];
Chris@82	1741 Tr = R0[WS(rs, 44)];
Chris@82	1742 Ts = Tq + Tr;
Chris@82	1743 T2c = Tq - Tr;
Chris@82	1744 }
Chris@82	1745 {
Chris@82	1746 E Tj, Tk, Tn, To;
Chris@82	1747 Tj = R0[WS(rs, 20)];
Chris@82	1748 Tk = R0[WS(rs, 52)];
Chris@82	1749 Tl = Tj + Tk;
Chris@82	1750 T29 = Tj - Tk;
Chris@82	1751 Tn = R0[WS(rs, 60)];
Chris@82	1752 To = R0[WS(rs, 28)];
Chris@82	1753 Tp = Tn + To;
Chris@82	1754 T2b = Tn - To;
Chris@82	1755 }
Chris@82	1756 {
Chris@82	1757 E Tm, Tt, TcE, TcF;
Chris@82	1758 Tm = Ti + Tl;
Chris@82	1759 Tt = Tp + Ts;
Chris@82	1760 Tu = Tm + Tt;
Chris@82	1761 Tbq = Tt - Tm;
Chris@82	1762 TcE = Ti - Tl;
Chris@82	1763 TcF = Tp - Ts;
Chris@82	1764 TcG = KP707106781 * (TcE + TcF);
Chris@82	1765 TdV = KP707106781 * (TcF - TcE);
Chris@82	1766 }
Chris@82	1767 {
Chris@82	1768 E T2a, T2d, T5T, T5U;
Chris@82	1769 T2a = FNMS(KP382683432, T29, KP923879532 * T28);
Chris@82	1770 T2d = FMA(KP923879532, T2b, KP382683432 * T2c);
Chris@82	1771 T2e = T2a + T2d;
Chris@82	1772 T8z = T2d - T2a;
Chris@82	1773 T5T = FNMS(KP923879532, T2c, KP382683432 * T2b);
Chris@82	1774 T5U = FMA(KP382683432, T28, KP923879532 * T29);
Chris@82	1775 T5V = T5T - T5U;
Chris@82	1776 T7s = T5U + T5T;
Chris@82	1777 }
Chris@82	1778 }
Chris@82	1779 {
Chris@82	1780 E Ty, T2g, TB, T2m, TF, T2l, TI, T2j;
Chris@82	1781 {
Chris@82	1782 E Tw, Tx, Tz, TA;
Chris@82	1783 Tw = R0[WS(rs, 2)];
Chris@82	1784 Tx = R0[WS(rs, 34)];
Chris@82	1785 Ty = Tw + Tx;
Chris@82	1786 T2g = Tw - Tx;
Chris@82	1787 Tz = R0[WS(rs, 18)];
Chris@82	1788 TA = R0[WS(rs, 50)];
Chris@82	1789 TB = Tz + TA;
Chris@82	1790 T2m = Tz - TA;
Chris@82	1791 {
Chris@82	1792 E TD, TE, T2h, TG, TH, T2i;
Chris@82	1793 TD = R0[WS(rs, 10)];
Chris@82	1794 TE = R0[WS(rs, 42)];
Chris@82	1795 T2h = TD - TE;
Chris@82	1796 TG = R0[WS(rs, 58)];
Chris@82	1797 TH = R0[WS(rs, 26)];
Chris@82	1798 T2i = TG - TH;
Chris@82	1799 TF = TD + TE;
Chris@82	1800 T2l = KP707106781 * (T2i - T2h);
Chris@82	1801 TI = TG + TH;
Chris@82	1802 T2j = KP707106781 * (T2h + T2i);
Chris@82	1803 }
Chris@82	1804 }
Chris@82	1805 {
Chris@82	1806 E TC, TJ, TcI, TcJ;
Chris@82	1807 TC = Ty + TB;
Chris@82	1808 TJ = TF + TI;
Chris@82	1809 TK = TC + TJ;
Chris@82	1810 Ta6 = TC - TJ;
Chris@82	1811 TcI = Ty - TB;
Chris@82	1812 TcJ = TI - TF;
Chris@82	1813 TcK = FMA(KP923879532, TcI, KP382683432 * TcJ);
Chris@82	1814 TdX = FNMS(KP382683432, TcI, KP923879532 * TcJ);
Chris@82	1815 }
Chris@82	1816 {
Chris@82	1817 E T2k, T2n, T7u, T7v;
Chris@82	1818 T2k = T2g + T2j;
Chris@82	1819 T2n = T2l - T2m;
Chris@82	1820 T2o = FMA(KP980785280, T2k, KP195090322 * T2n);
Chris@82	1821 T5X = FNMS(KP195090322, T2k, KP980785280 * T2n);
Chris@82	1822 T7u = T2g - T2j;
Chris@82	1823 T7v = T2m + T2l;
Chris@82	1824 T7w = FMA(KP831469612, T7u, KP555570233 * T7v);
Chris@82	1825 T8B = FNMS(KP555570233, T7u, KP831469612 * T7v);
Chris@82	1826 }
Chris@82	1827 }
Chris@82	1828 {
Chris@82	1829 E TN, T2p, TQ, T2v, TU, T2u, TX, T2s;
Chris@82	1830 {
Chris@82	1831 E TL, TM, TO, TP;
Chris@82	1832 TL = R0[WS(rs, 62)];
Chris@82	1833 TM = R0[WS(rs, 30)];
Chris@82	1834 TN = TL + TM;
Chris@82	1835 T2p = TL - TM;
Chris@82	1836 TO = R0[WS(rs, 14)];
Chris@82	1837 TP = R0[WS(rs, 46)];
Chris@82	1838 TQ = TO + TP;
Chris@82	1839 T2v = TO - TP;
Chris@82	1840 {
Chris@82	1841 E TS, TT, T2q, TV, TW, T2r;
Chris@82	1842 TS = R0[WS(rs, 6)];
Chris@82	1843 TT = R0[WS(rs, 38)];
Chris@82	1844 T2q = TS - TT;
Chris@82	1845 TV = R0[WS(rs, 54)];
Chris@82	1846 TW = R0[WS(rs, 22)];
Chris@82	1847 T2r = TV - TW;
Chris@82	1848 TU = TS + TT;
Chris@82	1849 T2u = KP707106781 * (T2r - T2q);
Chris@82	1850 TX = TV + TW;
Chris@82	1851 T2s = KP707106781 * (T2q + T2r);
Chris@82	1852 }
Chris@82	1853 }
Chris@82	1854 {
Chris@82	1855 E TR, TY, TcL, TcM;
Chris@82	1856 TR = TN + TQ;
Chris@82	1857 TY = TU + TX;
Chris@82	1858 TZ = TR + TY;
Chris@82	1859 Ta7 = TR - TY;
Chris@82	1860 TcL = TN - TQ;
Chris@82	1861 TcM = TX - TU;
Chris@82	1862 TcN = FNMS(KP382683432, TcM, KP923879532 * TcL);
Chris@82	1863 TdY = FMA(KP382683432, TcL, KP923879532 * TcM);
Chris@82	1864 }
Chris@82	1865 {
Chris@82	1866 E T2t, T2w, T7x, T7y;
Chris@82	1867 T2t = T2p + T2s;
Chris@82	1868 T2w = T2u - T2v;
Chris@82	1869 T2x = FNMS(KP195090322, T2w, KP980785280 * T2t);
Chris@82	1870 T5Y = FMA(KP195090322, T2t, KP980785280 * T2w);
Chris@82	1871 T7x = T2p - T2s;
Chris@82	1872 T7y = T2v + T2u;
Chris@82	1873 T7z = FNMS(KP555570233, T7y, KP831469612 * T7x);
Chris@82	1874 T8C = FMA(KP555570233, T7x, KP831469612 * T7y);
Chris@82	1875 }
Chris@82	1876 }
Chris@82	1877 {
Chris@82	1878 E T14, T2N, T17, T2D, T1b, T2O, T1e, T2C, T1j, T1m, T2K, TcR, T2Q, T1q, T1t;
Chris@82	1879 E T2H, TcS, T2R;
Chris@82	1880 {
Chris@82	1881 E T12, T13, T15, T16;
Chris@82	1882 T12 = R0[WS(rs, 1)];
Chris@82	1883 T13 = R0[WS(rs, 33)];
Chris@82	1884 T14 = T12 + T13;
Chris@82	1885 T2N = T12 - T13;
Chris@82	1886 T15 = R0[WS(rs, 17)];
Chris@82	1887 T16 = R0[WS(rs, 49)];
Chris@82	1888 T17 = T15 + T16;
Chris@82	1889 T2D = T15 - T16;
Chris@82	1890 }
Chris@82	1891 {
Chris@82	1892 E T19, T1a, T2B, T1c, T1d, T2A;
Chris@82	1893 T19 = R0[WS(rs, 9)];
Chris@82	1894 T1a = R0[WS(rs, 41)];
Chris@82	1895 T2B = T19 - T1a;
Chris@82	1896 T1c = R0[WS(rs, 57)];
Chris@82	1897 T1d = R0[WS(rs, 25)];
Chris@82	1898 T2A = T1c - T1d;
Chris@82	1899 T1b = T19 + T1a;
Chris@82	1900 T2O = KP707106781 * (T2B + T2A);
Chris@82	1901 T1e = T1c + T1d;
Chris@82	1902 T2C = KP707106781 * (T2A - T2B);
Chris@82	1903 }
Chris@82	1904 {
Chris@82	1905 E T2I, T2J, T2F, T2G;
Chris@82	1906 {
Chris@82	1907 E T1h, T1i, T1k, T1l;
Chris@82	1908 T1h = R0[WS(rs, 5)];
Chris@82	1909 T1i = R0[WS(rs, 37)];
Chris@82	1910 T1j = T1h + T1i;
Chris@82	1911 T2I = T1h - T1i;
Chris@82	1912 T1k = R0[WS(rs, 21)];
Chris@82	1913 T1l = R0[WS(rs, 53)];
Chris@82	1914 T1m = T1k + T1l;
Chris@82	1915 T2J = T1k - T1l;
Chris@82	1916 }
Chris@82	1917 T2K = FMA(KP382683432, T2I, KP923879532 * T2J);
Chris@82	1918 TcR = T1j - T1m;
Chris@82	1919 T2Q = FNMS(KP382683432, T2J, KP923879532 * T2I);
Chris@82	1920 {
Chris@82	1921 E T1o, T1p, T1r, T1s;
Chris@82	1922 T1o = R0[WS(rs, 61)];
Chris@82	1923 T1p = R0[WS(rs, 29)];
Chris@82	1924 T1q = T1o + T1p;
Chris@82	1925 T2F = T1o - T1p;
Chris@82	1926 T1r = R0[WS(rs, 13)];
Chris@82	1927 T1s = R0[WS(rs, 45)];
Chris@82	1928 T1t = T1r + T1s;
Chris@82	1929 T2G = T1r - T1s;
Chris@82	1930 }
Chris@82	1931 T2H = FNMS(KP923879532, T2G, KP382683432 * T2F);
Chris@82	1932 TcS = T1q - T1t;
Chris@82	1933 T2R = FMA(KP923879532, T2F, KP382683432 * T2G);
Chris@82	1934 }
Chris@82	1935 {
Chris@82	1936 E T18, T1f, TcQ, TcT;
Chris@82	1937 T18 = T14 + T17;
Chris@82	1938 T1f = T1b + T1e;
Chris@82	1939 T1g = T18 + T1f;
Chris@82	1940 Taa = T18 - T1f;
Chris@82	1941 TcQ = T14 - T17;
Chris@82	1942 TcT = KP707106781 * (TcR + TcS);
Chris@82	1943 TcU = TcQ + TcT;
Chris@82	1944 TeA = TcQ - TcT;
Chris@82	1945 }
Chris@82	1946 {
Chris@82	1947 E TcV, TcW, T1n, T1u;
Chris@82	1948 TcV = T1e - T1b;
Chris@82	1949 TcW = KP707106781 * (TcS - TcR);
Chris@82	1950 TcX = TcV + TcW;
Chris@82	1951 Tez = TcW - TcV;
Chris@82	1952 T1n = T1j + T1m;
Chris@82	1953 T1u = T1q + T1t;
Chris@82	1954 T1v = T1n + T1u;
Chris@82	1955 Tab = T1u - T1n;
Chris@82	1956 }
Chris@82	1957 {
Chris@82	1958 E T2E, T2L, T7C, T7D;
Chris@82	1959 T2E = T2C - T2D;
Chris@82	1960 T2L = T2H - T2K;
Chris@82	1961 T2M = T2E + T2L;
Chris@82	1962 T6z = T2L - T2E;
Chris@82	1963 T7C = T2N - T2O;
Chris@82	1964 T7D = T2K + T2H;
Chris@82	1965 T7E = T7C + T7D;
Chris@82	1966 T9e = T7C - T7D;
Chris@82	1967 }
Chris@82	1968 {
Chris@82	1969 E T7F, T7G, T2P, T2S;
Chris@82	1970 T7F = T2D + T2C;
Chris@82	1971 T7G = T2R - T2Q;
Chris@82	1972 T7H = T7F + T7G;
Chris@82	1973 T9d = T7G - T7F;
Chris@82	1974 T2P = T2N + T2O;
Chris@82	1975 T2S = T2Q + T2R;
Chris@82	1976 T2T = T2P + T2S;
Chris@82	1977 T6A = T2P - T2S;
Chris@82	1978 }
Chris@82	1979 }
Chris@82	1980 {
Chris@82	1981 E T4z, TaP, T5B, TaQ, T4G, TaT, T5y, TaS, Tbf, Tbg, T4O, Tdw, T5E, Tbc, Tbd;
Chris@82	1982 E T4V, Tdx, T5D;
Chris@82	1983 {
Chris@82	1984 E T4x, T4y, T5z, T5A;
Chris@82	1985 T4x = R1[WS(rs, 63)];
Chris@82	1986 T4y = R1[WS(rs, 31)];
Chris@82	1987 T4z = T4x - T4y;
Chris@82	1988 TaP = T4x + T4y;
Chris@82	1989 T5z = R1[WS(rs, 15)];
Chris@82	1990 T5A = R1[WS(rs, 47)];
Chris@82	1991 T5B = T5z - T5A;
Chris@82	1992 TaQ = T5z + T5A;
Chris@82	1993 }
Chris@82	1994 {
Chris@82	1995 E T4A, T4B, T4C, T4D, T4E, T4F;
Chris@82	1996 T4A = R1[WS(rs, 7)];
Chris@82	1997 T4B = R1[WS(rs, 39)];
Chris@82	1998 T4C = T4A - T4B;
Chris@82	1999 T4D = R1[WS(rs, 55)];
Chris@82	2000 T4E = R1[WS(rs, 23)];
Chris@82	2001 T4F = T4D - T4E;
Chris@82	2002 T4G = KP707106781 * (T4C + T4F);
Chris@82	2003 TaT = T4D + T4E;
Chris@82	2004 T5y = KP707106781 * (T4F - T4C);
Chris@82	2005 TaS = T4A + T4B;
Chris@82	2006 }
Chris@82	2007 {
Chris@82	2008 E T4K, T4N, T4R, T4U;
Chris@82	2009 {
Chris@82	2010 E T4I, T4J, T4L, T4M;
Chris@82	2011 T4I = R1[WS(rs, 3)];
Chris@82	2012 T4J = R1[WS(rs, 35)];
Chris@82	2013 T4K = T4I - T4J;
Chris@82	2014 Tbf = T4I + T4J;
Chris@82	2015 T4L = R1[WS(rs, 19)];
Chris@82	2016 T4M = R1[WS(rs, 51)];
Chris@82	2017 T4N = T4L - T4M;
Chris@82	2018 Tbg = T4L + T4M;
Chris@82	2019 }
Chris@82	2020 T4O = FNMS(KP382683432, T4N, KP923879532 * T4K);
Chris@82	2021 Tdw = Tbf - Tbg;
Chris@82	2022 T5E = FMA(KP382683432, T4K, KP923879532 * T4N);
Chris@82	2023 {
Chris@82	2024 E T4P, T4Q, T4S, T4T;
Chris@82	2025 T4P = R1[WS(rs, 59)];
Chris@82	2026 T4Q = R1[WS(rs, 27)];
Chris@82	2027 T4R = T4P - T4Q;
Chris@82	2028 Tbc = T4P + T4Q;
Chris@82	2029 T4S = R1[WS(rs, 11)];
Chris@82	2030 T4T = R1[WS(rs, 43)];
Chris@82	2031 T4U = T4S - T4T;
Chris@82	2032 Tbd = T4S + T4T;
Chris@82	2033 }
Chris@82	2034 T4V = FMA(KP923879532, T4R, KP382683432 * T4U);
Chris@82	2035 Tdx = Tbc - Tbd;
Chris@82	2036 T5D = FNMS(KP923879532, T4U, KP382683432 * T4R);
Chris@82	2037 }
Chris@82	2038 {
Chris@82	2039 E T4H, T4W, Tdv, Tdy;
Chris@82	2040 T4H = T4z + T4G;
Chris@82	2041 T4W = T4O + T4V;
Chris@82	2042 T4X = T4H + T4W;
Chris@82	2043 T6L = T4H - T4W;
Chris@82	2044 Tdv = TaP - TaQ;
Chris@82	2045 Tdy = KP707106781 * (Tdw + Tdx);
Chris@82	2046 Tdz = Tdv + Tdy;
Chris@82	2047 TeL = Tdv - Tdy;
Chris@82	2048 }
Chris@82	2049 {
Chris@82	2050 E TdI, TdJ, T5C, T5F;
Chris@82	2051 TdI = TaT - TaS;
Chris@82	2052 TdJ = KP707106781 * (Tdx - Tdw);
Chris@82	2053 TdK = TdI + TdJ;
Chris@82	2054 TeP = TdJ - TdI;
Chris@82	2055 T5C = T5y - T5B;
Chris@82	2056 T5F = T5D - T5E;
Chris@82	2057 T5G = T5C + T5F;
Chris@82	2058 T6P = T5F - T5C;
Chris@82	2059 }
Chris@82	2060 {
Chris@82	2061 E T8b, T8c, TaR, TaU;
Chris@82	2062 T8b = T4z - T4G;
Chris@82	2063 T8c = T5E + T5D;
Chris@82	2064 T8d = T8b + T8c;
Chris@82	2065 T9p = T8b - T8c;
Chris@82	2066 TaR = TaP + TaQ;
Chris@82	2067 TaU = TaS + TaT;
Chris@82	2068 TaV = TaR - TaU;
Chris@82	2069 Tc3 = TaR + TaU;
Chris@82	2070 }
Chris@82	2071 {
Chris@82	2072 E Tbe, Tbh, T8m, T8n;
Chris@82	2073 Tbe = Tbc + Tbd;
Chris@82	2074 Tbh = Tbf + Tbg;
Chris@82	2075 Tbi = Tbe - Tbh;
Chris@82	2076 Tc4 = Tbh + Tbe;
Chris@82	2077 T8m = T5B + T5y;
Chris@82	2078 T8n = T4V - T4O;
Chris@82	2079 T8o = T8m + T8n;
Chris@82	2080 T9t = T8n - T8m;
Chris@82	2081 }
Chris@82	2082 }
Chris@82	2083 {
Chris@82	2084 E T3k, Tai, T4m, Taj, T3r, Tam, T4j, Tal, TaI, TaJ, T3z, Tdb, T4p, TaF, TaG;
Chris@82	2085 E T3G, Tdc, T4o;
Chris@82	2086 {
Chris@82	2087 E T3i, T3j, T4k, T4l;
Chris@82	2088 T3i = R1[0];
Chris@82	2089 T3j = R1[WS(rs, 32)];
Chris@82	2090 T3k = T3i - T3j;
Chris@82	2091 Tai = T3i + T3j;
Chris@82	2092 T4k = R1[WS(rs, 16)];
Chris@82	2093 T4l = R1[WS(rs, 48)];
Chris@82	2094 T4m = T4k - T4l;
Chris@82	2095 Taj = T4k + T4l;
Chris@82	2096 }
Chris@82	2097 {
Chris@82	2098 E T3l, T3m, T3n, T3o, T3p, T3q;
Chris@82	2099 T3l = R1[WS(rs, 8)];
Chris@82	2100 T3m = R1[WS(rs, 40)];
Chris@82	2101 T3n = T3l - T3m;
Chris@82	2102 T3o = R1[WS(rs, 56)];
Chris@82	2103 T3p = R1[WS(rs, 24)];
Chris@82	2104 T3q = T3o - T3p;
Chris@82	2105 T3r = KP707106781 * (T3n + T3q);
Chris@82	2106 Tam = T3o + T3p;
Chris@82	2107 T4j = KP707106781 * (T3q - T3n);
Chris@82	2108 Tal = T3l + T3m;
Chris@82	2109 }
Chris@82	2110 {
Chris@82	2111 E T3v, T3y, T3C, T3F;
Chris@82	2112 {
Chris@82	2113 E T3t, T3u, T3w, T3x;
Chris@82	2114 T3t = R1[WS(rs, 4)];
Chris@82	2115 T3u = R1[WS(rs, 36)];
Chris@82	2116 T3v = T3t - T3u;
Chris@82	2117 TaI = T3t + T3u;
Chris@82	2118 T3w = R1[WS(rs, 20)];
Chris@82	2119 T3x = R1[WS(rs, 52)];
Chris@82	2120 T3y = T3w - T3x;
Chris@82	2121 TaJ = T3w + T3x;
Chris@82	2122 }
Chris@82	2123 T3z = FNMS(KP382683432, T3y, KP923879532 * T3v);
Chris@82	2124 Tdb = TaI - TaJ;
Chris@82	2125 T4p = FMA(KP382683432, T3v, KP923879532 * T3y);
Chris@82	2126 {
Chris@82	2127 E T3A, T3B, T3D, T3E;
Chris@82	2128 T3A = R1[WS(rs, 60)];
Chris@82	2129 T3B = R1[WS(rs, 28)];
Chris@82	2130 T3C = T3A - T3B;
Chris@82	2131 TaF = T3A + T3B;
Chris@82	2132 T3D = R1[WS(rs, 12)];
Chris@82	2133 T3E = R1[WS(rs, 44)];
Chris@82	2134 T3F = T3D - T3E;
Chris@82	2135 TaG = T3D + T3E;
Chris@82	2136 }
Chris@82	2137 T3G = FMA(KP923879532, T3C, KP382683432 * T3F);
Chris@82	2138 Tdc = TaF - TaG;
Chris@82	2139 T4o = FNMS(KP923879532, T3F, KP382683432 * T3C);
Chris@82	2140 }
Chris@82	2141 {
Chris@82	2142 E T3s, T3H, Tda, Tdd;
Chris@82	2143 T3s = T3k + T3r;
Chris@82	2144 T3H = T3z + T3G;
Chris@82	2145 T3I = T3s + T3H;
Chris@82	2146 T6H = T3s - T3H;
Chris@82	2147 Tda = Tai - Taj;
Chris@82	2148 Tdd = KP707106781 * (Tdb + Tdc);
Chris@82	2149 Tde = Tda + Tdd;
Chris@82	2150 TeH = Tda - Tdd;
Chris@82	2151 }
Chris@82	2152 {
Chris@82	2153 E Tdn, Tdo, T4n, T4q;
Chris@82	2154 Tdn = Tam - Tal;
Chris@82	2155 Tdo = KP707106781 * (Tdc - Tdb);
Chris@82	2156 Tdp = Tdn + Tdo;
Chris@82	2157 TeF = Tdo - Tdn;
Chris@82	2158 T4n = T4j - T4m;
Chris@82	2159 T4q = T4o - T4p;
Chris@82	2160 T4r = T4n + T4q;
Chris@82	2161 T6F = T4q - T4n;
Chris@82	2162 }
Chris@82	2163 {
Chris@82	2164 E T7S, T7T, Tak, Tan;
Chris@82	2165 T7S = T3k - T3r;
Chris@82	2166 T7T = T4p + T4o;
Chris@82	2167 T7U = T7S + T7T;
Chris@82	2168 T9l = T7S - T7T;
Chris@82	2169 Tak = Tai + Taj;
Chris@82	2170 Tan = Tal + Tam;
Chris@82	2171 Tao = Tak - Tan;
Chris@82	2172 TbW = Tak + Tan;
Chris@82	2173 }
Chris@82	2174 {
Chris@82	2175 E TaH, TaK, T83, T84;
Chris@82	2176 TaH = TaF + TaG;
Chris@82	2177 TaK = TaI + TaJ;
Chris@82	2178 TaL = TaH - TaK;
Chris@82	2179 TbX = TaK + TaH;
Chris@82	2180 T83 = T4m + T4j;
Chris@82	2181 T84 = T3G - T3z;
Chris@82	2182 T85 = T83 + T84;
Chris@82	2183 T9j = T84 - T83;
Chris@82	2184 }
Chris@82	2185 }
Chris@82	2186 {
Chris@82	2187 E T1z, T2V, T1C, T39, T1G, T38, T1J, T2Y, T1O, T1R, T32, Td0, T3c, T1V, T1Y;
Chris@82	2188 E T35, Td1, T3b;
Chris@82	2189 {
Chris@82	2190 E T1x, T1y, T1A, T1B;
Chris@82	2191 T1x = R0[WS(rs, 63)];
Chris@82	2192 T1y = R0[WS(rs, 31)];
Chris@82	2193 T1z = T1x + T1y;
Chris@82	2194 T2V = T1x - T1y;
Chris@82	2195 T1A = R0[WS(rs, 15)];
Chris@82	2196 T1B = R0[WS(rs, 47)];
Chris@82	2197 T1C = T1A + T1B;
Chris@82	2198 T39 = T1A - T1B;
Chris@82	2199 }
Chris@82	2200 {
Chris@82	2201 E T1E, T1F, T2W, T1H, T1I, T2X;
Chris@82	2202 T1E = R0[WS(rs, 7)];
Chris@82	2203 T1F = R0[WS(rs, 39)];
Chris@82	2204 T2W = T1E - T1F;
Chris@82	2205 T1H = R0[WS(rs, 55)];
Chris@82	2206 T1I = R0[WS(rs, 23)];
Chris@82	2207 T2X = T1H - T1I;
Chris@82	2208 T1G = T1E + T1F;
Chris@82	2209 T38 = KP707106781 * (T2X - T2W);
Chris@82	2210 T1J = T1H + T1I;
Chris@82	2211 T2Y = KP707106781 * (T2W + T2X);
Chris@82	2212 }
Chris@82	2213 {
Chris@82	2214 E T30, T31, T33, T34;
Chris@82	2215 {
Chris@82	2216 E T1M, T1N, T1P, T1Q;
Chris@82	2217 T1M = R0[WS(rs, 3)];
Chris@82	2218 T1N = R0[WS(rs, 35)];
Chris@82	2219 T1O = T1M + T1N;
Chris@82	2220 T30 = T1M - T1N;
Chris@82	2221 T1P = R0[WS(rs, 19)];
Chris@82	2222 T1Q = R0[WS(rs, 51)];
Chris@82	2223 T1R = T1P + T1Q;
Chris@82	2224 T31 = T1P - T1Q;
Chris@82	2225 }
Chris@82	2226 T32 = FNMS(KP382683432, T31, KP923879532 * T30);
Chris@82	2227 Td0 = T1O - T1R;
Chris@82	2228 T3c = FMA(KP382683432, T30, KP923879532 * T31);
Chris@82	2229 {
Chris@82	2230 E T1T, T1U, T1W, T1X;
Chris@82	2231 T1T = R0[WS(rs, 59)];
Chris@82	2232 T1U = R0[WS(rs, 27)];
Chris@82	2233 T1V = T1T + T1U;
Chris@82	2234 T33 = T1T - T1U;
Chris@82	2235 T1W = R0[WS(rs, 11)];
Chris@82	2236 T1X = R0[WS(rs, 43)];
Chris@82	2237 T1Y = T1W + T1X;
Chris@82	2238 T34 = T1W - T1X;
Chris@82	2239 }
Chris@82	2240 T35 = FMA(KP923879532, T33, KP382683432 * T34);
Chris@82	2241 Td1 = T1V - T1Y;
Chris@82	2242 T3b = FNMS(KP923879532, T34, KP382683432 * T33);
Chris@82	2243 }
Chris@82	2244 {
Chris@82	2245 E T1D, T1K, TcZ, Td2;
Chris@82	2246 T1D = T1z + T1C;
Chris@82	2247 T1K = T1G + T1J;
Chris@82	2248 T1L = T1D + T1K;
Chris@82	2249 Tad = T1D - T1K;
Chris@82	2250 TcZ = T1z - T1C;
Chris@82	2251 Td2 = KP707106781 * (Td0 + Td1);
Chris@82	2252 Td3 = TcZ + Td2;
Chris@82	2253 Tew = TcZ - Td2;
Chris@82	2254 }
Chris@82	2255 {
Chris@82	2256 E Td4, Td5, T1S, T1Z;
Chris@82	2257 Td4 = T1J - T1G;
Chris@82	2258 Td5 = KP707106781 * (Td1 - Td0);
Chris@82	2259 Td6 = Td4 + Td5;
Chris@82	2260 Tex = Td5 - Td4;
Chris@82	2261 T1S = T1O + T1R;
Chris@82	2262 T1Z = T1V + T1Y;
Chris@82	2263 T20 = T1S + T1Z;
Chris@82	2264 Tae = T1Z - T1S;
Chris@82	2265 }
Chris@82	2266 {
Chris@82	2267 E T2Z, T36, T7J, T7K;
Chris@82	2268 T2Z = T2V + T2Y;
Chris@82	2269 T36 = T32 + T35;
Chris@82	2270 T37 = T2Z + T36;
Chris@82	2271 T6x = T2Z - T36;
Chris@82	2272 T7J = T2V - T2Y;
Chris@82	2273 T7K = T3c + T3b;
Chris@82	2274 T7L = T7J + T7K;
Chris@82	2275 T9a = T7J - T7K;
Chris@82	2276 }
Chris@82	2277 {
Chris@82	2278 E T7M, T7N, T3a, T3d;
Chris@82	2279 T7M = T39 + T38;
Chris@82	2280 T7N = T35 - T32;
Chris@82	2281 T7O = T7M + T7N;
Chris@82	2282 T9b = T7N - T7M;
Chris@82	2283 T3a = T38 - T39;
Chris@82	2284 T3d = T3b - T3c;
Chris@82	2285 T3e = T3a + T3d;
Chris@82	2286 T6w = T3d - T3a;
Chris@82	2287 }
Chris@82	2288 }
Chris@82	2289 {
Chris@82	2290 E T3L, Tdf, T3X, Tar, T42, Tdi, T4e, Tay, T3S, Tdg, T3U, Tau, T49, Tdj, T4b;
Chris@82	2291 E TaB, Tdh, Tdk;
Chris@82	2292 {
Chris@82	2293 E T3J, T3K, Tap, T3V, T3W, Taq;
Chris@82	2294 T3J = R1[WS(rs, 2)];
Chris@82	2295 T3K = R1[WS(rs, 34)];
Chris@82	2296 Tap = T3J + T3K;
Chris@82	2297 T3V = R1[WS(rs, 18)];
Chris@82	2298 T3W = R1[WS(rs, 50)];
Chris@82	2299 Taq = T3V + T3W;
Chris@82	2300 T3L = T3J - T3K;
Chris@82	2301 Tdf = Tap - Taq;
Chris@82	2302 T3X = T3V - T3W;
Chris@82	2303 Tar = Tap + Taq;
Chris@82	2304 }
Chris@82	2305 {
Chris@82	2306 E T40, T41, Taw, T4c, T4d, Tax;
Chris@82	2307 T40 = R1[WS(rs, 62)];
Chris@82	2308 T41 = R1[WS(rs, 30)];
Chris@82	2309 Taw = T40 + T41;
Chris@82	2310 T4c = R1[WS(rs, 14)];
Chris@82	2311 T4d = R1[WS(rs, 46)];
Chris@82	2312 Tax = T4c + T4d;
Chris@82	2313 T42 = T40 - T41;
Chris@82	2314 Tdi = Taw - Tax;
Chris@82	2315 T4e = T4c - T4d;
Chris@82	2316 Tay = Taw + Tax;
Chris@82	2317 }
Chris@82	2318 {
Chris@82	2319 E T3O, Tas, T3R, Tat;
Chris@82	2320 {
Chris@82	2321 E T3M, T3N, T3P, T3Q;
Chris@82	2322 T3M = R1[WS(rs, 10)];
Chris@82	2323 T3N = R1[WS(rs, 42)];
Chris@82	2324 T3O = T3M - T3N;
Chris@82	2325 Tas = T3M + T3N;
Chris@82	2326 T3P = R1[WS(rs, 58)];
Chris@82	2327 T3Q = R1[WS(rs, 26)];
Chris@82	2328 T3R = T3P - T3Q;
Chris@82	2329 Tat = T3P + T3Q;
Chris@82	2330 }
Chris@82	2331 T3S = KP707106781 * (T3O + T3R);
Chris@82	2332 Tdg = Tat - Tas;
Chris@82	2333 T3U = KP707106781 * (T3R - T3O);
Chris@82	2334 Tau = Tas + Tat;
Chris@82	2335 }
Chris@82	2336 {
Chris@82	2337 E T45, Taz, T48, TaA;
Chris@82	2338 {
Chris@82	2339 E T43, T44, T46, T47;
Chris@82	2340 T43 = R1[WS(rs, 6)];
Chris@82	2341 T44 = R1[WS(rs, 38)];
Chris@82	2342 T45 = T43 - T44;
Chris@82	2343 Taz = T43 + T44;
Chris@82	2344 T46 = R1[WS(rs, 54)];
Chris@82	2345 T47 = R1[WS(rs, 22)];
Chris@82	2346 T48 = T46 - T47;
Chris@82	2347 TaA = T46 + T47;
Chris@82	2348 }
Chris@82	2349 T49 = KP707106781 * (T45 + T48);
Chris@82	2350 Tdj = TaA - Taz;
Chris@82	2351 T4b = KP707106781 * (T48 - T45);
Chris@82	2352 TaB = Taz + TaA;
Chris@82	2353 }
Chris@82	2354 TbZ = Tar + Tau;
Chris@82	2355 Tc0 = Tay + TaB;
Chris@82	2356 {
Chris@82	2357 E T3T, T3Y, Tdq, Tdr;
Chris@82	2358 T3T = T3L + T3S;
Chris@82	2359 T3Y = T3U - T3X;
Chris@82	2360 T3Z = FMA(KP980785280, T3T, KP195090322 * T3Y);
Chris@82	2361 T4s = FNMS(KP195090322, T3T, KP980785280 * T3Y);
Chris@82	2362 Tdq = FNMS(KP382683432, Tdf, KP923879532 * Tdg);
Chris@82	2363 Tdr = FMA(KP382683432, Tdi, KP923879532 * Tdj);
Chris@82	2364 Tds = Tdq + Tdr;
Chris@82	2365 TeI = Tdr - Tdq;
Chris@82	2366 }
Chris@82	2367 {
Chris@82	2368 E T4a, T4f, T7Y, T7Z;
Chris@82	2369 T4a = T42 + T49;
Chris@82	2370 T4f = T4b - T4e;
Chris@82	2371 T4g = FNMS(KP195090322, T4f, KP980785280 * T4a);
Chris@82	2372 T4t = FMA(KP195090322, T4a, KP980785280 * T4f);
Chris@82	2373 T7Y = T42 - T49;
Chris@82	2374 T7Z = T4e + T4b;
Chris@82	2375 T80 = FNMS(KP555570233, T7Z, KP831469612 * T7Y);
Chris@82	2376 T87 = FMA(KP555570233, T7Y, KP831469612 * T7Z);
Chris@82	2377 }
Chris@82	2378 Tdh = FMA(KP923879532, Tdf, KP382683432 * Tdg);
Chris@82	2379 Tdk = FNMS(KP382683432, Tdj, KP923879532 * Tdi);
Chris@82	2380 Tdl = Tdh + Tdk;
Chris@82	2381 TeE = Tdk - Tdh;
Chris@82	2382 {
Chris@82	2383 E T7V, T7W, Tav, TaC;
Chris@82	2384 T7V = T3L - T3S;
Chris@82	2385 T7W = T3X + T3U;
Chris@82	2386 T7X = FMA(KP831469612, T7V, KP555570233 * T7W);
Chris@82	2387 T86 = FNMS(KP555570233, T7V, KP831469612 * T7W);
Chris@82	2388 Tav = Tar - Tau;
Chris@82	2389 TaC = Tay - TaB;
Chris@82	2390 TaD = KP707106781 * (Tav + TaC);
Chris@82	2391 TaM = KP707106781 * (TaC - Tav);
Chris@82	2392 }
Chris@82	2393 }
Chris@82	2394 {
Chris@82	2395 E T50, TdA, T5c, TaY, T5h, TdD, T5t, Tb5, T57, TdB, T59, Tb1, T5o, TdE, T5q;
Chris@82	2396 E Tb8, TdC, TdF;
Chris@82	2397 {
Chris@82	2398 E T4Y, T4Z, TaW, T5a, T5b, TaX;
Chris@82	2399 T4Y = R1[WS(rs, 1)];
Chris@82	2400 T4Z = R1[WS(rs, 33)];
Chris@82	2401 TaW = T4Y + T4Z;
Chris@82	2402 T5a = R1[WS(rs, 17)];
Chris@82	2403 T5b = R1[WS(rs, 49)];
Chris@82	2404 TaX = T5a + T5b;
Chris@82	2405 T50 = T4Y - T4Z;
Chris@82	2406 TdA = TaW - TaX;
Chris@82	2407 T5c = T5a - T5b;
Chris@82	2408 TaY = TaW + TaX;
Chris@82	2409 }
Chris@82	2410 {
Chris@82	2411 E T5f, T5g, Tb3, T5r, T5s, Tb4;
Chris@82	2412 T5f = R1[WS(rs, 61)];
Chris@82	2413 T5g = R1[WS(rs, 29)];
Chris@82	2414 Tb3 = T5f + T5g;
Chris@82	2415 T5r = R1[WS(rs, 13)];
Chris@82	2416 T5s = R1[WS(rs, 45)];
Chris@82	2417 Tb4 = T5r + T5s;
Chris@82	2418 T5h = T5f - T5g;
Chris@82	2419 TdD = Tb3 - Tb4;
Chris@82	2420 T5t = T5r - T5s;
Chris@82	2421 Tb5 = Tb3 + Tb4;
Chris@82	2422 }
Chris@82	2423 {
Chris@82	2424 E T53, TaZ, T56, Tb0;
Chris@82	2425 {
Chris@82	2426 E T51, T52, T54, T55;
Chris@82	2427 T51 = R1[WS(rs, 9)];
Chris@82	2428 T52 = R1[WS(rs, 41)];
Chris@82	2429 T53 = T51 - T52;
Chris@82	2430 TaZ = T51 + T52;
Chris@82	2431 T54 = R1[WS(rs, 57)];
Chris@82	2432 T55 = R1[WS(rs, 25)];
Chris@82	2433 T56 = T54 - T55;
Chris@82	2434 Tb0 = T54 + T55;
Chris@82	2435 }
Chris@82	2436 T57 = KP707106781 * (T53 + T56);
Chris@82	2437 TdB = Tb0 - TaZ;
Chris@82	2438 T59 = KP707106781 * (T56 - T53);
Chris@82	2439 Tb1 = TaZ + Tb0;
Chris@82	2440 }
Chris@82	2441 {
Chris@82	2442 E T5k, Tb6, T5n, Tb7;
Chris@82	2443 {
Chris@82	2444 E T5i, T5j, T5l, T5m;
Chris@82	2445 T5i = R1[WS(rs, 5)];
Chris@82	2446 T5j = R1[WS(rs, 37)];
Chris@82	2447 T5k = T5i - T5j;
Chris@82	2448 Tb6 = T5i + T5j;
Chris@82	2449 T5l = R1[WS(rs, 53)];
Chris@82	2450 T5m = R1[WS(rs, 21)];
Chris@82	2451 T5n = T5l - T5m;
Chris@82	2452 Tb7 = T5l + T5m;
Chris@82	2453 }
Chris@82	2454 T5o = KP707106781 * (T5k + T5n);
Chris@82	2455 TdE = Tb7 - Tb6;
Chris@82	2456 T5q = KP707106781 * (T5n - T5k);
Chris@82	2457 Tb8 = Tb6 + Tb7;
Chris@82	2458 }
Chris@82	2459 Tc6 = TaY + Tb1;
Chris@82	2460 Tc7 = Tb5 + Tb8;
Chris@82	2461 {
Chris@82	2462 E T58, T5d, TdL, TdM;
Chris@82	2463 T58 = T50 + T57;
Chris@82	2464 T5d = T59 - T5c;
Chris@82	2465 T5e = FMA(KP980785280, T58, KP195090322 * T5d);
Chris@82	2466 T5H = FNMS(KP195090322, T58, KP980785280 * T5d);
Chris@82	2467 TdL = FNMS(KP382683432, TdA, KP923879532 * TdB);
Chris@82	2468 TdM = FMA(KP382683432, TdD, KP923879532 * TdE);
Chris@82	2469 TdN = TdL + TdM;
Chris@82	2470 TeM = TdM - TdL;
Chris@82	2471 }
Chris@82	2472 {
Chris@82	2473 E T5p, T5u, T8h, T8i;
Chris@82	2474 T5p = T5h + T5o;
Chris@82	2475 T5u = T5q - T5t;
Chris@82	2476 T5v = FNMS(KP195090322, T5u, KP980785280 * T5p);
Chris@82	2477 T5I = FMA(KP195090322, T5p, KP980785280 * T5u);
Chris@82	2478 T8h = T5h - T5o;
Chris@82	2479 T8i = T5t + T5q;
Chris@82	2480 T8j = FNMS(KP555570233, T8i, KP831469612 * T8h);
Chris@82	2481 T8q = FMA(KP555570233, T8h, KP831469612 * T8i);
Chris@82	2482 }
Chris@82	2483 TdC = FMA(KP923879532, TdA, KP382683432 * TdB);
Chris@82	2484 TdF = FNMS(KP382683432, TdE, KP923879532 * TdD);
Chris@82	2485 TdG = TdC + TdF;
Chris@82	2486 TeO = TdF - TdC;
Chris@82	2487 {
Chris@82	2488 E T8e, T8f, Tb2, Tb9;
Chris@82	2489 T8e = T50 - T57;
Chris@82	2490 T8f = T5c + T59;
Chris@82	2491 T8g = FMA(KP831469612, T8e, KP555570233 * T8f);
Chris@82	2492 T8p = FNMS(KP555570233, T8e, KP831469612 * T8f);
Chris@82	2493 Tb2 = TaY - Tb1;
Chris@82	2494 Tb9 = Tb5 - Tb8;
Chris@82	2495 Tba = KP707106781 * (Tb2 + Tb9);
Chris@82	2496 Tbj = KP707106781 * (Tb9 - Tb2);
Chris@82	2497 }
Chris@82	2498 }
Chris@82	2499 {
Chris@82	2500 E T11, TbV, Tc9, Tcf, T22, Tcb, Tc2, Tce;
Chris@82	2501 {
Chris@82	2502 E Tv, T10, Tc5, Tc8;
Chris@82	2503 Tv = Tf + Tu;
Chris@82	2504 T10 = TK + TZ;
Chris@82	2505 T11 = Tv + T10;
Chris@82	2506 TbV = Tv - T10;
Chris@82	2507 Tc5 = Tc3 + Tc4;
Chris@82	2508 Tc8 = Tc6 + Tc7;
Chris@82	2509 Tc9 = Tc5 - Tc8;
Chris@82	2510 Tcf = Tc5 + Tc8;
Chris@82	2511 }
Chris@82	2512 {
Chris@82	2513 E T1w, T21, TbY, Tc1;
Chris@82	2514 T1w = T1g + T1v;
Chris@82	2515 T21 = T1L + T20;
Chris@82	2516 T22 = T1w + T21;
Chris@82	2517 Tcb = T21 - T1w;
Chris@82	2518 TbY = TbW + TbX;
Chris@82	2519 Tc1 = TbZ + Tc0;
Chris@82	2520 Tc2 = TbY - Tc1;
Chris@82	2521 Tce = TbY + Tc1;
Chris@82	2522 }
Chris@82	2523 Cr[WS(csr, 32)] = T11 - T22;
Chris@82	2524 Ci[WS(csi, 32)] = Tcf - Tce;
Chris@82	2525 {
Chris@82	2526 E Tca, Tcc, Tcd, Tcg;
Chris@82	2527 Tca = KP707106781 * (Tc2 + Tc9);
Chris@82	2528 Cr[WS(csr, 48)] = TbV - Tca;
Chris@82	2529 Cr[WS(csr, 16)] = TbV + Tca;
Chris@82	2530 Tcc = KP707106781 * (Tc9 - Tc2);
Chris@82	2531 Ci[WS(csi, 16)] = Tcb + Tcc;
Chris@82	2532 Ci[WS(csi, 48)] = Tcc - Tcb;
Chris@82	2533 Tcd = T11 + T22;
Chris@82	2534 Tcg = Tce + Tcf;
Chris@82	2535 Cr[WS(csr, 64)] = Tcd - Tcg;
Chris@82	2536 Cr[0] = Tcd + Tcg;
Chris@82	2537 }
Chris@82	2538 }
Chris@82	2539 {
Chris@82	2540 E Tch, Tcu, Tck, Tct, Tco, Tcy, Tcr, Tcz, Tci, Tcj;
Chris@82	2541 Tch = Tf - Tu;
Chris@82	2542 Tcu = TZ - TK;
Chris@82	2543 Tci = T1g - T1v;
Chris@82	2544 Tcj = T1L - T20;
Chris@82	2545 Tck = KP707106781 * (Tci + Tcj);
Chris@82	2546 Tct = KP707106781 * (Tcj - Tci);
Chris@82	2547 {
Chris@82	2548 E Tcm, Tcn, Tcp, Tcq;
Chris@82	2549 Tcm = TbW - TbX;
Chris@82	2550 Tcn = Tc0 - TbZ;
Chris@82	2551 Tco = FMA(KP923879532, Tcm, KP382683432 * Tcn);
Chris@82	2552 Tcy = FNMS(KP382683432, Tcm, KP923879532 * Tcn);
Chris@82	2553 Tcp = Tc3 - Tc4;
Chris@82	2554 Tcq = Tc7 - Tc6;
Chris@82	2555 Tcr = FNMS(KP382683432, Tcq, KP923879532 * Tcp);
Chris@82	2556 Tcz = FMA(KP382683432, Tcp, KP923879532 * Tcq);
Chris@82	2557 }
Chris@82	2558 {
Chris@82	2559 E Tcl, Tcs, Tcx, TcA;
Chris@82	2560 Tcl = Tch + Tck;
Chris@82	2561 Tcs = Tco + Tcr;
Chris@82	2562 Cr[WS(csr, 56)] = Tcl - Tcs;
Chris@82	2563 Cr[WS(csr, 8)] = Tcl + Tcs;
Chris@82	2564 Tcx = Tcu + Tct;
Chris@82	2565 TcA = Tcy + Tcz;
Chris@82	2566 Ci[WS(csi, 8)] = Tcx + TcA;
Chris@82	2567 Ci[WS(csi, 56)] = TcA - Tcx;
Chris@82	2568 }
Chris@82	2569 {
Chris@82	2570 E Tcv, Tcw, TcB, TcC;
Chris@82	2571 Tcv = Tct - Tcu;
Chris@82	2572 Tcw = Tcr - Tco;
Chris@82	2573 Ci[WS(csi, 24)] = Tcv + Tcw;
Chris@82	2574 Ci[WS(csi, 40)] = Tcw - Tcv;
Chris@82	2575 TcB = Tch - Tck;
Chris@82	2576 TcC = Tcz - Tcy;
Chris@82	2577 Cr[WS(csr, 40)] = TcB - TcC;
Chris@82	2578 Cr[WS(csr, 24)] = TcB + TcC;
Chris@82	2579 }
Chris@82	2580 }
Chris@82	2581 {
Chris@82	2582 E Ta9, TbB, Tbs, TbM, Tag, TbL, TbJ, TbR, TaO, Tbw, Tbp, TbC, TbG, TbQ, Tbl;
Chris@82	2583 E Tbx, Ta8, Tbr;
Chris@82	2584 Ta8 = KP707106781 * (Ta6 + Ta7);
Chris@82	2585 Ta9 = Ta5 + Ta8;
Chris@82	2586 TbB = Ta5 - Ta8;
Chris@82	2587 Tbr = KP707106781 * (Ta7 - Ta6);
Chris@82	2588 Tbs = Tbq + Tbr;
Chris@82	2589 TbM = Tbr - Tbq;
Chris@82	2590 {
Chris@82	2591 E Tac, Taf, TbH, TbI;
Chris@82	2592 Tac = FMA(KP923879532, Taa, KP382683432 * Tab);
Chris@82	2593 Taf = FNMS(KP382683432, Tae, KP923879532 * Tad);
Chris@82	2594 Tag = Tac + Taf;
Chris@82	2595 TbL = Taf - Tac;
Chris@82	2596 TbH = TaV - Tba;
Chris@82	2597 TbI = Tbj - Tbi;
Chris@82	2598 TbJ = FNMS(KP555570233, TbI, KP831469612 * TbH);
Chris@82	2599 TbR = FMA(KP555570233, TbH, KP831469612 * TbI);
Chris@82	2600 }
Chris@82	2601 {
Chris@82	2602 E TaE, TaN, Tbn, Tbo;
Chris@82	2603 TaE = Tao + TaD;
Chris@82	2604 TaN = TaL + TaM;
Chris@82	2605 TaO = FMA(KP980785280, TaE, KP195090322 * TaN);
Chris@82	2606 Tbw = FNMS(KP195090322, TaE, KP980785280 * TaN);
Chris@82	2607 Tbn = FNMS(KP382683432, Taa, KP923879532 * Tab);
Chris@82	2608 Tbo = FMA(KP382683432, Tad, KP923879532 * Tae);
Chris@82	2609 Tbp = Tbn + Tbo;
Chris@82	2610 TbC = Tbo - Tbn;
Chris@82	2611 }
Chris@82	2612 {
Chris@82	2613 E TbE, TbF, Tbb, Tbk;
Chris@82	2614 TbE = Tao - TaD;
Chris@82	2615 TbF = TaM - TaL;
Chris@82	2616 TbG = FMA(KP831469612, TbE, KP555570233 * TbF);
Chris@82	2617 TbQ = FNMS(KP555570233, TbE, KP831469612 * TbF);
Chris@82	2618 Tbb = TaV + Tba;
Chris@82	2619 Tbk = Tbi + Tbj;
Chris@82	2620 Tbl = FNMS(KP195090322, Tbk, KP980785280 * Tbb);
Chris@82	2621 Tbx = FMA(KP195090322, Tbb, KP980785280 * Tbk);
Chris@82	2622 }
Chris@82	2623 {
Chris@82	2624 E Tah, Tbm, Tbv, Tby;
Chris@82	2625 Tah = Ta9 + Tag;
Chris@82	2626 Tbm = TaO + Tbl;
Chris@82	2627 Cr[WS(csr, 60)] = Tah - Tbm;
Chris@82	2628 Cr[WS(csr, 4)] = Tah + Tbm;
Chris@82	2629 Tbv = Tbs + Tbp;
Chris@82	2630 Tby = Tbw + Tbx;
Chris@82	2631 Ci[WS(csi, 4)] = Tbv + Tby;
Chris@82	2632 Ci[WS(csi, 60)] = Tby - Tbv;
Chris@82	2633 }
Chris@82	2634 {
Chris@82	2635 E Tbt, Tbu, Tbz, TbA;
Chris@82	2636 Tbt = Tbp - Tbs;
Chris@82	2637 Tbu = Tbl - TaO;
Chris@82	2638 Ci[WS(csi, 28)] = Tbt + Tbu;
Chris@82	2639 Ci[WS(csi, 36)] = Tbu - Tbt;
Chris@82	2640 Tbz = Ta9 - Tag;
Chris@82	2641 TbA = Tbx - Tbw;
Chris@82	2642 Cr[WS(csr, 36)] = Tbz - TbA;
Chris@82	2643 Cr[WS(csr, 28)] = Tbz + TbA;
Chris@82	2644 }
Chris@82	2645 {
Chris@82	2646 E TbD, TbK, TbP, TbS;
Chris@82	2647 TbD = TbB + TbC;
Chris@82	2648 TbK = TbG + TbJ;
Chris@82	2649 Cr[WS(csr, 52)] = TbD - TbK;
Chris@82	2650 Cr[WS(csr, 12)] = TbD + TbK;
Chris@82	2651 TbP = TbM + TbL;
Chris@82	2652 TbS = TbQ + TbR;
Chris@82	2653 Ci[WS(csi, 12)] = TbP + TbS;
Chris@82	2654 Ci[WS(csi, 52)] = TbS - TbP;
Chris@82	2655 }
Chris@82	2656 {
Chris@82	2657 E TbN, TbO, TbT, TbU;
Chris@82	2658 TbN = TbL - TbM;
Chris@82	2659 TbO = TbJ - TbG;
Chris@82	2660 Ci[WS(csi, 20)] = TbN + TbO;
Chris@82	2661 Ci[WS(csi, 44)] = TbO - TbN;
Chris@82	2662 TbT = TbB - TbC;
Chris@82	2663 TbU = TbR - TbQ;
Chris@82	2664 Cr[WS(csr, 44)] = TbT - TbU;
Chris@82	2665 Cr[WS(csr, 20)] = TbT + TbU;
Chris@82	2666 }
Chris@82	2667 }
Chris@82	2668 {
Chris@82	2669 E Tev, Tf7, Tfc, Tfm, Tff, Tfn, TeC, Tfh, TeK, Tf2, TeV, Tf8, TeY, Tfi, TeR;
Chris@82	2670 E Tf3;
Chris@82	2671 {
Chris@82	2672 E Tet, Teu, Tfa, Tfb;
Chris@82	2673 Tet = TcD - TcG;
Chris@82	2674 Teu = TdY - TdX;
Chris@82	2675 Tev = Tet - Teu;
Chris@82	2676 Tf7 = Tet + Teu;
Chris@82	2677 Tfa = TeF + TeE;
Chris@82	2678 Tfb = TeH + TeI;
Chris@82	2679 Tfc = FMA(KP290284677, Tfa, KP956940335 * Tfb);
Chris@82	2680 Tfm = FNMS(KP290284677, Tfb, KP956940335 * Tfa);
Chris@82	2681 }
Chris@82	2682 {
Chris@82	2683 E Tfd, Tfe, Tey, TeB;
Chris@82	2684 Tfd = TeL + TeM;
Chris@82	2685 Tfe = TeP + TeO;
Chris@82	2686 Tff = FNMS(KP290284677, Tfe, KP956940335 * Tfd);
Chris@82	2687 Tfn = FMA(KP956940335, Tfe, KP290284677 * Tfd);
Chris@82	2688 Tey = FMA(KP555570233, Tew, KP831469612 * Tex);
Chris@82	2689 TeB = FNMS(KP555570233, TeA, KP831469612 * Tez);
Chris@82	2690 TeC = Tey - TeB;
Chris@82	2691 Tfh = TeB + Tey;
Chris@82	2692 }
Chris@82	2693 {
Chris@82	2694 E TeG, TeJ, TeT, TeU;
Chris@82	2695 TeG = TeE - TeF;
Chris@82	2696 TeJ = TeH - TeI;
Chris@82	2697 TeK = FMA(KP471396736, TeG, KP881921264 * TeJ);
Chris@82	2698 Tf2 = FNMS(KP471396736, TeJ, KP881921264 * TeG);
Chris@82	2699 TeT = FNMS(KP555570233, Tex, KP831469612 * Tew);
Chris@82	2700 TeU = FMA(KP831469612, TeA, KP555570233 * Tez);
Chris@82	2701 TeV = TeT - TeU;
Chris@82	2702 Tf8 = TeU + TeT;
Chris@82	2703 }
Chris@82	2704 {
Chris@82	2705 E TeW, TeX, TeN, TeQ;
Chris@82	2706 TeW = TcN - TcK;
Chris@82	2707 TeX = TdV - TdU;
Chris@82	2708 TeY = TeW - TeX;
Chris@82	2709 Tfi = TeX + TeW;
Chris@82	2710 TeN = TeL - TeM;
Chris@82	2711 TeQ = TeO - TeP;
Chris@82	2712 TeR = FNMS(KP471396736, TeQ, KP881921264 * TeN);
Chris@82	2713 Tf3 = FMA(KP881921264, TeQ, KP471396736 * TeN);
Chris@82	2714 }
Chris@82	2715 {
Chris@82	2716 E TeD, TeS, Tf1, Tf4;
Chris@82	2717 TeD = Tev + TeC;
Chris@82	2718 TeS = TeK + TeR;
Chris@82	2719 Cr[WS(csr, 54)] = TeD - TeS;
Chris@82	2720 Cr[WS(csr, 10)] = TeD + TeS;
Chris@82	2721 Tf1 = TeY + TeV;
Chris@82	2722 Tf4 = Tf2 + Tf3;
Chris@82	2723 Ci[WS(csi, 10)] = Tf1 + Tf4;
Chris@82	2724 Ci[WS(csi, 54)] = Tf4 - Tf1;
Chris@82	2725 }
Chris@82	2726 {
Chris@82	2727 E TeZ, Tf0, Tf5, Tf6;
Chris@82	2728 TeZ = TeV - TeY;
Chris@82	2729 Tf0 = TeR - TeK;
Chris@82	2730 Ci[WS(csi, 22)] = TeZ + Tf0;
Chris@82	2731 Ci[WS(csi, 42)] = Tf0 - TeZ;
Chris@82	2732 Tf5 = Tev - TeC;
Chris@82	2733 Tf6 = Tf3 - Tf2;
Chris@82	2734 Cr[WS(csr, 42)] = Tf5 - Tf6;
Chris@82	2735 Cr[WS(csr, 22)] = Tf5 + Tf6;
Chris@82	2736 }
Chris@82	2737 {
Chris@82	2738 E Tf9, Tfg, Tfl, Tfo;
Chris@82	2739 Tf9 = Tf7 + Tf8;
Chris@82	2740 Tfg = Tfc + Tff;
Chris@82	2741 Cr[WS(csr, 58)] = Tf9 - Tfg;
Chris@82	2742 Cr[WS(csr, 6)] = Tf9 + Tfg;
Chris@82	2743 Tfl = Tfi + Tfh;
Chris@82	2744 Tfo = Tfm + Tfn;
Chris@82	2745 Ci[WS(csi, 6)] = Tfl + Tfo;
Chris@82	2746 Ci[WS(csi, 58)] = Tfo - Tfl;
Chris@82	2747 }
Chris@82	2748 {
Chris@82	2749 E Tfj, Tfk, Tfp, Tfq;
Chris@82	2750 Tfj = Tfh - Tfi;
Chris@82	2751 Tfk = Tff - Tfc;
Chris@82	2752 Ci[WS(csi, 26)] = Tfj + Tfk;
Chris@82	2753 Ci[WS(csi, 38)] = Tfk - Tfj;
Chris@82	2754 Tfp = Tf7 - Tf8;
Chris@82	2755 Tfq = Tfn - Tfm;
Chris@82	2756 Cr[WS(csr, 38)] = Tfp - Tfq;
Chris@82	2757 Cr[WS(csr, 26)] = Tfp + Tfq;
Chris@82	2758 }
Chris@82	2759 }
Chris@82	2760 {
Chris@82	2761 E TcP, Te9, Tee, Teo, Teh, Tep, Td8, Tej, Tdu, Te4, TdT, Tea, Te0, Tek, TdP;
Chris@82	2762 E Te5;
Chris@82	2763 {
Chris@82	2764 E TcH, TcO, Tec, Ted;
Chris@82	2765 TcH = TcD + TcG;
Chris@82	2766 TcO = TcK + TcN;
Chris@82	2767 TcP = TcH + TcO;
Chris@82	2768 Te9 = TcH - TcO;
Chris@82	2769 Tec = Tde - Tdl;
Chris@82	2770 Ted = Tds - Tdp;
Chris@82	2771 Tee = FMA(KP773010453, Tec, KP634393284 * Ted);
Chris@82	2772 Teo = FNMS(KP634393284, Tec, KP773010453 * Ted);
Chris@82	2773 }
Chris@82	2774 {
Chris@82	2775 E Tef, Teg, TcY, Td7;
Chris@82	2776 Tef = Tdz - TdG;
Chris@82	2777 Teg = TdN - TdK;
Chris@82	2778 Teh = FNMS(KP634393284, Teg, KP773010453 * Tef);
Chris@82	2779 Tep = FMA(KP634393284, Tef, KP773010453 * Teg);
Chris@82	2780 TcY = FMA(KP980785280, TcU, KP195090322 * TcX);
Chris@82	2781 Td7 = FNMS(KP195090322, Td6, KP980785280 * Td3);
Chris@82	2782 Td8 = TcY + Td7;
Chris@82	2783 Tej = Td7 - TcY;
Chris@82	2784 }
Chris@82	2785 {
Chris@82	2786 E Tdm, Tdt, TdR, TdS;
Chris@82	2787 Tdm = Tde + Tdl;
Chris@82	2788 Tdt = Tdp + Tds;
Chris@82	2789 Tdu = FMA(KP995184726, Tdm, KP098017140 * Tdt);
Chris@82	2790 Te4 = FNMS(KP098017140, Tdm, KP995184726 * Tdt);
Chris@82	2791 TdR = FNMS(KP195090322, TcU, KP980785280 * TcX);
Chris@82	2792 TdS = FMA(KP195090322, Td3, KP980785280 * Td6);
Chris@82	2793 TdT = TdR + TdS;
Chris@82	2794 Tea = TdS - TdR;
Chris@82	2795 }
Chris@82	2796 {
Chris@82	2797 E TdW, TdZ, TdH, TdO;
Chris@82	2798 TdW = TdU + TdV;
Chris@82	2799 TdZ = TdX + TdY;
Chris@82	2800 Te0 = TdW + TdZ;
Chris@82	2801 Tek = TdZ - TdW;
Chris@82	2802 TdH = Tdz + TdG;
Chris@82	2803 TdO = TdK + TdN;
Chris@82	2804 TdP = FNMS(KP098017140, TdO, KP995184726 * TdH);
Chris@82	2805 Te5 = FMA(KP098017140, TdH, KP995184726 * TdO);
Chris@82	2806 }
Chris@82	2807 {
Chris@82	2808 E Td9, TdQ, Te3, Te6;
Chris@82	2809 Td9 = TcP + Td8;
Chris@82	2810 TdQ = Tdu + TdP;
Chris@82	2811 Cr[WS(csr, 62)] = Td9 - TdQ;
Chris@82	2812 Cr[WS(csr, 2)] = Td9 + TdQ;
Chris@82	2813 Te3 = Te0 + TdT;
Chris@82	2814 Te6 = Te4 + Te5;
Chris@82	2815 Ci[WS(csi, 2)] = Te3 + Te6;
Chris@82	2816 Ci[WS(csi, 62)] = Te6 - Te3;
Chris@82	2817 }
Chris@82	2818 {
Chris@82	2819 E Te1, Te2, Te7, Te8;
Chris@82	2820 Te1 = TdT - Te0;
Chris@82	2821 Te2 = TdP - Tdu;
Chris@82	2822 Ci[WS(csi, 30)] = Te1 + Te2;
Chris@82	2823 Ci[WS(csi, 34)] = Te2 - Te1;
Chris@82	2824 Te7 = TcP - Td8;
Chris@82	2825 Te8 = Te5 - Te4;
Chris@82	2826 Cr[WS(csr, 34)] = Te7 - Te8;
Chris@82	2827 Cr[WS(csr, 30)] = Te7 + Te8;
Chris@82	2828 }
Chris@82	2829 {
Chris@82	2830 E Teb, Tei, Ten, Teq;
Chris@82	2831 Teb = Te9 + Tea;
Chris@82	2832 Tei = Tee + Teh;
Chris@82	2833 Cr[WS(csr, 50)] = Teb - Tei;
Chris@82	2834 Cr[WS(csr, 14)] = Teb + Tei;
Chris@82	2835 Ten = Tek + Tej;
Chris@82	2836 Teq = Teo + Tep;
Chris@82	2837 Ci[WS(csi, 14)] = Ten + Teq;
Chris@82	2838 Ci[WS(csi, 50)] = Teq - Ten;
Chris@82	2839 }
Chris@82	2840 {
Chris@82	2841 E Tel, Tem, Ter, Tes;
Chris@82	2842 Tel = Tej - Tek;
Chris@82	2843 Tem = Teh - Tee;
Chris@82	2844 Ci[WS(csi, 18)] = Tel + Tem;
Chris@82	2845 Ci[WS(csi, 46)] = Tem - Tel;
Chris@82	2846 Ter = Te9 - Tea;
Chris@82	2847 Tes = Tep - Teo;
Chris@82	2848 Cr[WS(csr, 46)] = Ter - Tes;
Chris@82	2849 Cr[WS(csr, 18)] = Ter + Tes;
Chris@82	2850 }
Chris@82	2851 }
Chris@82	2852 {
Chris@82	2853 E T6v, T77, T6C, T7h, T6Y, T7i, T6V, T78, T6R, T7n, T73, T7f, T6K, T7m, T72;
Chris@82	2854 E T7c;
Chris@82	2855 {
Chris@82	2856 E T6t, T6u, T6T, T6U;
Chris@82	2857 T6t = T27 - T2e;
Chris@82	2858 T6u = T5Y - T5X;
Chris@82	2859 T6v = T6t - T6u;
Chris@82	2860 T77 = T6t + T6u;
Chris@82	2861 {
Chris@82	2862 E T6y, T6B, T6W, T6X;
Chris@82	2863 T6y = FMA(KP773010453, T6w, KP634393284 * T6x);
Chris@82	2864 T6B = FNMS(KP634393284, T6A, KP773010453 * T6z);
Chris@82	2865 T6C = T6y - T6B;
Chris@82	2866 T7h = T6B + T6y;
Chris@82	2867 T6W = T2x - T2o;
Chris@82	2868 T6X = T5V - T5S;
Chris@82	2869 T6Y = T6W - T6X;
Chris@82	2870 T7i = T6X + T6W;
Chris@82	2871 }
Chris@82	2872 T6T = FNMS(KP634393284, T6w, KP773010453 * T6x);
Chris@82	2873 T6U = FMA(KP634393284, T6z, KP773010453 * T6A);
Chris@82	2874 T6V = T6T - T6U;
Chris@82	2875 T78 = T6U + T6T;
Chris@82	2876 {
Chris@82	2877 E T6N, T7d, T6Q, T7e, T6M, T6O;
Chris@82	2878 T6M = T5I - T5H;
Chris@82	2879 T6N = T6L - T6M;
Chris@82	2880 T7d = T6L + T6M;
Chris@82	2881 T6O = T5v - T5e;
Chris@82	2882 T6Q = T6O - T6P;
Chris@82	2883 T7e = T6P + T6O;
Chris@82	2884 T6R = FNMS(KP427555093, T6Q, KP903989293 * T6N);
Chris@82	2885 T7n = FMA(KP941544065, T7e, KP336889853 * T7d);
Chris@82	2886 T73 = FMA(KP903989293, T6Q, KP427555093 * T6N);
Chris@82	2887 T7f = FNMS(KP336889853, T7e, KP941544065 * T7d);
Chris@82	2888 }
Chris@82	2889 {
Chris@82	2890 E T6G, T7a, T6J, T7b, T6E, T6I;
Chris@82	2891 T6E = T4g - T3Z;
Chris@82	2892 T6G = T6E - T6F;
Chris@82	2893 T7a = T6F + T6E;
Chris@82	2894 T6I = T4t - T4s;
Chris@82	2895 T6J = T6H - T6I;
Chris@82	2896 T7b = T6H + T6I;
Chris@82	2897 T6K = FMA(KP427555093, T6G, KP903989293 * T6J);
Chris@82	2898 T7m = FNMS(KP336889853, T7b, KP941544065 * T7a);
Chris@82	2899 T72 = FNMS(KP427555093, T6J, KP903989293 * T6G);
Chris@82	2900 T7c = FMA(KP336889853, T7a, KP941544065 * T7b);
Chris@82	2901 }
Chris@82	2902 }
Chris@82	2903 {
Chris@82	2904 E T6D, T6S, T71, T74;
Chris@82	2905 T6D = T6v + T6C;
Chris@82	2906 T6S = T6K + T6R;
Chris@82	2907 Cr[WS(csr, 55)] = T6D - T6S;
Chris@82	2908 Cr[WS(csr, 9)] = T6D + T6S;
Chris@82	2909 T71 = T6Y + T6V;
Chris@82	2910 T74 = T72 + T73;
Chris@82	2911 Ci[WS(csi, 9)] = T71 + T74;
Chris@82	2912 Ci[WS(csi, 55)] = T74 - T71;
Chris@82	2913 }
Chris@82	2914 {
Chris@82	2915 E T6Z, T70, T75, T76;
Chris@82	2916 T6Z = T6V - T6Y;
Chris@82	2917 T70 = T6R - T6K;
Chris@82	2918 Ci[WS(csi, 23)] = T6Z + T70;
Chris@82	2919 Ci[WS(csi, 41)] = T70 - T6Z;
Chris@82	2920 T75 = T6v - T6C;
Chris@82	2921 T76 = T73 - T72;
Chris@82	2922 Cr[WS(csr, 41)] = T75 - T76;
Chris@82	2923 Cr[WS(csr, 23)] = T75 + T76;
Chris@82	2924 }
Chris@82	2925 {
Chris@82	2926 E T79, T7g, T7l, T7o;
Chris@82	2927 T79 = T77 + T78;
Chris@82	2928 T7g = T7c + T7f;
Chris@82	2929 Cr[WS(csr, 57)] = T79 - T7g;
Chris@82	2930 Cr[WS(csr, 7)] = T79 + T7g;
Chris@82	2931 T7l = T7i + T7h;
Chris@82	2932 T7o = T7m + T7n;
Chris@82	2933 Ci[WS(csi, 7)] = T7l + T7o;
Chris@82	2934 Ci[WS(csi, 57)] = T7o - T7l;
Chris@82	2935 }
Chris@82	2936 {
Chris@82	2937 E T7j, T7k, T7p, T7q;
Chris@82	2938 T7j = T7h - T7i;
Chris@82	2939 T7k = T7f - T7c;
Chris@82	2940 Ci[WS(csi, 25)] = T7j + T7k;
Chris@82	2941 Ci[WS(csi, 39)] = T7k - T7j;
Chris@82	2942 T7p = T77 - T78;
Chris@82	2943 T7q = T7n - T7m;
Chris@82	2944 Cr[WS(csr, 39)] = T7p - T7q;
Chris@82	2945 Cr[WS(csr, 25)] = T7p + T7q;
Chris@82	2946 }
Chris@82	2947 }
Chris@82	2948 {
Chris@82	2949 E T99, T9L, T9g, T9V, T9C, T9W, T9z, T9M, T9v, Ta1, T9H, T9T, T9o, Ta0, T9G;
Chris@82	2950 E T9Q;
Chris@82	2951 {
Chris@82	2952 E T97, T98, T9x, T9y;
Chris@82	2953 T97 = T7r - T7s;
Chris@82	2954 T98 = T8C - T8B;
Chris@82	2955 T99 = T97 - T98;
Chris@82	2956 T9L = T97 + T98;
Chris@82	2957 {
Chris@82	2958 E T9c, T9f, T9A, T9B;
Chris@82	2959 T9c = FMA(KP471396736, T9a, KP881921264 * T9b);
Chris@82	2960 T9f = FNMS(KP471396736, T9e, KP881921264 * T9d);
Chris@82	2961 T9g = T9c - T9f;
Chris@82	2962 T9V = T9f + T9c;
Chris@82	2963 T9A = T7z - T7w;
Chris@82	2964 T9B = T8z - T8y;
Chris@82	2965 T9C = T9A - T9B;
Chris@82	2966 T9W = T9B + T9A;
Chris@82	2967 }
Chris@82	2968 T9x = FNMS(KP471396736, T9b, KP881921264 * T9a);
Chris@82	2969 T9y = FMA(KP881921264, T9e, KP471396736 * T9d);
Chris@82	2970 T9z = T9x - T9y;
Chris@82	2971 T9M = T9y + T9x;
Chris@82	2972 {
Chris@82	2973 E T9r, T9R, T9u, T9S, T9q, T9s;
Chris@82	2974 T9q = T8q - T8p;
Chris@82	2975 T9r = T9p - T9q;
Chris@82	2976 T9R = T9p + T9q;
Chris@82	2977 T9s = T8j - T8g;
Chris@82	2978 T9u = T9s - T9t;
Chris@82	2979 T9S = T9t + T9s;
Chris@82	2980 T9v = FNMS(KP514102744, T9u, KP857728610 * T9r);
Chris@82	2981 Ta1 = FMA(KP970031253, T9S, KP242980179 * T9R);
Chris@82	2982 T9H = FMA(KP857728610, T9u, KP514102744 * T9r);
Chris@82	2983 T9T = FNMS(KP242980179, T9S, KP970031253 * T9R);
Chris@82	2984 }
Chris@82	2985 {
Chris@82	2986 E T9k, T9O, T9n, T9P, T9i, T9m;
Chris@82	2987 T9i = T80 - T7X;
Chris@82	2988 T9k = T9i - T9j;
Chris@82	2989 T9O = T9j + T9i;
Chris@82	2990 T9m = T87 - T86;
Chris@82	2991 T9n = T9l - T9m;
Chris@82	2992 T9P = T9l + T9m;
Chris@82	2993 T9o = FMA(KP514102744, T9k, KP857728610 * T9n);
Chris@82	2994 Ta0 = FNMS(KP242980179, T9P, KP970031253 * T9O);
Chris@82	2995 T9G = FNMS(KP514102744, T9n, KP857728610 * T9k);
Chris@82	2996 T9Q = FMA(KP242980179, T9O, KP970031253 * T9P);
Chris@82	2997 }
Chris@82	2998 }
Chris@82	2999 {
Chris@82	3000 E T9h, T9w, T9F, T9I;
Chris@82	3001 T9h = T99 + T9g;
Chris@82	3002 T9w = T9o + T9v;
Chris@82	3003 Cr[WS(csr, 53)] = T9h - T9w;
Chris@82	3004 Cr[WS(csr, 11)] = T9h + T9w;
Chris@82	3005 T9F = T9C + T9z;
Chris@82	3006 T9I = T9G + T9H;
Chris@82	3007 Ci[WS(csi, 11)] = T9F + T9I;
Chris@82	3008 Ci[WS(csi, 53)] = T9I - T9F;
Chris@82	3009 }
Chris@82	3010 {
Chris@82	3011 E T9D, T9E, T9J, T9K;
Chris@82	3012 T9D = T9z - T9C;
Chris@82	3013 T9E = T9v - T9o;
Chris@82	3014 Ci[WS(csi, 21)] = T9D + T9E;
Chris@82	3015 Ci[WS(csi, 43)] = T9E - T9D;
Chris@82	3016 T9J = T99 - T9g;
Chris@82	3017 T9K = T9H - T9G;
Chris@82	3018 Cr[WS(csr, 43)] = T9J - T9K;
Chris@82	3019 Cr[WS(csr, 21)] = T9J + T9K;
Chris@82	3020 }
Chris@82	3021 {
Chris@82	3022 E T9N, T9U, T9Z, Ta2;
Chris@82	3023 T9N = T9L + T9M;
Chris@82	3024 T9U = T9Q + T9T;
Chris@82	3025 Cr[WS(csr, 59)] = T9N - T9U;
Chris@82	3026 Cr[WS(csr, 5)] = T9N + T9U;
Chris@82	3027 T9Z = T9W + T9V;
Chris@82	3028 Ta2 = Ta0 + Ta1;
Chris@82	3029 Ci[WS(csi, 5)] = T9Z + Ta2;
Chris@82	3030 Ci[WS(csi, 59)] = Ta2 - T9Z;
Chris@82	3031 }
Chris@82	3032 {
Chris@82	3033 E T9X, T9Y, Ta3, Ta4;
Chris@82	3034 T9X = T9V - T9W;
Chris@82	3035 T9Y = T9T - T9Q;
Chris@82	3036 Ci[WS(csi, 27)] = T9X + T9Y;
Chris@82	3037 Ci[WS(csi, 37)] = T9Y - T9X;
Chris@82	3038 Ta3 = T9L - T9M;
Chris@82	3039 Ta4 = Ta1 - Ta0;
Chris@82	3040 Cr[WS(csr, 37)] = Ta3 - Ta4;
Chris@82	3041 Cr[WS(csr, 27)] = Ta3 + Ta4;
Chris@82	3042 }
Chris@82	3043 }
Chris@82	3044 {
Chris@82	3045 E T2z, T69, T3g, T6j, T60, T6k, T5P, T6a, T5L, T6p, T65, T6h, T4w, T6o, T64;
Chris@82	3046 E T6e;
Chris@82	3047 {
Chris@82	3048 E T2f, T2y, T5N, T5O;
Chris@82	3049 T2f = T27 + T2e;
Chris@82	3050 T2y = T2o + T2x;
Chris@82	3051 T2z = T2f + T2y;
Chris@82	3052 T69 = T2f - T2y;
Chris@82	3053 {
Chris@82	3054 E T2U, T3f, T5W, T5Z;
Chris@82	3055 T2U = FMA(KP098017140, T2M, KP995184726 * T2T);
Chris@82	3056 T3f = FNMS(KP098017140, T3e, KP995184726 * T37);
Chris@82	3057 T3g = T2U + T3f;
Chris@82	3058 T6j = T3f - T2U;
Chris@82	3059 T5W = T5S + T5V;
Chris@82	3060 T5Z = T5X + T5Y;
Chris@82	3061 T60 = T5W + T5Z;
Chris@82	3062 T6k = T5Z - T5W;
Chris@82	3063 }
Chris@82	3064 T5N = FNMS(KP098017140, T2T, KP995184726 * T2M);
Chris@82	3065 T5O = FMA(KP995184726, T3e, KP098017140 * T37);
Chris@82	3066 T5P = T5N + T5O;
Chris@82	3067 T6a = T5O - T5N;
Chris@82	3068 {
Chris@82	3069 E T5x, T6f, T5K, T6g, T5w, T5J;
Chris@82	3070 T5w = T5e + T5v;
Chris@82	3071 T5x = T4X + T5w;
Chris@82	3072 T6f = T4X - T5w;
Chris@82	3073 T5J = T5H + T5I;
Chris@82	3074 T5K = T5G + T5J;
Chris@82	3075 T6g = T5J - T5G;
Chris@82	3076 T5L = FNMS(KP049067674, T5K, KP998795456 * T5x);
Chris@82	3077 T6p = FMA(KP671558954, T6f, KP740951125 * T6g);
Chris@82	3078 T65 = FMA(KP049067674, T5x, KP998795456 * T5K);
Chris@82	3079 T6h = FNMS(KP671558954, T6g, KP740951125 * T6f);
Chris@82	3080 }
Chris@82	3081 {
Chris@82	3082 E T4i, T6c, T4v, T6d, T4h, T4u;
Chris@82	3083 T4h = T3Z + T4g;
Chris@82	3084 T4i = T3I + T4h;
Chris@82	3085 T6c = T3I - T4h;
Chris@82	3086 T4u = T4s + T4t;
Chris@82	3087 T4v = T4r + T4u;
Chris@82	3088 T6d = T4u - T4r;
Chris@82	3089 T4w = FMA(KP998795456, T4i, KP049067674 * T4v);
Chris@82	3090 T6o = FNMS(KP671558954, T6c, KP740951125 * T6d);
Chris@82	3091 T64 = FNMS(KP049067674, T4i, KP998795456 * T4v);
Chris@82	3092 T6e = FMA(KP740951125, T6c, KP671558954 * T6d);
Chris@82	3093 }
Chris@82	3094 }
Chris@82	3095 {
Chris@82	3096 E T3h, T5M, T63, T66;
Chris@82	3097 T3h = T2z + T3g;
Chris@82	3098 T5M = T4w + T5L;
Chris@82	3099 Cr[WS(csr, 63)] = T3h - T5M;
Chris@82	3100 Cr[WS(csr, 1)] = T3h + T5M;
Chris@82	3101 T63 = T60 + T5P;
Chris@82	3102 T66 = T64 + T65;
Chris@82	3103 Ci[WS(csi, 1)] = T63 + T66;
Chris@82	3104 Ci[WS(csi, 63)] = T66 - T63;
Chris@82	3105 }
Chris@82	3106 {
Chris@82	3107 E T61, T62, T67, T68;
Chris@82	3108 T61 = T5P - T60;
Chris@82	3109 T62 = T5L - T4w;
Chris@82	3110 Ci[WS(csi, 31)] = T61 + T62;
Chris@82	3111 Ci[WS(csi, 33)] = T62 - T61;
Chris@82	3112 T67 = T2z - T3g;
Chris@82	3113 T68 = T65 - T64;
Chris@82	3114 Cr[WS(csr, 33)] = T67 - T68;
Chris@82	3115 Cr[WS(csr, 31)] = T67 + T68;
Chris@82	3116 }
Chris@82	3117 {
Chris@82	3118 E T6b, T6i, T6n, T6q;
Chris@82	3119 T6b = T69 + T6a;
Chris@82	3120 T6i = T6e + T6h;
Chris@82	3121 Cr[WS(csr, 49)] = T6b - T6i;
Chris@82	3122 Cr[WS(csr, 15)] = T6b + T6i;
Chris@82	3123 T6n = T6k + T6j;
Chris@82	3124 T6q = T6o + T6p;
Chris@82	3125 Ci[WS(csi, 15)] = T6n + T6q;
Chris@82	3126 Ci[WS(csi, 49)] = T6q - T6n;
Chris@82	3127 }
Chris@82	3128 {
Chris@82	3129 E T6l, T6m, T6r, T6s;
Chris@82	3130 T6l = T6j - T6k;
Chris@82	3131 T6m = T6h - T6e;
Chris@82	3132 Ci[WS(csi, 17)] = T6l + T6m;
Chris@82	3133 Ci[WS(csi, 47)] = T6m - T6l;
Chris@82	3134 T6r = T69 - T6a;
Chris@82	3135 T6s = T6p - T6o;
Chris@82	3136 Cr[WS(csr, 47)] = T6r - T6s;
Chris@82	3137 Cr[WS(csr, 17)] = T6r + T6s;
Chris@82	3138 }
Chris@82	3139 }
Chris@82	3140 {
Chris@82	3141 E T7B, T8N, T7Q, T8X, T8E, T8Y, T8x, T8O, T8t, T93, T8J, T8V, T8a, T92, T8I;
Chris@82	3142 E T8S;
Chris@82	3143 {
Chris@82	3144 E T7t, T7A, T8v, T8w;
Chris@82	3145 T7t = T7r + T7s;
Chris@82	3146 T7A = T7w + T7z;
Chris@82	3147 T7B = T7t + T7A;
Chris@82	3148 T8N = T7t - T7A;
Chris@82	3149 {
Chris@82	3150 E T7I, T7P, T8A, T8D;
Chris@82	3151 T7I = FMA(KP956940335, T7E, KP290284677 * T7H);
Chris@82	3152 T7P = FNMS(KP290284677, T7O, KP956940335 * T7L);
Chris@82	3153 T7Q = T7I + T7P;
Chris@82	3154 T8X = T7P - T7I;
Chris@82	3155 T8A = T8y + T8z;
Chris@82	3156 T8D = T8B + T8C;
Chris@82	3157 T8E = T8A + T8D;
Chris@82	3158 T8Y = T8D - T8A;
Chris@82	3159 }
Chris@82	3160 T8v = FNMS(KP290284677, T7E, KP956940335 * T7H);
Chris@82	3161 T8w = FMA(KP290284677, T7L, KP956940335 * T7O);
Chris@82	3162 T8x = T8v + T8w;
Chris@82	3163 T8O = T8w - T8v;
Chris@82	3164 {
Chris@82	3165 E T8l, T8T, T8s, T8U, T8k, T8r;
Chris@82	3166 T8k = T8g + T8j;
Chris@82	3167 T8l = T8d + T8k;
Chris@82	3168 T8T = T8d - T8k;
Chris@82	3169 T8r = T8p + T8q;
Chris@82	3170 T8s = T8o + T8r;
Chris@82	3171 T8U = T8r - T8o;
Chris@82	3172 T8t = FNMS(KP146730474, T8s, KP989176509 * T8l);
Chris@82	3173 T93 = FMA(KP595699304, T8T, KP803207531 * T8U);
Chris@82	3174 T8J = FMA(KP146730474, T8l, KP989176509 * T8s);
Chris@82	3175 T8V = FNMS(KP595699304, T8U, KP803207531 * T8T);
Chris@82	3176 }
Chris@82	3177 {
Chris@82	3178 E T82, T8Q, T89, T8R, T81, T88;
Chris@82	3179 T81 = T7X + T80;
Chris@82	3180 T82 = T7U + T81;
Chris@82	3181 T8Q = T7U - T81;
Chris@82	3182 T88 = T86 + T87;
Chris@82	3183 T89 = T85 + T88;
Chris@82	3184 T8R = T88 - T85;
Chris@82	3185 T8a = FMA(KP989176509, T82, KP146730474 * T89);
Chris@82	3186 T92 = FNMS(KP595699304, T8Q, KP803207531 * T8R);
Chris@82	3187 T8I = FNMS(KP146730474, T82, KP989176509 * T89);
Chris@82	3188 T8S = FMA(KP803207531, T8Q, KP595699304 * T8R);
Chris@82	3189 }
Chris@82	3190 }
Chris@82	3191 {
Chris@82	3192 E T7R, T8u, T8H, T8K;
Chris@82	3193 T7R = T7B + T7Q;
Chris@82	3194 T8u = T8a + T8t;
Chris@82	3195 Cr[WS(csr, 61)] = T7R - T8u;
Chris@82	3196 Cr[WS(csr, 3)] = T7R + T8u;
Chris@82	3197 T8H = T8E + T8x;
Chris@82	3198 T8K = T8I + T8J;
Chris@82	3199 Ci[WS(csi, 3)] = T8H + T8K;
Chris@82	3200 Ci[WS(csi, 61)] = T8K - T8H;
Chris@82	3201 }
Chris@82	3202 {
Chris@82	3203 E T8F, T8G, T8L, T8M;
Chris@82	3204 T8F = T8x - T8E;
Chris@82	3205 T8G = T8t - T8a;
Chris@82	3206 Ci[WS(csi, 29)] = T8F + T8G;
Chris@82	3207 Ci[WS(csi, 35)] = T8G - T8F;
Chris@82	3208 T8L = T7B - T7Q;
Chris@82	3209 T8M = T8J - T8I;
Chris@82	3210 Cr[WS(csr, 35)] = T8L - T8M;
Chris@82	3211 Cr[WS(csr, 29)] = T8L + T8M;
Chris@82	3212 }
Chris@82	3213 {
Chris@82	3214 E T8P, T8W, T91, T94;
Chris@82	3215 T8P = T8N + T8O;
Chris@82	3216 T8W = T8S + T8V;
Chris@82	3217 Cr[WS(csr, 51)] = T8P - T8W;
Chris@82	3218 Cr[WS(csr, 13)] = T8P + T8W;
Chris@82	3219 T91 = T8Y + T8X;
Chris@82	3220 T94 = T92 + T93;
Chris@82	3221 Ci[WS(csi, 13)] = T91 + T94;
Chris@82	3222 Ci[WS(csi, 51)] = T94 - T91;
Chris@82	3223 }
Chris@82	3224 {
Chris@82	3225 E T8Z, T90, T95, T96;
Chris@82	3226 T8Z = T8X - T8Y;
Chris@82	3227 T90 = T8V - T8S;
Chris@82	3228 Ci[WS(csi, 19)] = T8Z + T90;
Chris@82	3229 Ci[WS(csi, 45)] = T90 - T8Z;
Chris@82	3230 T95 = T8N - T8O;
Chris@82	3231 T96 = T93 - T92;
Chris@82	3232 Cr[WS(csr, 45)] = T95 - T96;
Chris@82	3233 Cr[WS(csr, 19)] = T95 + T96;
Chris@82	3234 }
Chris@82	3235 }
Chris@82	3236 }
Chris@82	3237 }
Chris@82	3238 }
Chris@82	3239
Chris@82	3240 static const kr2c_desc desc = { 128, "r2cf_128", {812, 186, 144, 0}, &GENUS };
Chris@82	3241
Chris@82	3242 void X(codelet_r2cf_128) (planner *p) {
Chris@82	3243 X(kr2c_register) (p, r2cf_128, &desc);
Chris@82	3244 }
Chris@82	3245
Chris@82	3246 #endif

Mercurial > hg > sv-dependency-builds

annotate src/fftw-3.3.8/rdft/scalar/r2cf/r2cf_128.c @ 83:ae30d91d2ffe