WebSVN - shark - Blame - Rev 3 - /shark/trunk/ports/fftw/fftw/fn

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:54:31 EDT 1999 */
		22
		23	#include <ports/fftw-int.h>
		24	#include <ports/fftw.h>
		25
		26	/* Generated by: ./genfft -magic-alignment-check -magic-twiddle-load-all -magic-variables 4 -magic-loopi -notwiddle 32 */
		27
		28	/*
		29	* This function contains 372 FP additions, 84 FP multiplications,
		30	* (or, 340 additions, 52 multiplications, 32 fused multiply/add),
		31	* 92 stack variables, and 128 memory accesses
		32	*/
		33	static const fftw_real K831469612 = FFTW_KONST(+0.831469612302545237078788377617905756738560812);
		34	static const fftw_real K555570233 = FFTW_KONST(+0.555570233019602224742830813948532874374937191);
		35	static const fftw_real K195090322 = FFTW_KONST(+0.195090322016128267848284868477022240927691618);
		36	static const fftw_real K980785280 = FFTW_KONST(+0.980785280403230449126182236134239036973933731);
		37	static const fftw_real K923879532 = FFTW_KONST(+0.923879532511286756128183189396788286822416626);
		38	static const fftw_real K382683432 = FFTW_KONST(+0.382683432365089771728459984030398866761344562);
		39	static const fftw_real K707106781 = FFTW_KONST(+0.707106781186547524400844362104849039284835938);
		40
		41	/*
		42	* Generator Id's :
		43	* $Id: fn_32.c,v 1.1.1.1 2002-03-29 14:12:56 pj Exp $
		44	* $Id: fn_32.c,v 1.1.1.1 2002-03-29 14:12:56 pj Exp $
		45	* $Id: fn_32.c,v 1.1.1.1 2002-03-29 14:12:56 pj Exp $
		46	*/
		47
		48	void fftw_no_twiddle_32(const fftw_complex input, fftw_complex output, int istride, int ostride)
		49	{
		50	fftw_real tmp7;
		51	fftw_real tmp275;
		52	fftw_real tmp70;
		53	fftw_real tmp309;
		54	fftw_real tmp97;
		55	fftw_real tmp215;
		56	fftw_real tmp179;
		57	fftw_real tmp241;
		58	fftw_real tmp14;
		59	fftw_real tmp310;
		60	fftw_real tmp77;
		61	fftw_real tmp276;
		62	fftw_real tmp182;
		63	fftw_real tmp216;
		64	fftw_real tmp104;
		65	fftw_real tmp242;
		66	fftw_real tmp153;
		67	fftw_real tmp233;
		68	fftw_real tmp53;
		69	fftw_real tmp60;
		70	fftw_real tmp351;
		71	fftw_real tmp306;
		72	fftw_real tmp330;
		73	fftw_real tmp352;
		74	fftw_real tmp353;
		75	fftw_real tmp354;
		76	fftw_real tmp170;
		77	fftw_real tmp236;
		78	fftw_real tmp301;
		79	fftw_real tmp329;
		80	fftw_real tmp164;
		81	fftw_real tmp237;
		82	fftw_real tmp173;
		83	fftw_real tmp234;
		84	fftw_real tmp22;
		85	fftw_real tmp280;
		86	fftw_real tmp313;
		87	fftw_real tmp85;
		88	fftw_real tmp112;
		89	fftw_real tmp185;
		90	fftw_real tmp220;
		91	fftw_real tmp245;
		92	fftw_real tmp29;
		93	fftw_real tmp283;
		94	fftw_real tmp312;
		95	fftw_real tmp92;
		96	fftw_real tmp119;
		97	fftw_real tmp184;
		98	fftw_real tmp223;
		99	fftw_real tmp244;
		100	fftw_real tmp126;
		101	fftw_real tmp229;
		102	fftw_real tmp38;
		103	fftw_real tmp45;
		104	fftw_real tmp346;
		105	fftw_real tmp295;
		106	fftw_real tmp327;
		107	fftw_real tmp347;
		108	fftw_real tmp348;
		109	fftw_real tmp349;
		110	fftw_real tmp143;
		111	fftw_real tmp226;
		112	fftw_real tmp290;
		113	fftw_real tmp326;
		114	fftw_real tmp137;
		115	fftw_real tmp227;
		116	fftw_real tmp146;
		117	fftw_real tmp230;
		118	ASSERT_ALIGNED_DOUBLE();
		119	{
		120	fftw_real tmp3;
		121	fftw_real tmp95;
		122	fftw_real tmp66;
		123	fftw_real tmp178;
		124	fftw_real tmp6;
		125	fftw_real tmp177;
		126	fftw_real tmp69;
		127	fftw_real tmp96;
		128	ASSERT_ALIGNED_DOUBLE();
		129	{
		130	fftw_real tmp1;
		131	fftw_real tmp2;
		132	fftw_real tmp64;
		133	fftw_real tmp65;
		134	ASSERT_ALIGNED_DOUBLE();
		135	tmp1 = c_re(input[0]);
		136	tmp2 = c_re(input[16 * istride]);
		137	tmp3 = tmp1 + tmp2;
		138	tmp95 = tmp1 - tmp2;
		139	tmp64 = c_im(input[0]);
		140	tmp65 = c_im(input[16 * istride]);
		141	tmp66 = tmp64 + tmp65;
		142	tmp178 = tmp64 - tmp65;
		143	}
		144	{
		145	fftw_real tmp4;
		146	fftw_real tmp5;
		147	fftw_real tmp67;
		148	fftw_real tmp68;
		149	ASSERT_ALIGNED_DOUBLE();
		150	tmp4 = c_re(input[8 * istride]);
		151	tmp5 = c_re(input[24 * istride]);
		152	tmp6 = tmp4 + tmp5;
		153	tmp177 = tmp4 - tmp5;
		154	tmp67 = c_im(input[8 * istride]);
		155	tmp68 = c_im(input[24 * istride]);
		156	tmp69 = tmp67 + tmp68;
		157	tmp96 = tmp67 - tmp68;
		158	}
		159	tmp7 = tmp3 + tmp6;
		160	tmp275 = tmp3 - tmp6;
		161	tmp70 = tmp66 + tmp69;
		162	tmp309 = tmp66 - tmp69;
		163	tmp97 = tmp95 - tmp96;
		164	tmp215 = tmp95 + tmp96;
		165	tmp179 = tmp177 + tmp178;
		166	tmp241 = tmp178 - tmp177;
		167	}
		168	{
		169	fftw_real tmp10;
		170	fftw_real tmp99;
		171	fftw_real tmp73;
		172	fftw_real tmp98;
		173	fftw_real tmp13;
		174	fftw_real tmp101;
		175	fftw_real tmp76;
		176	fftw_real tmp102;
		177	ASSERT_ALIGNED_DOUBLE();
		178	{
		179	fftw_real tmp8;
		180	fftw_real tmp9;
		181	fftw_real tmp71;
		182	fftw_real tmp72;
		183	ASSERT_ALIGNED_DOUBLE();
		184	tmp8 = c_re(input[4 * istride]);
		185	tmp9 = c_re(input[20 * istride]);
		186	tmp10 = tmp8 + tmp9;
		187	tmp99 = tmp8 - tmp9;
		188	tmp71 = c_im(input[4 * istride]);
		189	tmp72 = c_im(input[20 * istride]);
		190	tmp73 = tmp71 + tmp72;
		191	tmp98 = tmp71 - tmp72;
		192	}
		193	{
		194	fftw_real tmp11;
		195	fftw_real tmp12;
		196	fftw_real tmp74;
		197	fftw_real tmp75;
		198	ASSERT_ALIGNED_DOUBLE();
		199	tmp11 = c_re(input[28 * istride]);
		200	tmp12 = c_re(input[12 * istride]);
		201	tmp13 = tmp11 + tmp12;
		202	tmp101 = tmp11 - tmp12;
		203	tmp74 = c_im(input[28 * istride]);
		204	tmp75 = c_im(input[12 * istride]);
		205	tmp76 = tmp74 + tmp75;
		206	tmp102 = tmp74 - tmp75;
		207	}
		208	tmp14 = tmp10 + tmp13;
		209	tmp310 = tmp13 - tmp10;
		210	tmp77 = tmp73 + tmp76;
		211	tmp276 = tmp73 - tmp76;
		212	{
		213	fftw_real tmp180;
		214	fftw_real tmp181;
		215	fftw_real tmp100;
		216	fftw_real tmp103;
		217	ASSERT_ALIGNED_DOUBLE();
		218	tmp180 = tmp101 - tmp102;
		219	tmp181 = tmp99 + tmp98;
		220	tmp182 = K707106781 * (tmp180 - tmp181);
		221	tmp216 = K707106781 * (tmp181 + tmp180);
		222	tmp100 = tmp98 - tmp99;
		223	tmp103 = tmp101 + tmp102;
		224	tmp104 = K707106781 * (tmp100 - tmp103);
		225	tmp242 = K707106781 * (tmp100 + tmp103);
		226	}
		227	}
		228	{
		229	fftw_real tmp49;
		230	fftw_real tmp149;
		231	fftw_real tmp169;
		232	fftw_real tmp302;
		233	fftw_real tmp52;
		234	fftw_real tmp166;
		235	fftw_real tmp152;
		236	fftw_real tmp303;
		237	fftw_real tmp56;
		238	fftw_real tmp157;
		239	fftw_real tmp156;
		240	fftw_real tmp298;
		241	fftw_real tmp59;
		242	fftw_real tmp159;
		243	fftw_real tmp162;
		244	fftw_real tmp299;
		245	ASSERT_ALIGNED_DOUBLE();
		246	{
		247	fftw_real tmp47;
		248	fftw_real tmp48;
		249	fftw_real tmp167;
		250	fftw_real tmp168;
		251	ASSERT_ALIGNED_DOUBLE();
		252	tmp47 = c_re(input[31 * istride]);
		253	tmp48 = c_re(input[15 * istride]);
		254	tmp49 = tmp47 + tmp48;
		255	tmp149 = tmp47 - tmp48;
		256	tmp167 = c_im(input[31 * istride]);
		257	tmp168 = c_im(input[15 * istride]);
		258	tmp169 = tmp167 - tmp168;
		259	tmp302 = tmp167 + tmp168;
		260	}
		261	{
		262	fftw_real tmp50;
		263	fftw_real tmp51;
		264	fftw_real tmp150;
		265	fftw_real tmp151;
		266	ASSERT_ALIGNED_DOUBLE();
		267	tmp50 = c_re(input[7 * istride]);
		268	tmp51 = c_re(input[23 * istride]);
		269	tmp52 = tmp50 + tmp51;
		270	tmp166 = tmp50 - tmp51;
		271	tmp150 = c_im(input[7 * istride]);
		272	tmp151 = c_im(input[23 * istride]);
		273	tmp152 = tmp150 - tmp151;
		274	tmp303 = tmp150 + tmp151;
		275	}
		276	{
		277	fftw_real tmp54;
		278	fftw_real tmp55;
		279	fftw_real tmp154;
		280	fftw_real tmp155;
		281	ASSERT_ALIGNED_DOUBLE();
		282	tmp54 = c_re(input[3 * istride]);
		283	tmp55 = c_re(input[19 * istride]);
		284	tmp56 = tmp54 + tmp55;
		285	tmp157 = tmp54 - tmp55;
		286	tmp154 = c_im(input[3 * istride]);
		287	tmp155 = c_im(input[19 * istride]);
		288	tmp156 = tmp154 - tmp155;
		289	tmp298 = tmp154 + tmp155;
		290	}
		291	{
		292	fftw_real tmp57;
		293	fftw_real tmp58;
		294	fftw_real tmp160;
		295	fftw_real tmp161;
		296	ASSERT_ALIGNED_DOUBLE();
		297	tmp57 = c_re(input[27 * istride]);
		298	tmp58 = c_re(input[11 * istride]);
		299	tmp59 = tmp57 + tmp58;
		300	tmp159 = tmp57 - tmp58;
		301	tmp160 = c_im(input[27 * istride]);
		302	tmp161 = c_im(input[11 * istride]);
		303	tmp162 = tmp160 - tmp161;
		304	tmp299 = tmp160 + tmp161;
		305	}
		306	{
		307	fftw_real tmp304;
		308	fftw_real tmp305;
		309	fftw_real tmp297;
		310	fftw_real tmp300;
		311	ASSERT_ALIGNED_DOUBLE();
		312	tmp153 = tmp149 - tmp152;
		313	tmp233 = tmp149 + tmp152;
		314	tmp53 = tmp49 + tmp52;
		315	tmp60 = tmp56 + tmp59;
		316	tmp351 = tmp53 - tmp60;
		317	tmp304 = tmp302 - tmp303;
		318	tmp305 = tmp59 - tmp56;
		319	tmp306 = tmp304 - tmp305;
		320	tmp330 = tmp305 + tmp304;
		321	tmp352 = tmp302 + tmp303;
		322	tmp353 = tmp298 + tmp299;
		323	tmp354 = tmp352 - tmp353;
		324	tmp170 = tmp166 + tmp169;
		325	tmp236 = tmp169 - tmp166;
		326	tmp297 = tmp49 - tmp52;
		327	tmp300 = tmp298 - tmp299;
		328	tmp301 = tmp297 - tmp300;
		329	tmp329 = tmp297 + tmp300;
		330	{
		331	fftw_real tmp158;
		332	fftw_real tmp163;
		333	fftw_real tmp171;
		334	fftw_real tmp172;
		335	ASSERT_ALIGNED_DOUBLE();
		336	tmp158 = tmp156 - tmp157;
		337	tmp163 = tmp159 + tmp162;
		338	tmp164 = K707106781 * (tmp158 - tmp163);
		339	tmp237 = K707106781 * (tmp158 + tmp163);
		340	tmp171 = tmp159 - tmp162;
		341	tmp172 = tmp157 + tmp156;
		342	tmp173 = K707106781 * (tmp171 - tmp172);
		343	tmp234 = K707106781 * (tmp172 + tmp171);
		344	}
		345	}
		346	}
		347	{
		348	fftw_real tmp18;
		349	fftw_real tmp109;
		350	fftw_real tmp81;
		351	fftw_real tmp107;
		352	fftw_real tmp21;
		353	fftw_real tmp106;
		354	fftw_real tmp84;
		355	fftw_real tmp110;
		356	fftw_real tmp278;
		357	fftw_real tmp279;
		358	ASSERT_ALIGNED_DOUBLE();
		359	{
		360	fftw_real tmp16;
		361	fftw_real tmp17;
		362	fftw_real tmp79;
		363	fftw_real tmp80;
		364	ASSERT_ALIGNED_DOUBLE();
		365	tmp16 = c_re(input[2 * istride]);
		366	tmp17 = c_re(input[18 * istride]);
		367	tmp18 = tmp16 + tmp17;
		368	tmp109 = tmp16 - tmp17;
		369	tmp79 = c_im(input[2 * istride]);
		370	tmp80 = c_im(input[18 * istride]);
		371	tmp81 = tmp79 + tmp80;
		372	tmp107 = tmp79 - tmp80;
		373	}
		374	{
		375	fftw_real tmp19;
		376	fftw_real tmp20;
		377	fftw_real tmp82;
		378	fftw_real tmp83;
		379	ASSERT_ALIGNED_DOUBLE();
		380	tmp19 = c_re(input[10 * istride]);
		381	tmp20 = c_re(input[26 * istride]);
		382	tmp21 = tmp19 + tmp20;
		383	tmp106 = tmp19 - tmp20;
		384	tmp82 = c_im(input[10 * istride]);
		385	tmp83 = c_im(input[26 * istride]);
		386	tmp84 = tmp82 + tmp83;
		387	tmp110 = tmp82 - tmp83;
		388	}
		389	tmp22 = tmp18 + tmp21;
		390	tmp278 = tmp81 - tmp84;
		391	tmp279 = tmp18 - tmp21;
		392	tmp280 = tmp278 - tmp279;
		393	tmp313 = tmp279 + tmp278;
		394	tmp85 = tmp81 + tmp84;
		395	{
		396	fftw_real tmp108;
		397	fftw_real tmp111;
		398	fftw_real tmp218;
		399	fftw_real tmp219;
		400	ASSERT_ALIGNED_DOUBLE();
		401	tmp108 = tmp106 + tmp107;
		402	tmp111 = tmp109 - tmp110;
		403	tmp112 = (K382683432 * tmp108) - (K923879532 * tmp111);
		404	tmp185 = (K923879532 * tmp108) + (K382683432 * tmp111);
		405	tmp218 = tmp107 - tmp106;
		406	tmp219 = tmp109 + tmp110;
		407	tmp220 = (K923879532 * tmp218) - (K382683432 * tmp219);
		408	tmp245 = (K382683432 * tmp218) + (K923879532 * tmp219);
		409	}
		410	}
		411	{
		412	fftw_real tmp25;
		413	fftw_real tmp116;
		414	fftw_real tmp88;
		415	fftw_real tmp114;
		416	fftw_real tmp28;
		417	fftw_real tmp113;
		418	fftw_real tmp91;
		419	fftw_real tmp117;
		420	fftw_real tmp281;
		421	fftw_real tmp282;
		422	ASSERT_ALIGNED_DOUBLE();
		423	{
		424	fftw_real tmp23;
		425	fftw_real tmp24;
		426	fftw_real tmp86;
		427	fftw_real tmp87;
		428	ASSERT_ALIGNED_DOUBLE();
		429	tmp23 = c_re(input[30 * istride]);
		430	tmp24 = c_re(input[14 * istride]);
		431	tmp25 = tmp23 + tmp24;
		432	tmp116 = tmp23 - tmp24;
		433	tmp86 = c_im(input[30 * istride]);
		434	tmp87 = c_im(input[14 * istride]);
		435	tmp88 = tmp86 + tmp87;
		436	tmp114 = tmp86 - tmp87;
		437	}
		438	{
		439	fftw_real tmp26;
		440	fftw_real tmp27;
		441	fftw_real tmp89;
		442	fftw_real tmp90;
		443	ASSERT_ALIGNED_DOUBLE();
		444	tmp26 = c_re(input[6 * istride]);
		445	tmp27 = c_re(input[22 * istride]);
		446	tmp28 = tmp26 + tmp27;
		447	tmp113 = tmp26 - tmp27;
		448	tmp89 = c_im(input[6 * istride]);
		449	tmp90 = c_im(input[22 * istride]);
		450	tmp91 = tmp89 + tmp90;
		451	tmp117 = tmp89 - tmp90;
		452	}
		453	tmp29 = tmp25 + tmp28;
		454	tmp281 = tmp25 - tmp28;
		455	tmp282 = tmp88 - tmp91;
		456	tmp283 = tmp281 + tmp282;
		457	tmp312 = tmp281 - tmp282;
		458	tmp92 = tmp88 + tmp91;
		459	{
		460	fftw_real tmp115;
		461	fftw_real tmp118;
		462	fftw_real tmp221;
		463	fftw_real tmp222;
		464	ASSERT_ALIGNED_DOUBLE();
		465	tmp115 = tmp113 + tmp114;
		466	tmp118 = tmp116 - tmp117;
		467	tmp119 = (K382683432 * tmp115) + (K923879532 * tmp118);
		468	tmp184 = (K382683432 * tmp118) - (K923879532 * tmp115);
		469	tmp221 = tmp114 - tmp113;
		470	tmp222 = tmp116 + tmp117;
		471	tmp223 = (K923879532 * tmp221) + (K382683432 * tmp222);
		472	tmp244 = (K923879532 * tmp222) - (K382683432 * tmp221);
		473	}
		474	}
		475	{
		476	fftw_real tmp34;
		477	fftw_real tmp139;
		478	fftw_real tmp125;
		479	fftw_real tmp286;
		480	fftw_real tmp37;
		481	fftw_real tmp122;
		482	fftw_real tmp142;
		483	fftw_real tmp287;
		484	fftw_real tmp41;
		485	fftw_real tmp132;
		486	fftw_real tmp135;
		487	fftw_real tmp292;
		488	fftw_real tmp44;
		489	fftw_real tmp127;
		490	fftw_real tmp130;
		491	fftw_real tmp293;
		492	ASSERT_ALIGNED_DOUBLE();
		493	{
		494	fftw_real tmp32;
		495	fftw_real tmp33;
		496	fftw_real tmp123;
		497	fftw_real tmp124;
		498	ASSERT_ALIGNED_DOUBLE();
		499	tmp32 = c_re(input[istride]);
		500	tmp33 = c_re(input[17 * istride]);
		501	tmp34 = tmp32 + tmp33;
		502	tmp139 = tmp32 - tmp33;
		503	tmp123 = c_im(input[istride]);
		504	tmp124 = c_im(input[17 * istride]);
		505	tmp125 = tmp123 - tmp124;
		506	tmp286 = tmp123 + tmp124;
		507	}
		508	{
		509	fftw_real tmp35;
		510	fftw_real tmp36;
		511	fftw_real tmp140;
		512	fftw_real tmp141;
		513	ASSERT_ALIGNED_DOUBLE();
		514	tmp35 = c_re(input[9 * istride]);
		515	tmp36 = c_re(input[25 * istride]);
		516	tmp37 = tmp35 + tmp36;
		517	tmp122 = tmp35 - tmp36;
		518	tmp140 = c_im(input[9 * istride]);
		519	tmp141 = c_im(input[25 * istride]);
		520	tmp142 = tmp140 - tmp141;
		521	tmp287 = tmp140 + tmp141;
		522	}
		523	{
		524	fftw_real tmp39;
		525	fftw_real tmp40;
		526	fftw_real tmp133;
		527	fftw_real tmp134;
		528	ASSERT_ALIGNED_DOUBLE();
		529	tmp39 = c_re(input[5 * istride]);
		530	tmp40 = c_re(input[21 * istride]);
		531	tmp41 = tmp39 + tmp40;
		532	tmp132 = tmp39 - tmp40;
		533	tmp133 = c_im(input[5 * istride]);
		534	tmp134 = c_im(input[21 * istride]);
		535	tmp135 = tmp133 - tmp134;
		536	tmp292 = tmp133 + tmp134;
		537	}
		538	{
		539	fftw_real tmp42;
		540	fftw_real tmp43;
		541	fftw_real tmp128;
		542	fftw_real tmp129;
		543	ASSERT_ALIGNED_DOUBLE();
		544	tmp42 = c_re(input[29 * istride]);
		545	tmp43 = c_re(input[13 * istride]);
		546	tmp44 = tmp42 + tmp43;
		547	tmp127 = tmp42 - tmp43;
		548	tmp128 = c_im(input[29 * istride]);
		549	tmp129 = c_im(input[13 * istride]);
		550	tmp130 = tmp128 - tmp129;
		551	tmp293 = tmp128 + tmp129;
		552	}
		553	{
		554	fftw_real tmp291;
		555	fftw_real tmp294;
		556	fftw_real tmp288;
		557	fftw_real tmp289;
		558	ASSERT_ALIGNED_DOUBLE();
		559	tmp126 = tmp122 + tmp125;
		560	tmp229 = tmp125 - tmp122;
		561	tmp38 = tmp34 + tmp37;
		562	tmp45 = tmp41 + tmp44;
		563	tmp346 = tmp38 - tmp45;
		564	tmp291 = tmp34 - tmp37;
		565	tmp294 = tmp292 - tmp293;
		566	tmp295 = tmp291 - tmp294;
		567	tmp327 = tmp291 + tmp294;
		568	tmp347 = tmp286 + tmp287;
		569	tmp348 = tmp292 + tmp293;
		570	tmp349 = tmp347 - tmp348;
		571	tmp143 = tmp139 - tmp142;
		572	tmp226 = tmp139 + tmp142;
		573	tmp288 = tmp286 - tmp287;
		574	tmp289 = tmp44 - tmp41;
		575	tmp290 = tmp288 - tmp289;
		576	tmp326 = tmp289 + tmp288;
		577	{
		578	fftw_real tmp131;
		579	fftw_real tmp136;
		580	fftw_real tmp144;
		581	fftw_real tmp145;
		582	ASSERT_ALIGNED_DOUBLE();
		583	tmp131 = tmp127 - tmp130;
		584	tmp136 = tmp132 + tmp135;
		585	tmp137 = K707106781 * (tmp131 - tmp136);
		586	tmp227 = K707106781 * (tmp136 + tmp131);
		587	tmp144 = tmp135 - tmp132;
		588	tmp145 = tmp127 + tmp130;
		589	tmp146 = K707106781 * (tmp144 - tmp145);
		590	tmp230 = K707106781 * (tmp144 + tmp145);
		591	}
		592	}
		593	}
		594	{
		595	fftw_real tmp285;
		596	fftw_real tmp317;
		597	fftw_real tmp320;
		598	fftw_real tmp322;
		599	fftw_real tmp308;
		600	fftw_real tmp316;
		601	fftw_real tmp315;
		602	fftw_real tmp321;
		603	ASSERT_ALIGNED_DOUBLE();
		604	{
		605	fftw_real tmp277;
		606	fftw_real tmp284;
		607	fftw_real tmp318;
		608	fftw_real tmp319;
		609	ASSERT_ALIGNED_DOUBLE();
		610	tmp277 = tmp275 - tmp276;
		611	tmp284 = K707106781 * (tmp280 - tmp283);
		612	tmp285 = tmp277 + tmp284;
		613	tmp317 = tmp277 - tmp284;
		614	tmp318 = (K382683432 * tmp290) - (K923879532 * tmp295);
		615	tmp319 = (K382683432 * tmp306) + (K923879532 * tmp301);
		616	tmp320 = tmp318 - tmp319;
		617	tmp322 = tmp318 + tmp319;
		618	}
		619	{
		620	fftw_real tmp296;
		621	fftw_real tmp307;
		622	fftw_real tmp311;
		623	fftw_real tmp314;
		624	ASSERT_ALIGNED_DOUBLE();
		625	tmp296 = (K923879532 * tmp290) + (K382683432 * tmp295);
		626	tmp307 = (K382683432 * tmp301) - (K923879532 * tmp306);
		627	tmp308 = tmp296 + tmp307;
		628	tmp316 = tmp307 - tmp296;
		629	tmp311 = tmp309 - tmp310;
		630	tmp314 = K707106781 * (tmp312 - tmp313);
		631	tmp315 = tmp311 - tmp314;
		632	tmp321 = tmp311 + tmp314;
		633	}
		634	c_re(output[22 * ostride]) = tmp285 - tmp308;
		635	c_re(output[6 * ostride]) = tmp285 + tmp308;
		636	c_im(output[30 * ostride]) = tmp315 - tmp316;
		637	c_im(output[14 * ostride]) = tmp315 + tmp316;
		638	c_re(output[30 * ostride]) = tmp317 - tmp320;
		639	c_re(output[14 * ostride]) = tmp317 + tmp320;
		640	c_im(output[22 * ostride]) = tmp321 - tmp322;
		641	c_im(output[6 * ostride]) = tmp321 + tmp322;
		642	}
		643	{
		644	fftw_real tmp325;
		645	fftw_real tmp337;
		646	fftw_real tmp340;
		647	fftw_real tmp342;
		648	fftw_real tmp332;
		649	fftw_real tmp336;
		650	fftw_real tmp335;
		651	fftw_real tmp341;
		652	ASSERT_ALIGNED_DOUBLE();
		653	{
		654	fftw_real tmp323;
		655	fftw_real tmp324;
		656	fftw_real tmp338;
		657	fftw_real tmp339;
		658	ASSERT_ALIGNED_DOUBLE();
		659	tmp323 = tmp275 + tmp276;
		660	tmp324 = K707106781 * (tmp313 + tmp312);
		661	tmp325 = tmp323 + tmp324;
		662	tmp337 = tmp323 - tmp324;
		663	tmp338 = (K923879532 * tmp326) - (K382683432 * tmp327);
		664	tmp339 = (K923879532 * tmp330) + (K382683432 * tmp329);
		665	tmp340 = tmp338 - tmp339;
		666	tmp342 = tmp338 + tmp339;
		667	}
		668	{
		669	fftw_real tmp328;
		670	fftw_real tmp331;
		671	fftw_real tmp333;
		672	fftw_real tmp334;
		673	ASSERT_ALIGNED_DOUBLE();
		674	tmp328 = (K382683432 * tmp326) + (K923879532 * tmp327);
		675	tmp331 = (K923879532 * tmp329) - (K382683432 * tmp330);
		676	tmp332 = tmp328 + tmp331;
		677	tmp336 = tmp331 - tmp328;
		678	tmp333 = tmp310 + tmp309;
		679	tmp334 = K707106781 * (tmp280 + tmp283);
		680	tmp335 = tmp333 - tmp334;
		681	tmp341 = tmp333 + tmp334;
		682	}
		683	c_re(output[18 * ostride]) = tmp325 - tmp332;
		684	c_re(output[2 * ostride]) = tmp325 + tmp332;
		685	c_im(output[26 * ostride]) = tmp335 - tmp336;
		686	c_im(output[10 * ostride]) = tmp335 + tmp336;
		687	c_re(output[26 * ostride]) = tmp337 - tmp340;
		688	c_re(output[10 * ostride]) = tmp337 + tmp340;
		689	c_im(output[18 * ostride]) = tmp341 - tmp342;
		690	c_im(output[2 * ostride]) = tmp341 + tmp342;
		691	}
		692	{
		693	fftw_real tmp345;
		694	fftw_real tmp361;
		695	fftw_real tmp364;
		696	fftw_real tmp366;
		697	fftw_real tmp356;
		698	fftw_real tmp360;
		699	fftw_real tmp359;
		700	fftw_real tmp365;
		701	ASSERT_ALIGNED_DOUBLE();
		702	{
		703	fftw_real tmp343;
		704	fftw_real tmp344;
		705	fftw_real tmp362;
		706	fftw_real tmp363;
		707	ASSERT_ALIGNED_DOUBLE();
		708	tmp343 = tmp7 - tmp14;
		709	tmp344 = tmp85 - tmp92;
		710	tmp345 = tmp343 + tmp344;
		711	tmp361 = tmp343 - tmp344;
		712	tmp362 = tmp349 - tmp346;
		713	tmp363 = tmp351 + tmp354;
		714	tmp364 = K707106781 * (tmp362 - tmp363);
		715	tmp366 = K707106781 * (tmp362 + tmp363);
		716	}
		717	{
		718	fftw_real tmp350;
		719	fftw_real tmp355;
		720	fftw_real tmp357;
		721	fftw_real tmp358;
		722	ASSERT_ALIGNED_DOUBLE();
		723	tmp350 = tmp346 + tmp349;
		724	tmp355 = tmp351 - tmp354;
		725	tmp356 = K707106781 * (tmp350 + tmp355);
		726	tmp360 = K707106781 * (tmp355 - tmp350);
		727	tmp357 = tmp70 - tmp77;
		728	tmp358 = tmp29 - tmp22;
		729	tmp359 = tmp357 - tmp358;
		730	tmp365 = tmp358 + tmp357;
		731	}
		732	c_re(output[20 * ostride]) = tmp345 - tmp356;
		733	c_re(output[4 * ostride]) = tmp345 + tmp356;
		734	c_im(output[28 * ostride]) = tmp359 - tmp360;
		735	c_im(output[12 * ostride]) = tmp359 + tmp360;
		736	c_re(output[28 * ostride]) = tmp361 - tmp364;
		737	c_re(output[12 * ostride]) = tmp361 + tmp364;
		738	c_im(output[20 * ostride]) = tmp365 - tmp366;
		739	c_im(output[4 * ostride]) = tmp365 + tmp366;
		740	}
		741	{
		742	fftw_real tmp31;
		743	fftw_real tmp367;
		744	fftw_real tmp370;
		745	fftw_real tmp372;
		746	fftw_real tmp62;
		747	fftw_real tmp63;
		748	fftw_real tmp94;
		749	fftw_real tmp371;
		750	ASSERT_ALIGNED_DOUBLE();
		751	{
		752	fftw_real tmp15;
		753	fftw_real tmp30;
		754	fftw_real tmp368;
		755	fftw_real tmp369;
		756	ASSERT_ALIGNED_DOUBLE();
		757	tmp15 = tmp7 + tmp14;
		758	tmp30 = tmp22 + tmp29;
		759	tmp31 = tmp15 + tmp30;
		760	tmp367 = tmp15 - tmp30;
		761	tmp368 = tmp347 + tmp348;
		762	tmp369 = tmp352 + tmp353;
		763	tmp370 = tmp368 - tmp369;
		764	tmp372 = tmp368 + tmp369;
		765	}
		766	{
		767	fftw_real tmp46;
		768	fftw_real tmp61;
		769	fftw_real tmp78;
		770	fftw_real tmp93;
		771	ASSERT_ALIGNED_DOUBLE();
		772	tmp46 = tmp38 + tmp45;
		773	tmp61 = tmp53 + tmp60;
		774	tmp62 = tmp46 + tmp61;
		775	tmp63 = tmp61 - tmp46;
		776	tmp78 = tmp70 + tmp77;
		777	tmp93 = tmp85 + tmp92;
		778	tmp94 = tmp78 - tmp93;
		779	tmp371 = tmp78 + tmp93;
		780	}
		781	c_re(output[16 * ostride]) = tmp31 - tmp62;
		782	c_re(output[0]) = tmp31 + tmp62;
		783	c_im(output[8 * ostride]) = tmp63 + tmp94;
		784	c_im(output[24 * ostride]) = tmp94 - tmp63;
		785	c_re(output[24 * ostride]) = tmp367 - tmp370;
		786	c_re(output[8 * ostride]) = tmp367 + tmp370;
		787	c_im(output[16 * ostride]) = tmp371 - tmp372;
		788	c_im(output[0]) = tmp371 + tmp372;
		789	}
		790	{
		791	fftw_real tmp121;
		792	fftw_real tmp189;
		793	fftw_real tmp187;
		794	fftw_real tmp193;
		795	fftw_real tmp148;
		796	fftw_real tmp190;
		797	fftw_real tmp175;
		798	fftw_real tmp191;
		799	ASSERT_ALIGNED_DOUBLE();
		800	{
		801	fftw_real tmp105;
		802	fftw_real tmp120;
		803	fftw_real tmp183;
		804	fftw_real tmp186;
		805	ASSERT_ALIGNED_DOUBLE();
		806	tmp105 = tmp97 - tmp104;
		807	tmp120 = tmp112 - tmp119;
		808	tmp121 = tmp105 + tmp120;
		809	tmp189 = tmp105 - tmp120;
		810	tmp183 = tmp179 - tmp182;
		811	tmp186 = tmp184 - tmp185;
		812	tmp187 = tmp183 - tmp186;
		813	tmp193 = tmp183 + tmp186;
		814	}
		815	{
		816	fftw_real tmp138;
		817	fftw_real tmp147;
		818	fftw_real tmp165;
		819	fftw_real tmp174;
		820	ASSERT_ALIGNED_DOUBLE();
		821	tmp138 = tmp126 - tmp137;
		822	tmp147 = tmp143 - tmp146;
		823	tmp148 = (K980785280 * tmp138) + (K195090322 * tmp147);
		824	tmp190 = (K195090322 * tmp138) - (K980785280 * tmp147);
		825	tmp165 = tmp153 - tmp164;
		826	tmp174 = tmp170 - tmp173;
		827	tmp175 = (K195090322 * tmp165) - (K980785280 * tmp174);
		828	tmp191 = (K195090322 * tmp174) + (K980785280 * tmp165);
		829	}
		830	{
		831	fftw_real tmp176;
		832	fftw_real tmp188;
		833	fftw_real tmp192;
		834	fftw_real tmp194;
		835	ASSERT_ALIGNED_DOUBLE();
		836	tmp176 = tmp148 + tmp175;
		837	c_re(output[23 * ostride]) = tmp121 - tmp176;
		838	c_re(output[7 * ostride]) = tmp121 + tmp176;
		839	tmp188 = tmp175 - tmp148;
		840	c_im(output[31 * ostride]) = tmp187 - tmp188;
		841	c_im(output[15 * ostride]) = tmp187 + tmp188;
		842	tmp192 = tmp190 - tmp191;
		843	c_re(output[31 * ostride]) = tmp189 - tmp192;
		844	c_re(output[15 * ostride]) = tmp189 + tmp192;
		845	tmp194 = tmp190 + tmp191;
		846	c_im(output[23 * ostride]) = tmp193 - tmp194;
		847	c_im(output[7 * ostride]) = tmp193 + tmp194;
		848	}
		849	}
		850	{
		851	fftw_real tmp197;
		852	fftw_real tmp209;
		853	fftw_real tmp207;
		854	fftw_real tmp213;
		855	fftw_real tmp200;
		856	fftw_real tmp210;
		857	fftw_real tmp203;
		858	fftw_real tmp211;
		859	ASSERT_ALIGNED_DOUBLE();
		860	{
		861	fftw_real tmp195;
		862	fftw_real tmp196;
		863	fftw_real tmp205;
		864	fftw_real tmp206;
		865	ASSERT_ALIGNED_DOUBLE();
		866	tmp195 = tmp97 + tmp104;
		867	tmp196 = tmp185 + tmp184;
		868	tmp197 = tmp195 + tmp196;
		869	tmp209 = tmp195 - tmp196;
		870	tmp205 = tmp179 + tmp182;
		871	tmp206 = tmp112 + tmp119;
		872	tmp207 = tmp205 - tmp206;
		873	tmp213 = tmp205 + tmp206;
		874	}
		875	{
		876	fftw_real tmp198;
		877	fftw_real tmp199;
		878	fftw_real tmp201;
		879	fftw_real tmp202;
		880	ASSERT_ALIGNED_DOUBLE();
		881	tmp198 = tmp126 + tmp137;
		882	tmp199 = tmp143 + tmp146;
		883	tmp200 = (K555570233 * tmp198) + (K831469612 * tmp199);
		884	tmp210 = (K831469612 * tmp198) - (K555570233 * tmp199);
		885	tmp201 = tmp153 + tmp164;
		886	tmp202 = tmp170 + tmp173;
		887	tmp203 = (K831469612 * tmp201) - (K555570233 * tmp202);
		888	tmp211 = (K831469612 * tmp202) + (K555570233 * tmp201);
		889	}
		890	{
		891	fftw_real tmp204;
		892	fftw_real tmp208;
		893	fftw_real tmp212;
		894	fftw_real tmp214;
		895	ASSERT_ALIGNED_DOUBLE();
		896	tmp204 = tmp200 + tmp203;
		897	c_re(output[19 * ostride]) = tmp197 - tmp204;
		898	c_re(output[3 * ostride]) = tmp197 + tmp204;
		899	tmp208 = tmp203 - tmp200;
		900	c_im(output[27 * ostride]) = tmp207 - tmp208;
		901	c_im(output[11 * ostride]) = tmp207 + tmp208;
		902	tmp212 = tmp210 - tmp211;
		903	c_re(output[27 * ostride]) = tmp209 - tmp212;
		904	c_re(output[11 * ostride]) = tmp209 + tmp212;
		905	tmp214 = tmp210 + tmp211;
		906	c_im(output[19 * ostride]) = tmp213 - tmp214;
		907	c_im(output[3 * ostride]) = tmp213 + tmp214;
		908	}
		909	}
		910	{
		911	fftw_real tmp225;
		912	fftw_real tmp249;
		913	fftw_real tmp247;
		914	fftw_real tmp253;
		915	fftw_real tmp232;
		916	fftw_real tmp250;
		917	fftw_real tmp239;
		918	fftw_real tmp251;
		919	ASSERT_ALIGNED_DOUBLE();
		920	{
		921	fftw_real tmp217;
		922	fftw_real tmp224;
		923	fftw_real tmp243;
		924	fftw_real tmp246;
		925	ASSERT_ALIGNED_DOUBLE();
		926	tmp217 = tmp215 - tmp216;
		927	tmp224 = tmp220 - tmp223;
		928	tmp225 = tmp217 + tmp224;
		929	tmp249 = tmp217 - tmp224;
		930	tmp243 = tmp241 - tmp242;
		931	tmp246 = tmp244 - tmp245;
		932	tmp247 = tmp243 - tmp246;
		933	tmp253 = tmp243 + tmp246;
		934	}
		935	{
		936	fftw_real tmp228;
		937	fftw_real tmp231;
		938	fftw_real tmp235;
		939	fftw_real tmp238;
		940	ASSERT_ALIGNED_DOUBLE();
		941	tmp228 = tmp226 - tmp227;
		942	tmp231 = tmp229 - tmp230;
		943	tmp232 = (K555570233 * tmp228) + (K831469612 * tmp231);
		944	tmp250 = (K555570233 * tmp231) - (K831469612 * tmp228);
		945	tmp235 = tmp233 - tmp234;
		946	tmp238 = tmp236 - tmp237;
		947	tmp239 = (K555570233 * tmp235) - (K831469612 * tmp238);
		948	tmp251 = (K831469612 * tmp235) + (K555570233 * tmp238);
		949	}
		950	{
		951	fftw_real tmp240;
		952	fftw_real tmp248;
		953	fftw_real tmp252;
		954	fftw_real tmp254;
		955	ASSERT_ALIGNED_DOUBLE();
		956	tmp240 = tmp232 + tmp239;
		957	c_re(output[21 * ostride]) = tmp225 - tmp240;
		958	c_re(output[5 * ostride]) = tmp225 + tmp240;
		959	tmp248 = tmp239 - tmp232;
		960	c_im(output[29 * ostride]) = tmp247 - tmp248;
		961	c_im(output[13 * ostride]) = tmp247 + tmp248;
		962	tmp252 = tmp250 - tmp251;
		963	c_re(output[29 * ostride]) = tmp249 - tmp252;
		964	c_re(output[13 * ostride]) = tmp249 + tmp252;
		965	tmp254 = tmp250 + tmp251;
		966	c_im(output[21 * ostride]) = tmp253 - tmp254;
		967	c_im(output[5 * ostride]) = tmp253 + tmp254;
		968	}
		969	}
		970	{
		971	fftw_real tmp257;
		972	fftw_real tmp269;
		973	fftw_real tmp267;
		974	fftw_real tmp273;
		975	fftw_real tmp260;
		976	fftw_real tmp270;
		977	fftw_real tmp263;
		978	fftw_real tmp271;
		979	ASSERT_ALIGNED_DOUBLE();
		980	{
		981	fftw_real tmp255;
		982	fftw_real tmp256;
		983	fftw_real tmp265;
		984	fftw_real tmp266;
		985	ASSERT_ALIGNED_DOUBLE();
		986	tmp255 = tmp215 + tmp216;
		987	tmp256 = tmp245 + tmp244;
		988	tmp257 = tmp255 + tmp256;
		989	tmp269 = tmp255 - tmp256;
		990	tmp265 = tmp241 + tmp242;
		991	tmp266 = tmp220 + tmp223;
		992	tmp267 = tmp265 - tmp266;
		993	tmp273 = tmp265 + tmp266;
		994	}
		995	{
		996	fftw_real tmp258;
		997	fftw_real tmp259;
		998	fftw_real tmp261;
		999	fftw_real tmp262;
		1000	ASSERT_ALIGNED_DOUBLE();
		1001	tmp258 = tmp226 + tmp227;
		1002	tmp259 = tmp229 + tmp230;
		1003	tmp260 = (K980785280 * tmp258) + (K195090322 * tmp259);
		1004	tmp270 = (K980785280 * tmp259) - (K195090322 * tmp258);
		1005	tmp261 = tmp233 + tmp234;
		1006	tmp262 = tmp236 + tmp237;
		1007	tmp263 = (K980785280 * tmp261) - (K195090322 * tmp262);
		1008	tmp271 = (K195090322 * tmp261) + (K980785280 * tmp262);
		1009	}
		1010	{
		1011	fftw_real tmp264;
		1012	fftw_real tmp268;
		1013	fftw_real tmp272;
		1014	fftw_real tmp274;
		1015	ASSERT_ALIGNED_DOUBLE();
		1016	tmp264 = tmp260 + tmp263;
		1017	c_re(output[17 * ostride]) = tmp257 - tmp264;
		1018	c_re(output[ostride]) = tmp257 + tmp264;
		1019	tmp268 = tmp263 - tmp260;
		1020	c_im(output[25 * ostride]) = tmp267 - tmp268;
		1021	c_im(output[9 * ostride]) = tmp267 + tmp268;
		1022	tmp272 = tmp270 - tmp271;
		1023	c_re(output[25 * ostride]) = tmp269 - tmp272;
		1024	c_re(output[9 * ostride]) = tmp269 + tmp272;
		1025	tmp274 = tmp270 + tmp271;
		1026	c_im(output[17 * ostride]) = tmp273 - tmp274;
		1027	c_im(output[ostride]) = tmp273 + tmp274;
		1028	}
		1029	}
		1030	}
		1031
		1032	fftw_codelet_desc fftw_no_twiddle_32_desc =
		1033	{
		1034	"fftw_no_twiddle_32",
		1035	(void (*)()) fftw_no_twiddle_32,
		1036	32,
		1037	FFTW_FORWARD,
		1038	FFTW_NOTW,
		1039	705,
		1040	0,
		1041	(const int *) 0,
		1042	};

Subversion Repositories shark

(root)/shark/trunk/ports/fftw/fftw/fn_32.c - Rev 3