n2fv_64.c

快速傅立叶变换库函数

						       T4D = VADD(T4B, T4C);
						       T13 = VADD(T11, T12);
						       T4F = VSUB(T11, T12);
						       T16 = VADD(T14, T15);
						       T4G = VSUB(T14, T15);
						       T1a = VADD(T18, T19);
						       T4I = VSUB(T18, T19);
						       T4J = VSUB(T1b, T1c);
						       T1d = VADD(T1b, T1c);
						  }
					     }
					}
				   }
				   {
					V T30, T10, T6k, T4E, T4Q, T4H, T17, T6n, T4P, T1e, T4K, T4R, T1n, T31;
					T30 = VADD(TW, TZ);
					T10 = VSUB(TW, TZ);
					T6k = VFNMS(LDK(KP707106781), T4D, T4A);
					T4E = VFMA(LDK(KP707106781), T4D, T4A);
					T4Q = VFMA(LDK(KP414213562), T4F, T4G);
					T4H = VFNMS(LDK(KP414213562), T4G, T4F);
					T33 = VADD(T13, T16);
					T17 = VSUB(T13, T16);
					T6n = VFNMS(LDK(KP707106781), T4O, T4N);
					T4P = VFMA(LDK(KP707106781), T4O, T4N);
					T34 = VADD(T1a, T1d);
					T1e = VSUB(T1a, T1d);
					T4K = VFMA(LDK(KP414213562), T4J, T4I);
					T4R = VFNMS(LDK(KP414213562), T4I, T4J);
					T1n = VSUB(T1j, T1m);
					T31 = VADD(T1j, T1m);
					{
					     V T1f, T1o, T6o, T4L, T4S, T6l;
					     T1f = VADD(T17, T1e);
					     T1o = VSUB(T17, T1e);
					     T6o = VSUB(T4H, T4K);
					     T4L = VADD(T4H, T4K);
					     T4S = VADD(T4Q, T4R);
					     T6l = VSUB(T4Q, T4R);
					     T3E = VSUB(T30, T31);
					     T32 = VADD(T30, T31);
					     T1p = VFMA(LDK(KP707106781), T1o, T1n);
					     T2v = VFNMS(LDK(KP707106781), T1o, T1n);
					     T1g = VFMA(LDK(KP707106781), T1f, T10);
					     T2u = VFNMS(LDK(KP707106781), T1f, T10);
					     T4M = VFMA(LDK(KP923879532), T4L, T4E);
					     T5K = VFNMS(LDK(KP923879532), T4L, T4E);
					     T6p = VFMA(LDK(KP923879532), T6o, T6n);
					     T6Z = VFNMS(LDK(KP923879532), T6o, T6n);
					     T6m = VFNMS(LDK(KP923879532), T6l, T6k);
					     T6Y = VFMA(LDK(KP923879532), T6l, T6k);
					     T5L = VFNMS(LDK(KP923879532), T4S, T4P);
					     T4T = VFMA(LDK(KP923879532), T4S, T4P);
					}
				   }
			      }
			 }
		    }
	       }
	       {
		    V T6b, T6F, T7n, T7o, T7p, T7q, T7v, T7w, T7x, T7y, T7z, T7A, T7B, T7C, T7f;
		    V T6X, T70, T79, T7a, T73, T6C, T76, T77, T6i;
		    {
			 V T2Z, T3r, T3s, T3m, T3d, T3v;
			 T2Z = VSUB(T2V, T2Y);
			 T3r = VADD(T2V, T2Y);
			 T3s = VADD(T3l, T3i);
			 T3m = VSUB(T3i, T3l);
			 T3d = VSUB(T39, T3c);
			 T3v = VADD(T39, T3c);
			 {
			      V T3x, T3t, T3P, T3J, T3D, T3V, T3Q, T3G, T36, T3u, T3Y, T3O, T6V, T6W;
			      {
				   V T3N, T3C, T3F, T35;
				   T3N = VSUB(T3B, T3A);
				   T3C = VADD(T3A, T3B);
				   T3F = VSUB(T33, T34);
				   T35 = VADD(T33, T34);
				   T3x = VSUB(T3r, T3s);
				   T3t = VADD(T3r, T3s);
				   T3P = VFMA(LDK(KP414213562), T3H, T3I);
				   T3J = VFNMS(LDK(KP414213562), T3I, T3H);
				   T3D = VFMA(LDK(KP707106781), T3C, T3z);
				   T3V = VFNMS(LDK(KP707106781), T3C, T3z);
				   T3Q = VFMA(LDK(KP414213562), T3E, T3F);
				   T3G = VFNMS(LDK(KP414213562), T3F, T3E);
				   T36 = VSUB(T32, T35);
				   T3u = VADD(T32, T35);
				   T3Y = VFNMS(LDK(KP707106781), T3N, T3M);
				   T3O = VFMA(LDK(KP707106781), T3N, T3M);
			      }
			      T6b = VFNMS(LDK(KP923879532), T6a, T69);
			      T6V = VFMA(LDK(KP923879532), T6a, T69);
			      T6W = VADD(T6E, T6D);
			      T6F = VSUB(T6D, T6E);
			      {
				   V T3K, T3Z, T3e, T3n;
				   T3K = VADD(T3G, T3J);
				   T3Z = VSUB(T3J, T3G);
				   T3e = VADD(T36, T3d);
				   T3n = VSUB(T3d, T36);
				   {
					V T3w, T3y, T3R, T3W;
					T3w = VADD(T3u, T3v);
					T3y = VSUB(T3v, T3u);
					T3R = VSUB(T3P, T3Q);
					T3W = VADD(T3Q, T3P);
					{
					     V T42, T40, T3L, T3T;
					     T42 = VFNMS(LDK(KP923879532), T3Z, T3Y);
					     T40 = VFMA(LDK(KP923879532), T3Z, T3Y);
					     T3L = VFNMS(LDK(KP923879532), T3K, T3D);
					     T3T = VFMA(LDK(KP923879532), T3K, T3D);
					     {
						  V T3o, T3q, T3f, T3p;
						  T3o = VFNMS(LDK(KP707106781), T3n, T3m);
						  T3q = VFMA(LDK(KP707106781), T3n, T3m);
						  T3f = VFNMS(LDK(KP707106781), T3e, T2Z);
						  T3p = VFMA(LDK(KP707106781), T3e, T2Z);
						  T7n = VFNMSI(T3y, T3x);
						  STM2(&(xo[96]), T7n, ovs, &(xo[0]));
						  T7o = VFMAI(T3y, T3x);
						  STM2(&(xo[32]), T7o, ovs, &(xo[0]));
						  T7p = VADD(T3t, T3w);
						  STM2(&(xo[0]), T7p, ovs, &(xo[0]));
						  T7q = VSUB(T3t, T3w);
						  STM2(&(xo[64]), T7q, ovs, &(xo[0]));
						  {
						       V T41, T3X, T3S, T3U;
						       T41 = VFMA(LDK(KP923879532), T3W, T3V);
						       T3X = VFNMS(LDK(KP923879532), T3W, T3V);
						       T3S = VFNMS(LDK(KP923879532), T3R, T3O);
						       T3U = VFMA(LDK(KP923879532), T3R, T3O);
						       T7r = VFMAI(T3q, T3p);
						       STM2(&(xo[16]), T7r, ovs, &(xo[0]));
						       T7s = VFNMSI(T3q, T3p);
						       STM2(&(xo[112]), T7s, ovs, &(xo[0]));
						       T7t = VFMAI(T3o, T3f);
						       STM2(&(xo[80]), T7t, ovs, &(xo[0]));
						       T7u = VFNMSI(T3o, T3f);
						       STM2(&(xo[48]), T7u, ovs, &(xo[0]));
						       T7v = VFNMSI(T40, T3X);
						       STM2(&(xo[88]), T7v, ovs, &(xo[0]));
						       T7w = VFMAI(T40, T3X);
						       STM2(&(xo[40]), T7w, ovs, &(xo[0]));
						       T7x = VFMAI(T42, T41);
						       STM2(&(xo[104]), T7x, ovs, &(xo[0]));
						       T7y = VFNMSI(T42, T41);
						       STM2(&(xo[24]), T7y, ovs, &(xo[0]));
						       T7z = VFMAI(T3U, T3T);
						       STM2(&(xo[8]), T7z, ovs, &(xo[0]));
						       T7A = VFNMSI(T3U, T3T);
						       STM2(&(xo[120]), T7A, ovs, &(xo[0]));
						       T7B = VFMAI(T3S, T3L);
						       STM2(&(xo[72]), T7B, ovs, &(xo[0]));
						       T7C = VFNMSI(T3S, T3L);
						       STM2(&(xo[56]), T7C, ovs, &(xo[0]));
						       T7f = VFNMS(LDK(KP831469612), T6W, T6V);
						       T6X = VFMA(LDK(KP831469612), T6W, T6V);
						  }
					     }
					}
				   }
			      }
			      T70 = VFMA(LDK(KP303346683), T6Z, T6Y);
			      T79 = VFNMS(LDK(KP303346683), T6Y, T6Z);
			      T7a = VFNMS(LDK(KP303346683), T71, T72);
			      T73 = VFMA(LDK(KP303346683), T72, T71);
			      T6C = VFNMS(LDK(KP923879532), T6B, T6A);
			      T76 = VFMA(LDK(KP923879532), T6B, T6A);
			      T77 = VSUB(T6e, T6h);
			      T6i = VADD(T6e, T6h);
			 }
		    }
		    {
			 V T2r, T2D, T2C, T2s, T5H, T5o, T5v, T5D, T7L, T7O, T7Q, T7S, T5r, T5I, T5x;
			 V T5h, T5F, T5B;
			 {
			      V TT, T2f, T7E, T7F, T7H, T7J, T2n, T1Y, T28, T2b, T2l, T2p, T2j, T2k;
			      {
				   V T1X, T2d, T7h, T7l, T2e, T1q, T75, T7d, T7m, T7k, T7c, T7e, Tn, TS;
				   T2r = VFNMS(LDK(KP707106781), Tm, T7);
				   Tn = VFMA(LDK(KP707106781), Tm, T7);
				   TS = VADD(TC, TR);
				   T2D = VSUB(TR, TC);
				   {
					V T7b, T7j, T74, T7i, T78, T7g;
					T1X = VFNMS(LDK(KP198912367), T1W, T1N);
					T2d = VFMA(LDK(KP198912367), T1N, T1W);
					T7g = VADD(T79, T7a);
					T7b = VSUB(T79, T7a);
					T7j = VSUB(T73, T70);
					T74 = VADD(T70, T73);
					T7i = VFNMS(LDK(KP831469612), T77, T76);
					T78 = VFMA(LDK(KP831469612), T77, T76);
					T2j = VFNMS(LDK(KP923879532), TS, Tn);
					TT = VFMA(LDK(KP923879532), TS, Tn);
					T7h = VFMA(LDK(KP956940335), T7g, T7f);
					T7l = VFNMS(LDK(KP956940335), T7g, T7f);
					T2e = VFMA(LDK(KP198912367), T1g, T1p);
					T1q = VFNMS(LDK(KP198912367), T1p, T1g);
					T75 = VFNMS(LDK(KP956940335), T74, T6X);
					T7d = VFMA(LDK(KP956940335), T74, T6X);
					T7m = VFNMS(LDK(KP956940335), T7j, T7i);
					T7k = VFMA(LDK(KP956940335), T7j, T7i);
					T7c = VFNMS(LDK(KP956940335), T7b, T78);
					T7e = VFMA(LDK(KP956940335), T7b, T78);
				   }
				   T2k = VADD(T2e, T2d);
				   T2f = VSUB(T2d, T2e);
				   {
					V T7D, T7G, T7I, T7K;
					T7D = VFNMSI(T7k, T7h);
					STM2(&(xo[90]), T7D, ovs, &(xo[2]));
					STN2(&(xo[88]), T7v, T7D, ovs);
					T7E = VFMAI(T7k, T7h);
					STM2(&(xo[38]), T7E, ovs, &(xo[2]));
					T7F = VFMAI(T7m, T7l);
					STM2(&(xo[102]), T7F, ovs, &(xo[2]));
					T7G = VFNMSI(T7m, T7l);
					STM2(&(xo[26]), T7G, ovs, &(xo[2]));
					STN2(&(xo[24]), T7y, T7G, ovs);
					T7H = VFMAI(T7e, T7d);
					STM2(&(xo[6]), T7H, ovs, &(xo[2]));
					T7I = VFNMSI(T7e, T7d);
					STM2(&(xo[122]), T7I, ovs, &(xo[2]));
					STN2(&(xo[120]), T7A, T7I, ovs);
					T7J = VFMAI(T7c, T75);
					STM2(&(xo[70]), T7J, ovs, &(xo[2]));
					T7K = VFNMSI(T7c, T75);
					STM2(&(xo[58]), T7K, ovs, &(xo[2]));
					STN2(&(xo[56]), T7C, T7K, ovs);
					T2n = VSUB(T1X, T1q);
					T1Y = VADD(T1q, T1X);
				   }
				   T2C = VFNMS(LDK(KP707106781), T27, T26);
				   T28 = VFMA(LDK(KP707106781), T27, T26);
				   T2b = VSUB(T29, T2a);
				   T2s = VADD(T2a, T29);
			      }
			      T2l = VFNMS(LDK(KP980785280), T2k, T2j);
			      T2p = VFMA(LDK(KP980785280), T2k, T2j);
			      {
				   V T5z, T4z, T5A, T5g;
				   {
					V T4f, T4y, T1Z, T2h, T4U, T5t, T2m, T2c, T5u, T5f;
					T5H = VFNMS(LDK(KP923879532), T4e, T47);
					T4f = VFMA(LDK(KP923879532), T4e, T47);
					T4y = VADD(T4o, T4x);
					T5T = VSUB(T4x, T4o);
					T1Z = VFNMS(LDK(KP980785280), T1Y, TT);
					T2h = VFMA(LDK(KP980785280), T1Y, TT);
					T4U = VFNMS(LDK(KP098491403), T4T, T4M);
					T5t = VFMA(LDK(KP098491403), T4M, T4T);
					T2m = VFNMS(LDK(KP923879532), T2b, T28);
					T2c = VFMA(LDK(KP923879532), T2b, T28);
					T5u = VFMA(LDK(KP098491403), T57, T5e);
					T5f = VFNMS(LDK(KP098491403), T5e, T57);
					T5z = VFNMS(LDK(KP980785280), T4y, T4f);
					T4z = VFMA(LDK(KP980785280), T4y, T4f);
					T5S = VFNMS(LDK(KP923879532), T5n, T5k);
					T5o = VFMA(LDK(KP923879532), T5n, T5k);
					{
					     V T2o, T2q, T2i, T2g;
					     T2o = VFMA(LDK(KP980785280), T2n, T2m);
					     T2q = VFNMS(LDK(KP980785280), T2n, T2m);
					     T2i = VFMA(LDK(KP980785280), T2f, T2c);
					     T2g = VFNMS(LDK(KP980785280), T2f, T2c);
					     T5A = VADD(T5t, T5u);
					     T5v = VSUB(T5t, T5u);
					     T5D = VSUB(T5f, T4U);
					     T5g = VADD(T4U, T5f);
					     T7L = VFNMSI(T2o, T2l);
					     STM2(&(xo[92]), T7L, ovs, &(xo[0]));
					     {
						  V T7M, T7N, T7P, T7R;
						  T7M = VFMAI(T2o, T2l);
						  STM2(&(xo[36]), T7M, ovs, &(xo[0]));
						  STN2(&(xo[36]), T7M, T7E, ovs);
						  T7N = VFMAI(T2q, T2p);
						  STM2(&(xo[100]), T7N, ovs, &(xo[0]));
						  STN2(&(xo[100]), T7N, T7F, ovs);
						  T7O = VFNMSI(T2q, T2p);
						  STM2(&(xo[28]), T7O, ovs, &(xo[0]));
						  T7P = VFMAI(T2i, T2h);
						  STM2(&(xo[4]), T7P, ovs, &(xo[0]));
						  STN2(&(xo[4]), T7P, T7H, ovs);
						  T7Q = VFNMSI(T2i, T2h);
						  STM2(&(xo[124]), T7Q, ovs, &(xo[0]));
						  T7R = VFMAI(T2g, T1Z);
						  STM2(&(xo[68]), T7R, ovs, &(xo[0]));
						  STN2(&(xo[68]), T7R, T7J, ovs);
						  T7S = VFNMSI(T2g, T1Z);
						  STM2(&(xo[60]), T7S, ovs, &(xo[0]));
						  T5r = VSUB(T5p, T5q);
						  T5I = VADD(T5p, T5q);
					     }
					}
				   }
				   T5x = VFMA(LDK(KP995184726), T5g, T4z);
				   T5h = VFNMS(LDK(KP995184726), T5g, T4z);
				   T5F = VFMA(LDK(KP995184726), T5A, T5z);
				   T5B = VFNMS(LDK(KP995184726), T5A, T5z);
			      }
			 }
			 {
			      V T6J, T6R, T6L, T6z, T6T, T6P;
			      {
				   V T6N, T6j, T6O, T6y;
				   {
					V T6q, T6H, T5C, T5s, T6I, T6x;
					T6q = VFNMS(LDK(KP534511135), T6p, T6m);
					T6H = VFMA(LDK(KP534511135), T6m, T6p);
					T5C = VFNMS(LDK(KP980785280), T5r, T5o);
					T5s = VFMA(LDK(KP980785280), T5r, T5o);
					T6I = VFMA(LDK(KP534511135), T6t, T6w);
					T6x = VFNMS(LDK(KP534511135), T6w, T6t);
					T6N = VFMA(LDK(KP831469612), T6i, T6b);
					T6j = VFNMS(LDK(KP831469612), T6i, T6b);
					{
					     V T5E, T5G, T5y, T5w;
					     T5E = VFNMS(LDK(KP995184726), T5D, T5C);
					     T5G = VFMA(LDK(KP995184726), T5D, T5C);
					     T5y = VFMA(LDK(KP995184726), T5v, T5s);
					     T5w = VFNMS(LDK(KP995184726), T5v, T5s);
					     T6O = VADD(T6H, T6I);
					     T6J = VSUB(T6H, T6I);
					     T6R = VSUB(T6x, T6q);
					     T6y = VADD(T6q, T6x);
					     {
						  V T7T, T7U, T7V, T7W;
						  T7T = VFMAI(T5E, T5B);
						  STM2(&(xo[94]), T7T, ovs, &(xo[2]));
						  STN2(&(xo[92]), T7L, T7T, ovs);
						  T7U = VFNMSI(T5E, T5B);
						  STM2(&(xo[34]), T7U, ovs, &(xo[2]));
						  STN2(&(xo[32]), T7o, T7U, ovs);
						  T7V = VFNMSI(T5G, T5F);
						  STM2(&(xo[98]), T7V, ovs, &(xo[2]));
						  STN2(&(xo[96]), T7n, T7V, ovs);
						  T7W = VFMAI(T5G, T5F);
						  STM2(&(xo[30]), T7W, ovs, &(xo[2]));
						  STN2(&(xo[28]), T7O, T7W, ovs);
						  {
						       V T7X, T7Y, T7Z, T80;
						       T7X = VFMAI(T5y, T5x);
						       STM2(&(xo[126]), T7X, ovs, &(xo[2]));
						       STN2(&(xo[124]), T7Q, T7X, ovs);
						       T7Y = VFNMSI(T5y, T5x);
						       STM2(&(xo[2]), T7Y, ovs, &(xo[2]));
						       STN2(&(xo[0]), T7p, T7Y, ovs);
						       T7Z = VFMAI(T5w, T5h);
						       STM2(&(xo[62]), T7Z, ovs, &(xo[2]));
						       STN2(&(xo[60]), T7S, T7Z, ovs);
						       T80 = VFNMSI(T5w, T5h);
						       STM2(&(xo[66]), T80, ovs, &(xo[2]));
						       STN2(&(xo[64]), T7q, T80, ovs);
						  }
					     }
					}
				   }
				   T6L = VFMA(LDK(KP881921264), T6y, T6j);
				   T6z = VFNMS(LDK(KP881921264), T6y, T6j);
				   T6T = VFMA(LDK(KP881921264), T6O, T6N);
				   T6P = VFNMS(LDK(KP881921264), T6O, T6N);