hb_64.c

fftw-3.0.1

	       }
	       {
		    E T4p, T4q, T4s, T4t;
		    T4p = iio[-WS(ios, 31)];
		    T4q = rio[WS(ios, 63)];
		    T4r = T4p - T4q;
		    T4s = iio[-WS(ios, 15)];
		    T4t = rio[WS(ios, 47)];
		    T4u = T4s - T4t;
		    T7Q = T4r + T4u;
		    Tbh = T4s + T4t;
		    TaR = T4q + T4p;
	       }
	       {
		    E TaJ, TaK, TaM, TaN;
		    {
			 E T1E, T1F, T3T, T3U;
			 T1E = rio[WS(ios, 7)];
			 T1F = iio[-WS(ios, 39)];
			 T1G = T1E + T1F;
			 TaJ = T1E - T1F;
			 T3T = iio[-WS(ios, 7)];
			 T3U = rio[WS(ios, 39)];
			 T3V = T3T - T3U;
			 TaK = T3T + T3U;
		    }
		    {
			 E T1H, T1I, T3W, T3X;
			 T1H = iio[-WS(ios, 55)];
			 T1I = rio[WS(ios, 23)];
			 T1J = T1H + T1I;
			 TaM = T1H - T1I;
			 T3W = iio[-WS(ios, 23)];
			 T3X = rio[WS(ios, 55)];
			 T3Y = T3W - T3X;
			 TaN = T3X + T3W;
		    }
		    T1K = T1G + T1J;
		    T7R = T3V + T3Y;
		    Tbe = TaK - TaJ;
		    Tbd = TaM - TaN;
		    TaO = TaM + TaN;
		    TaL = TaJ + TaK;
	       }
	       {
		    E T1O, TaX, T44, TaV, T1R, TaU, T47, TaY;
		    {
			 E T1M, T1N, T42, T43;
			 T1M = rio[WS(ios, 3)];
			 T1N = iio[-WS(ios, 35)];
			 T1O = T1M + T1N;
			 TaX = T1M - T1N;
			 T42 = iio[-WS(ios, 3)];
			 T43 = rio[WS(ios, 35)];
			 T44 = T42 - T43;
			 TaV = T42 + T43;
		    }
		    {
			 E T1P, T1Q, T45, T46;
			 T1P = rio[WS(ios, 19)];
			 T1Q = iio[-WS(ios, 51)];
			 T1R = T1P + T1Q;
			 TaU = T1P - T1Q;
			 T45 = iio[-WS(ios, 19)];
			 T46 = rio[WS(ios, 51)];
			 T47 = T45 - T46;
			 TaY = T45 + T46;
		    }
		    T1S = T1O + T1R;
		    TfV = TaV - TaU;
		    TfW = TaX + TaY;
		    T41 = T1O - T1R;
		    T48 = T44 - T47;
		    TaW = TaU + TaV;
		    TaZ = TaX - TaY;
		    T7L = T44 + T47;
	       }
	       {
		    E T1V, Tb4, T4d, Tb2, T1Y, Tb1, T4g, Tb5;
		    {
			 E T1T, T1U, T4b, T4c;
			 T1T = iio[-WS(ios, 59)];
			 T1U = rio[WS(ios, 27)];
			 T1V = T1T + T1U;
			 Tb4 = T1T - T1U;
			 T4b = iio[-WS(ios, 27)];
			 T4c = rio[WS(ios, 59)];
			 T4d = T4b - T4c;
			 Tb2 = T4c + T4b;
		    }
		    {
			 E T1W, T1X, T4e, T4f;
			 T1W = rio[WS(ios, 11)];
			 T1X = iio[-WS(ios, 43)];
			 T1Y = T1W + T1X;
			 Tb1 = T1W - T1X;
			 T4e = iio[-WS(ios, 11)];
			 T4f = rio[WS(ios, 43)];
			 T4g = T4e - T4f;
			 Tb5 = T4e + T4f;
		    }
		    T1Z = T1V + T1Y;
		    TfY = Tb4 + Tb5;
		    TfZ = Tb1 + Tb2;
		    T4a = T1V - T1Y;
		    T4h = T4d - T4g;
		    Tb3 = Tb1 - Tb2;
		    Tb6 = Tb4 - Tb5;
		    T7M = T4g + T4d;
	       }
	       T1L = T1D + T1K;
	       T20 = T1S + T1Z;
	       T9c = T1L - T20;
	       T9d = T7R + T7Q;
	       T9e = T7L + T7M;
	       T9f = T9d - T9e;
	       {
		    E T3S, T3Z, TfX, Tg0;
		    T3S = T1z - T1C;
		    T3Z = T3V - T3Y;
		    T40 = T3S + T3Z;
		    T66 = T3S - T3Z;
		    TfX = FNMS(KP382683432, TfW, KP923879532 * TfV);
		    Tg0 = FNMS(KP923879532, TfZ, KP382683432 * TfY);
		    Tg1 = TfX + Tg0;
		    Thu = Tg0 - TfX;
	       }
	       {
		    E Tg6, Tg7, Tg3, Tg4;
		    Tg6 = KP707106781 * (TaL + TaO);
		    Tg7 = Tbg + Tbh;
		    Tg8 = Tg6 + Tg7;
		    Thv = Tg7 - Tg6;
		    Tg3 = FMA(KP382683432, TfV, KP923879532 * TfW);
		    Tg4 = FMA(KP923879532, TfY, KP382683432 * TfZ);
		    Tg5 = Tg3 + Tg4;
		    Thr = Tg3 - Tg4;
	       }
	       {
		    E T4l, T4m, T49, T4i;
		    T4l = T48 - T41;
		    T4m = T4a + T4h;
		    T4n = KP707106781 * (T4l + T4m);
		    T67 = KP707106781 * (T4m - T4l);
		    T49 = T41 + T48;
		    T4i = T4a - T4h;
		    T4j = KP707106781 * (T49 + T4i);
		    T69 = KP707106781 * (T49 - T4i);
	       }
	       {
		    E T4o, T4v, TaP, TaS;
		    T4o = T1J - T1G;
		    T4v = T4r - T4u;
		    T4w = T4o + T4v;
		    T6a = T4v - T4o;
		    TaP = KP707106781 * (TaL - TaO);
		    TaS = TaQ - TaR;
		    TaT = TaP + TaS;
		    TdW = TaS - TaP;
	       }
	       {
		    E Tb0, Tb7, TfS, TfT;
		    Tb0 = FMA(KP923879532, TaW, KP382683432 * TaZ);
		    Tb7 = FNMS(KP382683432, Tb6, KP923879532 * Tb3);
		    Tb8 = Tb0 + Tb7;
		    TdZ = Tb0 - Tb7;
		    TfS = KP707106781 * (Tbe + Tbd);
		    TfT = TaQ + TaR;
		    TfU = TfS - TfT;
		    Ths = TfS + TfT;
	       }
	       {
		    E T7K, T7N, T7P, T7S;
		    T7K = T1D - T1K;
		    T7N = T7L - T7M;
		    T7O = T7K + T7N;
		    T8y = T7K - T7N;
		    T7P = T1Z - T1S;
		    T7S = T7Q - T7R;
		    T7T = T7P + T7S;
		    T8z = T7S - T7P;
	       }
	       {
		    E Tba, Tbb, Tbf, Tbi;
		    Tba = FNMS(KP382683432, TaW, KP923879532 * TaZ);
		    Tbb = FMA(KP923879532, Tb6, KP382683432 * Tb3);
		    Tbc = Tba + Tbb;
		    TdX = Tbb - Tba;
		    Tbf = KP707106781 * (Tbd - Tbe);
		    Tbi = Tbg - Tbh;
		    Tbj = Tbf + Tbi;
		    Te0 = Tbi - Tbf;
	       }
	  }
	  {
	       E T14, T17, T18, Tax, Tas, T3K, T3N, T7F, Tay, Tat, T1b, T3e, T1e, T3h, T1f;
	       E T7G, TaB, TaA, Taq, Tan, T1n, TfC, TfD, T3k, T3r, Ta8, Tab, T7A, T1u, TfF;
	       E TfG, T3t, T3A, Taf, Tai, T7B;
	       {
		    E T12, T13, T15, T16;
		    T12 = rio[WS(ios, 1)];
		    T13 = iio[-WS(ios, 33)];
		    T14 = T12 + T13;
		    T15 = rio[WS(ios, 17)];
		    T16 = iio[-WS(ios, 49)];
		    T17 = T15 + T16;
		    T18 = T14 + T17;
		    Tax = T15 - T16;
		    Tas = T12 - T13;
	       }
	       {
		    E T3I, T3J, T3L, T3M;
		    T3I = iio[-WS(ios, 1)];
		    T3J = rio[WS(ios, 33)];
		    T3K = T3I - T3J;
		    T3L = iio[-WS(ios, 17)];
		    T3M = rio[WS(ios, 49)];
		    T3N = T3L - T3M;
		    T7F = T3K + T3N;
		    Tay = T3I + T3J;
		    Tat = T3L + T3M;
	       }
	       {
		    E Tap, Tao, Tal, Tam;
		    {
			 E T19, T1a, T3c, T3d;
			 T19 = rio[WS(ios, 9)];
			 T1a = iio[-WS(ios, 41)];
			 T1b = T19 + T1a;
			 Tap = T19 - T1a;
			 T3c = iio[-WS(ios, 9)];
			 T3d = rio[WS(ios, 41)];
			 T3e = T3c - T3d;
			 Tao = T3c + T3d;
		    }
		    {
			 E T1c, T1d, T3f, T3g;
			 T1c = iio[-WS(ios, 57)];
			 T1d = rio[WS(ios, 25)];
			 T1e = T1c + T1d;
			 Tal = T1c - T1d;
			 T3f = iio[-WS(ios, 25)];
			 T3g = rio[WS(ios, 57)];
			 T3h = T3f - T3g;
			 Tam = T3g + T3f;
		    }
		    T1f = T1b + T1e;
		    T7G = T3e + T3h;
		    TaB = Tal + Tam;
		    TaA = Tap + Tao;
		    Taq = Tao - Tap;
		    Tan = Tal - Tam;
	       }
	       {
		    E T1j, Ta6, T3n, Taa, T1m, Ta9, T3q, Ta7;
		    {
			 E T1h, T1i, T3l, T3m;
			 T1h = rio[WS(ios, 5)];
			 T1i = iio[-WS(ios, 37)];
			 T1j = T1h + T1i;
			 Ta6 = T1h - T1i;
			 T3l = iio[-WS(ios, 5)];
			 T3m = rio[WS(ios, 37)];
			 T3n = T3l - T3m;
			 Taa = T3l + T3m;
		    }
		    {
			 E T1k, T1l, T3o, T3p;
			 T1k = rio[WS(ios, 21)];
			 T1l = iio[-WS(ios, 53)];
			 T1m = T1k + T1l;
			 Ta9 = T1k - T1l;
			 T3o = iio[-WS(ios, 21)];
			 T3p = rio[WS(ios, 53)];
			 T3q = T3o - T3p;
			 Ta7 = T3o + T3p;
		    }
		    T1n = T1j + T1m;
		    TfC = Taa - Ta9;
		    TfD = Ta6 + Ta7;
		    T3k = T1j - T1m;
		    T3r = T3n - T3q;
		    Ta8 = Ta6 - Ta7;
		    Tab = Ta9 + Taa;
		    T7A = T3n + T3q;
	       }
	       {
		    E T1q, Tad, T3w, Tah, T1t, Tag, T3z, Tae;
		    {
			 E T1o, T1p, T3u, T3v;
			 T1o = iio[-WS(ios, 61)];
			 T1p = rio[WS(ios, 29)];
			 T1q = T1o + T1p;
			 Tad = T1o - T1p;
			 T3u = iio[-WS(ios, 29)];
			 T3v = rio[WS(ios, 61)];
			 T3w = T3u - T3v;
			 Tah = T3v + T3u;
		    }
		    {
			 E T1r, T1s, T3x, T3y;
			 T1r = rio[WS(ios, 13)];
			 T1s = iio[-WS(ios, 45)];
			 T1t = T1r + T1s;
			 Tag = T1r - T1s;
			 T3x = iio[-WS(ios, 13)];
			 T3y = rio[WS(ios, 45)];
			 T3z = T3x - T3y;
			 Tae = T3x + T3y;
		    }
		    T1u = T1q + T1t;
		    TfF = Tad + Tae;
		    TfG = Tag + Tah;
		    T3t = T1q - T1t;
		    T3A = T3w - T3z;
		    Taf = Tad - Tae;
		    Tai = Tag - Tah;
		    T7B = T3z + T3w;
	       }
	       T1g = T18 + T1f;
	       T1v = T1n + T1u;
	       T97 = T1g - T1v;
	       T98 = T7G + T7F;
	       T99 = T7A + T7B;
	       T9a = T98 - T99;
	       {
		    E T3b, T3i, TfE, TfH;
		    T3b = T14 - T17;
		    T3i = T3e - T3h;
		    T3j = T3b + T3i;
		    T5Z = T3b - T3i;
		    TfE = FNMS(KP382683432, TfD, KP923879532 * TfC);
		    TfH = FNMS(KP923879532, TfG, KP382683432 * TfF);
		    TfI = TfE + TfH;
		    Thk = TfH - TfE;
	       }
	       {
		    E TfN, TfO, TfK, TfL;
		    TfN = KP707106781 * (TaA + TaB);
		    TfO = Tas + Tat;
		    TfP = TfN + TfO;
		    Thl = TfO - TfN;
		    TfK = FMA(KP382683432, TfC, KP923879532 * TfD);
		    TfL = FMA(KP923879532, TfF, KP382683432 * TfG);
		    TfM = TfK + TfL;
		    Tho = TfK - TfL;
	       }
	       {
		    E T3E, T3F, T3s, T3B;
		    T3E = T3r - T3k;
		    T3F = T3t + T3A;
		    T3G = KP707106781 * (T3E + T3F);
		    T60 = KP707106781 * (T3F - T3E);
		    T3s = T3k + T3r;
		    T3B = T3t - T3A;
		    T3C = KP707106781 * (T3s + T3B);
		    T62 = KP707106781 * (T3s - T3B);
	       }
	       {
		    E T3H, T3O, Tac, Taj;
		    T3H = T1e - T1b;
		    T3O = T3K - T3N;
		    T3P = T3H + T3O;
		    T63 = T3O - T3H;
		    Tac = FNMS(KP382683432, Tab, KP923879532 * Ta8);
		    Taj = FMA(KP923879532, Taf, KP382683432 * Tai);
		    Tak = Tac + Taj;
		    TdQ = Taj - Tac;
	       }
	       {
		    E Tar, Tau, Tfz, TfA;
		    Tar = KP707106781 * (Tan - Taq);
		    Tau = Tas - Tat;
		    Tav = Tar + Tau;
		    TdT = Tau - Tar;
		    Tfz = Tay - Tax;
		    TfA = KP707106781 * (Taq + Tan);
		    TfB = Tfz + TfA;
		    Thn = Tfz - TfA;
	       }
	       {
		    E T7z, T7C, T7E, T7H;
		    T7z = T18 - T1f;
		    T7C = T7A - T7B;
		    T7D = T7z + T7C;
		    T8v = T7z - T7C;
		    T7E = T1u - T1n;
		    T7H = T7F - T7G;
		    T7I = T7E + T7H;
		    T8w = T7H - T7E;
	       }
	       {
		    E Taz, TaC, TaE, TaF;
		    Taz = Tax + Tay;
		    TaC = KP707106781 * (TaA - TaB);
		    TaD = Taz + TaC;
		    TdP = Taz - TaC;
		    TaE = FMA(KP923879532, Tab, KP382683432 * Ta8);
		    TaF = FNMS(KP382683432, Taf, KP923879532 * Tai);
		    TaG = TaE + TaF;
		    TdS = TaE - TaF;
	       }
	  }
	  {
	       E T11, T9K, T9T, Ta2, T22, T9Q, T9N, Ta3;
	       {
		    E Tv, T10, T9R, T9S;
		    Tv = Tf + Tu;
		    T10 = TK + TZ;
		    T11 = Tv + T10;
		    T9K = Tv - T10;
		    T9R = T9p + T9o;
		    T9S = T93 + T94;
		    T9T = T9R - T9S;
		    Ta2 = T9S + T9R;
	       }
	       {
		    E T1w, T21, T9L, T9M;
		    T1w = T1g + T1v;
		    T21 = T1L + T20;
		    T22 = T1w + T21;
		    T9Q = T21 - T1w;
		    T9L = T99 + T98;
		    T9M = T9e + T9d;
		    T9N = T9L - T9M;
		    Ta3 = T9L + T9M;
	       }
	       rio[0] = T11 + T22;
	       iio[-WS(ios, 63)] = Ta3 + Ta2;
	       {
		    E T9O, T9U, T9J, T9P;
		    T9O = T9K + T9N;
		    T9U = T9Q + T9T;
		    T9J = W[94];
		    T9P = W[95];
		    rio[WS(ios, 48)] = FNMS(T9P, T9U, T9J * T9O);
		    iio[-WS(ios, 15)] = FMA(T9P, T9O, T9J * T9U);
	       }
	       {
		    E T9W, T9Y, T9V, T9X;
		    T9W = T9K - T9N;
		    T9Y = T9T - T9Q;
		    T9V = W[30];
		    T9X = W[31];
		    rio[WS(ios, 16)] = FNMS(T9X, T9Y, T9V * T9W);
		    iio[-WS(ios, 47)] = FMA(T9X, T9W, T9V * T9Y);
	       }
	       {
		    E Ta0, Ta4, T9Z, Ta1;
		    Ta0 = T11 - T22;
		    Ta4 = Ta2 - Ta3;
		    T9Z = W[62];
		    Ta1 = W[63];
		    rio[WS(ios, 32)] = FNMS(Ta1, Ta4, T9Z * Ta0);
		    iio[-WS(ios, 31)] = FMA(Ta1, Ta0, T9Z * Ta4);
	       }
	  }
	  {
	       E T96, T9y, T9r, T9D, T9h, T9C, T9m, T9z;
	       {
		    E T92, T95, T9n, T9q;
		    T92 = Tf - Tu;
		    T95 = T93 - T94;
		    T96 = T92 + T95;
		    T9y = T92 - T95;
		    T9n = TZ - TK;
		    T9q = T9o - T9p;
		    T9r = T9n + T9q;
		    T9D = T9q - T9n;
	       }
	       {
		    E T9b, T9g, T9k, T9l;
		    T9b = T97 + T9a;
		    T9g = T9c - T9f;
		    T9h = KP707106781 * (T9b + T9g);
		    T9C = KP707106781 * (T9b - T9g);
		    T9k = T9a - T97;
		    T9l = T9c + T9f;
		    T9m = KP707106781 * (T9k + T9l);
		    T9z = KP707106781 * (T9l - T9k);
	       }
	       {
		    E T9i, T9s, T91, T9j;
		    T9i = T96 + T9h;
		    T9s = T9m + T9r;
		    T91 = W[110];
		    T9j = W[111];
		    rio[WS(ios, 56)] = FNMS(T9j, T9s, T91 * T9i);
		    iio[-WS(ios, 7)] = FMA(T9j, T9i, T91 * T9s);
	       }
	       {
		    E T9G, T9I, T9F, T9H;
		    T9G = T9y - T9z;
		    T9I = T9D - T9C;