WebSVN - shark - Blame - Rev 760 - /shark/trunk/ports/fftw/fftw/ftw

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:22 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 -twiddle 10 */
		27
		28	/*
		29	* This function contains 102 FP additions, 60 FP multiplications,
		30	* (or, 72 additions, 30 multiplications, 30 fused multiply/add),
		31	* 42 stack variables, and 40 memory accesses
		32	*/
		33	static const fftw_real K587785252 = FFTW_KONST(+0.587785252292473129168705954639072768597652438);
		34	static const fftw_real K951056516 = FFTW_KONST(+0.951056516295153572116439333379382143405698634);
		35	static const fftw_real K250000000 = FFTW_KONST(+0.250000000000000000000000000000000000000000000);
		36	static const fftw_real K559016994 = FFTW_KONST(+0.559016994374947424102293417182819058860154590);
		37
		38	/*
		39	* Generator Id's :
107	pj	40	* $Id: ftw_10.c,v 1.2 2003-03-24 11:14:52 pj Exp $
		41	* $Id: ftw_10.c,v 1.2 2003-03-24 11:14:52 pj Exp $
		42	* $Id: ftw_10.c,v 1.2 2003-03-24 11:14:52 pj Exp $
2	pj	43	*/
		44
		45	void fftw_twiddle_10(fftw_complex A, const fftw_complex W, int iostride, int m, int dist)
		46	{
		47	int i;
		48	fftw_complex *inout;
		49	inout = A;
		50	for (i = m; i > 0; i = i - 1, inout = inout + dist, W = W + 9) {
		51	fftw_real tmp7;
		52	fftw_real tmp55;
		53	fftw_real tmp100;
		54	fftw_real tmp112;
		55	fftw_real tmp41;
		56	fftw_real tmp52;
		57	fftw_real tmp53;
		58	fftw_real tmp59;
		59	fftw_real tmp60;
		60	fftw_real tmp61;
		61	fftw_real tmp75;
		62	fftw_real tmp78;
		63	fftw_real tmp110;
		64	fftw_real tmp86;
		65	fftw_real tmp87;
		66	fftw_real tmp96;
		67	fftw_real tmp18;
		68	fftw_real tmp29;
		69	fftw_real tmp30;
		70	fftw_real tmp56;
		71	fftw_real tmp57;
		72	fftw_real tmp58;
		73	fftw_real tmp68;
		74	fftw_real tmp71;
		75	fftw_real tmp109;
		76	fftw_real tmp89;
		77	fftw_real tmp90;
		78	fftw_real tmp95;
		79	ASSERT_ALIGNED_DOUBLE();
		80	{
		81	fftw_real tmp1;
		82	fftw_real tmp99;
		83	fftw_real tmp6;
		84	fftw_real tmp98;
		85	ASSERT_ALIGNED_DOUBLE();
		86	tmp1 = c_re(inout[0]);
		87	tmp99 = c_im(inout[0]);
		88	{
		89	fftw_real tmp3;
		90	fftw_real tmp5;
		91	fftw_real tmp2;
		92	fftw_real tmp4;
		93	ASSERT_ALIGNED_DOUBLE();
		94	tmp3 = c_re(inout[5 * iostride]);
		95	tmp5 = c_im(inout[5 * iostride]);
		96	tmp2 = c_re(W[4]);
		97	tmp4 = c_im(W[4]);
		98	tmp6 = (tmp2 * tmp3) - (tmp4 * tmp5);
		99	tmp98 = (tmp4 * tmp3) + (tmp2 * tmp5);
		100	}
		101	tmp7 = tmp1 - tmp6;
		102	tmp55 = tmp1 + tmp6;
		103	tmp100 = tmp98 + tmp99;
		104	tmp112 = tmp99 - tmp98;
		105	}
		106	{
		107	fftw_real tmp35;
		108	fftw_real tmp73;
		109	fftw_real tmp51;
		110	fftw_real tmp77;
		111	fftw_real tmp40;
		112	fftw_real tmp74;
		113	fftw_real tmp46;
		114	fftw_real tmp76;
		115	ASSERT_ALIGNED_DOUBLE();
		116	{
		117	fftw_real tmp32;
		118	fftw_real tmp34;
		119	fftw_real tmp31;
		120	fftw_real tmp33;
		121	ASSERT_ALIGNED_DOUBLE();
		122	tmp32 = c_re(inout[4 * iostride]);
		123	tmp34 = c_im(inout[4 * iostride]);
		124	tmp31 = c_re(W[3]);
		125	tmp33 = c_im(W[3]);
		126	tmp35 = (tmp31 * tmp32) - (tmp33 * tmp34);
		127	tmp73 = (tmp33 * tmp32) + (tmp31 * tmp34);
		128	}
		129	{
		130	fftw_real tmp48;
		131	fftw_real tmp50;
		132	fftw_real tmp47;
		133	fftw_real tmp49;
		134	ASSERT_ALIGNED_DOUBLE();
		135	tmp48 = c_re(inout[iostride]);
		136	tmp50 = c_im(inout[iostride]);
		137	tmp47 = c_re(W[0]);
		138	tmp49 = c_im(W[0]);
		139	tmp51 = (tmp47 * tmp48) - (tmp49 * tmp50);
		140	tmp77 = (tmp49 * tmp48) + (tmp47 * tmp50);
		141	}
		142	{
		143	fftw_real tmp37;
		144	fftw_real tmp39;
		145	fftw_real tmp36;
		146	fftw_real tmp38;
		147	ASSERT_ALIGNED_DOUBLE();
		148	tmp37 = c_re(inout[9 * iostride]);
		149	tmp39 = c_im(inout[9 * iostride]);
		150	tmp36 = c_re(W[8]);
		151	tmp38 = c_im(W[8]);
		152	tmp40 = (tmp36 * tmp37) - (tmp38 * tmp39);
		153	tmp74 = (tmp38 * tmp37) + (tmp36 * tmp39);
		154	}
		155	{
		156	fftw_real tmp43;
		157	fftw_real tmp45;
		158	fftw_real tmp42;
		159	fftw_real tmp44;
		160	ASSERT_ALIGNED_DOUBLE();
		161	tmp43 = c_re(inout[6 * iostride]);
		162	tmp45 = c_im(inout[6 * iostride]);
		163	tmp42 = c_re(W[5]);
		164	tmp44 = c_im(W[5]);
		165	tmp46 = (tmp42 * tmp43) - (tmp44 * tmp45);
		166	tmp76 = (tmp44 * tmp43) + (tmp42 * tmp45);
		167	}
		168	tmp41 = tmp35 - tmp40;
		169	tmp52 = tmp46 - tmp51;
		170	tmp53 = tmp41 + tmp52;
		171	tmp59 = tmp35 + tmp40;
		172	tmp60 = tmp46 + tmp51;
		173	tmp61 = tmp59 + tmp60;
		174	tmp75 = tmp73 - tmp74;
		175	tmp78 = tmp76 - tmp77;
		176	tmp110 = tmp75 + tmp78;
		177	tmp86 = tmp73 + tmp74;
		178	tmp87 = tmp76 + tmp77;
		179	tmp96 = tmp86 + tmp87;
		180	}
		181	{
		182	fftw_real tmp12;
		183	fftw_real tmp66;
		184	fftw_real tmp28;
		185	fftw_real tmp70;
		186	fftw_real tmp17;
		187	fftw_real tmp67;
		188	fftw_real tmp23;
		189	fftw_real tmp69;
		190	ASSERT_ALIGNED_DOUBLE();
		191	{
		192	fftw_real tmp9;
		193	fftw_real tmp11;
		194	fftw_real tmp8;
		195	fftw_real tmp10;
		196	ASSERT_ALIGNED_DOUBLE();
		197	tmp9 = c_re(inout[2 * iostride]);
		198	tmp11 = c_im(inout[2 * iostride]);
		199	tmp8 = c_re(W[1]);
		200	tmp10 = c_im(W[1]);
		201	tmp12 = (tmp8 * tmp9) - (tmp10 * tmp11);
		202	tmp66 = (tmp10 * tmp9) + (tmp8 * tmp11);
		203	}
		204	{
		205	fftw_real tmp25;
		206	fftw_real tmp27;
		207	fftw_real tmp24;
		208	fftw_real tmp26;
		209	ASSERT_ALIGNED_DOUBLE();
		210	tmp25 = c_re(inout[3 * iostride]);
		211	tmp27 = c_im(inout[3 * iostride]);
		212	tmp24 = c_re(W[2]);
		213	tmp26 = c_im(W[2]);
		214	tmp28 = (tmp24 * tmp25) - (tmp26 * tmp27);
		215	tmp70 = (tmp26 * tmp25) + (tmp24 * tmp27);
		216	}
		217	{
		218	fftw_real tmp14;
		219	fftw_real tmp16;
		220	fftw_real tmp13;
		221	fftw_real tmp15;
		222	ASSERT_ALIGNED_DOUBLE();
		223	tmp14 = c_re(inout[7 * iostride]);
		224	tmp16 = c_im(inout[7 * iostride]);
		225	tmp13 = c_re(W[6]);
		226	tmp15 = c_im(W[6]);
		227	tmp17 = (tmp13 * tmp14) - (tmp15 * tmp16);
		228	tmp67 = (tmp15 * tmp14) + (tmp13 * tmp16);
		229	}
		230	{
		231	fftw_real tmp20;
		232	fftw_real tmp22;
		233	fftw_real tmp19;
		234	fftw_real tmp21;
		235	ASSERT_ALIGNED_DOUBLE();
		236	tmp20 = c_re(inout[8 * iostride]);
		237	tmp22 = c_im(inout[8 * iostride]);
		238	tmp19 = c_re(W[7]);
		239	tmp21 = c_im(W[7]);
		240	tmp23 = (tmp19 * tmp20) - (tmp21 * tmp22);
		241	tmp69 = (tmp21 * tmp20) + (tmp19 * tmp22);
		242	}
		243	tmp18 = tmp12 - tmp17;
		244	tmp29 = tmp23 - tmp28;
		245	tmp30 = tmp18 + tmp29;
		246	tmp56 = tmp12 + tmp17;
		247	tmp57 = tmp23 + tmp28;
		248	tmp58 = tmp56 + tmp57;
		249	tmp68 = tmp66 - tmp67;
		250	tmp71 = tmp69 - tmp70;
		251	tmp109 = tmp68 + tmp71;
		252	tmp89 = tmp66 + tmp67;
		253	tmp90 = tmp69 + tmp70;
		254	tmp95 = tmp89 + tmp90;
		255	}
		256	{
		257	fftw_real tmp63;
		258	fftw_real tmp54;
		259	fftw_real tmp64;
		260	fftw_real tmp80;
		261	fftw_real tmp82;
		262	fftw_real tmp72;
		263	fftw_real tmp79;
		264	fftw_real tmp81;
		265	fftw_real tmp65;
		266	ASSERT_ALIGNED_DOUBLE();
		267	tmp63 = K559016994 * (tmp30 - tmp53);
		268	tmp54 = tmp30 + tmp53;
		269	tmp64 = tmp7 - (K250000000 * tmp54);
		270	tmp72 = tmp68 - tmp71;
		271	tmp79 = tmp75 - tmp78;
		272	tmp80 = (K951056516 * tmp72) + (K587785252 * tmp79);
		273	tmp82 = (K951056516 * tmp79) - (K587785252 * tmp72);
		274	c_re(inout[5 * iostride]) = tmp7 + tmp54;
		275	tmp81 = tmp64 - tmp63;
		276	c_re(inout[7 * iostride]) = tmp81 - tmp82;
		277	c_re(inout[3 * iostride]) = tmp81 + tmp82;
		278	tmp65 = tmp63 + tmp64;
		279	c_re(inout[9 * iostride]) = tmp65 - tmp80;
		280	c_re(inout[iostride]) = tmp65 + tmp80;
		281	}
		282	{
		283	fftw_real tmp111;
		284	fftw_real tmp113;
		285	fftw_real tmp114;
		286	fftw_real tmp118;
		287	fftw_real tmp120;
		288	fftw_real tmp116;
		289	fftw_real tmp117;
		290	fftw_real tmp119;
		291	fftw_real tmp115;
		292	ASSERT_ALIGNED_DOUBLE();
		293	tmp111 = K559016994 * (tmp109 - tmp110);
		294	tmp113 = tmp109 + tmp110;
		295	tmp114 = tmp112 - (K250000000 * tmp113);
		296	tmp116 = tmp18 - tmp29;
		297	tmp117 = tmp41 - tmp52;
		298	tmp118 = (K951056516 * tmp116) + (K587785252 * tmp117);
		299	tmp120 = (K951056516 * tmp117) - (K587785252 * tmp116);
		300	c_im(inout[5 * iostride]) = tmp113 + tmp112;
		301	tmp119 = tmp114 - tmp111;
		302	c_im(inout[3 * iostride]) = tmp119 - tmp120;
		303	c_im(inout[7 * iostride]) = tmp120 + tmp119;
		304	tmp115 = tmp111 + tmp114;
		305	c_im(inout[iostride]) = tmp115 - tmp118;
		306	c_im(inout[9 * iostride]) = tmp118 + tmp115;
		307	}
		308	{
		309	fftw_real tmp84;
		310	fftw_real tmp62;
		311	fftw_real tmp83;
		312	fftw_real tmp92;
		313	fftw_real tmp94;
		314	fftw_real tmp88;
		315	fftw_real tmp91;
		316	fftw_real tmp93;
		317	fftw_real tmp85;
		318	ASSERT_ALIGNED_DOUBLE();
		319	tmp84 = K559016994 * (tmp58 - tmp61);
		320	tmp62 = tmp58 + tmp61;
		321	tmp83 = tmp55 - (K250000000 * tmp62);
		322	tmp88 = tmp86 - tmp87;
		323	tmp91 = tmp89 - tmp90;
		324	tmp92 = (K951056516 * tmp88) - (K587785252 * tmp91);
		325	tmp94 = (K951056516 * tmp91) + (K587785252 * tmp88);
		326	c_re(inout[0]) = tmp55 + tmp62;
		327	tmp93 = tmp84 + tmp83;
		328	c_re(inout[4 * iostride]) = tmp93 - tmp94;
		329	c_re(inout[6 * iostride]) = tmp93 + tmp94;
		330	tmp85 = tmp83 - tmp84;
		331	c_re(inout[2 * iostride]) = tmp85 - tmp92;
		332	c_re(inout[8 * iostride]) = tmp85 + tmp92;
		333	}
		334	{
		335	fftw_real tmp105;
		336	fftw_real tmp97;
		337	fftw_real tmp104;
		338	fftw_real tmp103;
		339	fftw_real tmp107;
		340	fftw_real tmp101;
		341	fftw_real tmp102;
		342	fftw_real tmp108;
		343	fftw_real tmp106;
		344	ASSERT_ALIGNED_DOUBLE();
		345	tmp105 = K559016994 * (tmp95 - tmp96);
		346	tmp97 = tmp95 + tmp96;
		347	tmp104 = tmp100 - (K250000000 * tmp97);
		348	tmp101 = tmp59 - tmp60;
		349	tmp102 = tmp56 - tmp57;
		350	tmp103 = (K951056516 * tmp101) - (K587785252 * tmp102);
		351	tmp107 = (K951056516 * tmp102) + (K587785252 * tmp101);
		352	c_im(inout[0]) = tmp97 + tmp100;
		353	tmp108 = tmp105 + tmp104;
		354	c_im(inout[4 * iostride]) = tmp107 + tmp108;
		355	c_im(inout[6 * iostride]) = tmp108 - tmp107;
		356	tmp106 = tmp104 - tmp105;
		357	c_im(inout[2 * iostride]) = tmp103 + tmp106;
		358	c_im(inout[8 * iostride]) = tmp106 - tmp103;
		359	}
		360	}
		361	}
		362
		363	static const int twiddle_order[] =
		364	{1, 2, 3, 4, 5, 6, 7, 8, 9};
		365	fftw_codelet_desc fftw_twiddle_10_desc =
		366	{
		367	"fftw_twiddle_10",
		368	(void (*)()) fftw_twiddle_10,
		369	10,
		370	FFTW_FORWARD,
		371	FFTW_TWIDDLE,
		372	220,
		373	9,
		374	twiddle_order,
		375	};

Subversion Repositories shark

(root)/shark/trunk/ports/fftw/fftw/ftw_10.c - Rev 760