n1_13.c

快速fft变换

/*
 * Copyright (c) 2003, 2007-8 Matteo Frigo
 * Copyright (c) 2003, 2007-8 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 Nov 15 20:36:54 EST 2008 */

#include "codelet-dft.h"

#ifdef HAVE_FMA

/* Generated by: ../../../genfft/gen_notw -fma -reorder-insns -schedule-for-pipeline -compact -variables 4 -pipeline-latency 4 -n 13 -name n1_13 -include n.h */

/*
 * This function contains 176 FP additions, 114 FP multiplications,
 * (or, 62 additions, 0 multiplications, 114 fused multiply/add),
 * 87 stack variables, 25 constants, and 52 memory accesses
 */
#include "n.h"

static void n1_13(const R *ri, const R *ii, R *ro, R *io, stride is, stride os, INT v, INT ivs, INT ovs)
{
     DK(KP875502302, +0.875502302409147941146295545768755143177842006);
     DK(KP520028571, +0.520028571888864619117130500499232802493238139);
     DK(KP575140729, +0.575140729474003121368385547455453388461001608);
     DK(KP600477271, +0.600477271932665282925769253334763009352012849);
     DK(KP300462606, +0.300462606288665774426601772289207995520941381);
     DK(KP516520780, +0.516520780623489722840901288569017135705033622);
     DK(KP968287244, +0.968287244361984016049539446938120421179794516);
     DK(KP503537032, +0.503537032863766627246873853868466977093348562);
     DK(KP251768516, +0.251768516431883313623436926934233488546674281);
     DK(KP581704778, +0.581704778510515730456870384989698884939833902);
     DK(KP859542535, +0.859542535098774820163672132761689612766401925);
     DK(KP083333333, +0.083333333333333333333333333333333333333333333);
     DK(KP957805992, +0.957805992594665126462521754605754580515587217);
     DK(KP522026385, +0.522026385161275033714027226654165028300441940);
     DK(KP853480001, +0.853480001859823990758994934970528322872359049);
     DK(KP769338817, +0.769338817572980603471413688209101117038278899);
     DK(KP612264650, +0.612264650376756543746494474777125408779395514);
     DK(KP038632954, +0.038632954644348171955506895830342264440241080);
     DK(KP302775637, +0.302775637731994646559610633735247973125648287);
     DK(KP514918778, +0.514918778086315755491789696138117261566051239);
     DK(KP686558370, +0.686558370781754340655719594850823015421401653);
     DK(KP226109445, +0.226109445035782405468510155372505010481906348);
     DK(KP301479260, +0.301479260047709873958013540496673347309208464);
     DK(KP866025403, +0.866025403784438646763723170752936183471402627);
     DK(KP500000000, +0.500000000000000000000000000000000000000000000);
     INT i;
     for (i = v; i > 0; i = i - 1, ri = ri + ivs, ii = ii + ivs, ro = ro + ovs, io = io + ovs, MAKE_VOLATILE_STRIDE(is), MAKE_VOLATILE_STRIDE(os)) {
	  E T2B, T2H, T2I, T2G;
	  {
	       E T1, T1P, T2n, T2o, To, TH, T2h, T2k, TE, TB, TF, Tw, T2j, T2c, T1m;
	       E T1W, T1X, T1c, T19, T1j, T12, T1f, T21, T24, T27, T1U;
	       T1 = ri[0];
	       T1P = ii[0];
	       {
		    E T2b, Tv, Ts, T2a;
		    {
			 E T2d, Tf, Tq, Ty, Tb, Tr, T6, Tx, Ti, Tt, Tu, Tl;
			 {
			      E T7, T8, T9, Td, Te;
			      Td = ri[WS(is, 8)];
			      Te = ri[WS(is, 5)];
			      T7 = ri[WS(is, 12)];
			      T8 = ri[WS(is, 10)];
			      T9 = ri[WS(is, 4)];
			      T2d = Td - Te;
			      Tf = Td + Te;
			      {
				   E T2, Ta, T3, T4;
				   T2 = ri[WS(is, 1)];
				   Ta = T8 + T9;
				   Tq = T8 - T9;
				   T3 = ri[WS(is, 3)];
				   T4 = ri[WS(is, 9)];
				   {
					E Tg, T5, Th, Tj, Tk;
					Tg = ri[WS(is, 11)];
					Ty = FMS(KP500000000, Ta, T7);
					Tb = T7 + Ta;
					Tr = T4 - T3;
					T5 = T3 + T4;
					Th = ri[WS(is, 6)];
					Tj = ri[WS(is, 7)];
					Tk = ri[WS(is, 2)];
					T6 = T2 + T5;
					Tx = FNMS(KP500000000, T5, T2);
					Ti = Tg + Th;
					Tt = Tg - Th;
					Tu = Tj - Tk;
					Tl = Tj + Tk;
				   }
			      }
			 }
			 {
			      E Tc, Tm, T2e, T2g;
			      Tc = T6 + Tb;
			      T2n = T6 - Tb;
			      T2b = Ti - Tl;
			      Tm = Ti + Tl;
			      T2e = Tt + Tu;
			      Tv = Tt - Tu;
			      Ts = Tq - Tr;
			      T2g = Tr + Tq;
			      {
				   E Tz, TA, Tn, T2f;
				   Tz = Tx - Ty;
				   T2a = Tx + Ty;
				   TA = FNMS(KP500000000, Tm, Tf);
				   Tn = Tf + Tm;
				   T2f = FNMS(KP500000000, T2e, T2d);
				   T2o = T2d + T2e;
				   To = Tc + Tn;
				   TH = Tc - Tn;
				   T2h = FMA(KP866025403, T2g, T2f);
				   T2k = FNMS(KP866025403, T2g, T2f);
				   TE = Tz - TA;
				   TB = Tz + TA;
			      }
			 }
		    }
		    {
			 E T1R, TM, T10, T18, T1l, TX, T1k, T15, TP, T1a, T1b, TS;
			 {
			      E T16, TY, TZ, TK, TL;
			      TK = ii[WS(is, 8)];
			      TF = Ts - Tv;
			      Tw = Ts + Tv;
			      T2j = FNMS(KP866025403, T2b, T2a);
			      T2c = FMA(KP866025403, T2b, T2a);
			      TL = ii[WS(is, 5)];
			      T16 = ii[WS(is, 12)];
			      TY = ii[WS(is, 10)];
			      TZ = ii[WS(is, 4)];
			      T1R = TK + TL;
			      TM = TK - TL;
			      {
				   E T13, T17, TV, TW;
				   T13 = ii[WS(is, 1)];
				   T17 = TY + TZ;
				   T10 = TY - TZ;
				   TV = ii[WS(is, 9)];
				   TW = ii[WS(is, 3)];
				   {
					E TN, T14, TO, TQ, TR;
					TN = ii[WS(is, 11)];
					T18 = FMS(KP500000000, T17, T16);
					T1l = T16 + T17;
					TX = TV - TW;
					T14 = TW + TV;
					TO = ii[WS(is, 6)];
					TQ = ii[WS(is, 7)];
					TR = ii[WS(is, 2)];
					T1k = T13 + T14;
					T15 = FNMS(KP500000000, T14, T13);
					TP = TN - TO;
					T1a = TN + TO;
					T1b = TQ + TR;
					TS = TQ - TR;
				   }
			      }
			 }
			 {
			      E T1Q, T11, TT, T1S;
			      T1Q = T1k + T1l;
			      T1m = T1k - T1l;
			      T11 = TX + T10;
			      T1W = T10 - TX;
			      T1X = TP - TS;
			      TT = TP + TS;
			      T1S = T1a + T1b;
			      T1c = T1a - T1b;
			      {
				   E T1Z, TU, T1T, T20;
				   T19 = T15 + T18;
				   T1Z = T15 - T18;
				   T1j = TM + TT;
				   TU = FNMS(KP500000000, TT, TM);
				   T1T = T1R + T1S;
				   T20 = FNMS(KP500000000, T1S, T1R);
				   T12 = FMA(KP866025403, T11, TU);
				   T1f = FNMS(KP866025403, T11, TU);
				   T21 = T1Z + T20;
				   T24 = T1Z - T20;
				   T27 = T1Q - T1T;
				   T1U = T1Q + T1T;
			      }
			 }
		    }
	       }
	       {
		    E T1g, T1d, T25, T1Y;
		    ro[0] = T1 + To;
		    T1g = FNMS(KP866025403, T1c, T19);
		    T1d = FMA(KP866025403, T1c, T19);
		    T25 = T1W - T1X;
		    T1Y = T1W + T1X;
		    io[0] = T1P + T1U;
		    {
			 E T1C, T1B, T1F, T1K;
			 {
			      E TC, T1J, T1z, T1w, T1I, T1O, Tp, T1E, T1q, TI, T1o, T1s;
			      {
				   E TG, T1n, T1G, T1u, T1e, T1h, T1v, T1x, T1y, T1H, T1i;
				   TC = FMA(KP301479260, TB, Tw);
				   T1x = FNMS(KP226109445, Tw, TB);
				   T1y = FMA(KP686558370, TE, TF);
				   TG = FNMS(KP514918778, TF, TE);
				   T1n = FNMS(KP302775637, T1m, T1j);
				   T1G = FMA(KP302775637, T1j, T1m);
				   T1u = FNMS(KP038632954, T12, T1d);
				   T1e = FMA(KP038632954, T1d, T12);
				   T1h = FMA(KP612264650, T1g, T1f);
				   T1v = FNMS(KP612264650, T1f, T1g);
				   T1J = FMA(KP769338817, T1y, T1x);
				   T1z = FNMS(KP769338817, T1y, T1x);
				   T1H = FNMS(KP853480001, T1v, T1u);
				   T1w = FMA(KP853480001, T1v, T1u);
				   T1I = FNMS(KP522026385, T1H, T1G);
				   T1O = FMA(KP957805992, T1G, T1H);
				   Tp = FNMS(KP083333333, To, T1);
				   T1E = FMA(KP853480001, T1h, T1e);
				   T1i = FNMS(KP853480001, T1h, T1e);
				   T1q = FNMS(KP859542535, TG, TH);
				   TI = FMA(KP581704778, TH, TG);
				   T1o = FMA(KP957805992, T1n, T1i);
				   T1s = FNMS(KP522026385, T1i, T1n);
			      }
			      {
				   E T1A, T1D, T1t, T1L, T1M;
				   {
					E T1p, TD, TJ, T1N, T1r;
					T1p = FNMS(KP251768516, TC, Tp);
					TD = FMA(KP503537032, TC, Tp);
					T1C = FNMS(KP968287244, T1z, T1w);
					T1A = FMA(KP968287244, T1z, T1w);
					TJ = FMA(KP516520780, TI, TD);
					T1N = FNMS(KP516520780, TI, TD);
					T1D = FNMS(KP300462606, T1q, T1p);
					T1r = FMA(KP300462606, T1q, T1p);
					ro[WS(os, 8)] = FNMS(KP600477271, T1O, T1N);
					ro[WS(os, 12)] = FMA(KP600477271, T1o, TJ);
					ro[WS(os, 1)] = FNMS(KP600477271, T1o, TJ);
					T1t = FNMS(KP575140729, T1s, T1r);
					T1B = FMA(KP575140729, T1s, T1r);
					ro[WS(os, 5)] = FMA(KP600477271, T1O, T1N);
				   }
				   T1L = FNMS(KP520028571, T1E, T1D);
				   T1F = FMA(KP520028571, T1E, T1D);
				   T1K = FMA(KP875502302, T1J, T1I);
				   T1M = FNMS(KP875502302, T1J, T1I);
				   ro[WS(os, 3)] = FMA(KP520028571, T1A, T1t);
				   ro[WS(os, 9)] = FNMS(KP520028571, T1A, T1t);
				   ro[WS(os, 6)] = FMA(KP575140729, T1M, T1L);
				   ro[WS(os, 11)] = FNMS(KP575140729, T1M, T1L);
			      }
			 }
			 {
			      E T22, T2F, T2N, T2K, T2w, T2A, T1V, T2C, T28, T2y, T2M, T2q;
			      {
				   E T26, T2v, T2p, T2i, T2s, T2t, T2l, T2D, T2E, T2u, T2m;
				   T2D = FNMS(KP226109445, T1Y, T21);
				   T22 = FMA(KP301479260, T21, T1Y);
				   ro[WS(os, 2)] = FMA(KP575140729, T1K, T1F);
				   ro[WS(os, 7)] = FNMS(KP575140729, T1K, T1F);
				   ro[WS(os, 4)] = FMA(KP520028571, T1C, T1B);
				   ro[WS(os, 10)] = FNMS(KP520028571, T1C, T1B);
				   T26 = FNMS(KP514918778, T25, T24);
				   T2E = FMA(KP686558370, T24, T25);
				   T2v = FNMS(KP302775637, T2n, T2o);
				   T2p = FMA(KP302775637, T2o, T2n);
				   T2i = FNMS(KP038632954, T2h, T2c);
				   T2s = FMA(KP038632954, T2c, T2h);
				   T2t = FMA(KP612264650, T2j, T2k);
				   T2l = FNMS(KP612264650, T2k, T2j);
				   T2F = FNMS(KP769338817, T2E, T2D);
				   T2N = FMA(KP769338817, T2E, T2D);
				   T2K = FMA(KP853480001, T2t, T2s);
				   T2u = FNMS(KP853480001, T2t, T2s);
				   T2w = FMA(KP957805992, T2v, T2u);
				   T2A = FNMS(KP522026385, T2u, T2v);
				   T1V = FNMS(KP083333333, T1U, T1P);
				   T2m = FNMS(KP853480001, T2l, T2i);
				   T2C = FMA(KP853480001, T2l, T2i);
				   T28 = FMA(KP581704778, T27, T26);
				   T2y = FNMS(KP859542535, T26, T27);
				   T2M = FNMS(KP522026385, T2m, T2p);
				   T2q = FMA(KP957805992, T2p, T2m);
			      }
			      {
				   E T2O, T2Q, T2z, T2P, T2L;
				   {
					E T23, T2x, T2r, T29, T2J;
					T23 = FMA(KP503537032, T22, T1V);
					T2x = FNMS(KP251768516, T22, T1V);
					T2O = FNMS(KP875502302, T2N, T2M);
					T2Q = FMA(KP875502302, T2N, T2M);
					T2r = FMA(KP516520780, T28, T23);
					T29 = FNMS(KP516520780, T28, T23);
					T2z = FMA(KP300462606, T2y, T2x);
					T2J = FNMS(KP300462606, T2y, T2x);
					io[WS(os, 12)] = FNMS(KP600477271, T2w, T2r);
					io[WS(os, 1)] = FMA(KP600477271, T2w, T2r);
					io[WS(os, 8)] = FMA(KP600477271, T2q, T29);
					io[WS(os, 5)] = FNMS(KP600477271, T2q, T29);
					T2P = FMA(KP520028571, T2K, T2J);
					T2L = FNMS(KP520028571, T2K, T2J);
				   }
				   T2B = FMA(KP575140729, T2A, T2z);
				   T2H = FNMS(KP575140729, T2A, T2z);
				   io[WS(os, 11)] = FMA(KP575140729, T2Q, T2P);
				   io[WS(os, 6)] = FNMS(KP575140729, T2Q, T2P);
				   io[WS(os, 7)] = FMA(KP575140729, T2O, T2L);
				   io[WS(os, 2)] = FNMS(KP575140729, T2O, T2L);
				   T2I = FMA(KP968287244, T2F, T2C);
				   T2G = FNMS(KP968287244, T2F, T2C);
			      }
			 }
		    }
	       }
	  }