r2cfii_64.c

快速fft变换

					T1E = R1[WS(rs, 26)];
				   }
			      }
			 }
		    }
	       }
	       {
		    E T6A, T4r, T4y, T3h, T3k, T36, T6k, T40, T5X, T6c, T6b, T60;
		    {
			 E T5w, T5t, T2Z, T6U, T6T, T32;
			 {
			      E Tx, T2N, T2v, T6V, T6P, T6Q, T1e, T2G, T31, T2X, T2L, T1Y, T6W, T2Q, T30;
			      E T2U;
			      {
				   E T1W, T1L, T2O, T2P, T2V, T2W, T6O, TS, T1d;
				   {
					E T4q, T4w, T1V, T1B, T1J, T4m, T4l, T1G, Tw, T1A, T4p;
					T6A = Tv + Tm;
					Tw = Tm - Tv;
					T1A = FMA(KP707106781, T1z, T1y);
					T4p = FMS(KP707106781, T1z, T1y);
					{
					     E T4o, T1x, T1I, T1F;
					     T4o = FMA(KP707106781, T1w, T1t);
					     T1x = FNMS(KP707106781, T1w, T1t);
					     T1I = T1D - T1E;
					     T1F = T1D + T1E;
					     T4q = FNMS(KP198912367, T4p, T4o);
					     T4w = FMA(KP198912367, T4o, T4p);
					     T1V = FMA(KP668178637, T1x, T1A);
					     T1B = FNMS(KP668178637, T1A, T1x);
					     T1J = FNMS(KP707106781, T1I, T1H);
					     T4m = FMA(KP707106781, T1I, T1H);
					     T4l = FMA(KP707106781, T1F, T1C);
					     T1G = FNMS(KP707106781, T1F, T1C);
					     Tx = FNMS(KP831469612, Tw, Td);
					     T2N = FMA(KP831469612, Tw, Td);
					}
					{
					     E T4n, T4x, T1U, T1K;
					     T4n = FMA(KP198912367, T4m, T4l);
					     T4x = FNMS(KP198912367, T4l, T4m);
					     T1U = FMA(KP668178637, T1G, T1J);
					     T1K = FNMS(KP668178637, T1J, T1G);
					     T5w = T4n + T4q;
					     T4r = T4n - T4q;
					     T5t = T4x + T4w;
					     T4y = T4w - T4x;
					     T3h = T1U - T1V;
					     T1W = T1U + T1V;
					     T3k = T1K + T1B;
					     T1L = T1B - T1K;
					     T6O = T34 + T35;
					     T36 = T34 - T35;
					}
				   }
				   T2O = FNMS(KP534511135, TK, TR);
				   TS = FMA(KP534511135, TR, TK);
				   T1d = FMA(KP534511135, T1c, T15);
				   T2P = FNMS(KP534511135, T15, T1c);
				   T2v = FMA(KP831469612, T2u, T2b);
				   T2V = FNMS(KP831469612, T2u, T2b);
				   T6V = FNMS(KP831469612, T6O, T6N);
				   T6P = FMA(KP831469612, T6O, T6N);
				   T6Q = TS + T1d;
				   T1e = TS - T1d;
				   T2W = FMA(KP831469612, T2F, T2C);
				   T2G = FNMS(KP831469612, T2F, T2C);
				   {
					E T2S, T2T, T1M, T1X;
					T2S = FMA(KP831469612, T1L, T1s);
					T1M = FNMS(KP831469612, T1L, T1s);
					T1X = FNMS(KP831469612, T1W, T1T);
					T2T = FMA(KP831469612, T1W, T1T);
					T31 = FMA(KP250486960, T2V, T2W);
					T2X = FNMS(KP250486960, T2W, T2V);
					T2L = FNMS(KP599376933, T1M, T1X);
					T1Y = FMA(KP599376933, T1X, T1M);
					T6W = T2O + T2P;
					T2Q = T2O - T2P;
					T30 = FMA(KP250486960, T2S, T2T);
					T2U = FNMS(KP250486960, T2T, T2S);
				   }
			      }
			      {
				   E T2J, T1f, T6X, T6Z, T2K, T2H;
				   T2J = FNMS(KP881921264, T1e, Tx);
				   T1f = FMA(KP881921264, T1e, Tx);
				   T6X = FNMS(KP881921264, T6W, T6V);
				   T6Z = FMA(KP881921264, T6W, T6V);
				   T2K = FNMS(KP599376933, T2v, T2G);
				   T2H = FMA(KP599376933, T2G, T2v);
				   {
					E T2R, T2Y, T6R, T6S;
					T2Z = FNMS(KP881921264, T2Q, T2N);
					T2R = FMA(KP881921264, T2Q, T2N);
					{
					     E T2M, T6Y, T70, T2I;
					     T2M = T2K - T2L;
					     T6Y = T2L + T2K;
					     T70 = T1Y + T2H;
					     T2I = T1Y - T2H;
					     Cr[WS(csr, 10)] = FMA(KP857728610, T2M, T2J);
					     Cr[WS(csr, 21)] = FNMS(KP857728610, T2M, T2J);
					     Ci[WS(csi, 5)] = FMA(KP857728610, T6Y, T6X);
					     Ci[WS(csi, 26)] = FMS(KP857728610, T6Y, T6X);
					     Ci[WS(csi, 21)] = FNMS(KP857728610, T70, T6Z);
					     Ci[WS(csi, 10)] = -(FMA(KP857728610, T70, T6Z));
					     Cr[WS(csr, 5)] = FMA(KP857728610, T2I, T1f);
					     Cr[WS(csr, 26)] = FNMS(KP857728610, T2I, T1f);
					     T2Y = T2U - T2X;
					     T6U = T2U + T2X;
					}
					T6T = FNMS(KP881921264, T6Q, T6P);
					T6R = FMA(KP881921264, T6Q, T6P);
					T6S = T30 + T31;
					T32 = T30 - T31;
					Cr[WS(csr, 2)] = FMA(KP970031253, T2Y, T2R);
					Cr[WS(csr, 29)] = FNMS(KP970031253, T2Y, T2R);
					Ci[WS(csi, 29)] = FNMS(KP970031253, T6S, T6R);
					Ci[WS(csi, 2)] = -(FMA(KP970031253, T6S, T6R));
				   }
			      }
			 }
			 {
			      E T5j, T5L, T5B, T6d, T67, T68, T5q, T5E, T5Z, T5V, T5J, T5y, T6e, T5O, T5Y;
			      E T5S;
			      {
				   E T5M, T5N, T5T, T5U;
				   {
					E T66, T5i, T5m, T5p;
					T6k = T5g + T5h;
					T5i = T5g - T5h;
					Cr[WS(csr, 13)] = FMA(KP970031253, T32, T2Z);
					Cr[WS(csr, 18)] = FNMS(KP970031253, T32, T2Z);
					Ci[WS(csi, 13)] = FNMS(KP970031253, T6U, T6T);
					Ci[WS(csi, 18)] = -(FMA(KP970031253, T6U, T6T));
					T5j = FNMS(KP980785280, T5i, T5f);
					T5L = FMA(KP980785280, T5i, T5f);
					T66 = T3W + T3Z;
					T40 = T3W - T3Z;
					T5M = FNMS(KP098491403, T5k, T5l);
					T5m = FMA(KP098491403, T5l, T5k);
					T5p = FMA(KP098491403, T5o, T5n);
					T5N = FNMS(KP098491403, T5n, T5o);
					T5B = FNMS(KP980785280, T5A, T5z);
					T5T = FMA(KP980785280, T5A, T5z);
					T6d = FNMS(KP980785280, T66, T65);
					T67 = FMA(KP980785280, T66, T65);
					T68 = T5m + T5p;
					T5q = T5m - T5p;
					T5U = FMA(KP980785280, T5D, T5C);
					T5E = FNMS(KP980785280, T5D, T5C);
				   }
				   {
					E T5Q, T5R, T5u, T5x;
					T5Q = FMA(KP980785280, T5t, T5s);
					T5u = FNMS(KP980785280, T5t, T5s);
					T5x = FNMS(KP980785280, T5w, T5v);
					T5R = FMA(KP980785280, T5w, T5v);
					T5Z = FNMS(KP049126849, T5T, T5U);
					T5V = FMA(KP049126849, T5U, T5T);
					T5J = FNMS(KP906347169, T5u, T5x);
					T5y = FMA(KP906347169, T5x, T5u);
					T6e = T5M + T5N;
					T5O = T5M - T5N;
					T5Y = FMA(KP049126849, T5Q, T5R);
					T5S = FNMS(KP049126849, T5R, T5Q);
				   }
			      }
			      {
				   E T5H, T5r, T6f, T6h, T5I, T5F;
				   T5H = FNMS(KP995184726, T5q, T5j);
				   T5r = FMA(KP995184726, T5q, T5j);
				   T6f = FNMS(KP995184726, T6e, T6d);
				   T6h = FMA(KP995184726, T6e, T6d);
				   T5I = FMA(KP906347169, T5B, T5E);
				   T5F = FNMS(KP906347169, T5E, T5B);
				   {
					E T5P, T5W, T69, T6a;
					T5X = FNMS(KP995184726, T5O, T5L);
					T5P = FMA(KP995184726, T5O, T5L);
					{
					     E T5K, T6g, T6i, T5G;
					     T5K = T5I - T5J;
					     T6g = T5J + T5I;
					     T6i = T5F - T5y;
					     T5G = T5y + T5F;
					     Cr[WS(csr, 8)] = FMA(KP740951125, T5K, T5H);
					     Cr[WS(csr, 23)] = FNMS(KP740951125, T5K, T5H);
					     Ci[WS(csi, 7)] = FMA(KP740951125, T6g, T6f);
					     Ci[WS(csi, 24)] = FMS(KP740951125, T6g, T6f);
					     Ci[WS(csi, 23)] = FMA(KP740951125, T6i, T6h);
					     Ci[WS(csi, 8)] = FMS(KP740951125, T6i, T6h);
					     Cr[WS(csr, 7)] = FMA(KP740951125, T5G, T5r);
					     Cr[WS(csr, 24)] = FNMS(KP740951125, T5G, T5r);
					     T5W = T5S + T5V;
					     T6c = T5V - T5S;
					}
					T6b = FNMS(KP995184726, T68, T67);
					T69 = FMA(KP995184726, T68, T67);
					T6a = T5Y + T5Z;
					T60 = T5Y - T5Z;
					Cr[0] = FMA(KP998795456, T5W, T5P);
					Cr[WS(csr, 31)] = FNMS(KP998795456, T5W, T5P);
					Ci[WS(csi, 31)] = FNMS(KP998795456, T6a, T69);
					Ci[0] = -(FMA(KP998795456, T6a, T69));
				   }
			      }
			 }
		    }
		    {
			 E T3L, T6G, T6F, T3O;
			 {
			      E T37, T3z, T3p, T6H, T6B, T6C, T3e, T3s, T3M, T3J, T3w, T3m, T6I, T3C, T3N;
			      E T3G;
			      {
				   E T3B, T3A, T3H, T3I, T3a, T3d;
				   Cr[WS(csr, 15)] = FMA(KP998795456, T60, T5X);
				   Cr[WS(csr, 16)] = FNMS(KP998795456, T60, T5X);
				   Ci[WS(csi, 15)] = FMA(KP998795456, T6c, T6b);
				   Ci[WS(csi, 16)] = FMS(KP998795456, T6c, T6b);
				   T37 = FNMS(KP831469612, T36, T33);
				   T3z = FMA(KP831469612, T36, T33);
				   T3B = FMA(KP303346683, T38, T39);
				   T3a = FNMS(KP303346683, T39, T38);
				   T3d = FNMS(KP303346683, T3c, T3b);
				   T3A = FMA(KP303346683, T3b, T3c);
				   T3p = FMA(KP831469612, T3o, T3n);
				   T3H = FNMS(KP831469612, T3o, T3n);
				   T6H = FNMS(KP831469612, T6A, T6z);
				   T6B = FMA(KP831469612, T6A, T6z);
				   T6C = T3d + T3a;
				   T3e = T3a - T3d;
				   T3I = FMA(KP831469612, T3r, T3q);
				   T3s = FNMS(KP831469612, T3r, T3q);
				   {
					E T3E, T3F, T3i, T3l;
					T3E = FMA(KP831469612, T3h, T3g);
					T3i = FNMS(KP831469612, T3h, T3g);
					T3l = FNMS(KP831469612, T3k, T3j);
					T3F = FMA(KP831469612, T3k, T3j);
					T3M = FNMS(KP148335987, T3H, T3I);
					T3J = FMA(KP148335987, T3I, T3H);
					T3w = FMA(KP741650546, T3i, T3l);
					T3m = FNMS(KP741650546, T3l, T3i);
					T6I = T3A + T3B;
					T3C = T3A - T3B;
					T3N = FNMS(KP148335987, T3E, T3F);
					T3G = FMA(KP148335987, T3F, T3E);
				   }
			      }
			      {
				   E T3v, T3f, T6J, T6L, T3x, T3t;
				   T3v = FNMS(KP956940335, T3e, T37);
				   T3f = FMA(KP956940335, T3e, T37);
				   T6J = FMA(KP956940335, T6I, T6H);
				   T6L = FNMS(KP956940335, T6I, T6H);
				   T3x = FMA(KP741650546, T3p, T3s);
				   T3t = FNMS(KP741650546, T3s, T3p);
				   {
					E T3D, T3K, T6D, T6E;
					T3L = FNMS(KP956940335, T3C, T3z);
					T3D = FMA(KP956940335, T3C, T3z);
					{
					     E T3y, T6K, T6M, T3u;
					     T3y = T3w - T3x;
					     T6K = T3w + T3x;
					     T6M = T3m + T3t;
					     T3u = T3m - T3t;
					     Cr[WS(csr, 9)] = FMA(KP803207531, T3y, T3v);
					     Cr[WS(csr, 22)] = FNMS(KP803207531, T3y, T3v);
					     Ci[WS(csi, 25)] = FNMS(KP803207531, T6K, T6J);
					     Ci[WS(csi, 6)] = -(FMA(KP803207531, T6K, T6J));
					     Ci[WS(csi, 9)] = FNMS(KP803207531, T6M, T6L);
					     Ci[WS(csi, 22)] = -(FMA(KP803207531, T6M, T6L));
					     Cr[WS(csr, 6)] = FMA(KP803207531, T3u, T3f);
					     Cr[WS(csr, 25)] = FNMS(KP803207531, T3u, T3f);
					     T3K = T3G - T3J;
					     T6G = T3G + T3J;
					}
					T6F = FNMS(KP956940335, T6C, T6B);
					T6D = FMA(KP956940335, T6C, T6B);
					T6E = T3N + T3M;
					T3O = T3M - T3N;
					Cr[WS(csr, 1)] = FMA(KP989176509, T3K, T3D);
					Cr[WS(csr, 30)] = FNMS(KP989176509, T3K, T3D);
					Ci[WS(csi, 1)] = FMA(KP989176509, T6E, T6D);
					Ci[WS(csi, 30)] = FMS(KP989176509, T6E, T6D);
				   }
			      }
			 }
			 {
			      E T41, T4Z, T4L, T6r, T6l, T6m, T4g, T4S, T5c, T59, T4W, T4A, T6s, T52, T5d;
			      E T56;
			      {
				   E T51, T50, T57, T58, T48, T4f;
				   Cr[WS(csr, 14)] = FMA(KP989176509, T3O, T3L);
				   Cr[WS(csr, 17)] = FNMS(KP989176509, T3O, T3L);
				   Ci[WS(csi, 17)] = FNMS(KP989176509, T6G, T6F);
				   Ci[WS(csi, 14)] = -(FMA(KP989176509, T6G, T6F));
				   T41 = FNMS(KP980785280, T40, T3T);
				   T4Z = FMA(KP980785280, T40, T3T);
				   T51 = FMA(KP820678790, T44, T47);
				   T48 = FNMS(KP820678790, T47, T44);
				   T4f = FNMS(KP820678790, T4e, T4b);
				   T50 = FMA(KP820678790, T4b, T4e);
				   T4L = FNMS(KP980785280, T4K, T4D);
				   T57 = FMA(KP980785280, T4K, T4D);
				   T6r = FMA(KP980785280, T6k, T6j);
				   T6l = FNMS(KP980785280, T6k, T6j);
				   T6m = T4f + T48;
				   T4g = T48 - T4f;
				   T58 = FMA(KP980785280, T4R, T4O);
				   T4S = FNMS(KP980785280, T4R, T4O);
				   {
					E T54, T55, T4s, T4z;
					T54 = FMA(KP980785280, T4r, T4k);
					T4s = FNMS(KP980785280, T4r, T4k);
					T4z = FNMS(KP980785280, T4y, T4v);
					T55 = FMA(KP980785280, T4y, T4v);
					T5c = FMA(KP357805721, T57, T58);
					T59 = FNMS(KP357805721, T58, T57);
					T4W = FMA(KP472964775, T4s, T4z);
					T4A = FNMS(KP472964775, T4z, T4s);
					T6s = T50 + T51;
					T52 = T50 - T51;
					T5d = FNMS(KP357805721, T54, T55);
					T56 = FMA(KP357805721, T55, T54);
				   }
			      }
			      {
				   E T4V, T4h, T6t, T6v, T4X, T4T;
				   T4V = FNMS(KP773010453, T4g, T41);
				   T4h = FMA(KP773010453, T4g, T41);
				   T6t = FMA(KP773010453, T6s, T6r);
				   T6v = FNMS(KP773010453, T6s, T6r);
				   T4X = FNMS(KP472964775, T4L, T4S);
				   T4T = FMA(KP472964775, T4S, T4L);
				   {
					E T53, T5a, T6n, T6o;
					T5b = FNMS(KP773010453, T52, T4Z);
					T53 = FMA(KP773010453, T52, T4Z);
					{
					     E T4Y, T6u, T6w, T4U;
					     T4Y = T4W - T4X;
					     T6u = T4W + T4X;
					     T6w = T4T - T4A;
					     T4U = T4A + T4T;
					     Cr[WS(csr, 11)] = FMA(KP903989293, T4Y, T4V);
					     Cr[WS(csr, 20)] = FNMS(KP903989293, T4Y, T4V);
					     Ci[WS(csi, 27)] = FNMS(KP903989293, T6u, T6t);
					     Ci[WS(csi, 4)] = -(FMA(KP903989293, T6u, T6t));
					     Ci[WS(csi, 11)] = FMA(KP903989293, T6w, T6v);
					     Ci[WS(csi, 20)] = FMS(KP903989293, T6w, T6v);
					     Cr[WS(csr, 4)] = FMA(KP903989293, T4U, T4h);
					     Cr[WS(csr, 27)] = FNMS(KP903989293, T4U, T4h);
					     T5a = T56 + T59;
					     T6q = T59 - T56;
					}
					T6p = FNMS(KP773010453, T6m, T6l);
					T6n = FMA(KP773010453, T6m, T6l);
					T6o = T5d + T5c;
					T5e = T5c - T5d;
					Cr[WS(csr, 3)] = FMA(KP941544065, T5a, T53);
					Cr[WS(csr, 28)] = FNMS(KP941544065, T5a, T53);
					Ci[WS(csi, 3)] = FMA(KP941544065, T6o, T6n);
					Ci[WS(csi, 28)] = FMS(KP941544065, T6o, T6n);
				   }
			      }
			 }
		    }
	       }
	  }
	  Cr[WS(csr, 12)] = FMA(KP941544065, T5e, T5b);
	  Cr[WS(csr, 19)] = FNMS(KP941544065, T5e, T5b);
	  Ci[WS(csi, 19)] = FMA(KP941544065, T6q, T6p);
	  Ci[WS(csi, 12)] = FMS(KP941544065, T6q, T6p);
     }
}

static const kr2c_desc desc = { 64, "r2cfII_64", {114, 0, 320, 0}, &GENUS };