WebSVN - shark - Blame - Rev 1688 - /shark/trunk/ports/fftw/fftw/ftwi

Rev	Author	Line No.	Line
2	pj	1	/*
		2	* Copyright (c) 1997-1999 Massachusetts Institute of Technology
		3	*
		4	* This program is free software; you can redistribute it and/or modify
		5	* it under the terms of the GNU General Public License as published by
		6	* the Free Software Foundation; either version 2 of the License, or
		7	* (at your option) any later version.
		8	*
		9	* This program is distributed in the hope that it will be useful,
		10	* but WITHOUT ANY WARRANTY; without even the implied warranty of
		11	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
		12	* GNU General Public License for more details.
		13	*
		14	* You should have received a copy of the GNU General Public License
		15	* along with this program; if not, write to the Free Software
		16	* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
		17	*
		18	*/
		19
		20	/* This file was automatically generated --- DO NOT EDIT */
		21	/* Generated on Tue May 18 13:55:50 EDT 1999 */
		22
107	pj	23	#include <fftw-int.h>
		24	#include <fftw.h>
2	pj	25
		26	/* Generated by: ./genfft -magic-alignment-check -magic-twiddle-load-all -magic-variables 4 -magic-loopi -twiddleinv 9 */
		27
		28	/*
		29	* This function contains 96 FP additions, 72 FP multiplications,
		30	* (or, 60 additions, 36 multiplications, 36 fused multiply/add),
		31	* 34 stack variables, and 36 memory accesses
		32	*/
		33	static const fftw_real K642787609 = FFTW_KONST(+0.642787609686539326322643409907263432907559884);
		34	static const fftw_real K766044443 = FFTW_KONST(+0.766044443118978035202392650555416673935832457);
		35	static const fftw_real K939692620 = FFTW_KONST(+0.939692620785908384054109277324731469936208134);
		36	static const fftw_real K342020143 = FFTW_KONST(+0.342020143325668733044099614682259580763083368);
		37	static const fftw_real K984807753 = FFTW_KONST(+0.984807753012208059366743024589523013670643252);
		38	static const fftw_real K173648177 = FFTW_KONST(+0.173648177666930348851716626769314796000375677);
		39	static const fftw_real K500000000 = FFTW_KONST(+0.500000000000000000000000000000000000000000000);
		40	static const fftw_real K866025403 = FFTW_KONST(+0.866025403784438646763723170752936183471402627);
		41
		42	/*
		43	* Generator Id's :
107	pj	44	* $Id: ftwi_9.c,v 1.2 2003-03-24 11:14:54 pj Exp $
		45	* $Id: ftwi_9.c,v 1.2 2003-03-24 11:14:54 pj Exp $
		46	* $Id: ftwi_9.c,v 1.2 2003-03-24 11:14:54 pj Exp $
2	pj	47	*/
		48
		49	void fftwi_twiddle_9(fftw_complex A, const fftw_complex W, int iostride, int m, int dist)
		50	{
		51	int i;
		52	fftw_complex *inout;
		53	inout = A;
		54	for (i = m; i > 0; i = i - 1, inout = inout + dist, W = W + 8) {
		55	fftw_real tmp1;
		56	fftw_real tmp99;
		57	fftw_real tmp64;
		58	fftw_real tmp98;
		59	fftw_real tmp105;
		60	fftw_real tmp104;
		61	fftw_real tmp12;
		62	fftw_real tmp61;
		63	fftw_real tmp47;
		64	fftw_real tmp78;
		65	fftw_real tmp89;
		66	fftw_real tmp54;
		67	fftw_real tmp75;
		68	fftw_real tmp90;
		69	fftw_real tmp30;
		70	fftw_real tmp68;
		71	fftw_real tmp86;
		72	fftw_real tmp59;
		73	fftw_real tmp71;
		74	fftw_real tmp87;
		75	ASSERT_ALIGNED_DOUBLE();
		76	{
		77	fftw_real tmp6;
		78	fftw_real tmp63;
		79	fftw_real tmp11;
		80	fftw_real tmp62;
		81	ASSERT_ALIGNED_DOUBLE();
		82	tmp1 = c_re(inout[0]);
		83	tmp99 = c_im(inout[0]);
		84	{
		85	fftw_real tmp3;
		86	fftw_real tmp5;
		87	fftw_real tmp2;
		88	fftw_real tmp4;
		89	ASSERT_ALIGNED_DOUBLE();
		90	tmp3 = c_re(inout[3 * iostride]);
		91	tmp5 = c_im(inout[3 * iostride]);
		92	tmp2 = c_re(W[2]);
		93	tmp4 = c_im(W[2]);
		94	tmp6 = (tmp2 * tmp3) + (tmp4 * tmp5);
		95	tmp63 = (tmp2 * tmp5) - (tmp4 * tmp3);
		96	}
		97	{
		98	fftw_real tmp8;
		99	fftw_real tmp10;
		100	fftw_real tmp7;
		101	fftw_real tmp9;
		102	ASSERT_ALIGNED_DOUBLE();
		103	tmp8 = c_re(inout[6 * iostride]);
		104	tmp10 = c_im(inout[6 * iostride]);
		105	tmp7 = c_re(W[5]);
		106	tmp9 = c_im(W[5]);
		107	tmp11 = (tmp7 * tmp8) + (tmp9 * tmp10);
		108	tmp62 = (tmp7 * tmp10) - (tmp9 * tmp8);
		109	}
		110	tmp64 = K866025403 * (tmp62 - tmp63);
		111	tmp98 = tmp63 + tmp62;
		112	tmp105 = tmp99 - (K500000000 * tmp98);
		113	tmp104 = K866025403 * (tmp6 - tmp11);
		114	tmp12 = tmp6 + tmp11;
		115	tmp61 = tmp1 - (K500000000 * tmp12);
		116	}
		117	{
		118	fftw_real tmp35;
		119	fftw_real tmp50;
		120	fftw_real tmp40;
		121	fftw_real tmp51;
		122	fftw_real tmp45;
		123	fftw_real tmp52;
		124	fftw_real tmp46;
		125	fftw_real tmp53;
		126	ASSERT_ALIGNED_DOUBLE();
		127	{
		128	fftw_real tmp32;
		129	fftw_real tmp34;
		130	fftw_real tmp31;
		131	fftw_real tmp33;
		132	ASSERT_ALIGNED_DOUBLE();
		133	tmp32 = c_re(inout[2 * iostride]);
		134	tmp34 = c_im(inout[2 * iostride]);
		135	tmp31 = c_re(W[1]);
		136	tmp33 = c_im(W[1]);
		137	tmp35 = (tmp31 * tmp32) + (tmp33 * tmp34);
		138	tmp50 = (tmp31 * tmp34) - (tmp33 * tmp32);
		139	}
		140	{
		141	fftw_real tmp37;
		142	fftw_real tmp39;
		143	fftw_real tmp36;
		144	fftw_real tmp38;
		145	ASSERT_ALIGNED_DOUBLE();
		146	tmp37 = c_re(inout[5 * iostride]);
		147	tmp39 = c_im(inout[5 * iostride]);
		148	tmp36 = c_re(W[4]);
		149	tmp38 = c_im(W[4]);
		150	tmp40 = (tmp36 * tmp37) + (tmp38 * tmp39);
		151	tmp51 = (tmp36 * tmp39) - (tmp38 * tmp37);
		152	}
		153	{
		154	fftw_real tmp42;
		155	fftw_real tmp44;
		156	fftw_real tmp41;
		157	fftw_real tmp43;
		158	ASSERT_ALIGNED_DOUBLE();
		159	tmp42 = c_re(inout[8 * iostride]);
		160	tmp44 = c_im(inout[8 * iostride]);
		161	tmp41 = c_re(W[7]);
		162	tmp43 = c_im(W[7]);
		163	tmp45 = (tmp41 * tmp42) + (tmp43 * tmp44);
		164	tmp52 = (tmp41 * tmp44) - (tmp43 * tmp42);
		165	}
		166	tmp46 = tmp40 + tmp45;
		167	tmp53 = tmp51 + tmp52;
		168	{
		169	fftw_real tmp76;
		170	fftw_real tmp77;
		171	fftw_real tmp73;
		172	fftw_real tmp74;
		173	ASSERT_ALIGNED_DOUBLE();
		174	tmp47 = tmp35 + tmp46;
		175	tmp76 = tmp35 - (K500000000 * tmp46);
		176	tmp77 = K866025403 * (tmp52 - tmp51);
		177	tmp78 = tmp76 - tmp77;
		178	tmp89 = tmp76 + tmp77;
		179	tmp54 = tmp50 + tmp53;
		180	tmp73 = tmp50 - (K500000000 * tmp53);
		181	tmp74 = K866025403 * (tmp40 - tmp45);
		182	tmp75 = tmp73 - tmp74;
		183	tmp90 = tmp74 + tmp73;
		184	}
		185	}
		186	{
		187	fftw_real tmp18;
		188	fftw_real tmp55;
		189	fftw_real tmp23;
		190	fftw_real tmp56;
		191	fftw_real tmp28;
		192	fftw_real tmp57;
		193	fftw_real tmp29;
		194	fftw_real tmp58;
		195	ASSERT_ALIGNED_DOUBLE();
		196	{
		197	fftw_real tmp15;
		198	fftw_real tmp17;
		199	fftw_real tmp14;
		200	fftw_real tmp16;
		201	ASSERT_ALIGNED_DOUBLE();
		202	tmp15 = c_re(inout[iostride]);
		203	tmp17 = c_im(inout[iostride]);
		204	tmp14 = c_re(W[0]);
		205	tmp16 = c_im(W[0]);
		206	tmp18 = (tmp14 * tmp15) + (tmp16 * tmp17);
		207	tmp55 = (tmp14 * tmp17) - (tmp16 * tmp15);
		208	}
		209	{
		210	fftw_real tmp20;
		211	fftw_real tmp22;
		212	fftw_real tmp19;
		213	fftw_real tmp21;
		214	ASSERT_ALIGNED_DOUBLE();
		215	tmp20 = c_re(inout[4 * iostride]);
		216	tmp22 = c_im(inout[4 * iostride]);
		217	tmp19 = c_re(W[3]);
		218	tmp21 = c_im(W[3]);
		219	tmp23 = (tmp19 * tmp20) + (tmp21 * tmp22);
		220	tmp56 = (tmp19 * tmp22) - (tmp21 * tmp20);
		221	}
		222	{
		223	fftw_real tmp25;
		224	fftw_real tmp27;
		225	fftw_real tmp24;
		226	fftw_real tmp26;
		227	ASSERT_ALIGNED_DOUBLE();
		228	tmp25 = c_re(inout[7 * iostride]);
		229	tmp27 = c_im(inout[7 * iostride]);
		230	tmp24 = c_re(W[6]);
		231	tmp26 = c_im(W[6]);
		232	tmp28 = (tmp24 * tmp25) + (tmp26 * tmp27);
		233	tmp57 = (tmp24 * tmp27) - (tmp26 * tmp25);
		234	}
		235	tmp29 = tmp23 + tmp28;
		236	tmp58 = tmp56 + tmp57;
		237	{
		238	fftw_real tmp66;
		239	fftw_real tmp67;
		240	fftw_real tmp69;
		241	fftw_real tmp70;
		242	ASSERT_ALIGNED_DOUBLE();
		243	tmp30 = tmp18 + tmp29;
		244	tmp66 = tmp18 - (K500000000 * tmp29);
		245	tmp67 = K866025403 * (tmp57 - tmp56);
		246	tmp68 = tmp66 - tmp67;
		247	tmp86 = tmp66 + tmp67;
		248	tmp59 = tmp55 + tmp58;
		249	tmp69 = tmp55 - (K500000000 * tmp58);
		250	tmp70 = K866025403 * (tmp23 - tmp28);
		251	tmp71 = tmp69 - tmp70;
		252	tmp87 = tmp70 + tmp69;
		253	}
		254	}
		255	{
		256	fftw_real tmp60;
		257	fftw_real tmp13;
		258	fftw_real tmp48;
		259	fftw_real tmp49;
		260	ASSERT_ALIGNED_DOUBLE();
		261	tmp60 = K866025403 * (tmp54 - tmp59);
		262	tmp13 = tmp1 + tmp12;
		263	tmp48 = tmp30 + tmp47;
		264	tmp49 = tmp13 - (K500000000 * tmp48);
		265	c_re(inout[0]) = tmp13 + tmp48;
		266	c_re(inout[3 * iostride]) = tmp49 + tmp60;
		267	c_re(inout[6 * iostride]) = tmp49 - tmp60;
		268	}
		269	{
		270	fftw_real tmp101;
		271	fftw_real tmp97;
		272	fftw_real tmp100;
		273	fftw_real tmp102;
		274	ASSERT_ALIGNED_DOUBLE();
		275	tmp101 = K866025403 * (tmp30 - tmp47);
		276	tmp97 = tmp59 + tmp54;
		277	tmp100 = tmp98 + tmp99;
		278	tmp102 = tmp100 - (K500000000 * tmp97);
		279	c_im(inout[0]) = tmp97 + tmp100;
		280	c_im(inout[6 * iostride]) = tmp102 - tmp101;
		281	c_im(inout[3 * iostride]) = tmp101 + tmp102;
		282	}
		283	{
		284	fftw_real tmp65;
		285	fftw_real tmp110;
		286	fftw_real tmp80;
		287	fftw_real tmp111;
		288	fftw_real tmp84;
		289	fftw_real tmp109;
		290	fftw_real tmp81;
		291	fftw_real tmp112;
		292	ASSERT_ALIGNED_DOUBLE();
		293	tmp65 = tmp61 - tmp64;
		294	tmp110 = tmp105 - tmp104;
		295	{
		296	fftw_real tmp72;
		297	fftw_real tmp79;
		298	fftw_real tmp82;
		299	fftw_real tmp83;
		300	ASSERT_ALIGNED_DOUBLE();
		301	tmp72 = (K173648177 * tmp68) - (K984807753 * tmp71);
		302	tmp79 = (K342020143 * tmp75) + (K939692620 * tmp78);
		303	tmp80 = tmp72 - tmp79;
		304	tmp111 = K866025403 * (tmp72 + tmp79);
		305	tmp82 = (K342020143 * tmp78) - (K939692620 * tmp75);
		306	tmp83 = (K173648177 * tmp71) + (K984807753 * tmp68);
		307	tmp84 = K866025403 * (tmp82 - tmp83);
		308	tmp109 = tmp83 + tmp82;
		309	}
		310	c_re(inout[2 * iostride]) = tmp65 + tmp80;
		311	tmp81 = tmp65 - (K500000000 * tmp80);
		312	c_re(inout[8 * iostride]) = tmp81 - tmp84;
		313	c_re(inout[5 * iostride]) = tmp81 + tmp84;
		314	c_im(inout[2 * iostride]) = tmp109 + tmp110;
		315	tmp112 = tmp110 - (K500000000 * tmp109);
		316	c_im(inout[5 * iostride]) = tmp111 + tmp112;
		317	c_im(inout[8 * iostride]) = tmp112 - tmp111;
		318	}
		319	{
		320	fftw_real tmp85;
		321	fftw_real tmp106;
		322	fftw_real tmp92;
		323	fftw_real tmp107;
		324	fftw_real tmp96;
		325	fftw_real tmp103;
		326	fftw_real tmp93;
		327	fftw_real tmp108;
		328	ASSERT_ALIGNED_DOUBLE();
		329	tmp85 = tmp61 + tmp64;
		330	tmp106 = tmp104 + tmp105;
		331	{
		332	fftw_real tmp88;
		333	fftw_real tmp91;
		334	fftw_real tmp94;
		335	fftw_real tmp95;
		336	ASSERT_ALIGNED_DOUBLE();
		337	tmp88 = (K766044443 * tmp86) - (K642787609 * tmp87);
		338	tmp91 = (K173648177 * tmp89) - (K984807753 * tmp90);
		339	tmp92 = tmp88 + tmp91;
		340	tmp107 = K866025403 * (tmp88 - tmp91);
		341	tmp94 = (K173648177 * tmp90) + (K984807753 * tmp89);
		342	tmp95 = (K766044443 * tmp87) + (K642787609 * tmp86);
		343	tmp96 = K866025403 * (tmp94 - tmp95);
		344	tmp103 = tmp95 + tmp94;
		345	}
		346	c_re(inout[iostride]) = tmp85 + tmp92;
		347	tmp93 = tmp85 - (K500000000 * tmp92);
		348	c_re(inout[7 * iostride]) = tmp93 - tmp96;
		349	c_re(inout[4 * iostride]) = tmp93 + tmp96;
		350	c_im(inout[iostride]) = tmp103 + tmp106;
		351	tmp108 = tmp106 - (K500000000 * tmp103);
		352	c_im(inout[4 * iostride]) = tmp107 + tmp108;
		353	c_im(inout[7 * iostride]) = tmp108 - tmp107;
		354	}
		355	}
		356	}
		357
		358	static const int twiddle_order[] =
		359	{1, 2, 3, 4, 5, 6, 7, 8};
		360	fftw_codelet_desc fftwi_twiddle_9_desc =
		361	{
		362	"fftwi_twiddle_9",
		363	(void (*)()) fftwi_twiddle_9,
		364	9,
		365	FFTW_BACKWARD,
		366	FFTW_TWIDDLE,
		367	209,
		368	8,
		369	twiddle_order,
		370	};

Subversion Repositories shark

(root)/shark/trunk/ports/fftw/fftw/ftwi_9.c - Rev 1688