mhc2r_64.c

fftw-3.0.1

/*
 * 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 22:11:18 EDT 2003 */

#include "codelet-rdft.h"

/* Generated by: /homee/stevenj/cvs/fftw3.0.1/genfft/gen_hc2r_noinline -compact -variables 4 -sign 1 -n 64 -name mhc2r_64 -include hc2r.h */

/*
 * This function contains 394 FP additions, 134 FP multiplications,
 * (or, 342 additions, 82 multiplications, 52 fused multiply/add),
 * 109 stack variables, and 128 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_hc2r_noinline.ml,v 1.1 2003/04/17 19:25:50 athena Exp $
 */

#include "hc2r.h"

static void mhc2r_64_0(const R *ri, const R *ii, R *O, stride ris, stride iis, stride os)
{
     DK(KP1_268786568, +1.268786568327290996430343226450986741351374190);
     DK(KP1_546020906, +1.546020906725473921621813219516939601942082586);
     DK(KP196034280, +0.196034280659121203988391127777283691722273346);
     DK(KP1_990369453, +1.990369453344393772489673906218959843150949737);
     DK(KP942793473, +0.942793473651995297112775251810508755314920638);
     DK(KP1_763842528, +1.763842528696710059425513727320776699016885241);
     DK(KP580569354, +0.580569354508924735272384751634790549382952557);
     DK(KP1_913880671, +1.913880671464417729871595773960539938965698411);
     DK(KP1_111140466, +1.111140466039204449485661627897065748749874382);
     DK(KP1_662939224, +1.662939224605090474157576755235811513477121624);
     DK(KP390180644, +0.390180644032256535696569736954044481855383236);
     DK(KP1_961570560, +1.961570560806460898252364472268478073947867462);
     DK(KP923879532, +0.923879532511286756128183189396788286822416626);
     DK(KP382683432, +0.382683432365089771728459984030398866761344562);
     DK(KP707106781, +0.707106781186547524400844362104849039284835938);
     DK(KP765366864, +0.765366864730179543456919968060797733522689125);
     DK(KP1_847759065, +1.847759065022573512256366378793576573644833252);
     DK(KP1_414213562, +1.414213562373095048801688724209698078569671875);
     DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
     {
	  E Ta, T2S, T18, T2u, T3F, T4V, T5l, T61, Th, T2T, T1h, T2v, T3M, T4W, T5o;
	  E T62, T3Q, T5q, T5u, T44, Tp, Tw, T2V, T2W, T2X, T2Y, T3X, T5t, T1r, T2x;
	  E T41, T5r, T1A, T2y, T4a, T5y, T5N, T4H, TN, T31, T4E, T5z, T39, T3q, T1L;
	  E T2B, T4h, T5M, T2h, T2F, T12, T36, T5D, T5J, T5G, T5K, T1U, T26, T23, T27;
	  E T4p, T4z, T4w, T4A, T34, T3r;
	  {
	       E T5, T3A, T3, T3y, T9, T3C, T17, T3D, T6, T14;
	       {
		    E T4, T3z, T1, T2;
		    T4 = ri[WS(ris, 16)];
		    T5 = KP2_000000000 * T4;
		    T3z = ii[WS(iis, 16)];
		    T3A = KP2_000000000 * T3z;
		    T1 = ri[0];
		    T2 = ri[WS(ris, 32)];
		    T3 = T1 + T2;
		    T3y = T1 - T2;
		    {
			 E T7, T8, T15, T16;
			 T7 = ri[WS(ris, 8)];
			 T8 = ri[WS(ris, 24)];
			 T9 = KP2_000000000 * (T7 + T8);
			 T3C = T7 - T8;
			 T15 = ii[WS(iis, 8)];
			 T16 = ii[WS(iis, 24)];
			 T17 = KP2_000000000 * (T15 - T16);
			 T3D = T15 + T16;
		    }
	       }
	       T6 = T3 + T5;
	       Ta = T6 + T9;
	       T2S = T6 - T9;
	       T14 = T3 - T5;
	       T18 = T14 - T17;
	       T2u = T14 + T17;
	       {
		    E T3B, T3E, T5j, T5k;
		    T3B = T3y - T3A;
		    T3E = KP1_414213562 * (T3C - T3D);
		    T3F = T3B + T3E;
		    T4V = T3B - T3E;
		    T5j = T3y + T3A;
		    T5k = KP1_414213562 * (T3C + T3D);
		    T5l = T5j - T5k;
		    T61 = T5j + T5k;
	       }
	  }
	  {
	       E Td, T3G, T1c, T3K, Tg, T3J, T1f, T3H, T19, T1g;
	       {
		    E Tb, Tc, T1a, T1b;
		    Tb = ri[WS(ris, 4)];
		    Tc = ri[WS(ris, 28)];
		    Td = Tb + Tc;
		    T3G = Tb - Tc;
		    T1a = ii[WS(iis, 4)];
		    T1b = ii[WS(iis, 28)];
		    T1c = T1a - T1b;
		    T3K = T1a + T1b;
	       }
	       {
		    E Te, Tf, T1d, T1e;
		    Te = ri[WS(ris, 20)];
		    Tf = ri[WS(ris, 12)];
		    Tg = Te + Tf;
		    T3J = Te - Tf;
		    T1d = ii[WS(iis, 20)];
		    T1e = ii[WS(iis, 12)];
		    T1f = T1d - T1e;
		    T3H = T1d + T1e;
	       }
	       Th = KP2_000000000 * (Td + Tg);
	       T2T = KP2_000000000 * (T1f + T1c);
	       T19 = Td - Tg;
	       T1g = T1c - T1f;
	       T1h = KP1_414213562 * (T19 - T1g);
	       T2v = KP1_414213562 * (T19 + T1g);
	       {
		    E T3I, T3L, T5m, T5n;
		    T3I = T3G - T3H;
		    T3L = T3J + T3K;
		    T3M = FNMS(KP765366864, T3L, KP1_847759065 * T3I);
		    T4W = FMA(KP765366864, T3I, KP1_847759065 * T3L);
		    T5m = T3G + T3H;
		    T5n = T3K - T3J;
		    T5o = FNMS(KP1_847759065, T5n, KP765366864 * T5m);
		    T62 = FMA(KP1_847759065, T5m, KP765366864 * T5n);
	       }
	  }
	  {
	       E Tl, T3O, T1v, T43, To, T42, T1y, T3P, Ts, T3R, T1p, T3S, Tv, T3U, T1m;
	       E T3V;
	       {
		    E Tj, Tk, T1t, T1u;
		    Tj = ri[WS(ris, 2)];
		    Tk = ri[WS(ris, 30)];
		    Tl = Tj + Tk;
		    T3O = Tj - Tk;
		    T1t = ii[WS(iis, 2)];
		    T1u = ii[WS(iis, 30)];
		    T1v = T1t - T1u;
		    T43 = T1t + T1u;
	       }
	       {
		    E Tm, Tn, T1w, T1x;
		    Tm = ri[WS(ris, 18)];
		    Tn = ri[WS(ris, 14)];
		    To = Tm + Tn;
		    T42 = Tm - Tn;
		    T1w = ii[WS(iis, 18)];
		    T1x = ii[WS(iis, 14)];
		    T1y = T1w - T1x;
		    T3P = T1w + T1x;
	       }
	       {
		    E Tq, Tr, T1n, T1o;
		    Tq = ri[WS(ris, 10)];
		    Tr = ri[WS(ris, 22)];
		    Ts = Tq + Tr;
		    T3R = Tq - Tr;
		    T1n = ii[WS(iis, 10)];
		    T1o = ii[WS(iis, 22)];
		    T1p = T1n - T1o;
		    T3S = T1n + T1o;
	       }
	       {
		    E Tt, Tu, T1k, T1l;
		    Tt = ri[WS(ris, 6)];
		    Tu = ri[WS(ris, 26)];
		    Tv = Tt + Tu;
		    T3U = Tt - Tu;
		    T1k = ii[WS(iis, 26)];
		    T1l = ii[WS(iis, 6)];
		    T1m = T1k - T1l;
		    T3V = T1l + T1k;
	       }
	       T3Q = T3O - T3P;
	       T5q = T3O + T3P;
	       T5u = T43 - T42;
	       T44 = T42 + T43;
	       Tp = Tl + To;
	       Tw = Ts + Tv;
	       T2V = Tp - Tw;
	       {
		    E T3T, T3W, T1j, T1q;
		    T2W = T1y + T1v;
		    T2X = T1p + T1m;
		    T2Y = T2W - T2X;
		    T3T = T3R - T3S;
		    T3W = T3U - T3V;
		    T3X = KP707106781 * (T3T + T3W);
		    T5t = KP707106781 * (T3T - T3W);
		    T1j = Tl - To;
		    T1q = T1m - T1p;
		    T1r = T1j + T1q;
		    T2x = T1j - T1q;
		    {
			 E T3Z, T40, T1s, T1z;
			 T3Z = T3R + T3S;
			 T40 = T3U + T3V;
			 T41 = KP707106781 * (T3Z - T40);
			 T5r = KP707106781 * (T3Z + T40);
			 T1s = Ts - Tv;
			 T1z = T1v - T1y;
			 T1A = T1s + T1z;
			 T2y = T1z - T1s;
		    }
	       }
	  }
	  {
	       E TB, T48, T2c, T4G, TE, T4F, T2f, T49, TI, T4b, T1J, T4c, TL, T4e, T1G;
	       E T4f;
	       {
		    E Tz, TA, T2a, T2b;
		    Tz = ri[WS(ris, 1)];
		    TA = ri[WS(ris, 31)];
		    TB = Tz + TA;
		    T48 = Tz - TA;
		    T2a = ii[WS(iis, 1)];
		    T2b = ii[WS(iis, 31)];
		    T2c = T2a - T2b;
		    T4G = T2a + T2b;
	       }
	       {
		    E TC, TD, T2d, T2e;
		    TC = ri[WS(ris, 17)];
		    TD = ri[WS(ris, 15)];
		    TE = TC + TD;
		    T4F = TC - TD;
		    T2d = ii[WS(iis, 17)];
		    T2e = ii[WS(iis, 15)];
		    T2f = T2d - T2e;
		    T49 = T2d + T2e;
	       }
	       {
		    E TG, TH, T1H, T1I;
		    TG = ri[WS(ris, 9)];
		    TH = ri[WS(ris, 23)];
		    TI = TG + TH;
		    T4b = TG - TH;
		    T1H = ii[WS(iis, 9)];
		    T1I = ii[WS(iis, 23)];
		    T1J = T1H - T1I;
		    T4c = T1H + T1I;
	       }
	       {
		    E TJ, TK, T1E, T1F;
		    TJ = ri[WS(ris, 7)];
		    TK = ri[WS(ris, 25)];
		    TL = TJ + TK;
		    T4e = TJ - TK;
		    T1E = ii[WS(iis, 25)];
		    T1F = ii[WS(iis, 7)];
		    T1G = T1E - T1F;
		    T4f = T1F + T1E;
	       }
	       {
		    E TF, TM, T1D, T1K;
		    T4a = T48 - T49;
		    T5y = T48 + T49;
		    T5N = T4G - T4F;
		    T4H = T4F + T4G;
		    TF = TB + TE;
		    TM = TI + TL;
		    TN = TF + TM;
		    T31 = TF - TM;
		    {
			 E T4C, T4D, T37, T38;
			 T4C = T4b + T4c;
			 T4D = T4e + T4f;
			 T4E = KP707106781 * (T4C - T4D);
			 T5z = KP707106781 * (T4C + T4D);
			 T37 = T2f + T2c;
			 T38 = T1J + T1G;
			 T39 = T37 - T38;
			 T3q = T38 + T37;
		    }
		    T1D = TB - TE;
		    T1K = T1G - T1J;
		    T1L = T1D + T1K;
		    T2B = T1D - T1K;
		    {
			 E T4d, T4g, T29, T2g;
			 T4d = T4b - T4c;
			 T4g = T4e - T4f;
			 T4h = KP707106781 * (T4d + T4g);
			 T5M = KP707106781 * (T4d - T4g);
			 T29 = TI - TL;
			 T2g = T2c - T2f;
			 T2h = T29 + T2g;
			 T2F = T2g - T29;
		    }
	       }
	  }
	  {
	       E TQ, T4j, T1P, T4n, TT, T4m, T1S, T4k, TX, T4q, T1Y, T4u, T10, T4t, T21;
	       E T4r;
	       {
		    E TO, TP, T1N, T1O;
		    TO = ri[WS(ris, 5)];
		    TP = ri[WS(ris, 27)];
		    TQ = TO + TP;
		    T4j = TO - TP;
		    T1N = ii[WS(iis, 5)];
		    T1O = ii[WS(iis, 27)];
		    T1P = T1N - T1O;
		    T4n = T1N + T1O;
	       }
	       {
		    E TR, TS, T1Q, T1R;
		    TR = ri[WS(ris, 21)];
		    TS = ri[WS(ris, 11)];
		    TT = TR + TS;
		    T4m = TR - TS;
		    T1Q = ii[WS(iis, 21)];
		    T1R = ii[WS(iis, 11)];
		    T1S = T1Q - T1R;
		    T4k = T1Q + T1R;
	       }
	       {
		    E TV, TW, T1W, T1X;
		    TV = ri[WS(ris, 3)];
		    TW = ri[WS(ris, 29)];
		    TX = TV + TW;
		    T4q = TV - TW;
		    T1W = ii[WS(iis, 29)];
		    T1X = ii[WS(iis, 3)];
		    T1Y = T1W - T1X;
		    T4u = T1X + T1W;
	       }
	       {
		    E TY, TZ, T1Z, T20;
		    TY = ri[WS(ris, 13)];
		    TZ = ri[WS(ris, 19)];
		    T10 = TY + TZ;
		    T4t = TY - TZ;
		    T1Z = ii[WS(iis, 13)];
		    T20 = ii[WS(iis, 19)];
		    T21 = T1Z - T20;
		    T4r = T1Z + T20;
	       }
	       {
		    E TU, T11, T5B, T5C;
		    TU = TQ + TT;