text: /*
   * Copyright (c) 2003 Matteo Frigo
   * Copyright (c) 2003 Massachusetts Institute of Technology
   *
   * This program is free software; you can redistribute it and/or modify
   * it under the terms of the GNU General Public License as published by
   * the Free Software Foundation; either version 2 of the License, or
   * (at your option) any later version.
   *
   * This program is distributed in the hope that it will be useful,
   * but WITHOUT ANY WARRANTY; without even the implied warranty of
   * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   * GNU General Public License for more details.
   *
   * You should have received a copy of the GNU General Public License
   * along with this program; if not, write to the Free Software
   * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
   *
   */
  
  /* This file was automatically generated --- DO NOT EDIT */
  /* Generated on Sat Jul  5 21:30:32 EDT 2003 */
  
  #include "codelet-dft.h"
  
  /* Generated by: /homee/stevenj/cvs/fftw3.0.1/genfft/gen_twiddle -compact -variables 4 -twiddle-log3 -n 64 -name t2_64 -include t.h */
  
  /*
   * This function contains 1154 FP additions, 660 FP multiplications,
   * (or, 880 additions, 386 multiplications, 274 fused multiply/add),
   * 382 stack variables, and 256 memory accesses
   */
  /*
   * Generator Id's : 
   * $Id: algsimp.ml,v 1.7 2003/03/15 20:29:42 stevenj Exp $
   * $Id: fft.ml,v 1.2 2003/03/15 20:29:42 stevenj Exp $
   * $Id: gen_twiddle.ml,v 1.16 2003/04/16 19:51:27 athena Exp $
   */
  
  #include "t.h"
  
  static const R *t2_64(R *ri, R *ii, const R *W, stride ios, int m, int dist)
  {
       DK(KP290284677, +0.290284677254462367636192375817395274691476278);
       DK(KP956940335, +0.956940335732208864935797886980269969482849206);
       DK(KP881921264, +0.881921264348355029712756863660388349508442621);
       DK(KP471396736, +0.471396736825997648556387625905254377657460319);
       DK(KP098017140, +0.098017140329560601994195563888641845861136673);
       DK(KP995184726, +0.995184726672196886244836953109479921575474869);
       DK(KP773010453, +0.773010453362736960810906609758469800971041293);
       DK(KP634393284, +0.634393284163645498215171613225493370675687095);
       DK(KP555570233, +0.555570233019602224742830813948532874374937191);
       DK(KP831469612, +0.831469612302545237078788377617905756738560812);
       DK(KP980785280, +0.980785280403230449126182236134239036973933731);
       DK(KP195090322, +0.195090322016128267848284868477022240927691618);
       DK(KP707106781, +0.707106781186547524400844362104849039284835938);
       DK(KP923879532, +0.923879532511286756128183189396788286822416626);
       DK(KP382683432, +0.382683432365089771728459984030398866761344562);
       int i;
       for (i = m; i > 0; i = i - 1, ri = ri + dist, ii = ii + dist, W = W + 10) {
  	  E T1, Ti1, Tp, Tt, TH, TL, T6a, T6c, T4J, T4H, T1g, T91, T7W, T7m, T2O;
  	  E T4j, T7P, T4P, T8y, T2w, T8t, T2Z, T8e, T48, T1z, T7s, T1I, T7t, T8p, Ten;
  	  E T1Y, T7D, T2t, T7O, T7L, Te6, T3N, T8E, T7A, Te0, T4C, TeA, T8S, T9v, T65;
  	  E Tfi, T9J, Taq, T6K, Tf6, Ta2, Ta5, T73, Tfc, Tad, Tag, T3z, T83, T3u, T82;
  	  E T81, T84, T15, T9K, T68, T7j, T43, T9w, T4F, T8G, T5l, TeL, T9k, T9n, T6o;
  	  E Tf2, T9Q, T9R, T6z, Tf3, T9T, T9W, To, Ts, T4o, T8u, T4U, T92, T5a, TeT;
  	  E T8V, T8Y, T5G, TeG, T97, T9e, T27, T7X, T2T, T7E, T7b, Tai, T6T, Ta3, Tf7;
  	  E Ta8, T7Q, T2H, T2c, T76, Tah, T7F, T4d, T8z, TG, TK, T69, T6b, T3b, T87;
  	  E T5u, T9l, TeM, T9q, T88, T89, T3o, T86, T5P, T9f, TeH, T9a, T34, T8f, T1r;
  	  E T7n, T3S, T8F, T4G, T4I;
  	  T1 = ri[0];
  	  Ti1 = ii[0];
  	  Tp = ri[WS(ios, 32)];
  	  Tt = ii[WS(ios, 32)];
  	  TH = ri[WS(ios, 16)];
  	  TL = ii[WS(ios, 16)];
  	  T6a = ri[WS(ios, 47)];
  	  T6c = ii[WS(ios, 47)];
  	  T4J = ii[WS(ios, 49)];
  	  T4H = ri[WS(ios, 49)];
  	  {
  	       E T12, T14, T1b, T1f, T1q, T1m, T1w, T1y, T1D, T1H, T1S, T1M, T1N, T1W, T2M;
  	       E T2g, T2b, T29, T26, T22, T2i, T2E, T2G, T2v, T2u, T2r, T2n, T2N, T3a, T38;
  	       E T3l, T3n, T3r, T3t, T33, T31, T2Y, T2W, T4g, T2S, T2Q, T3w, T3y, T3E, T3G;
  	       E T3P, T3J, T3L, T3R, T4a, T4c, T47, T46, T42, T40, T4i, T6P, T6R, T6M, T6L;
  	       E T6I, T6G, T6W, T6Y, T74, T75, T5I, T78, T7a, T6x, T6v, T6s, T6q, T6h, T6m;
  	       E T6k, T6g, T5N, T5L, T5Z, T63, T66, T67, T5H, T54, T4D, T4A, T4y, T4n, T4l;
  	       E T4E, T4X, T4Z, T4T, T4R, T4O, T4N, T58, T5s, T5q, T5x, T5z, T5C, T5E, T5n;
  	       E T5m, T5j, T5h, T5d, T5e, Ta, Ty, Tf, Tw, T2, Tj, T3, Tc, T1E, T1B;
  	       E T1F, T1A, T1R, T3x, T2m, T3K, T61, T1V, T60, T3I, T51, T52, T2V, T56, T5X;
  	       E T3v, T55, T2X, T2q, T5W, T4w, T6E, Ta0, T8Q, Tac, T72, Tb, Tg, Th, T3e;
  	       E T3f, T3h, T1a, T2x, T2B, TU, TV, TY, T1e, T2y, T2A, TC, TD, T1u, Tk;
  	       E Tl, Tm, T39, T3U, T3W, T37, T3T, T3X, TQ, TR, TZ, T3c, T3d, T3i, Tx;
  	       E Tz, T1t, TN, TX, T2f, T5V, Tao, T2h, T3D, T4f, T4h, T3F, T3q, T3s;
  	       T12 = ri[WS(ios, 48)];
  	       T14 = ii[WS(ios, 48)];
  	       T1b = ri[WS(ios, 8)];
  	       T1f = ii[WS(ios, 8)];
  	       T1q = ii[WS(ios, 40)];
  	       T1m = ri[WS(ios, 40)];
  	       T1w = ri[WS(ios, 56)];
  	       T1y = ii[WS(ios, 56)];
  	       T1D = ri[WS(ios, 24)];
  	       T1H = ii[WS(ios, 24)];
  	       T1S = ri[WS(ios, 36)];
  	       T1M = ri[WS(ios, 4)];
  	       T1N = ii[WS(ios, 4)];
  	       T1W = ii[WS(ios, 36)];
  	       T2M = ri[WS(ios, 2)];
  	       T2g = ri[WS(ios, 60)];
  	       T2b = ii[WS(ios, 52)];
  	       T29 = ri[WS(ios, 52)];
  	       T26 = ii[WS(ios, 20)];
  	       T22 = ri[WS(ios, 20)];
  	       T2i = ii[WS(ios, 60)];
  	       T2E = ri[WS(ios, 44)];
  	       T2G = ii[WS(ios, 44)];
  	       T2v = ii[WS(ios, 12)];
  	       T2u = ri[WS(ios, 12)];
  	       T2r = ii[WS(ios, 28)];
  	       T2n = ri[WS(ios, 28)];
  	       T2N = ii[WS(ios, 2)];
  	       T3a = ii[WS(ios, 10)];
  	       T38 = ri[WS(ios, 10)];
  	       T3l = ri[WS(ios, 42)];
  	       T3n = ii[WS(ios, 42)];
  	       T3r = ri[WS(ios, 58)];
  	       T3t = ii[WS(ios, 58)];
  	       T33 = ii[WS(ios, 50)];
  	       T31 = ri[WS(ios, 50)];
  	       T2Y = ii[WS(ios, 18)];
  	       T2W = ri[WS(ios, 18)];
  	       T4g = ri[WS(ios, 54)];
  	       T2S = ii[WS(ios, 34)];
  	       T2Q = ri[WS(ios, 34)];
  	       T3w = ri[WS(ios, 26)];
  	       T3y = ii[WS(ios, 26)];
  	       T3E = ri[WS(ios, 62)];
  	       T3G = ii[WS(ios, 62)];
  	       T3P = ri[WS(ios, 14)];
  	       T3J = ri[WS(ios, 30)];
  	       T3L = ii[WS(ios, 30)];
  	       T3R = ii[WS(ios, 14)];
  	       T4a = ri[WS(ios, 38)];
  	       T4c = ii[WS(ios, 38)];
  	       T47 = ii[WS(ios, 6)];
  	       T46 = ri[WS(ios, 6)];
  	       T42 = ii[WS(ios, 46)];
  	       T40 = ri[WS(ios, 46)];
  	       T4i = ii[WS(ios, 54)];
  	       T6P = ri[WS(ios, 51)];
  	       T6R = ii[WS(ios, 51)];
  	       T6M = ii[WS(ios, 19)];
  	       T6L = ri[WS(ios, 19)];
  	       T6I = ii[WS(ios, 35)];
  	       T6G = ri[WS(ios, 35)];
  	       T6W = ri[WS(ios, 59)];
  	       T6Y = ii[WS(ios, 59)];
  	       T74 = ri[WS(ios, 11)];
  	       T75 = ii[WS(ios, 11)];
  	       T5I = ii[WS(ios, 13)];
  	       T78 = ri[WS(ios, 43)];
  	       T7a = ii[WS(ios, 43)];
  	       T6x = ii[WS(ios, 23)];
  	       T6v = ri[WS(ios, 23)];
  	       T6s = ii[WS(ios, 55)];
  	       T6q = ri[WS(ios, 55)];
  	       T6h = ii[WS(ios, 7)];
  	       T6m = ii[WS(ios, 39)];
  	       T6k = ri[WS(ios, 39)];
  	       T6g = ri[WS(ios, 7)];
  	       T5N = ii[WS(ios, 45)];
  	       T5L = ri[WS(ios, 45)];
  	       T5Z = ri[WS(ios, 31)];
  	       T63 = ii[WS(ios, 31)];
  	       T66 = ri[WS(ios, 15)];
  	       T67 = ii[WS(ios, 15)];
  	       T5H = ri[WS(ios, 13)];
  	       T54 = ri[WS(ios, 25)];
  	       T4D = ri[WS(ios, 17)];
  	       T4A = ii[WS(ios, 33)];
  	       T4y = ri[WS(ios, 33)];
  	       T4n = ii[WS(ios, 22)];
  	       T4l = ri[WS(ios, 22)];
  	       T4E = ii[WS(ios, 17)];
  	       T4X = ri[WS(ios, 57)];
  	       T4Z = ii[WS(ios, 57)];
  	       T4T = ii[WS(ios, 41)];
  	       T4R = ri[WS(ios, 41)];
  	       T4O = ii[WS(ios, 9)];
  	       T4N = ri[WS(ios, 9)];
  	       T58 = ii[WS(ios, 25)];
  	       T5s = ii[WS(ios, 53)];
  	       T5q = ri[WS(ios, 53)];
  	       T5x = ri[WS(ios, 61)];
  	       T5z = ii[WS(ios, 61)];
  	       T5C = ri[WS(ios, 29)];
  	       T5E = ii[WS(ios, 29)];
  	       T5n = ii[WS(ios, 21)];
  	       T5m = ri[WS(ios, 21)];
  	       T5j = ii[WS(ios, 37)];
  	       T5h = ri[WS(ios, 37)];
  	       T5d = ri[WS(ios, 5)];
  	       T5e = ii[WS(ios, 5)];
  	       {
  		    E T4u, T4v, T5T, T5U, T6C, T6D, T70, T71, T4, T7, T5, T8, TO, TP, T1U;
  		    E T2p, T18, T2k, T2l, T2o, TT, TS, T19, T1c, T1T, T1P, T1Q, T1d;
  		    T4u = ri[WS(ios, 1)];
  		    T4v = ii[WS(ios, 1)];
  		    T5T = ri[WS(ios, 63)];
  		    T5U = ii[WS(ios, 63)];
  		    T6C = ri[WS(ios, 3)];
  		    T6D = ii[WS(ios, 3)];
  		    T70 = ri[WS(ios, 27)];
  		    T71 = ii[WS(ios, 27)];
  		    {
  			 E T6, Te, T9, Td;
  			 T4 = W[2];
  			 T7 = W[3];
  			 T5 = W[0];
  			 T8 = W[1];
  			 T6 = T4 * T5;
  			 Te = T7 * T5;
  			 T9 = T7 * T8;
  			 Td = T4 * T8;
  			 Ta = T6 - T9;
  			 Ty = Td - Te;
  			 Tf = Td + Te;
  			 Tw = T6 + T9;
  			 T2 = W[6];
  			 Tj = W[7];
  			 T3 = W[4];
  			 Tc = W[5];
  			 TO = T3 * T4;
  			 TP = Tc * T7;
  			 T1U = Tj * T3;
  			 T2p = Tj * T5;
  			 T18 = T3 * T5;
  			 T2k = T2 * T5;
  			 T2l = Tj * T8;
  			 T2o = T2 * T8;
  			 TT = Tc * T4;
  			 TS = T3 * T7;
  			 T19 = Tc * T8;
  			 T1c = T3 * T8;
  			 T1E = T2 * T7;
  			 T1T = T2 * Tc;
  			 T1B = Tj * T7;
  			 T1F = Tj * T4;
  			 T1P = T2 * T3;
  			 T1Q = Tj * Tc;
  			 T1A = T2 * T4;
  			 T1d = Tc * T5;
  		    }
  		    T1R = T1P - T1Q;
  		    T3x = T2o - T2p;
  		    T2m = T2k - T2l;
  		    T3K = T1E + T1F;
  		    T61 = Tj * Ta;
  		    T1V = T1T + T1U;
  		    T60 = T2 * Tf;
  		    T3I = T1A - T1B;
  		    T51 = T2 * Tw;
  		    T52 = Tj * Ty;
  		    T2V = T1P + T1Q;
  		    T56 = Tj * Tw;
  		    T5X = Tj * Tf;
  		    T3v = T2k + T2l;
  		    T55 = T2 * Ty;
  		    T2X = T1T - T1U;
  		    T2q = T2o + T2p;
  		    T5W = T2 * Ta;
  		    T4w = FMA(T5, T4u, T8 * T4v);
  		    T6E = FMA(T4, T6C, T7 * T6D);
  		    Ta0 = FNMS(T7, T6C, T4 * T6D);
  		    T8Q = FNMS(T8, T4u, T5 * T4v);
  		    Tac = FNMS(Tj, T70, T2 * T71);
  		    T72 = FMA(T2, T70, Tj * T71);
  		    Tb = T3 * Ta;
  		    Tg = Tc * Tf;
  		    Th = Tb + Tg;
  		    T3e = TS - TT;
  		    T3f = Tj * T3e;
  		    T3h = T2 * T3e;
  		    T1a = T18 + T19;
  		    T2x = T2 * T1a;
  		    T2B = Tj * T1a;
  		    TU = TS + TT;
  		    TV = Tj * TU;
  		    TY = T2 * TU;
  		    T1e = T1c - T1d;
  		    T2y = Tj * T1e;
  		    T2A = T2 * T1e;
  		    TC = T3 * Ty;
  		    TD = Tc * Tw;
  		    T1u = TC + TD;
  		    Tk = T3 * Tf;
  		    Tl = Tc * Ta;
  		    Tm = Tk - Tl;
  		    T39 = T1c + T1d;
  		    T3U = Tj * T39;
  		    T3W = T2 * T39;
  		    T37 = T18 - T19;
  		    T3T = T2 * T37;
  		    T3X = Tj * T37;
  		    TQ = TO - TP;
  		    TR = T2 * TQ;
  		    TZ = Tj * TQ;
  		    T3c = TO + TP;
  		    T3d = T2 * T3c;
  		    T3i = Tj * T3c;
  		    Tx = T3 * Tw;
  		    Tz = Tc * Ty;
  		    T1t = Tx - Tz;
  		    TN = W[8];
  		    TX = W[9];
  		    T2f = FMA(TN, T4, TX * T7);
  		    T5V = FMA(TN, T5T, TX * T5U);
  		    Tao = FNMS(TX, T5T, TN * T5U);
  		    T2h = FNMS(TX, T4, TN * T7);
  		    T3D = FMA(TN, T5, TX * T8);
  		    T4f = FMA(TN, T3, TX * Tc);
  		    T4h = FNMS(TX, T3, TN * Tc);
  		    T3F = FNMS(TX, T5, TN * T8);
  	       }
  	       T1g = FNMS(T1e, T1f, T1a * T1b);
  	       T91 = FNMS(Tc, T4N, T3 * T4O);
  	       T7W = FMA(Ty, T2M, Tw * T2N);
  	       T7m = FMA(T1e, T1b, T1a * T1f);
  	       T2O = FNMS(Ty, T2N, Tw * T2M);
  	       T4j = FNMS(T4h, T4i, T4f * T4g);
  	       T7P = FNMS(TU, T2u, TQ * T2v);
  	       T4P = FMA(T3, T4N, Tc * T4O);
  	       T8y = FMA(T3e, T46, T3c * T47);
  	       T2w = FMA(TQ, T2u, TU * T2v);
  	       {
  		    E T1v, T1x, T1O, T1X;
  		    T8t = FMA(T4h, T4g, T4f * T4i);
  		    T2Z = FNMS(T2X, T2Y, T2V * T2W);
  		    T8e = FMA(T2X, T2W, T2V * T2Y);
  		    T48 = FNMS(T3e, T47, T3c * T46);
  		    T1v = FMA(TN, T1t, TX * T1u);
  		    T1x = FNMS(TX, T1t, TN * T1u);
  		    T1z = FNMS(T1x, T1y, T1v * T1w);
  		    T7s = FMA(T1x, T1w, T1v * T1y);
  		    {
  			 E T1C, T1G, T8n, T8o;
  			 T1C = T1A + T1B;
  			 T1G = T1E - T1F;
  			 T1I = FNMS(T1G, T1H, T1C * T1D);
  			 T7t = FMA(T1G, T1D, T1C * T1H);
  			 T8n = FMA(T3F, T3E, T3D * T3G);
  			 T8o = FNMS(T3K, T3J, T3I * T3L);
  			 T8p = T8n - T8o;
  			 Ten = T8n + T8o;
  		    }
  		    T1O = FMA(Ta, T1M, Tf * T1N);
  		    T1X = FMA(T1R, T1S, T1V * T1W);
  		    T1Y = T1O + T1X;
  		    T7D = T1O - T1X;
  		    {
  			 E T2j, T2s, T7J, T7K;
  			 T2j = FNMS(T2h, T2i, T2f * T2g);
  			 T2s = FMA(T2m, T2n, T2q * T2r);
  			 T2t = T2j + T2s;
  			 T7O = T2j - T2s;
  			 T7J = FMA(T2h, T2g, T2f * T2i);
  			 T7K = FNMS(T2q, T2n, T2m * T2r);
  			 T7L = T7J - T7K;
  			 Te6 = T7J + T7K;
  		    }
  	       }
  	       {
  		    E T3H, T3M, T7y, T7z;
  		    T3H = FNMS(T3F, T3G, T3D * T3E);
  		    T3M = FMA(T3I, T3J, T3K * T3L);
  		    T3N = T3H + T3M;
  		    T8E = T3H - T3M;
  		    T7y = FNMS(Tf, T1M, Ta * T1N);
  		    T7z = FNMS(T1V, T1S, T1R * T1W);
  		    T7A = T7y - T7z;
  		    Te0 = T7y + T7z;
  	       }
  	       {
  		    E T4B, T8R, T4x, T4z;
  		    T4x = T3d + T3f;
  		    T4z = T3h - T3i;
  		    T4B = FNMS(T4z, T4A, T4x * T4y);
  		    T8R = FMA(T4z, T4y, T4x * T4A);
  		    T4C = T4w + T4B;
  		    TeA = T8Q + T8R;
  		    T8S = T8Q - T8R;
  		    T9v = T4w - T4B;
  	       }
  	       {
  		    E T64, Tap, T5Y, T62;
  		    T5Y = T5W - T5X;
  		    T62 = T60 + T61;
  		    T64 = FMA(T5Y, T5Z, T62 * T63);
  		    Tap = FNMS(T62, T5Z, T5Y * T63);
  		    T65 = T5V + T64;
  		    Tfi = Tao + Tap;
  		    T9J = T5V - T64;
  		    Taq = Tao - Tap;
  	       }
  	       {
  		    E T6J, Ta1, T6F, T6H;
  		    T6F = T2x + T2y;
  		    T6H = T2A - T2B;
  		    T6J = FNMS(T6H, T6I, T6F * T6G);
  		    Ta1 = FMA(T6H, T6G, T6F * T6I);
  		    T6K = T6E + T6J;
  		    Tf6 = Ta0 + Ta1;
  		    Ta2 = Ta0 - Ta1;
  		    Ta5 = T6E - T6J;
  	       }
  	       {
  		    E T6Z, Tab, T6V, T6X;
  		    T6V = FMA(TN, Ta, TX * Tf);
  		    T6X = FNMS(TX, Ta, TN * Tf);
  		    T6Z = FNMS(T6X, T6Y, T6V * T6W);
  		    Tab = FMA(T6X, T6W, T6V * T6Y);
  		    T73 = T6Z + T72;
  		    Tfc = Tab + Tac;
  		    Tad = Tab - Tac;
  		    Tag = T6Z - T72;
  	       }
  	       T3z = FNMS(T3x, T3y, T3v * T3w);
  	       T83 = FMA(T3x, T3w, T3v * T3y);
  	       T3q = FNMS(TX, Tm, TN * Th);
  	       T3s = FMA(TN, Tm, TX * Th);
  	       T3u = FMA(T3q, T3r, T3s * T3t);
  	       T82 = FNMS(T3s, T3r, T3q * T3t);
  	       T81 = T3u - T3z;
  	       T84 = T82 - T83;
  	       {
  		    E TW, T10, T11, T13;
  		    TW = TR + TV;
  		    T10 = TY - TZ;
  		    T11 = FNMS(TX, T10, TN * TW);
  		    T13 = FMA(TN, T10, TX * TW);
  		    T15 = FMA(T11, T12, T13 * T14);
  		    T9K = FMA(T10, T66, TW * T67);
  		    T68 = FNMS(T10, T67, TW * T66);
  		    T7j = FNMS(T13, T12, T11 * T14);
  	       }
  	       {
  		    E T3V, T3Y, T3Z, T41;
  		    T3V = T3T + T3U;
  		    T3Y = T3W - T3X;
  		    T3Z = FNMS(TX, T3Y, TN * T3V);
  		    T41 = FMA(TN, T3Y, TX * T3V);
  		    T43 = FMA(T3Z, T40, T41 * T42);
  		    T9w = FMA(T3Y, T4D, T3V * T4E);
  		    T4F = FNMS(T3Y, T4E, T3V * T4D);
  		    T8G = FNMS(T41, T40, T3Z * T42);
  	       }
  	       {
  		    E T5f, T9i, T5k, T9j, T5g, T5i;
  		    T5f = FNMS(Tm, T5e, Th * T5d);
  		    T9i = FMA(Tm, T5d, Th * T5e);
  		    T5g = T3T - T3U;
  		    T5i = T3W + T3X;
  		    T5k = FMA(T5g, T5h, T5i * T5j);
  		    T9j = FNMS(T5i, T5h, T5g * T5j);
  		    T5l = T5f + T5k;
  		    TeL = T9i + T9j;
  		    T9k = T9i - T9j;
  		    T9n = T5f - T5k;
  	       }
  	       {
  		    E T6i, T9O, T6n, T9P, T6j, T6l;
  		    T6i = FMA(T1t, T6g, T1u * T6h);
  		    T9O = FNMS(T1u, T6g, T1t * T6h);
  		    T6j = TR - TV;
  		    T6l = TY + TZ;
  		    T6n = FMA(T6j, T6k, T6l * T6m);
  		    T9P = FNMS(T6l, T6k, T6j * T6m);
  		    T6o = T6i + T6n;
  		    Tf2 = T9O + T9P;
  		    T9Q = T9O - T9P;
  		    T9R = T6i - T6n;
  	       }
  	       {
  		    E T6t, T9U, T6y, T9V;
  		    {
  			 E T6p, T6r, T6u, T6w;
  			 T6p = FNMS(TX, T1e, TN * T1a);
  			 T6r = FMA(TN, T1e, TX * T1a);
  			 T6t = FMA(T6p, T6q, T6r * T6s);
  			 T9U = FNMS(T6r, T6q, T6p * T6s);
  			 T6u = T5W + T5X;
  			 T6w = T60 - T61;
  			 T6y = FNMS(T6w, T6x, T6u * T6v);
  			 T9V = FMA(T6w, T6v, T6u * T6x);
  		    }
  		    T6z = T6t + T6y;
  		    Tf3 = T9U + T9V;
  		    T9T = T6t - T6y;
  		    T9W = T9U - T9V;
  	       }
  	       {
  		    E Ti, Tn, T4k, Tq, Tr, T4m, T4Q, T4S;
  		    Ti = T2 * Th;
  		    Tn = Tj * Tm;
  		    T4k = Ti - Tn;
  		    Tq = T2 * Tm;
  		    Tr = Tj * Th;
  		    T4m = Tq + Tr;
  		    To = Ti + Tn;
  		    Ts = Tq - Tr;
  		    T4o = FMA(T4k, T4l, T4m * T4n);
  		    T8u = FNMS(T4m, T4l, T4k * T4n);
  		    T4Q = FMA(TN, T4k, TX * T4m);
  		    T4S = FNMS(TX, T4k, TN * T4m);
  		    T4U = FNMS(T4S, T4T, T4Q * T4R);
  		    T92 = FMA(T4S, T4R, T4Q * T4T);
  	       }
  	       {
  		    E T50, T8W, T59, T8X;
  		    {
  			 E T4W, T4Y, T53, T57;
  			 T4W = FNMS(TX, T3e, TN * T3c);
  			 T4Y = FMA(TN, T3e, TX * T3c);
  			 T50 = FMA(T4W, T4X, T4Y * T4Z);
  			 T8W = FNMS(T4Y, T4X, T4W * T4Z);
  			 T53 = T51 - T52;
  			 T57 = T55 + T56;
  			 T59 = FMA(T53, T54, T57 * T58);
  			 T8X = FNMS(T57, T54, T53 * T58);
  		    }
  		    T5a = T50 + T59;
  		    TeT = T8W + T8X;
  		    T8V = T50 - T59;
  		    T8Y = T8W - T8X;
  	       }
  	       {
  		    E T5A, T9c, T5F, T9d;
  		    {
  			 E T5w, T5y, T5B, T5D;
  			 T5w = FNMS(TX, Ty, TN * Tw);
  			 T5y = FMA(TN, Ty, TX * Tw);
  			 T5A = FMA(T5w, T5x, T5y * T5z);
  			 T9c = FNMS(T5y, T5x, T5w * T5z);
  			 T5B = T51 + T52;
  			 T5D = T55 - T56;
  			 T5F = FNMS(T5D, T5E, T5B * T5C);
  			 T9d = FMA(T5D, T5C, T5B * T5E);
  		    }
  		    T5G = T5A + T5F;
  		    TeG = T9c + T9d;
  		    T97 = T5A - T5F;
  		    T9e = T9c - T9d;
  	       }
  	       {
  		    E T21, T2P, T25, T2R, T77, T79;
  		    {
  			 E T1Z, T20, T23, T24;
  			 T1Z = T2 * T1t;
  			 T20 = Tj * T1u;
  			 T21 = T1Z + T20;
  			 T2P = T1Z - T20;
  			 T23 = T2 * T1u;
  			 T24 = Tj * T1t;
  			 T25 = T23 - T24;
  			 T2R = T23 + T24;
  		    }
  		    T27 = FNMS(T25, T26, T21 * T22);
  		    T7X = FNMS(T2R, T2Q, T2P * T2S);
  		    T2T = FMA(T2P, T2Q, T2R * T2S);
  		    T7E = FMA(T25, T22, T21 * T26);
  		    T77 = FNMS(TX, T25, TN * T21);
  		    T79 = FMA(TN, T25, TX * T21);
  		    T7b = FMA(T77, T78, T79 * T7a);
  		    Tai = FNMS(T79, T78, T77 * T7a);
  	       }
  	       {
  		    E T6S, Ta7, T2D, Ta6, T2F, T6N;
  		    {
  			 E T6O, T6Q, T2z, T2C;
  			 T6O = FMA(TN, TQ, TX * TU);
  			 T6Q = FNMS(TX, TQ, TN * TU);
  			 T6S = FNMS(T6Q, T6R, T6O * T6P);
  			 Ta7 = FMA(T6Q, T6P, T6O * T6R);
  			 T2z = T2x - T2y;
  			 T2C = T2A + T2B;
  			 T2D = FMA(TN, T2z, TX * T2C);
  			 Ta6 = FNMS(T2C, T6L, T2z * T6M);
  			 T2F = FNMS(TX, T2z, TN * T2C);
  			 T6N = FMA(T2z, T6L, T2C * T6M);
  		    }
  		    T6T = T6N + T6S;
  		    Ta3 = T6N - T6S;
  		    Tf7 = Ta6 + Ta7;
  		    Ta8 = Ta6 - Ta7;
  		    T7Q = FMA(T2F, T2E, T2D * T2G);
  		    T2H = FNMS(T2F, T2G, T2D * T2E);
  	       }
  	       {
  		    E TA, TE, TB, TF, TJ, TI, T2a, T28, T49, T4b;
  		    TA = Tx + Tz;
  		    TE = TC - TD;
  		    TB = T2 * TA;
  		    TF = Tj * TE;
  		    TJ = Tj * TA;
  		    TI = T2 * TE;
  		    T2a = FMA(TN, TE, TX * TA);
  		    T28 = FNMS(TX, TE, TN * TA);
  		    T2c = FMA(T28, T29, T2a * T2b);
  		    T76 = FNMS(TE, T75, TA * T74);
  		    Tah = FMA(TE, T74, TA * T75);
  		    T7F = FNMS(T2a, T29, T28 * T2b);
  		    T49 = TB + TF;
  		    T4b = TI - TJ;
  		    T4d = FNMS(T4b, T4c, T49 * T4a);
  		    T8z = FMA(T4b, T4a, T49 * T4c);
  		    TG = TB - TF;
  		    TK = TI + TJ;
  		    T69 = FMA(TN, TG, TX * TK);
  		    T6b = FNMS(TX, TG, TN * TK);
  	       }
  	       {
  		    E T5t, T9p, T3k, T9o, T3m, T5o;
  		    T3b = FMA(T37, T38, T39 * T3a);
  		    T87 = FNMS(T39, T38, T37 * T3a);
  		    {
  			 E T5p, T5r, T3g, T3j;
  			 T5p = FMA(TN, T37, TX * T39);
  			 T5r = FNMS(TX, T37, TN * T39);
  			 T5t = FNMS(T5r, T5s, T5p * T5q);
  			 T9p = FMA(T5r, T5q, T5p * T5s);
  			 T3g = T3d - T3f;
  			 T3j = T3h + T3i;
  			 T3k = FMA(TN, T3g, TX * T3j);
  			 T9o = FNMS(T3j, T5m, T3g * T5n);
  			 T3m = FNMS(TX, T3g, TN * T3j);
  			 T5o = FMA(T3g, T5m, T3j * T5n);
  		    }
  		    T5u = T5o + T5t;
  		    T9l = T5o - T5t;
  		    TeM = T9o + T9p;
  		    T9q = T9o - T9p;
  		    T88 = FMA(T3m, T3l, T3k * T3n);
  		    T89 = T87 - T88;
  		    T3o = FNMS(T3m, T3n, T3k * T3l);
  		    T86 = T3b - T3o;
  	       }
  	       {
  		    E T5O, T99, T1i, T1n, T1o, T1k, T30, T5J, T98, T32;
  		    {
  			 E T5K, T5M, T1h, T1j;
  			 T5K = FNMS(TX, T2X, TN * T2V);
  			 T5M = FMA(TN, T2X, TX * T2V);
  			 T5O = FMA(T5K, T5L, T5M * T5N);
  			 T99 = FNMS(T5M, T5L, T5K * T5N);
  			 T1h = Tb - Tg;
  			 T1j = Tk + Tl;
  			 T1i = T2 * T1h;
  			 T1n = T2 * T1j;
  			 T1o = Tj * T1h;
  			 T1k = Tj * T1j;
  			 T30 = FMA(TN, T1h, TX * T1j);
  			 T5J = FMA(T1h, T5H, T1j * T5I);
  			 T98 = FNMS(T1j, T5H, T1h * T5I);
  			 T32 = FNMS(TX, T1h, TN * T1j);
  		    }
  		    T5P = T5J + T5O;
  		    T9f = T5J - T5O;
  		    TeH = T98 + T99;
  		    T9a = T98 - T99;
  		    T34 = FNMS(T32, T33, T30 * T31);
  		    T8f = FMA(T32, T31, T30 * T33);
  		    {
  			 E T1l, T1p, T3O, T3Q;
  			 T1l = T1i - T1k;
  			 T1p = T1n + T1o;
  			 T1r = FMA(T1l, T1m, T1p * T1q);
  			 T7n = FNMS(T1p, T1m, T1l * T1q);
  			 T3O = T1i + T1k;
  			 T3Q = T1n - T1o;
  			 T3S = FNMS(T3Q, T3R, T3O * T3P);
  			 T8F = FMA(T3Q, T3P, T3O * T3R);
  			 T4G = FNMS(TX, T3Q, TN * T3O);
  			 T4I = FMA(TN, T3Q, TX * T3O);
  		    }
  	       }
  	  }
  	  {
  	       E T5R, TgT, TgY, ThE, T9t, Tbe, T9G, Tbb, Tcl, Tdq, Tcs, Tdn, TeP, Tg4, TeY;
  	       E Tg1, T7e, Th4, ThJ, Th9, Tfp, Tg8, Tfg, Tgb, T2K, TgC, Tih, ThX, TfQ, TiL;
  	       E Tea, Tiv, Tam, Tbl, TcL, Tdu, Taz, Tbi, TcE, Tdx, T7U, Tjv, Tdc, Tjh, Tb0;
  	       E TjL, TbU, TiZ, T8D, Tb5, Tc8, Tdi, T8M, Tb6, Tc5, Tdh, T4r, Thz, Tex, Tfz;
  	       E TfX, Tgl, TgN, Thj, T8m, TaI, Tdg, TdG, Tb4, Tbu, Tc2, TcU, T3C, Thy, Tem;
  	       E Tfy, TfU, Tgk, TgI, Thi, T6B, Th1, Tfm, Tga, Th8, ThI, T9Z, Tbh, Taw, Tbk;
  	       E TcI, Tdw, Tf5, Tg7, Tcx, Tdt, T5c, TgV, TeV, Tg0, TgS, ThD, TeE, Tg3, T96;
  	       E Tbd, Tce, Tdp, Tcp, Tdm, T9D, Tba, T1L, Tgz, Ti4, Tii, Tiy, TiM, TdZ, TfN;
  	       E T7x, TaX, Tj4, Tji, Tjy, TjM, TbN, Td9;
  	       {
  		    E T5v, T5Q, TgW, TgX;
  		    T5v = T5l + T5u;
  		    T5Q = T5G + T5P;
  		    T5R = T5v + T5Q;
  		    TgT = T5Q - T5v;
  		    TgW = TeL + TeM;
  		    TgX = TeG + TeH;
  		    TgY = TgW - TgX;
  		    ThE = TgW + TgX;
  	       }
  	       {
  		    E T9h, T9F, T9s, T9E;
  		    {
  			 E T9b, T9g, T9m, T9r;
  			 T9b = T97 - T9a;
  			 T9g = T9e + T9f;
  			 T9h = FNMS(KP923879532, T9g, KP382683432 * T9b);
  			 T9F = FMA(KP382683432, T9g, KP923879532 * T9b);
  			 T9m = T9k + T9l;
  			 T9r = T9n - T9q;
  			 T9s = FMA(KP923879532, T9m, KP382683432 * T9r);
  			 T9E = FNMS(KP923879532, T9r, KP382683432 * T9m);
  		    }
  		    T9t = T9h - T9s;
  		    Tbe = T9E + T9F;
  		    T9G = T9E - T9F;
  		    Tbb = T9s + T9h;
  	       }
  	       {
  		    E Tch, Tcr, Tck, Tcq;
  		    {
  			 E Tcf, Tcg, Tci, Tcj;
  			 Tcf = T97 + T9a;
  			 Tcg = T9e - T9f;
  			 Tch = FNMS(KP382683432, Tcg, KP923879532 * Tcf);
  			 Tcr = FMA(KP923879532, Tcg, KP382683432 * Tcf);
  			 Tci = T9k - T9l;
  			 Tcj = T9n + T9q;
  			 Tck = FMA(KP382683432, Tci, KP923879532 * Tcj);
  			 Tcq = FNMS(KP382683432, Tcj, KP923879532 * Tci);
  		    }
  		    Tcl = Tch - Tck;
  		    Tdq = Tcq + Tcr;
  		    Tcs = Tcq - Tcr;
  		    Tdn = Tck + Tch;
  	       }
  	       {
  		    E TeJ, TeX, TeO, TeW;
  		    {
  			 E TeF, TeI, TeK, TeN;
  			 TeF = T5G - T5P;
  			 TeI = TeG - TeH;
  			 TeJ = TeF - TeI;
  			 TeX = TeF + TeI;
  			 TeK = T5l - T5u;
  			 TeN = TeL - TeM;
  			 TeO = TeK + TeN;
  			 TeW = TeN - TeK;
  		    }
  		    TeP = KP707106781 * (TeJ - TeO);
  		    Tg4 = KP707106781 * (TeW + TeX);
  		    TeY = KP707106781 * (TeW - TeX);
  		    Tg1 = KP707106781 * (TeO + TeJ);
  	       }
  	       {
  		    E T6U, Th2, T7d, Tfb, Tfe, Th3, Tfa, Tfo, Tfn, Tff;
  		    T6U = T6K + T6T;
  		    Th2 = Tf6 + Tf7;
  		    {
  			 E T7c, Tfd, Tf8, Tf9;
  			 T7c = T76 + T7b;
  			 T7d = T73 + T7c;
  			 Tfb = T73 - T7c;
  			 Tfd = Tah + Tai;
  			 Tfe = Tfc - Tfd;
  			 Th3 = Tfc + Tfd;
  			 Tf8 = Tf6 - Tf7;
  			 Tf9 = T6K - T6T;
  			 Tfa = Tf8 - Tf9;
  			 Tfo = Tf9 + Tf8;
  		    }
  		    T7e = T6U + T7d;
  		    Th4 = Th2 - Th3;
  		    ThJ = Th2 + Th3;
  		    Th9 = T7d - T6U;
  		    Tfn = Tfb - Tfe;
  		    Tfp = KP707106781 * (Tfn - Tfo);
  		    Tg8 = KP707106781 * (Tfo + Tfn);
  		    Tff = Tfb + Tfe;
  		    Tfg = KP707106781 * (Tfa - Tff);
  		    Tgb = KP707106781 * (Tfa + Tff);
  	       }
  	       {
  		    E T2e, Te3, Te8, TgB, T2J, Te5, Te2, TgA;
  		    {
  			 E T2d, Te7, T2I, Te1;
  			 T2d = T27 + T2c;
  			 T2e = T1Y + T2d;
  			 Te3 = T1Y - T2d;
  			 Te7 = T7P + T7Q;
  			 Te8 = Te6 - Te7;
  			 TgB = Te6 + Te7;
  			 T2I = T2w + T2H;
  			 T2J = T2t + T2I;
  			 Te5 = T2t - T2I;
  			 Te1 = T7E + T7F;
  			 Te2 = Te0 - Te1;
  			 TgA = Te0 + Te1;
  		    }
  		    T2K = T2e + T2J;
  		    TgC = TgA - TgB;
  		    Tih = T2J - T2e;
  		    ThX = TgA + TgB;
  		    {
  			 E TfO, TfP, Te4, Te9;
  			 TfO = Te3 + Te2;
  			 TfP = Te5 - Te8;
  			 TfQ = KP707106781 * (TfO + TfP);
  			 TiL = KP707106781 * (TfP - TfO);
  			 Te4 = Te2 - Te3;
  			 Te9 = Te5 + Te8;
  			 Tea = KP707106781 * (Te4 - Te9);
  			 Tiv = KP707106781 * (Te4 + Te9);
  		    }
  	       }
  	       {
  		    E Taf, TcB, Tak, TcC, Taa, Tay, TcA, TcK, Tae, Taj;
  		    Tae = T76 - T7b;
  		    Taf = Tad + Tae;
  		    TcB = Tad - Tae;
  		    Taj = Tah - Tai;
  		    Tak = Tag - Taj;
  		    TcC = Tag + Taj;
  		    {
  			 E Ta4, Ta9, Tcy, Tcz;
  			 Ta4 = Ta2 + Ta3;
  			 Ta9 = Ta5 - Ta8;
  			 Taa = FNMS(KP923879532, Ta9, KP382683432 * Ta4);
  			 Tay = FMA(KP923879532, Ta4, KP382683432 * Ta9);
  			 Tcy = Ta2 - Ta3;
  			 Tcz = Ta5 + Ta8;
  			 TcA = FNMS(KP382683432, Tcz, KP923879532 * Tcy);
  			 TcK = FMA(KP382683432, Tcy, KP923879532 * Tcz);
  		    }
  		    {
  			 E Tal, TcJ, Tax, TcD;
  			 Tal = FMA(KP382683432, Taf, KP923879532 * Tak);
  			 Tam = Taa - Tal;
  			 Tbl = Taa + Tal;
  			 TcJ = FNMS(KP382683432, TcB, KP923879532 * TcC);
  			 TcL = TcJ - TcK;
  			 Tdu = TcK + TcJ;
  			 Tax = FNMS(KP923879532, Taf, KP382683432 * Tak);
  			 Taz = Tax - Tay;
  			 Tbi = Tay + Tax;
  			 TcD = FMA(KP923879532, TcB, KP382683432 * TcC);
  			 TcE = TcA - TcD;
  			 Tdx = TcA + TcD;
  		    }
  	       }
  	       {
  		    E T7C, TbO, T7S, TbS, T7H, TbP, T7N, TbR;
  		    {
  			 E T7B, T7R, T7G, T7M;
  			 T7B = T27 - T2c;
  			 T7C = T7A + T7B;
  			 TbO = T7A - T7B;
  			 T7R = T7P - T7Q;
  			 T7S = T7O - T7R;
  			 TbS = T7O + T7R;
  			 T7G = T7E - T7F;
  			 T7H = T7D - T7G;
  			 TbP = T7D + T7G;
  			 T7M = T2w - T2H;
  			 T7N = T7L + T7M;
  			 TbR = T7L - T7M;
  		    }
  		    {
  			 E T7I, T7T, Tda, Tdb;
  			 T7I = FNMS(KP923879532, T7H, KP382683432 * T7C);
  			 T7T = FMA(KP382683432, T7N, KP923879532 * T7S);
  			 T7U = T7I - T7T;
  			 Tjv = T7I + T7T;
  			 Tda = FMA(KP382683432, TbO, KP923879532 * TbP);
  			 Tdb = FNMS(KP382683432, TbR, KP923879532 * TbS);
  			 Tdc = Tda + Tdb;
  			 Tjh = Tdb - Tda;
  		    }
  		    {
  			 E TaY, TaZ, TbQ, TbT;
  			 TaY = FMA(KP923879532, T7C, KP382683432 * T7H);
  			 TaZ = FNMS(KP923879532, T7N, KP382683432 * T7S);
  			 Tb0 = TaY + TaZ;
  			 TjL = TaZ - TaY;
  			 TbQ = FNMS(KP382683432, TbP, KP923879532 * TbO);
  			 TbT = FMA(KP923879532, TbR, KP382683432 * TbS);
  			 TbU = TbQ - TbT;
  			 TiZ = TbQ + TbT;
  		    }
  	       }
  	       {
  		    E T8r, Tc6, T8I, Tc3, T8w, T8K, T8B, T8J, T8q, T8H;
  		    T8q = T3S - T43;
  		    T8r = T8p + T8q;
  		    Tc6 = T8p - T8q;
  		    T8H = T8F - T8G;
  		    T8I = T8E - T8H;
  		    Tc3 = T8E + T8H;
  		    {
  			 E T8s, T8v, T8x, T8A;
  			 T8s = T4j - T4o;
  			 T8v = T8t - T8u;
  			 T8w = T8s - T8v;
  			 T8K = T8s + T8v;
  			 T8x = T48 - T4d;
  			 T8A = T8y - T8z;
  			 T8B = T8x + T8A;
  			 T8J = T8A - T8x;
  		    }
  		    {
  			 E T8C, Tc7, T8L, Tc4;
  			 T8C = KP707106781 * (T8w - T8B);
  			 T8D = T8r - T8C;
  			 Tb5 = T8r + T8C;
  			 Tc7 = KP707106781 * (T8J + T8K);
  			 Tc8 = Tc6 - Tc7;
  			 Tdi = Tc6 + Tc7;
  			 T8L = KP707106781 * (T8J - T8K);
  			 T8M = T8I - T8L;
  			 Tb6 = T8I + T8L;
  			 Tc4 = KP707106781 * (T8B + T8w);
  			 Tc5 = Tc3 - Tc4;
  			 Tdh = Tc3 + Tc4;
  		    }
  	       }
  	       {
  		    E T45, Tes, Tep, TgK, T4q, Teq, Tev, TgL, T44, Teo, Ter, Tew;
  		    T44 = T3S + T43;
  		    T45 = T3N + T44;
  		    Tes = T3N - T44;
  		    Teo = T8F + T8G;
  		    Tep = Ten - Teo;
  		    TgK = Ten + Teo;
  		    {
  			 E T4e, T4p, Tet, Teu;
  			 T4e = T48 + T4d;
  			 T4p = T4j + T4o;
  			 T4q = T4e + T4p;
  			 Teq = T4p - T4e;
  			 Tet = T8y + T8z;
  			 Teu = T8t + T8u;
  			 Tev = Tet - Teu;
  			 TgL = Tet + Teu;
  		    }
  		    T4r = T45 + T4q;
  		    Thz = TgK + TgL;
  		    Ter = Tep - Teq;
  		    Tew = Tes - Tev;
  		    Tex = FMA(KP382683432, Ter, KP923879532 * Tew);
  		    Tfz = FNMS(KP923879532, Ter, KP382683432 * Tew);
  		    {
  			 E TfV, TfW, TgJ, TgM;
  			 TfV = Tep + Teq;
  			 TfW = Tes + Tev;
  			 TfX = FMA(KP923879532, TfV, KP382683432 * TfW);
  			 Tgl = FNMS(KP382683432, TfV, KP923879532 * TfW);
  			 TgJ = T45 - T4q;
  			 TgM = TgK - TgL;
  			 TgN = TgJ + TgM;
  			 Thj = TgJ - TgM;
  		    }
  	       }
  	       {
  		    E T80, TbW, T8k, TbX, T8b, Tc0, T8h, TbZ;
  		    {
  			 E T7Y, T7Z, T8i, T8j;
  			 T7Y = T7W - T7X;
  			 T7Z = T2Z - T34;
  			 T80 = T7Y + T7Z;
  			 TbW = T7Y - T7Z;
  			 T8i = T89 - T86;
  			 T8j = T81 + T84;
  			 T8k = KP707106781 * (T8i - T8j);
  			 TbX = KP707106781 * (T8i + T8j);
  		    }
  		    {
  			 E T85, T8a, T8d, T8g;
  			 T85 = T81 - T84;
  			 T8a = T86 + T89;
  			 T8b = KP707106781 * (T85 - T8a);
  			 Tc0 = KP707106781 * (T8a + T85);
  			 T8d = T2O - T2T;
  			 T8g = T8e - T8f;
  			 T8h = T8d - T8g;
  			 TbZ = T8d + T8g;
  		    }
  		    {
  			 E T8c, T8l, Tde, Tdf;
  			 T8c = T80 - T8b;
  			 T8l = T8h - T8k;
  			 T8m = FNMS(KP980785280, T8l, KP195090322 * T8c);
  			 TaI = FMA(KP980785280, T8c, KP195090322 * T8l);
  			 Tde = TbW + TbX;
  			 Tdf = TbZ + Tc0;
  			 Tdg = FNMS(KP195090322, Tdf, KP980785280 * Tde);
  			 TdG = FMA(KP980785280, Tdf, KP195090322 * Tde);
  		    }
  		    {
  			 E Tb2, Tb3, TbY, Tc1;
  			 Tb2 = T80 + T8b;
  			 Tb3 = T8h + T8k;
  			 Tb4 = FNMS(KP555570233, Tb3, KP831469612 * Tb2);
  			 Tbu = FMA(KP555570233, Tb2, KP831469612 * Tb3);
  			 TbY = TbW - TbX;
  			 Tc1 = TbZ - Tc0;
  			 Tc2 = FNMS(KP831469612, Tc1, KP555570233 * TbY);
  			 TcU = FMA(KP555570233, Tc1, KP831469612 * TbY);
  		    }
  	       }
  	       {
  		    E T36, Teh, Tek, TgF, T3B, Tef, Tee, TgE, Teg, Tel;
  		    {
  			 E T2U, T35, Tei, Tej;
  			 T2U = T2O + T2T;
  			 T35 = T2Z + T34;
  			 T36 = T2U + T35;
  			 Teh = T2U - T35;
  			 Tei = T87 + T88;
  			 Tej = T82 + T83;
  			 Tek = Tei - Tej;
  			 TgF = Tei + Tej;
  		    }
  		    {
  			 E T3p, T3A, Tec, Ted;
  			 T3p = T3b + T3o;
  			 T3A = T3u + T3z;
  			 T3B = T3p + T3A;
  			 Tef = T3A - T3p;
  			 Tec = T7W + T7X;
  			 Ted = T8e + T8f;
  			 Tee = Tec - Ted;
  			 TgE = Tec + Ted;
  		    }
  		    T3C = T36 + T3B;
  		    Thy = TgE + TgF;
  		    Teg = Tee - Tef;
  		    Tel = Teh - Tek;
  		    Tem = FNMS(KP923879532, Tel, KP382683432 * Teg);
  		    Tfy = FMA(KP923879532, Teg, KP382683432 * Tel);
  		    {
  			 E TfS, TfT, TgG, TgH;
  			 TfS = Tee + Tef;
  			 TfT = Teh + Tek;
  			 TfU = FNMS(KP382683432, TfT, KP923879532 * TfS);
  			 Tgk = FMA(KP382683432, TfS, KP923879532 * TfT);
  			 TgG = TgE - TgF;
  			 TgH = T36 - T3B;
  			 TgI = TgG - TgH;
  			 Thi = TgH + TgG;
  		    }
  	       }
  	       {
  		    E T6A, Tfl, Th7, Tf4, T6e, Tar, T9Y, TcH, Tav, Tcw, T9M, Tfj;
  		    T6A = T6o + T6z;
  		    Tfl = T6z - T6o;
  		    Th7 = Tf2 + Tf3;
  		    Tf4 = Tf2 - Tf3;
  		    {
  			 E T6d, T9S, T9X, Tat, Tau, T9L;
  			 T6d = FNMS(T6b, T6c, T69 * T6a);
  			 T6e = T68 + T6d;
  			 Tar = T68 - T6d;
  			 T9S = T9Q - T9R;
  			 T9X = T9T + T9W;
  			 T9Y = KP707106781 * (T9S - T9X);
  			 TcH = KP707106781 * (T9S + T9X);
  			 Tat = T9T - T9W;
  			 Tau = T9R + T9Q;
  			 Tav = KP707106781 * (Tat - Tau);
  			 Tcw = KP707106781 * (Tau + Tat);
  			 T9L = FMA(T6b, T6a, T69 * T6c);
  			 T9M = T9K - T9L;
  			 Tfj = T9K + T9L;
  		    }
  		    {
  			 E T6f, Tfk, Th6, T9N;
  			 T6f = T65 + T6e;
  			 T6B = T6f + T6A;
  			 Th1 = T6f - T6A;
  			 Tfk = Tfi - Tfj;
  			 Tfm = Tfk - Tfl;
  			 Tga = Tfk + Tfl;
  			 Th6 = Tfi + Tfj;
  			 Th8 = Th6 - Th7;
  			 ThI = Th6 + Th7;
  			 T9N = T9J - T9M;
  			 T9Z = T9N - T9Y;
  			 Tbh = T9N + T9Y;
  		    }
  		    {
  			 E Tas, TcG, Tf1, Tcv;
  			 Tas = Taq + Tar;
  			 Taw = Tas - Tav;
  			 Tbk = Tas + Tav;
  			 TcG = Taq - Tar;
  			 TcI = TcG - TcH;
  			 Tdw = TcG + TcH;
  			 Tf1 = T65 - T6e;
  			 Tf5 = Tf1 - Tf4;
  			 Tg7 = Tf1 + Tf4;
  			 Tcv = T9J + T9M;
  			 Tcx = Tcv - Tcw;
  			 Tdt = Tcv + Tcw;
  		    }
  	       }
  	       {
  		    E T8Z, T9B, T5b, TeD, TeU, TgR, T94, T9A, T4L, T8T, T9y, TeB, T4V;
  		    T8Z = T8V - T8Y;
  		    T9B = T8V + T8Y;
  		    T4V = T4P + T4U;
  		    T5b = T4V + T5a;
  		    TeD = T5a - T4V;
  		    {
  			 E TeS, T90, T93, T4K, T9x;
  			 TeS = T91 + T92;
  			 TeU = TeS - TeT;
  			 TgR = TeS + TeT;
  			 T90 = T4P - T4U;
  			 T93 = T91 - T92;
  			 T94 = T90 + T93;
  			 T9A = T93 - T90;
  			 T4K = FMA(T4G, T4H, T4I * T4J);
  			 T4L = T4F + T4K;
  			 T8T = T4F - T4K;
  			 T9x = FNMS(T4I, T4H, T4G * T4J);
  			 T9y = T9w - T9x;
  			 TeB = T9w + T9x;
  		    }
  		    {
  			 E T4M, TeR, TgQ, TeC;
  			 T4M = T4C + T4L;
  			 T5c = T4M + T5b;
  			 TgV = T4M - T5b;
  			 TeR = T4C - T4L;
  			 TeV = TeR - TeU;
  			 Tg0 = TeR + TeU;
  			 TgQ = TeA + TeB;
  			 TgS = TgQ - TgR;
  			 ThD = TgQ + TgR;
  			 TeC = TeA - TeB;
  			 TeE = TeC - TeD;
  			 Tg3 = TeC + TeD;
  		    }
  		    {
  			 E T8U, T95, Tcc, Tcd;
  			 T8U = T8S + T8T;
  			 T95 = KP707106781 * (T8Z - T94);
  			 T96 = T8U - T95;
  			 Tbd = T8U + T95;
  			 Tcc = T8S - T8T;
  			 Tcd = KP707106781 * (T9A + T9B);
  			 Tce = Tcc - Tcd;
  			 Tdp = Tcc + Tcd;
  		    }
  		    {
  			 E Tcn, Tco, T9z, T9C;
  			 Tcn = T9v + T9y;
  			 Tco = KP707106781 * (T94 + T8Z);
  			 Tcp = Tcn - Tco;
  			 Tdm = Tcn + Tco;
  			 T9z = T9v - T9y;
  			 T9C = KP707106781 * (T9A - T9B);
  			 T9D = T9z - T9C;
  			 Tba = T9z + T9C;
  		    }
  	       }
  	       {
  		    E Tv, T7h, TdY, ThY, Ti2, Tj1, T16, Tj2, T1K, Tiw, T7q, TbK, T7v, TbL, T7k;
  		    E ThZ, T7r, T7u, T7i;
  		    {
  			 E Tu, TdW, TdX, Ti0, TM;
  			 Tu = FNMS(Ts, Tt, To * Tp);
  			 Tv = T1 + Tu;
  			 T7h = T1 - Tu;
  			 TdW = T7m + T7n;
  			 TdX = T7s + T7t;
  			 TdY = TdW - TdX;
  			 ThY = TdW + TdX;
  			 Ti0 = FMA(Ts, Tp, To * Tt);
  			 Ti2 = Ti0 + Ti1;
  			 Tj1 = Ti1 - Ti0;
  			 TM = FMA(TG, TH, TK * TL);
  			 T16 = TM + T15;
  			 Tj2 = TM - T15;
  		    }
  		    {
  			 E T1s, T1J, T7o, T7p;
  			 T1s = T1g + T1r;
  			 T1J = T1z + T1I;
  			 T1K = T1s + T1J;
  			 Tiw = T1J - T1s;
  			 T7o = T7m - T7n;
  			 T7p = T1g - T1r;
  			 T7q = T7o - T7p;
  			 TbK = T7p + T7o;
  		    }
  		    T7r = T1z - T1I;
  		    T7u = T7s - T7t;
  		    T7v = T7r + T7u;
  		    TbL = T7r - T7u;
  		    T7i = FNMS(TK, TH, TG * TL);
  		    T7k = T7i - T7j;
  		    ThZ = T7i + T7j;
  		    {
  			 E T17, Ti3, Tix, TdV;
  			 T17 = Tv + T16;
  			 T1L = T17 + T1K;
  			 Tgz = T17 - T1K;
  			 Ti3 = ThZ + Ti2;
  			 Ti4 = ThY + Ti3;
  			 Tii = Ti3 - ThY;
  			 Tix = Ti2 - ThZ;
  			 Tiy = Tiw + Tix;
  			 TiM = Tix - Tiw;
  			 TdV = Tv - T16;
  			 TdZ = TdV - TdY;
  			 TfN = TdV + TdY;
  		    }
  		    {
  			 E T7l, T7w, Tj0, Tj3;
  			 T7l = T7h - T7k;
  			 T7w = KP707106781 * (T7q - T7v);
  			 T7x = T7l - T7w;
  			 TaX = T7l + T7w;
  			 Tj0 = KP707106781 * (T7q + T7v);
  			 Tj3 = Tj1 - Tj2;
  			 Tj4 = Tj0 + Tj3;
  			 Tji = Tj3 - Tj0;
  		    }
  		    {
  			 E Tjw, Tjx, TbJ, TbM;
  			 Tjw = KP707106781 * (TbL - TbK);
  			 Tjx = Tj2 + Tj1;
  			 Tjy = Tjw + Tjx;
  			 TjM = Tjx - Tjw;
  			 TbJ = T7h + T7k;
  			 TbM = KP707106781 * (TbK + TbL);
  			 TbN = TbJ - TbM;
  			 Td9 = TbJ + TbM;
  		    }
  	       }
  	       {
  		    E T4t, ThR, Ti6, Ti8, T7g, Ti7, ThU, ThV;
  		    {
  			 E T2L, T4s, ThW, Ti5;
  			 T2L = T1L + T2K;
  			 T4s = T3C + T4r;
  			 T4t = T2L + T4s;
  			 ThR = T2L - T4s;
  			 ThW = Thy + Thz;
  			 Ti5 = ThX + Ti4;
  			 Ti6 = ThW + Ti5;
  			 Ti8 = Ti5 - ThW;
  		    }
  		    {
  			 E T5S, T7f, ThS, ThT;
  			 T5S = T5c + T5R;
  			 T7f = T6B + T7e;
  			 T7g = T5S + T7f;
  			 Ti7 = T7f - T5S;
  			 ThS = ThD + ThE;
  			 ThT = ThI + ThJ;
  			 ThU = ThS - ThT;
  			 ThV = ThS + ThT;
  		    }
  		    ri[WS(ios, 32)] = T4t - T7g;
  		    ii[WS(ios, 32)] = Ti6 - ThV;
  		    ri[0] = T4t + T7g;
  		    ii[0] = ThV + Ti6;
  		    ri[WS(ios, 48)] = ThR - ThU;
  		    ii[WS(ios, 48)] = Ti8 - Ti7;
  		    ri[WS(ios, 16)] = ThR + ThU;
  		    ii[WS(ios, 16)] = Ti7 + Ti8;
  	       }
  	       {
  		    E ThB, ThN, Tic, Tie, ThG, ThO, ThL, ThP;
  		    {
  			 E Thx, ThA, Tia, Tib;
  			 Thx = T1L - T2K;
  			 ThA = Thy - Thz;
  			 ThB = Thx + ThA;
  			 ThN = Thx - ThA;
  			 Tia = T4r - T3C;
  			 Tib = Ti4 - ThX;
  			 Tic = Tia + Tib;
  			 Tie = Tib - Tia;
  		    }
  		    {
  			 E ThC, ThF, ThH, ThK;
  			 ThC = T5c - T5R;
  			 ThF = ThD - ThE;
  			 ThG = ThC + ThF;
  			 ThO = ThF - ThC;
  			 ThH = T6B - T7e;
  			 ThK = ThI - ThJ;
  			 ThL = ThH - ThK;
  			 ThP = ThH + ThK;
  		    }
  		    {
  			 E ThM, Ti9, ThQ, Tid;
  			 ThM = KP707106781 * (ThG + ThL);
  			 ri[WS(ios, 40)] = ThB - ThM;
  			 ri[WS(ios, 8)] = ThB + ThM;
  			 Ti9 = KP707106781 * (ThO + ThP);
  			 ii[WS(ios, 8)] = Ti9 + Tic;
  			 ii[WS(ios, 40)] = Tic - Ti9;
  			 ThQ = KP707106781 * (ThO - ThP);
  			 ri[WS(ios, 56)] = ThN - ThQ;
  			 ri[WS(ios, 24)] = ThN + ThQ;
  			 Tid = KP707106781 * (ThL - ThG);
  			 ii[WS(ios, 24)] = Tid + Tie;
  			 ii[WS(ios, 56)] = Tie - Tid;
  		    }
  	       }
  	       {
  		    E TgP, Thd, Tiq, Tis, Th0, The, Thb, Thf;
  		    {
  			 E TgD, TgO, Tio, Tip;
  			 TgD = Tgz - TgC;
  			 TgO = KP707106781 * (TgI - TgN);
  			 TgP = TgD + TgO;
  			 Thd = TgD - TgO;
  			 Tio = KP707106781 * (Thj - Thi);
  			 Tip = Tii - Tih;
  			 Tiq = Tio + Tip;
  			 Tis = Tip - Tio;
  		    }
  		    {
  			 E TgU, TgZ, Th5, Tha;
  			 TgU = TgS - TgT;
  			 TgZ = TgV - TgY;
  			 Th0 = FMA(KP923879532, TgU, KP382683432 * TgZ);
  			 The = FNMS(KP923879532, TgZ, KP382683432 * TgU);
  			 Th5 = Th1 - Th4;
  			 Tha = Th8 - Th9;
  			 Thb = FNMS(KP923879532, Tha, KP382683432 * Th5);
  			 Thf = FMA(KP382683432, Tha, KP923879532 * Th5);
  		    }
  		    {
  			 E Thc, Tin, Thg, Tir;
  			 Thc = Th0 + Thb;
  			 ri[WS(ios, 44)] = TgP - Thc;
  			 ri[WS(ios, 12)] = TgP + Thc;
  			 Tin = The + Thf;
  			 ii[WS(ios, 12)] = Tin + Tiq;
  			 ii[WS(ios, 44)] = Tiq - Tin;
  			 Thg = The - Thf;
  			 ri[WS(ios, 60)] = Thd - Thg;
  			 ri[WS(ios, 28)] = Thd + Thg;
  			 Tir = Thb - Th0;
  			 ii[WS(ios, 28)] = Tir + Tis;
  			 ii[WS(ios, 60)] = Tis - Tir;
  		    }
  	       }
  	       {
  		    E TfB, TfJ, TiO, TiQ, TfE, TfK, TfH, TfL;
  		    {
  			 E Tfx, TfA, TiK, TiN;
  			 Tfx = TdZ + Tea;
  			 TfA = Tfy + Tfz;
  			 TfB = Tfx + TfA;
  			 TfJ = Tfx - TfA;
  			 TiK = Tem + Tex;
  			 TiN = TiL + TiM;
  			 TiO = TiK + TiN;
  			 TiQ = TiN - TiK;
  		    }
  		    {
  			 E TfC, TfD, TfF, TfG;
  			 TfC = TeE + TeP;
  			 TfD = TeV + TeY;
  			 TfE = FMA(KP555570233, TfC, KP831469612 * TfD);
  			 TfK = FNMS(KP555570233, TfD, KP831469612 * TfC);
  			 TfF = Tf5 + Tfg;
  			 TfG = Tfm + Tfp;
  			 TfH = FNMS(KP555570233, TfG, KP831469612 * TfF);
  			 TfL = FMA(KP831469612, TfG, KP555570233 * TfF);
  		    }
  		    {
  			 E TfI, TiJ, TfM, TiP;
  			 TfI = TfE + TfH;
  			 ri[WS(ios, 38)] = TfB - TfI;
  			 ri[WS(ios, 6)] = TfB + TfI;
  			 TiJ = TfK + TfL;
  			 ii[WS(ios, 6)] = TiJ + TiO;
  			 ii[WS(ios, 38)] = TiO - TiJ;
  			 TfM = TfK - TfL;
  			 ri[WS(ios, 54)] = TfJ - TfM;
  			 ri[WS(ios, 22)] = TfJ + TfM;
  			 TiP = TfH - TfE;
  			 ii[WS(ios, 22)] = TiP + TiQ;
  			 ii[WS(ios, 54)] = TiQ - TiP;
  		    }
  	       }
  	       {
  		    E Thl, Tht, Tik, Tim, Tho, Thu, Thr, Thv;
  		    {
  			 E Thh, Thk, Tig, Tij;
  			 Thh = Tgz + TgC;
  			 Thk = KP707106781 * (Thi + Thj);
  			 Thl = Thh + Thk;
  			 Tht = Thh - Thk;
  			 Tig = KP707106781 * (TgI + TgN);
  			 Tij = Tih + Tii;
  			 Tik = Tig + Tij;
  			 Tim = Tij - Tig;
  		    }
  		    {
  			 E Thm, Thn, Thp, Thq;
  			 Thm = TgS + TgT;
  			 Thn = TgV + TgY;
  			 Tho = FMA(KP382683432, Thm, KP923879532 * Thn);
  			 Thu = FNMS(KP382683432, Thn, KP923879532 * Thm);
  			 Thp = Th1 + Th4;
  			 Thq = Th8 + Th9;
  			 Thr = FNMS(KP382683432, Thq, KP923879532 * Thp);
  			 Thv = FMA(KP923879532, Thq, KP382683432 * Thp);
  		    }
  		    {
  			 E Ths, Tif, Thw, Til;
  			 Ths = Tho + Thr;
  			 ri[WS(ios, 36)] = Thl - Ths;
  			 ri[WS(ios, 4)] = Thl + Ths;
  			 Tif = Thu + Thv;
  			 ii[WS(ios, 4)] = Tif + Tik;
  			 ii[WS(ios, 36)] = Tik - Tif;
  			 Thw = Thu - Thv;
  			 ri[WS(ios, 52)] = Tht - Thw;
  			 ri[WS(ios, 20)] = Tht + Thw;
  			 Til = Thr - Tho;
  			 ii[WS(ios, 20)] = Til + Tim;
  			 ii[WS(ios, 52)] = Tim - Til;
  		    }
  	       }
  	       {
  		    E Tez, Tft, TiU, TiW, Tf0, Tfu, Tfr, Tfv;
  		    {
  			 E Teb, Tey, TiS, TiT;
  			 Teb = TdZ - Tea;
  			 Tey = Tem - Tex;
  			 Tez = Teb + Tey;
  			 Tft = Teb - Tey;
  			 TiS = Tfz - Tfy;
  			 TiT = TiM - TiL;
  			 TiU = TiS + TiT;
  			 TiW = TiT - TiS;
  		    }
  		    {
  			 E TeQ, TeZ, Tfh, Tfq;
  			 TeQ = TeE - TeP;
  			 TeZ = TeV - TeY;
  			 Tf0 = FMA(KP980785280, TeQ, KP195090322 * TeZ);
  			 Tfu = FNMS(KP980785280, TeZ, KP195090322 * TeQ);
  			 Tfh = Tf5 - Tfg;
  			 Tfq = Tfm - Tfp;
  			 Tfr = FNMS(KP980785280, Tfq, KP195090322 * Tfh);
  			 Tfv = FMA(KP195090322, Tfq, KP980785280 * Tfh);
  		    }
  		    {
  			 E Tfs, TiR, Tfw, TiV;
  			 Tfs = Tf0 + Tfr;
  			 ri[WS(ios, 46)] = Tez - Tfs;
  			 ri[WS(ios, 14)] = Tez + Tfs;
  			 TiR = Tfu + Tfv;
  			 ii[WS(ios, 14)] = TiR + TiU;
  			 ii[WS(ios, 46)] = TiU - TiR;
  			 Tfw = Tfu - Tfv;
  			 ri[WS(ios, 62)] = Tft - Tfw;
  			 ri[WS(ios, 30)] = Tft + Tfw;
  			 TiV = Tfr - Tf0;
  			 ii[WS(ios, 30)] = TiV + TiW;
  			 ii[WS(ios, 62)] = TiW - TiV;
  		    }
  	       }
  	       {
  		    E TfZ, Tgf, TiG, TiI, Tg6, Tgg, Tgd, Tgh;
  		    {
  			 E TfR, TfY, TiE, TiF;
  			 TfR = TfN - TfQ;
  			 TfY = TfU - TfX;
  			 TfZ = TfR + TfY;
  			 Tgf = TfR - TfY;
  			 TiE = Tgl - Tgk;
  			 TiF = Tiy - Tiv;
  			 TiG = TiE + TiF;
  			 TiI = TiF - TiE;
  		    }
  		    {
  			 E Tg2, Tg5, Tg9, Tgc;
  			 Tg2 = Tg0 - Tg1;
  			 Tg5 = Tg3 - Tg4;
  			 Tg6 = FMA(KP555570233, Tg2, KP831469612 * Tg5);
  			 Tgg = FNMS(KP831469612, Tg2, KP555570233 * Tg5);
  			 Tg9 = Tg7 - Tg8;
  			 Tgc = Tga - Tgb;
  			 Tgd = FNMS(KP831469612, Tgc, KP555570233 * Tg9);
  			 Tgh = FMA(KP831469612, Tg9, KP555570233 * Tgc);
  		    }
  		    {
  			 E Tge, TiD, Tgi, TiH;
  			 Tge = Tg6 + Tgd;
  			 ri[WS(ios, 42)] = TfZ - Tge;
  			 ri[WS(ios, 10)] = TfZ + Tge;
  			 TiD = Tgg + Tgh;
  			 ii[WS(ios, 10)] = TiD + TiG;
  			 ii[WS(ios, 42)] = TiG - TiD;
  			 Tgi = Tgg - Tgh;
  			 ri[WS(ios, 58)] = Tgf - Tgi;
  			 ri[WS(ios, 26)] = Tgf + Tgi;
  			 TiH = Tgd - Tg6;
  			 ii[WS(ios, 26)] = TiH + TiI;
  			 ii[WS(ios, 58)] = TiI - TiH;
  		    }
  	       }
  	       {
  		    E Tgn, Tgv, TiA, TiC, Tgq, Tgw, Tgt, Tgx;
  		    {
  			 E Tgj, Tgm, Tiu, Tiz;
  			 Tgj = TfN + TfQ;
  			 Tgm = Tgk + Tgl;
  			 Tgn = Tgj + Tgm;
  			 Tgv = Tgj - Tgm;
  			 Tiu = TfU + TfX;
  			 Tiz = Tiv + Tiy;
  			 TiA = Tiu + Tiz;
  			 TiC = Tiz - Tiu;
  		    }
  		    {
  			 E Tgo, Tgp, Tgr, Tgs;
  			 Tgo = Tg0 + Tg1;
  			 Tgp = Tg3 + Tg4;
  			 Tgq = FMA(KP980785280, Tgo, KP195090322 * Tgp);
  			 Tgw = FNMS(KP195090322, Tgo, KP980785280 * Tgp);
  			 Tgr = Tg7 + Tg8;
  			 Tgs = Tga + Tgb;
  			 Tgt = FNMS(KP195090322, Tgs, KP980785280 * Tgr);
  			 Tgx = FMA(KP195090322, Tgr, KP980785280 * Tgs);
  		    }
  		    {
  			 E Tgu, Tit, Tgy, TiB;
  			 Tgu = Tgq + Tgt;
  			 ri[WS(ios, 34)] = Tgn - Tgu;
  			 ri[WS(ios, 2)] = Tgn + Tgu;
  			 Tit = Tgw + Tgx;
  			 ii[WS(ios, 2)] = Tit + TiA;
  			 ii[WS(ios, 34)] = TiA - Tit;
  			 Tgy = Tgw - Tgx;
  			 ri[WS(ios, 50)] = Tgv - Tgy;
  			 ri[WS(ios, 18)] = Tgv + Tgy;
  			 TiB = Tgt - Tgq;
  			 ii[WS(ios, 18)] = TiB + TiC;
  			 ii[WS(ios, 50)] = TiC - TiB;
  		    }
  	       }
  	       {
  		    E T7V, TjN, TjT, TaH, T8O, TjK, TaK, TjS, TaO, TaU, T9I, TaE, TaR, TaV, TaB;
  		    E TaF, T8N;
  		    T7V = T7x - T7U;
  		    TjN = TjL + TjM;
  		    TjT = TjM - TjL;
  		    TaH = T7x + T7U;
  		    T8N = FMA(KP195090322, T8D, KP980785280 * T8M);
  		    T8O = T8m - T8N;
  		    TjK = T8m + T8N;
  		    {
  			 E TaJ, TaM, TaN, T9u, T9H;
  			 TaJ = FNMS(KP980785280, T8D, KP195090322 * T8M);
  			 TaK = TaI + TaJ;
  			 TjS = TaJ - TaI;
  			 TaM = T96 + T9t;
  			 TaN = T9D + T9G;
  			 TaO = FMA(KP634393284, TaM, KP773010453 * TaN);
  			 TaU = FNMS(KP634393284, TaN, KP773010453 * TaM);
  			 T9u = T96 - T9t;
  			 T9H = T9D - T9G;
  			 T9I = FMA(KP995184726, T9u, KP098017140 * T9H);
  			 TaE = FNMS(KP995184726, T9H, KP098017140 * T9u);
  			 {
  			      E TaP, TaQ, Tan, TaA;
  			      TaP = T9Z + Tam;
  			      TaQ = Taw + Taz;
  			      TaR = FNMS(KP634393284, TaQ, KP773010453 * TaP);
  			      TaV = FMA(KP773010453, TaQ, KP634393284 * TaP);
  			      Tan = T9Z - Tam;
  			      TaA = Taw - Taz;
  			      TaB = FNMS(KP995184726, TaA, KP098017140 * Tan);
  			      TaF = FMA(KP098017140, TaA, KP995184726 * Tan);
  			 }
  		    }
  		    {
  			 E T8P, TaC, TjR, TjU;
  			 T8P = T7V + T8O;
  			 TaC = T9I + TaB;
  			 ri[WS(ios, 47)] = T8P - TaC;
  			 ri[WS(ios, 15)] = T8P + TaC;
  			 TjR = TaE + TaF;
  			 TjU = TjS + TjT;
  			 ii[WS(ios, 15)] = TjR + TjU;
  			 ii[WS(ios, 47)] = TjU - TjR;
  		    }
  		    {
  			 E TaD, TaG, TjV, TjW;
  			 TaD = T7V - T8O;
  			 TaG = TaE - TaF;
  			 ri[WS(ios, 63)] = TaD - TaG;
  			 ri[WS(ios, 31)] = TaD + TaG;
  			 TjV = TaB - T9I;
  			 TjW = TjT - TjS;
  			 ii[WS(ios, 31)] = TjV + TjW;
  			 ii[WS(ios, 63)] = TjW - TjV;
  		    }
  		    {
  			 E TaL, TaS, TjJ, TjO;
  			 TaL = TaH + TaK;
  			 TaS = TaO + TaR;
  			 ri[WS(ios, 39)] = TaL - TaS;
  			 ri[WS(ios, 7)] = TaL + TaS;
  			 TjJ = TaU + TaV;
  			 TjO = TjK + TjN;
  			 ii[WS(ios, 7)] = TjJ + TjO;
  			 ii[WS(ios, 39)] = TjO - TjJ;
  		    }
  		    {
  			 E TaT, TaW, TjP, TjQ;
  			 TaT = TaH - TaK;
  			 TaW = TaU - TaV;
  			 ri[WS(ios, 55)] = TaT - TaW;
  			 ri[WS(ios, 23)] = TaT + TaW;
  			 TjP = TaR - TaO;
  			 TjQ = TjN - TjK;
  			 ii[WS(ios, 23)] = TjP + TjQ;
  			 ii[WS(ios, 55)] = TjQ - TjP;
  		    }
  	       }
  	       {
  		    E TbV, Tjj, Tjp, TcT, Tca, Tjg, TcW, Tjo, Td0, Td6, Tcu, TcQ, Td3, Td7, TcN;
  		    E TcR, Tc9;
  		    TbV = TbN - TbU;
  		    Tjj = Tjh + Tji;
  		    Tjp = Tji - Tjh;
  		    TcT = TbN + TbU;
  		    Tc9 = FMA(KP831469612, Tc5, KP555570233 * Tc8);
  		    Tca = Tc2 - Tc9;
  		    Tjg = Tc2 + Tc9;
  		    {
  			 E TcV, TcY, TcZ, Tcm, Tct;
  			 TcV = FNMS(KP831469612, Tc8, KP555570233 * Tc5);
  			 TcW = TcU + TcV;
  			 Tjo = TcV - TcU;
  			 TcY = Tce + Tcl;
  			 TcZ = Tcp + Tcs;
  			 Td0 = FMA(KP471396736, TcY, KP881921264 * TcZ);
  			 Td6 = FNMS(KP471396736, TcZ, KP881921264 * TcY);
  			 Tcm = Tce - Tcl;
  			 Tct = Tcp - Tcs;
  			 Tcu = FMA(KP956940335, Tcm, KP290284677 * Tct);
  			 TcQ = FNMS(KP956940335, Tct, KP290284677 * Tcm);
  			 {
  			      E Td1, Td2, TcF, TcM;
  			      Td1 = Tcx + TcE;
  			      Td2 = TcI + TcL;
  			      Td3 = FNMS(KP471396736, Td2, KP881921264 * Td1);
  			      Td7 = FMA(KP881921264, Td2, KP471396736 * Td1);
  			      TcF = Tcx - TcE;
  			      TcM = TcI - TcL;
  			      TcN = FNMS(KP956940335, TcM, KP290284677 * TcF);
  			      TcR = FMA(KP290284677, TcM, KP956940335 * TcF);
  			 }
  		    }
  		    {
  			 E Tcb, TcO, Tjn, Tjq;
  			 Tcb = TbV + Tca;
  			 TcO = Tcu + TcN;
  			 ri[WS(ios, 45)] = Tcb - TcO;
  			 ri[WS(ios, 13)] = Tcb + TcO;
  			 Tjn = TcQ + TcR;
  			 Tjq = Tjo + Tjp;
  			 ii[WS(ios, 13)] = Tjn + Tjq;
  			 ii[WS(ios, 45)] = Tjq - Tjn;
  		    }
  		    {
  			 E TcP, TcS, Tjr, Tjs;
  			 TcP = TbV - Tca;
  			 TcS = TcQ - TcR;
  			 ri[WS(ios, 61)] = TcP - TcS;
  			 ri[WS(ios, 29)] = TcP + TcS;
  			 Tjr = TcN - Tcu;
  			 Tjs = Tjp - Tjo;
  			 ii[WS(ios, 29)] = Tjr + Tjs;
  			 ii[WS(ios, 61)] = Tjs - Tjr;
  		    }
  		    {
  			 E TcX, Td4, Tjf, Tjk;
  			 TcX = TcT + TcW;
  			 Td4 = Td0 + Td3;
  			 ri[WS(ios, 37)] = TcX - Td4;
  			 ri[WS(ios, 5)] = TcX + Td4;
  			 Tjf = Td6 + Td7;
  			 Tjk = Tjg + Tjj;
  			 ii[WS(ios, 5)] = Tjf + Tjk;
  			 ii[WS(ios, 37)] = Tjk - Tjf;
  		    }
  		    {
  			 E Td5, Td8, Tjl, Tjm;
  			 Td5 = TcT - TcW;
  			 Td8 = Td6 - Td7;
  			 ri[WS(ios, 53)] = Td5 - Td8;
  			 ri[WS(ios, 21)] = Td5 + Td8;
  			 Tjl = Td3 - Td0;
  			 Tjm = Tjj - Tjg;
  			 ii[WS(ios, 21)] = Tjl + Tjm;
  			 ii[WS(ios, 53)] = Tjm - Tjl;
  		    }
  	       }
  	       {
  		    E Tb1, Tjz, TjF, Tbt, Tb8, Tju, Tbw, TjE, TbA, TbG, Tbg, Tbq, TbD, TbH, Tbn;
  		    E Tbr, Tb7;
  		    Tb1 = TaX - Tb0;
  		    Tjz = Tjv + Tjy;
  		    TjF = Tjy - Tjv;
  		    Tbt = TaX + Tb0;
  		    Tb7 = FMA(KP831469612, Tb5, KP555570233 * Tb6);
  		    Tb8 = Tb4 - Tb7;
  		    Tju = Tb4 + Tb7;
  		    {
  			 E Tbv, Tby, Tbz, Tbc, Tbf;
  			 Tbv = FNMS(KP555570233, Tb5, KP831469612 * Tb6);
  			 Tbw = Tbu + Tbv;
  			 TjE = Tbv - Tbu;
  			 Tby = Tba + Tbb;
  			 Tbz = Tbd + Tbe;
  			 TbA = FMA(KP956940335, Tby, KP290284677 * Tbz);
  			 TbG = FNMS(KP290284677, Tby, KP956940335 * Tbz);
  			 Tbc = Tba - Tbb;
  			 Tbf = Tbd - Tbe;
  			 Tbg = FMA(KP471396736, Tbc, KP881921264 * Tbf);
  			 Tbq = FNMS(KP881921264, Tbc, KP471396736 * Tbf);
  			 {
  			      E TbB, TbC, Tbj, Tbm;
  			      TbB = Tbh + Tbi;
  			      TbC = Tbk + Tbl;
  			      TbD = FNMS(KP290284677, TbC, KP956940335 * TbB);
  			      TbH = FMA(KP290284677, TbB, KP956940335 * TbC);
  			      Tbj = Tbh - Tbi;
  			      Tbm = Tbk - Tbl;
  			      Tbn = FNMS(KP881921264, Tbm, KP471396736 * Tbj);
  			      Tbr = FMA(KP881921264, Tbj, KP471396736 * Tbm);
  			 }
  		    }
  		    {
  			 E Tb9, Tbo, TjD, TjG;
  			 Tb9 = Tb1 + Tb8;
  			 Tbo = Tbg + Tbn;
  			 ri[WS(ios, 43)] = Tb9 - Tbo;
  			 ri[WS(ios, 11)] = Tb9 + Tbo;
  			 TjD = Tbq + Tbr;
  			 TjG = TjE + TjF;
  			 ii[WS(ios, 11)] = TjD + TjG;
  			 ii[WS(ios, 43)] = TjG - TjD;
  		    }
  		    {
  			 E Tbp, Tbs, TjH, TjI;
  			 Tbp = Tb1 - Tb8;
  			 Tbs = Tbq - Tbr;
  			 ri[WS(ios, 59)] = Tbp - Tbs;
  			 ri[WS(ios, 27)] = Tbp + Tbs;
  			 TjH = Tbn - Tbg;
  			 TjI = TjF - TjE;
  			 ii[WS(ios, 27)] = TjH + TjI;
  			 ii[WS(ios, 59)] = TjI - TjH;
  		    }
  		    {
  			 E Tbx, TbE, Tjt, TjA;
  			 Tbx = Tbt + Tbw;
  			 TbE = TbA + TbD;
  			 ri[WS(ios, 35)] = Tbx - TbE;
  			 ri[WS(ios, 3)] = Tbx + TbE;
  			 Tjt = TbG + TbH;
  			 TjA = Tju + Tjz;
  			 ii[WS(ios, 3)] = Tjt + TjA;
  			 ii[WS(ios, 35)] = TjA - Tjt;
  		    }
  		    {
  			 E TbF, TbI, TjB, TjC;
  			 TbF = Tbt - Tbw;
  			 TbI = TbG - TbH;
  			 ri[WS(ios, 51)] = TbF - TbI;
  			 ri[WS(ios, 19)] = TbF + TbI;
  			 TjB = TbD - TbA;
  			 TjC = Tjz - Tju;
  			 ii[WS(ios, 19)] = TjB + TjC;
  			 ii[WS(ios, 51)] = TjC - TjB;
  		    }
  	       }
  	       {
  		    E Tdd, Tj5, Tjb, TdF, Tdk, TiY, TdI, Tja, TdM, TdS, Tds, TdC, TdP, TdT, Tdz;
  		    E TdD, Tdj;
  		    Tdd = Td9 - Tdc;
  		    Tj5 = TiZ + Tj4;
  		    Tjb = Tj4 - TiZ;
  		    TdF = Td9 + Tdc;
  		    Tdj = FMA(KP195090322, Tdh, KP980785280 * Tdi);
  		    Tdk = Tdg - Tdj;
  		    TiY = Tdg + Tdj;
  		    {
  			 E TdH, TdK, TdL, Tdo, Tdr;
  			 TdH = FNMS(KP195090322, Tdi, KP980785280 * Tdh);
  			 TdI = TdG + TdH;
  			 Tja = TdH - TdG;
  			 TdK = Tdm + Tdn;
  			 TdL = Tdp + Tdq;
  			 TdM = FMA(KP995184726, TdK, KP098017140 * TdL);
  			 TdS = FNMS(KP098017140, TdK, KP995184726 * TdL);
  			 Tdo = Tdm - Tdn;
  			 Tdr = Tdp - Tdq;
  			 Tds = FMA(KP634393284, Tdo, KP773010453 * Tdr);
  			 TdC = FNMS(KP773010453, Tdo, KP634393284 * Tdr);
  			 {
  			      E TdN, TdO, Tdv, Tdy;
  			      TdN = Tdt + Tdu;
  			      TdO = Tdw + Tdx;
  			      TdP = FNMS(KP098017140, TdO, KP995184726 * TdN);
  			      TdT = FMA(KP098017140, TdN, KP995184726 * TdO);
  			      Tdv = Tdt - Tdu;
  			      Tdy = Tdw - Tdx;
  			      Tdz = FNMS(KP773010453, Tdy, KP634393284 * Tdv);
  			      TdD = FMA(KP773010453, Tdv, KP634393284 * Tdy);
  			 }
  		    }
  		    {
  			 E Tdl, TdA, Tj9, Tjc;
  			 Tdl = Tdd + Tdk;
  			 TdA = Tds + Tdz;
  			 ri[WS(ios, 41)] = Tdl - TdA;
  			 ri[WS(ios, 9)] = Tdl + TdA;
  			 Tj9 = TdC + TdD;
  			 Tjc = Tja + Tjb;
  			 ii[WS(ios, 9)] = Tj9 + Tjc;
  			 ii[WS(ios, 41)] = Tjc - Tj9;
  		    }
  		    {
  			 E TdB, TdE, Tjd, Tje;
  			 TdB = Tdd - Tdk;
  			 TdE = TdC - TdD;
  			 ri[WS(ios, 57)] = TdB - TdE;
  			 ri[WS(ios, 25)] = TdB + TdE;
  			 Tjd = Tdz - Tds;
  			 Tje = Tjb - Tja;
  			 ii[WS(ios, 25)] = Tjd + Tje;
  			 ii[WS(ios, 57)] = Tje - Tjd;
  		    }
  		    {
  			 E TdJ, TdQ, TiX, Tj6;
  			 TdJ = TdF + TdI;
  			 TdQ = TdM + TdP;
  			 ri[WS(ios, 33)] = TdJ - TdQ;
  			 ri[WS(ios, 1)] = TdJ + TdQ;
  			 TiX = TdS + TdT;
  			 Tj6 = TiY + Tj5;
  			 ii[WS(ios, 1)] = TiX + Tj6;
  			 ii[WS(ios, 33)] = Tj6 - TiX;
  		    }
  		    {
  			 E TdR, TdU, Tj7, Tj8;
  			 TdR = TdF - TdI;
  			 TdU = TdS - TdT;
  			 ri[WS(ios, 49)] = TdR - TdU;
  			 ri[WS(ios, 17)] = TdR + TdU;
  			 Tj7 = TdP - TdM;
  			 Tj8 = Tj5 - TiY;
  			 ii[WS(ios, 17)] = Tj7 + Tj8;
  			 ii[WS(ios, 49)] = Tj8 - Tj7;
  		    }
  	       }
  	  }
       }
       return W;
  }
  
  static const tw_instr twinstr[] = {
       {TW_COS, 0, 1},
       {TW_SIN, 0, 1},
       {TW_COS, 0, 3},
       {TW_SIN, 0, 3},
       {TW_COS, 0, 9},
       {TW_SIN, 0, 9},
       {TW_COS, 0, 27},
       {TW_SIN, 0, 27},
       {TW_COS, 0, 63},
       {TW_SIN, 0, 63},
       {TW_NEXT, 1, 0}
  };
  
  static const ct_desc desc = { 64, "t2_64", twinstr, {880, 386, 274, 0}, &GENUS, 0, 0, 0 };
  
  void X(codelet_t2_64) (planner *p) {
       X(kdft_dit_register) (p, t2_64, &desc);
  }
  