/* * 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 21:09:05 EST 2008 */#include "codelet-rdft.h"#ifdef HAVE_FMA/* Generated by: ../../../genfft/gen_r2cb -fma -reorder-insns -schedule-for-pipeline -compact -variables 4 -pipeline-latency 4 -sign 1 -n 32 -name r2cbIII_32 -dft-III -include r2cbIII.h *//* * This function contains 174 FP additions, 100 FP multiplications, * (or, 106 additions, 32 multiplications, 68 fused multiply/add), * 101 stack variables, 18 constants, and 64 memory accesses */#include "r2cbIII.h"static void r2cbIII_32(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs){     DK(KP534511135, +0.534511135950791641089685961295362908582039528);     DK(KP1_763842528, +1.763842528696710059425513727320776699016885241);     DK(KP303346683, +0.303346683607342391675883946941299872384187453);     DK(KP1_913880671, +1.913880671464417729871595773960539938965698411);     DK(KP098491403, +0.098491403357164253077197521291327432293052451);     DK(KP1_990369453, +1.990369453344393772489673906218959843150949737);     DK(KP820678790, +0.820678790828660330972281985331011598767386482);     DK(KP1_546020906, +1.546020906725473921621813219516939601942082586);     DK(KP1_847759065, +1.847759065022573512256366378793576573644833252);     DK(KP923879532, +0.923879532511286756128183189396788286822416626);     DK(KP668178637, +0.668178637919298919997757686523080761552472251);     DK(KP1_662939224, +1.662939224605090474157576755235811513477121624);     DK(KP198912367, +0.198912367379658006911597622644676228597850501);     DK(KP1_961570560, +1.961570560806460898252364472268478073947867462);     DK(KP707106781, +0.707106781186547524400844362104849039284835938);     DK(KP1_414213562, +1.414213562373095048801688724209698078569671875);     DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);     DK(KP414213562, +0.414213562373095048801688724209698078569671875);     INT i;     for (i = v; i > 0; i = i - 1, R0 = R0 + ovs, R1 = R1 + ovs, Cr = Cr + ivs, Ci = Ci + ivs, MAKE_VOLATILE_STRIDE(rs), MAKE_VOLATILE_STRIDE(csr), MAKE_VOLATILE_STRIDE(csi)) {	  E T1N, T1K, T1Q, T1H, T1O, T1P;	  {	       E T1I, T1e, T1Z, T7, T2E, T2i, T1x, Tz, Te, T2j, T22, T2F, T1h, T1y, TK;	       E T1J, Tm, T2B, TX, Tp, T2m, T28, T1M, T1C, T1k, TW, TY, T2a, T14, T15;	       E Ts, TZ;	       {		    E TE, T1g, TJ, T1f;		    {			 E T4, Tv, T3, T2g, T1d, T5, Tw, Tx;			 {			      E T1, T2, T1b, T1c;			      T1 = Cr[0];			      T2 = Cr[WS(csr, 15)];			      T1b = Ci[0];			      T1c = Ci[WS(csi, 15)];			      T4 = Cr[WS(csr, 8)];			      Tv = T1 - T2;			      T3 = T1 + T2;			      T2g = T1c - T1b;			      T1d = T1b + T1c;			      T5 = Cr[WS(csr, 7)];			      Tw = Ci[WS(csi, 8)];			      Tx = Ci[WS(csi, 7)];			 }			 {			      E Tb, TA, Ta, T20, TD, Tc, TG, TH;			      {				   E T8, T9, TB, TC;				   T8 = Cr[WS(csr, 4)];				   {					E T1a, T6, T2h, Ty;					T1a = T4 - T5;					T6 = T4 + T5;					T2h = Tx - Tw;					Ty = Tw + Tx;					T1I = T1a - T1d;					T1e = T1a + T1d;					T1Z = T3 - T6;					T7 = T3 + T6;					T2E = T2h + T2g;					T2i = T2g - T2h;					T1x = Tv + Ty;					Tz = Tv - Ty;					T9 = Cr[WS(csr, 11)];				   }				   TB = Ci[WS(csi, 4)];				   TC = Ci[WS(csi, 11)];				   Tb = Cr[WS(csr, 3)];				   TA = T8 - T9;				   Ta = T8 + T9;				   T20 = TC - TB;				   TD = TB + TC;				   Tc = Cr[WS(csr, 12)];				   TG = Ci[WS(csi, 3)];				   TH = Ci[WS(csi, 12)];			      }			      {				   E TF, Td, T21, TI;				   TE = TA - TD;				   T1g = TA + TD;				   TF = Tb - Tc;				   Td = Tb + Tc;				   T21 = TG - TH;				   TI = TG + TH;				   Te = Ta + Td;				   T2j = Ta - Td;				   T22 = T20 - T21;				   T2F = T20 + T21;				   TJ = TF - TI;				   T1f = TF + TI;			      }			 }		    }		    {			 E TM, Ti, TN, T25, TU, TR, Tl, TO;			 {			      E TS, TT, Tg, Th, Tj, Tk;			      Tg = Cr[WS(csr, 2)];			      Th = Cr[WS(csr, 13)];			      T1h = T1f - T1g;			      T1y = T1g + T1f;			      TK = TE + TJ;			      T1J = TE - TJ;			      TM = Tg - Th;			      Ti = Tg + Th;			      TS = Ci[WS(csi, 2)];			      TT = Ci[WS(csi, 13)];			      Tj = Cr[WS(csr, 10)];			      Tk = Cr[WS(csr, 5)];			      TN = Ci[WS(csi, 10)];			      T25 = TS - TT;			      TU = TS + TT;			      TR = Tj - Tk;			      Tl = Tj + Tk;			      TO = Ci[WS(csi, 5)];			 }			 {			      E T12, T13, Tq, Tr;			      {				   E Tn, T1A, TV, T24, T26, TP, To, T27, T1B, TQ;				   Tn = Cr[WS(csr, 1)];				   T1A = TR - TU;				   TV = TR + TU;				   T24 = Ti - Tl;				   Tm = Ti + Tl;				   T26 = TN - TO;				   TP = TN + TO;				   To = Cr[WS(csr, 14)];				   T12 = Ci[WS(csi, 1)];				   T27 = T25 - T26;				   T2B = T26 + T25;				   T1B = TM + TP;				   TQ = TM - TP;				   TX = Tn - To;				   Tp = Tn + To;				   T2m = T24 + T27;				   T28 = T24 - T27;				   T1M = FNMS(KP414213562, T1A, T1B);				   T1C = FMA(KP414213562, T1B, T1A);				   T1k = FMA(KP414213562, TQ, TV);				   TW = FNMS(KP414213562, TV, TQ);				   T13 = Ci[WS(csi, 14)];			      }			      Tq = Cr[WS(csr, 6)];			      Tr = Cr[WS(csr, 9)];			      TY = Ci[WS(csi, 6)];			      T2a = T13 - T12;			      T14 = T12 + T13;			      T15 = Tq - Tr;			      Ts = Tq + Tr;			      TZ = Ci[WS(csi, 9)];			 }		    }	       }	       {		    E T1L, T1F, T23, T2n, T2k, T2e, T1p, T1t, T1s, T1i, T1o, T19, T1l, T1q;		    {			 E T2z, T2G, T2H, T2C, T1j, T17, T2r, T2s, T2u, T2v, T2K, T2D;			 {			      E T2L, T2d, T2l, T2O;			      {				   E Tf, T2N, Tu, T2M;				   {					E T1D, T16, T29, Tt, T2b, T10;					T2z = T7 - Te;					Tf = T7 + Te;					T1D = T15 + T14;					T16 = T14 - T15;					T29 = Tp - Ts;					Tt = Tp + Ts;					T2b = TY - TZ;					T10 = TY + TZ;					T2N = T2F + T2E;					T2G = T2E - T2F;					T2H = Tm - Tt;					Tu = Tm + Tt;					{					     E T2c, T2A, T1E, T11;					     T2c = T2a - T2b;					     T2A = T2b + T2a;					     T1E = TX + T10;					     T11 = TX - T10;					     T2L = Tf - Tu;					     T2d = T29 + T2c;					     T2l = T29 - T2c;					     T2C = T2A - T2B;					     T2M = T2B + T2A;					     T1L = FMA(KP414213562, T1D, T1E);					     T1F = FNMS(KP414213562, T1E, T1D);					     T1j = FMA(KP414213562, T11, T16);					     T17 = FNMS(KP414213562, T16, T11);					     T2O = T2M + T2N;					}				   }				   R0[0] = KP2_000000000 * (Tf + Tu);				   R0[WS(rs, 8)] = KP2_000000000 * (T2N - T2M);			      }			      T23 = T1Z + T22;			      T2r = T1Z - T22;			      R0[WS(rs, 12)] = KP1_414213562 * (T2O - T2L);			      R0[WS(rs, 4)] = KP1_414213562 * (T2L + T2O);			      T2s = T2m + T2l;			      T2n = T2l - T2m;			      T2k = T2i - T2j;			      T2u = T2j + T2i;			      T2v = T28 - T2d;			      T2e = T28 + T2d;			 }			 {			      E T2y, T2t, T2x, T2w;			      T2y = FMA(KP707106781, T2s, T2r);			      T2t = FNMS(KP707106781, T2s, T2r);			      T2x = FMA(KP707106781, T2v, T2u);			      T2w = FNMS(KP707106781, T2v, T2u);			      R0[WS(rs, 7)] = KP1_961570560 * (FMA(KP198912367, T2y, T2x));			      R0[WS(rs, 15)] = -(KP1_961570560 * (FNMS(KP198912367, T2x, T2y)));			      R0[WS(rs, 11)] = KP1_662939224 * (FNMS(KP668178637, T2t, T2w));			      R0[WS(rs, 3)] = KP1_662939224 * (FMA(KP668178637, T2w, T2t));			      T2K = T2z - T2C;			      T2D = T2z + T2C;			 }			 {			      E TL, T18, T2J, T2I;			      T1p = FNMS(KP707106781, TK, Tz);			      TL = FMA(KP707106781, TK, Tz);			      T18 = TW + T17;			      T1t = TW - T17;			      T1s = FMA(KP707106781, T1h, T1e);			      T1i = FNMS(KP707106781, T1h, T1e);			      T2J = T2H + T2G;			      T2I = T2G - T2H;			      T1o = FNMS(KP923879532, T18, TL);			      T19 = FMA(KP923879532, T18, TL);			      R0[WS(rs, 6)] = KP1_847759065 * (FMA(KP414213562, T2K, T2J));			      R0[WS(rs, 14)] = -(KP1_847759065 * (FNMS(KP414213562, T2J, T2K)));			      R0[WS(rs, 10)] = KP1_847759065 * (FNMS(KP414213562, T2D, T2I));			      R0[WS(rs, 2)] = KP1_847759065 * (FMA(KP414213562, T2I, T2D));			      T1l = T1j - T1k;			      T1q = T1k + T1j;			 }		    }		    {			 E T1z, T1U, T1Y, T1T, T1V, T1G;			 {			      E T1w, T1r, T1n, T1m;			      T1n = FMA(KP923879532, T1l, T1i);			      T1m = FNMS(KP923879532, T1l, T1i);			      T1w = FMA(KP923879532, T1q, T1p);			      T1r = FNMS(KP923879532, T1q, T1p);			      R1[WS(rs, 4)] = -(KP1_546020906 * (FNMS(KP820678790, T1o, T1n)));			      R1[WS(rs, 12)] = -(KP1_546020906 * (FMA(KP820678790, T1n, T1o)));			      R1[WS(rs, 8)] = -(KP1_990369453 * (FMA(KP098491403, T19, T1m)));			      R1[0] = KP1_990369453 * (FNMS(KP098491403, T1m, T19));			      {				   E T1R, T1S, T1v, T1u;				   T1z = FNMS(KP707106781, T1y, T1x);				   T1R = FMA(KP707106781, T1y, T1x);				   T1S = T1M + T1L;				   T1N = T1L - T1M;				   T1K = FNMS(KP707106781, T1J, T1I);				   T1U = FMA(KP707106781, T1J, T1I);				   T1v = FNMS(KP923879532, T1t, T1s);				   T1u = FMA(KP923879532, T1t, T1s);				   T1Y = FMA(KP923879532, T1S, T1R);				   T1T = FNMS(KP923879532, T1S, T1R);				   R1[WS(rs, 6)] = -(KP1_913880671 * (FNMS(KP303346683, T1w, T1v)));				   R1[WS(rs, 14)] = -(KP1_913880671 * (FMA(KP303346683, T1v, T1w)));				   R1[WS(rs, 10)] = -(KP1_763842528 * (FMA(KP534511135, T1r, T1u)));				   R1[WS(rs, 2)] = KP1_763842528 * (FNMS(KP534511135, T1u, T1r));				   T1V = T1C + T1F;				   T1G = T1C - T1F;			      }			 }			 {			      E T2q, T2f, T1X, T1W, T2p, T2o;			      T1X = FMA(KP923879532, T1V, T1U);			      T1W = FNMS(KP923879532, T1V, T1U);			      T2q = FNMS(KP707106781, T2e, T23);			      T2f = FMA(KP707106781, T2e, T23);			      R1[WS(rs, 7)] = KP1_990369453 * (FMA(KP098491403, T1Y, T1X));			      R1[WS(rs, 15)] = -(KP1_990369453 * (FNMS(KP098491403, T1X, T1Y)));			      R1[WS(rs, 11)] = KP1_546020906 * (FNMS(KP820678790, T1T, T1W));			      R1[WS(rs, 3)] = KP1_546020906 * (FMA(KP820678790, T1W, T1T));			      T2p = FNMS(KP707106781, T2n, T2k);			      T2o = FMA(KP707106781, T2n, T2k);			      T1Q = FNMS(KP923879532, T1G, T1z);			      T1H = FMA(KP923879532, T1G, T1z);			      R0[WS(rs, 5)] = KP1_662939224 * (FMA(KP668178637, T2q, T2p));			      R0[WS(rs, 13)] = -(KP1_662939224 * (FNMS(KP668178637, T2p, T2q)));			      R0[WS(rs, 9)] = KP1_961570560 * (FNMS(KP198912367, T2f, T2o));			      R0[WS(rs, 1)] = KP1_961570560 * (FMA(KP198912367, T2o, T2f));			 }		    }	       }	  }	  T1O = FMA(KP923879532, T1N, T1K);	  T1P = FNMS(KP923879532, T1N, T1K);	  R1[WS(rs, 5)] = KP1_763842528 * (FMA(KP534511135, T1Q, T1P));	  R1[WS(rs, 13)] = -(KP1_763842528 * (FNMS(KP534511135, T1P, T1Q)));	  R1[WS(rs, 9)] = KP1_913880671 * (FNMS(KP303346683, T1H, T1O));	  R1[WS(rs, 1)] = KP1_913880671 * (FMA(KP303346683, T1O, T1H));     }}static const kr2c_desc desc = { 32, "r2cbIII_32", {106, 32, 68, 0}, &GENUS };void X(codelet_r2cbIII_32) (planner *p) {     X(kr2c_register) (p, r2cbIII_32, &desc);