WebSVN - shark - Blame - Rev 1618 - /shark/branches/xen/ports/fftw/fftw/fni

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:57 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 -notwiddleinv 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 K195090322 = FFTW_KONST(+0.195090322016128267848284868477022240927691618);
		34	static const fftw_real K980785280 = FFTW_KONST(+0.980785280403230449126182236134239036973933731);
		35	static const fftw_real K831469612 = FFTW_KONST(+0.831469612302545237078788377617905756738560812);
		36	static const fftw_real K555570233 = FFTW_KONST(+0.555570233019602224742830813948532874374937191);
		37	static const fftw_real K382683432 = FFTW_KONST(+0.382683432365089771728459984030398866761344562);
		38	static const fftw_real K923879532 = FFTW_KONST(+0.923879532511286756128183189396788286822416626);
		39	static const fftw_real K707106781 = FFTW_KONST(+0.707106781186547524400844362104849039284835938);
		40
		41	/*
		42	* Generator Id's :
107	pj	43	* $Id: fni_32.c,v 1.2 2003-03-24 11:14:51 pj Exp $
		44	* $Id: fni_32.c,v 1.2 2003-03-24 11:14:51 pj Exp $
		45	* $Id: fni_32.c,v 1.2 2003-03-24 11:14:51 pj Exp $
2	pj	46	*/
		47
		48	void fftwi_no_twiddle_32(const fftw_complex input, fftw_complex output, int istride, int ostride)
		49	{
		50	fftw_real tmp7;
		51	fftw_real tmp339;
		52	fftw_real tmp70;
		53	fftw_real tmp313;
		54	fftw_real tmp97;
		55	fftw_real tmp215;
		56	fftw_real tmp179;
		57	fftw_real tmp241;
		58	fftw_real tmp14;
		59	fftw_real tmp314;
		60	fftw_real tmp77;
		61	fftw_real tmp340;
		62	fftw_real tmp182;
		63	fftw_real tmp216;
		64	fftw_real tmp104;
		65	fftw_real tmp242;
		66	fftw_real tmp153;
		67	fftw_real tmp236;
		68	fftw_real tmp53;
		69	fftw_real tmp60;
		70	fftw_real tmp287;
		71	fftw_real tmp336;
		72	fftw_real tmp360;
		73	fftw_real tmp290;
		74	fftw_real tmp293;
		75	fftw_real tmp294;
		76	fftw_real tmp170;
		77	fftw_real tmp233;
		78	fftw_real tmp333;
		79	fftw_real tmp359;
		80	fftw_real tmp164;
		81	fftw_real tmp234;
		82	fftw_real tmp173;
		83	fftw_real tmp237;
		84	fftw_real tmp22;
		85	fftw_real tmp318;
		86	fftw_real tmp343;
		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 tmp321;
		94	fftw_real tmp342;
		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 tmp278;
		105	fftw_real tmp329;
		106	fftw_real tmp357;
		107	fftw_real tmp281;
		108	fftw_real tmp284;
		109	fftw_real tmp285;
		110	fftw_real tmp143;
		111	fftw_real tmp226;
		112	fftw_real tmp326;
		113	fftw_real tmp356;
		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 tmp177;
		122	fftw_real tmp66;
		123	fftw_real tmp96;
		124	fftw_real tmp6;
		125	fftw_real tmp95;
		126	fftw_real tmp69;
		127	fftw_real tmp178;
		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	tmp177 = tmp1 - tmp2;
		139	tmp64 = c_im(input[0]);
		140	tmp65 = c_im(input[16 * istride]);
		141	tmp66 = tmp64 + tmp65;
		142	tmp96 = 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	tmp95 = tmp4 - tmp5;
		154	tmp67 = c_im(input[8 * istride]);
		155	tmp68 = c_im(input[24 * istride]);
		156	tmp69 = tmp67 + tmp68;
		157	tmp178 = tmp67 - tmp68;
		158	}
		159	tmp7 = tmp3 + tmp6;
		160	tmp339 = tmp3 - tmp6;
		161	tmp70 = tmp66 + tmp69;
		162	tmp313 = tmp66 - tmp69;
		163	tmp97 = tmp95 + tmp96;
		164	tmp215 = tmp96 - tmp95;
		165	tmp179 = tmp177 - tmp178;
		166	tmp241 = tmp177 + tmp178;
		167	}
		168	{
		169	fftw_real tmp10;
		170	fftw_real tmp98;
		171	fftw_real tmp73;
		172	fftw_real tmp99;
		173	fftw_real tmp13;
		174	fftw_real tmp102;
		175	fftw_real tmp76;
		176	fftw_real tmp101;
		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	tmp98 = tmp8 - tmp9;
		188	tmp71 = c_im(input[4 * istride]);
		189	tmp72 = c_im(input[20 * istride]);
		190	tmp73 = tmp71 + tmp72;
		191	tmp99 = 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	tmp102 = tmp11 - tmp12;
		203	tmp74 = c_im(input[28 * istride]);
		204	tmp75 = c_im(input[12 * istride]);
		205	tmp76 = tmp74 + tmp75;
		206	tmp101 = tmp74 - tmp75;
		207	}
		208	tmp14 = tmp10 + tmp13;
		209	tmp314 = tmp10 - tmp13;
		210	tmp77 = tmp73 + tmp76;
		211	tmp340 = tmp76 - tmp73;
		212	{
		213	fftw_real tmp180;
		214	fftw_real tmp181;
		215	fftw_real tmp100;
		216	fftw_real tmp103;
		217	ASSERT_ALIGNED_DOUBLE();
		218	tmp180 = tmp98 - tmp99;
		219	tmp181 = tmp102 + tmp101;
		220	tmp182 = K707106781 * (tmp180 + tmp181);
		221	tmp216 = K707106781 * (tmp180 - tmp181);
		222	tmp100 = tmp98 + tmp99;
		223	tmp103 = tmp101 - tmp102;
		224	tmp104 = K707106781 * (tmp100 + tmp103);
		225	tmp242 = K707106781 * (tmp103 - tmp100);
		226	}
		227	}
		228	{
		229	fftw_real tmp49;
		230	fftw_real tmp149;
		231	fftw_real tmp169;
		232	fftw_real tmp288;
		233	fftw_real tmp52;
		234	fftw_real tmp166;
		235	fftw_real tmp152;
		236	fftw_real tmp289;
		237	fftw_real tmp56;
		238	fftw_real tmp154;
		239	fftw_real tmp157;
		240	fftw_real tmp291;
		241	fftw_real tmp59;
		242	fftw_real tmp159;
		243	fftw_real tmp162;
		244	fftw_real tmp292;
		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	tmp288 = 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	tmp289 = tmp150 + tmp151;
		275	}
		276	{
		277	fftw_real tmp54;
		278	fftw_real tmp55;
		279	fftw_real tmp155;
		280	fftw_real tmp156;
		281	ASSERT_ALIGNED_DOUBLE();
		282	tmp54 = c_re(input[3 * istride]);
		283	tmp55 = c_re(input[19 * istride]);
		284	tmp56 = tmp54 + tmp55;
		285	tmp154 = tmp54 - tmp55;
		286	tmp155 = c_im(input[3 * istride]);
		287	tmp156 = c_im(input[19 * istride]);
		288	tmp157 = tmp155 - tmp156;
		289	tmp291 = tmp155 + tmp156;
		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	tmp292 = tmp160 + tmp161;
		305	}
		306	{
		307	fftw_real tmp334;
		308	fftw_real tmp335;
		309	fftw_real tmp331;
		310	fftw_real tmp332;
		311	ASSERT_ALIGNED_DOUBLE();
		312	tmp153 = tmp149 - tmp152;
		313	tmp236 = tmp149 + tmp152;
		314	tmp53 = tmp49 + tmp52;
		315	tmp60 = tmp56 + tmp59;
		316	tmp287 = tmp53 - tmp60;
		317	tmp334 = tmp49 - tmp52;
		318	tmp335 = tmp292 - tmp291;
		319	tmp336 = tmp334 - tmp335;
		320	tmp360 = tmp334 + tmp335;
		321	tmp290 = tmp288 + tmp289;
		322	tmp293 = tmp291 + tmp292;
		323	tmp294 = tmp290 - tmp293;
		324	tmp170 = tmp166 + tmp169;
		325	tmp233 = tmp169 - tmp166;
		326	tmp331 = tmp288 - tmp289;
		327	tmp332 = tmp56 - tmp59;
		328	tmp333 = tmp331 - tmp332;
		329	tmp359 = tmp332 + tmp331;
		330	{
		331	fftw_real tmp158;
		332	fftw_real tmp163;
		333	fftw_real tmp171;
		334	fftw_real tmp172;
		335	ASSERT_ALIGNED_DOUBLE();
		336	tmp158 = tmp154 - tmp157;
		337	tmp163 = tmp159 + tmp162;
		338	tmp164 = K707106781 * (tmp158 + tmp163);
		339	tmp234 = K707106781 * (tmp158 - tmp163);
		340	tmp171 = tmp154 + tmp157;
		341	tmp172 = tmp162 - tmp159;
		342	tmp173 = K707106781 * (tmp171 + tmp172);
		343	tmp237 = K707106781 * (tmp172 - tmp171);
		344	}
		345	}
		346	}
		347	{
		348	fftw_real tmp18;
		349	fftw_real tmp106;
		350	fftw_real tmp81;
		351	fftw_real tmp110;
		352	fftw_real tmp21;
		353	fftw_real tmp109;
		354	fftw_real tmp84;
		355	fftw_real tmp107;
		356	fftw_real tmp316;
		357	fftw_real tmp317;
		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	tmp106 = tmp16 - tmp17;
		369	tmp79 = c_im(input[2 * istride]);
		370	tmp80 = c_im(input[18 * istride]);
		371	tmp81 = tmp79 + tmp80;
		372	tmp110 = 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	tmp109 = tmp19 - tmp20;
		384	tmp82 = c_im(input[10 * istride]);
		385	tmp83 = c_im(input[26 * istride]);
		386	tmp84 = tmp82 + tmp83;
		387	tmp107 = tmp82 - tmp83;
		388	}
		389	tmp22 = tmp18 + tmp21;
		390	tmp316 = tmp18 - tmp21;
		391	tmp317 = tmp81 - tmp84;
		392	tmp318 = tmp316 - tmp317;
		393	tmp343 = tmp316 + tmp317;
		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 = (K923879532 * tmp108) - (K382683432 * tmp111);
		404	tmp185 = (K923879532 * tmp111) + (K382683432 * tmp108);
		405	tmp218 = tmp106 + tmp107;
		406	tmp219 = tmp110 - tmp109;
		407	tmp220 = (K382683432 * tmp218) - (K923879532 * tmp219);
		408	tmp245 = (K382683432 * tmp219) + (K923879532 * tmp218);
		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 tmp319;
		421	fftw_real tmp320;
		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	tmp319 = tmp25 - tmp28;
		455	tmp320 = tmp88 - tmp91;
		456	tmp321 = tmp319 + tmp320;
		457	tmp342 = tmp320 - tmp319;
		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 = (K923879532 * tmp115) - (K382683432 * tmp118);
		469	tmp221 = tmp114 - tmp113;
		470	tmp222 = tmp116 + tmp117;
		471	tmp223 = (K923879532 * tmp221) + (K382683432 * tmp222);
		472	tmp244 = (K382683432 * tmp221) - (K923879532 * tmp222);
		473	}
		474	}
		475	{
		476	fftw_real tmp34;
		477	fftw_real tmp122;
		478	fftw_real tmp142;
		479	fftw_real tmp279;
		480	fftw_real tmp37;
		481	fftw_real tmp139;
		482	fftw_real tmp125;
		483	fftw_real tmp280;
		484	fftw_real tmp41;
		485	fftw_real tmp127;
		486	fftw_real tmp130;
		487	fftw_real tmp282;
		488	fftw_real tmp44;
		489	fftw_real tmp132;
		490	fftw_real tmp135;
		491	fftw_real tmp283;
		492	ASSERT_ALIGNED_DOUBLE();
		493	{
		494	fftw_real tmp32;
		495	fftw_real tmp33;
		496	fftw_real tmp140;
		497	fftw_real tmp141;
		498	ASSERT_ALIGNED_DOUBLE();
		499	tmp32 = c_re(input[istride]);
		500	tmp33 = c_re(input[17 * istride]);
		501	tmp34 = tmp32 + tmp33;
		502	tmp122 = tmp32 - tmp33;
		503	tmp140 = c_im(input[istride]);
		504	tmp141 = c_im(input[17 * istride]);
		505	tmp142 = tmp140 - tmp141;
		506	tmp279 = tmp140 + tmp141;
		507	}
		508	{
		509	fftw_real tmp35;
		510	fftw_real tmp36;
		511	fftw_real tmp123;
		512	fftw_real tmp124;
		513	ASSERT_ALIGNED_DOUBLE();
		514	tmp35 = c_re(input[9 * istride]);
		515	tmp36 = c_re(input[25 * istride]);
		516	tmp37 = tmp35 + tmp36;
		517	tmp139 = tmp35 - tmp36;
		518	tmp123 = c_im(input[9 * istride]);
		519	tmp124 = c_im(input[25 * istride]);
		520	tmp125 = tmp123 - tmp124;
		521	tmp280 = tmp123 + tmp124;
		522	}
		523	{
		524	fftw_real tmp39;
		525	fftw_real tmp40;
		526	fftw_real tmp128;
		527	fftw_real tmp129;
		528	ASSERT_ALIGNED_DOUBLE();
		529	tmp39 = c_re(input[5 * istride]);
		530	tmp40 = c_re(input[21 * istride]);
		531	tmp41 = tmp39 + tmp40;
		532	tmp127 = tmp39 - tmp40;
		533	tmp128 = c_im(input[5 * istride]);
		534	tmp129 = c_im(input[21 * istride]);
		535	tmp130 = tmp128 - tmp129;
		536	tmp282 = tmp128 + tmp129;
		537	}
		538	{
		539	fftw_real tmp42;
		540	fftw_real tmp43;
		541	fftw_real tmp133;
		542	fftw_real tmp134;
		543	ASSERT_ALIGNED_DOUBLE();
		544	tmp42 = c_re(input[29 * istride]);
		545	tmp43 = c_re(input[13 * istride]);
		546	tmp44 = tmp42 + tmp43;
		547	tmp132 = tmp42 - tmp43;
		548	tmp133 = c_im(input[29 * istride]);
		549	tmp134 = c_im(input[13 * istride]);
		550	tmp135 = tmp133 - tmp134;
		551	tmp283 = tmp133 + tmp134;
		552	}
		553	{
		554	fftw_real tmp327;
		555	fftw_real tmp328;
		556	fftw_real tmp324;
		557	fftw_real tmp325;
		558	ASSERT_ALIGNED_DOUBLE();
		559	tmp126 = tmp122 - tmp125;
		560	tmp229 = tmp122 + tmp125;
		561	tmp38 = tmp34 + tmp37;
		562	tmp45 = tmp41 + tmp44;
		563	tmp278 = tmp38 - tmp45;
		564	tmp327 = tmp34 - tmp37;
		565	tmp328 = tmp283 - tmp282;
		566	tmp329 = tmp327 - tmp328;
		567	tmp357 = tmp327 + tmp328;
		568	tmp281 = tmp279 + tmp280;
		569	tmp284 = tmp282 + tmp283;
		570	tmp285 = tmp281 - tmp284;
		571	tmp143 = tmp139 + tmp142;
		572	tmp226 = tmp142 - tmp139;
		573	tmp324 = tmp279 - tmp280;
		574	tmp325 = tmp41 - tmp44;
		575	tmp326 = tmp324 - tmp325;
		576	tmp356 = tmp325 + tmp324;
		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 * (tmp131 - tmp136);
		587	tmp144 = tmp127 + tmp130;
		588	tmp145 = tmp135 - tmp132;
		589	tmp146 = K707106781 * (tmp144 + tmp145);
		590	tmp230 = K707106781 * (tmp145 - tmp144);
		591	}
		592	}
		593	}
		594	{
		595	fftw_real tmp277;
		596	fftw_real tmp301;
		597	fftw_real tmp304;
		598	fftw_real tmp306;
		599	fftw_real tmp296;
		600	fftw_real tmp300;
		601	fftw_real tmp299;
		602	fftw_real tmp305;
		603	ASSERT_ALIGNED_DOUBLE();
		604	{
		605	fftw_real tmp275;
		606	fftw_real tmp276;
		607	fftw_real tmp302;
		608	fftw_real tmp303;
		609	ASSERT_ALIGNED_DOUBLE();
		610	tmp275 = tmp70 - tmp77;
		611	tmp276 = tmp22 - tmp29;
		612	tmp277 = tmp275 - tmp276;
		613	tmp301 = tmp276 + tmp275;
		614	tmp302 = tmp278 + tmp285;
		615	tmp303 = tmp294 - tmp287;
		616	tmp304 = K707106781 * (tmp302 + tmp303);
		617	tmp306 = K707106781 * (tmp303 - tmp302);
		618	}
		619	{
		620	fftw_real tmp286;
		621	fftw_real tmp295;
		622	fftw_real tmp297;
		623	fftw_real tmp298;
		624	ASSERT_ALIGNED_DOUBLE();
		625	tmp286 = tmp278 - tmp285;
		626	tmp295 = tmp287 + tmp294;
		627	tmp296 = K707106781 * (tmp286 - tmp295);
		628	tmp300 = K707106781 * (tmp286 + tmp295);
		629	tmp297 = tmp7 - tmp14;
		630	tmp298 = tmp92 - tmp85;
		631	tmp299 = tmp297 + tmp298;
		632	tmp305 = tmp297 - tmp298;
		633	}
		634	c_im(output[28 * ostride]) = tmp277 - tmp296;
		635	c_im(output[12 * ostride]) = tmp277 + tmp296;
		636	c_re(output[20 * ostride]) = tmp299 - tmp300;
		637	c_re(output[4 * ostride]) = tmp299 + tmp300;
		638	c_im(output[20 * ostride]) = tmp301 - tmp304;
		639	c_im(output[4 * ostride]) = tmp301 + tmp304;
		640	c_re(output[28 * ostride]) = tmp305 - tmp306;
		641	c_re(output[12 * ostride]) = tmp305 + tmp306;
		642	}
		643	{
		644	fftw_real tmp31;
		645	fftw_real tmp311;
		646	fftw_real tmp310;
		647	fftw_real tmp312;
		648	fftw_real tmp62;
		649	fftw_real tmp63;
		650	fftw_real tmp94;
		651	fftw_real tmp307;
		652	ASSERT_ALIGNED_DOUBLE();
		653	{
		654	fftw_real tmp15;
		655	fftw_real tmp30;
		656	fftw_real tmp308;
		657	fftw_real tmp309;
		658	ASSERT_ALIGNED_DOUBLE();
		659	tmp15 = tmp7 + tmp14;
		660	tmp30 = tmp22 + tmp29;
		661	tmp31 = tmp15 + tmp30;
		662	tmp311 = tmp15 - tmp30;
		663	tmp308 = tmp281 + tmp284;
		664	tmp309 = tmp290 + tmp293;
		665	tmp310 = tmp308 + tmp309;
		666	tmp312 = tmp309 - tmp308;
		667	}
		668	{
		669	fftw_real tmp46;
		670	fftw_real tmp61;
		671	fftw_real tmp78;
		672	fftw_real tmp93;
		673	ASSERT_ALIGNED_DOUBLE();
		674	tmp46 = tmp38 + tmp45;
		675	tmp61 = tmp53 + tmp60;
		676	tmp62 = tmp46 + tmp61;
		677	tmp63 = tmp46 - tmp61;
		678	tmp78 = tmp70 + tmp77;
		679	tmp93 = tmp85 + tmp92;
		680	tmp94 = tmp78 - tmp93;
		681	tmp307 = tmp78 + tmp93;
		682	}
		683	c_re(output[16 * ostride]) = tmp31 - tmp62;
		684	c_re(output[0]) = tmp31 + tmp62;
		685	c_im(output[8 * ostride]) = tmp63 + tmp94;
		686	c_im(output[24 * ostride]) = tmp94 - tmp63;
		687	c_im(output[16 * ostride]) = tmp307 - tmp310;
		688	c_im(output[0]) = tmp307 + tmp310;
		689	c_re(output[24 * ostride]) = tmp311 - tmp312;
		690	c_re(output[8 * ostride]) = tmp311 + tmp312;
		691	}
		692	{
		693	fftw_real tmp121;
		694	fftw_real tmp189;
		695	fftw_real tmp187;
		696	fftw_real tmp193;
		697	fftw_real tmp148;
		698	fftw_real tmp190;
		699	fftw_real tmp175;
		700	fftw_real tmp191;
		701	ASSERT_ALIGNED_DOUBLE();
		702	{
		703	fftw_real tmp105;
		704	fftw_real tmp120;
		705	fftw_real tmp183;
		706	fftw_real tmp186;
		707	ASSERT_ALIGNED_DOUBLE();
		708	tmp105 = tmp97 - tmp104;
		709	tmp120 = tmp112 - tmp119;
		710	tmp121 = tmp105 - tmp120;
		711	tmp189 = tmp105 + tmp120;
		712	tmp183 = tmp179 - tmp182;
		713	tmp186 = tmp184 - tmp185;
		714	tmp187 = tmp183 + tmp186;
		715	tmp193 = tmp183 - tmp186;
		716	}
		717	{
		718	fftw_real tmp138;
		719	fftw_real tmp147;
		720	fftw_real tmp165;
		721	fftw_real tmp174;
		722	ASSERT_ALIGNED_DOUBLE();
		723	tmp138 = tmp126 - tmp137;
		724	tmp147 = tmp143 - tmp146;
		725	tmp148 = (K555570233 * tmp138) - (K831469612 * tmp147);
		726	tmp190 = (K831469612 * tmp138) + (K555570233 * tmp147);
		727	tmp165 = tmp153 - tmp164;
		728	tmp174 = tmp170 - tmp173;
		729	tmp175 = (K555570233 * tmp165) + (K831469612 * tmp174);
		730	tmp191 = (K555570233 * tmp174) - (K831469612 * tmp165);
		731	}
		732	{
		733	fftw_real tmp176;
		734	fftw_real tmp188;
		735	fftw_real tmp192;
		736	fftw_real tmp194;
		737	ASSERT_ALIGNED_DOUBLE();
		738	tmp176 = tmp148 - tmp175;
		739	c_im(output[29 * ostride]) = tmp121 - tmp176;
		740	c_im(output[13 * ostride]) = tmp121 + tmp176;
		741	tmp188 = tmp148 + tmp175;
		742	c_re(output[21 * ostride]) = tmp187 - tmp188;
		743	c_re(output[5 * ostride]) = tmp187 + tmp188;
		744	tmp192 = tmp190 + tmp191;
		745	c_im(output[21 * ostride]) = tmp189 - tmp192;
		746	c_im(output[5 * ostride]) = tmp189 + tmp192;
		747	tmp194 = tmp191 - tmp190;
		748	c_re(output[29 * ostride]) = tmp193 - tmp194;
		749	c_re(output[13 * ostride]) = tmp193 + tmp194;
		750	}
		751	}
		752	{
		753	fftw_real tmp197;
		754	fftw_real tmp209;
		755	fftw_real tmp207;
		756	fftw_real tmp213;
		757	fftw_real tmp200;
		758	fftw_real tmp210;
		759	fftw_real tmp203;
		760	fftw_real tmp211;
		761	ASSERT_ALIGNED_DOUBLE();
		762	{
		763	fftw_real tmp195;
		764	fftw_real tmp196;
		765	fftw_real tmp205;
		766	fftw_real tmp206;
		767	ASSERT_ALIGNED_DOUBLE();
		768	tmp195 = tmp97 + tmp104;
		769	tmp196 = tmp185 + tmp184;
		770	tmp197 = tmp195 - tmp196;
		771	tmp209 = tmp195 + tmp196;
		772	tmp205 = tmp179 + tmp182;
		773	tmp206 = tmp112 + tmp119;
		774	tmp207 = tmp205 + tmp206;
		775	tmp213 = tmp205 - tmp206;
		776	}
		777	{
		778	fftw_real tmp198;
		779	fftw_real tmp199;
		780	fftw_real tmp201;
		781	fftw_real tmp202;
		782	ASSERT_ALIGNED_DOUBLE();
		783	tmp198 = tmp126 + tmp137;
		784	tmp199 = tmp143 + tmp146;
		785	tmp200 = (K980785280 * tmp198) - (K195090322 * tmp199);
		786	tmp210 = (K195090322 * tmp198) + (K980785280 * tmp199);
		787	tmp201 = tmp153 + tmp164;
		788	tmp202 = tmp170 + tmp173;
		789	tmp203 = (K980785280 * tmp201) + (K195090322 * tmp202);
		790	tmp211 = (K980785280 * tmp202) - (K195090322 * tmp201);
		791	}
		792	{
		793	fftw_real tmp204;
		794	fftw_real tmp208;
		795	fftw_real tmp212;
		796	fftw_real tmp214;
		797	ASSERT_ALIGNED_DOUBLE();
		798	tmp204 = tmp200 - tmp203;
		799	c_im(output[25 * ostride]) = tmp197 - tmp204;
		800	c_im(output[9 * ostride]) = tmp197 + tmp204;
		801	tmp208 = tmp200 + tmp203;
		802	c_re(output[17 * ostride]) = tmp207 - tmp208;
		803	c_re(output[ostride]) = tmp207 + tmp208;
		804	tmp212 = tmp210 + tmp211;
		805	c_im(output[17 * ostride]) = tmp209 - tmp212;
		806	c_im(output[ostride]) = tmp209 + tmp212;
		807	tmp214 = tmp211 - tmp210;
		808	c_re(output[25 * ostride]) = tmp213 - tmp214;
		809	c_re(output[9 * ostride]) = tmp213 + tmp214;
		810	}
		811	}
		812	{
		813	fftw_real tmp323;
		814	fftw_real tmp347;
		815	fftw_real tmp350;
		816	fftw_real tmp352;
		817	fftw_real tmp338;
		818	fftw_real tmp346;
		819	fftw_real tmp345;
		820	fftw_real tmp351;
		821	ASSERT_ALIGNED_DOUBLE();
		822	{
		823	fftw_real tmp315;
		824	fftw_real tmp322;
		825	fftw_real tmp348;
		826	fftw_real tmp349;
		827	ASSERT_ALIGNED_DOUBLE();
		828	tmp315 = tmp313 - tmp314;
		829	tmp322 = K707106781 * (tmp318 - tmp321);
		830	tmp323 = tmp315 + tmp322;
		831	tmp347 = tmp315 - tmp322;
		832	tmp348 = (K382683432 * tmp329) - (K923879532 * tmp326);
		833	tmp349 = (K923879532 * tmp333) + (K382683432 * tmp336);
		834	tmp350 = tmp348 - tmp349;
		835	tmp352 = tmp348 + tmp349;
		836	}
		837	{
		838	fftw_real tmp330;
		839	fftw_real tmp337;
		840	fftw_real tmp341;
		841	fftw_real tmp344;
		842	ASSERT_ALIGNED_DOUBLE();
		843	tmp330 = (K382683432 * tmp326) + (K923879532 * tmp329);
		844	tmp337 = (K382683432 * tmp333) - (K923879532 * tmp336);
		845	tmp338 = tmp330 + tmp337;
		846	tmp346 = tmp337 - tmp330;
		847	tmp341 = tmp339 - tmp340;
		848	tmp344 = K707106781 * (tmp342 - tmp343);
		849	tmp345 = tmp341 - tmp344;
		850	tmp351 = tmp341 + tmp344;
		851	}
		852	c_im(output[22 * ostride]) = tmp323 - tmp338;
		853	c_im(output[6 * ostride]) = tmp323 + tmp338;
		854	c_re(output[30 * ostride]) = tmp345 - tmp346;
		855	c_re(output[14 * ostride]) = tmp345 + tmp346;
		856	c_im(output[30 * ostride]) = tmp347 - tmp350;
		857	c_im(output[14 * ostride]) = tmp347 + tmp350;
		858	c_re(output[22 * ostride]) = tmp351 - tmp352;
		859	c_re(output[6 * ostride]) = tmp351 + tmp352;
		860	}
		861	{
		862	fftw_real tmp355;
		863	fftw_real tmp367;
		864	fftw_real tmp370;
		865	fftw_real tmp372;
		866	fftw_real tmp362;
		867	fftw_real tmp366;
		868	fftw_real tmp365;
		869	fftw_real tmp371;
		870	ASSERT_ALIGNED_DOUBLE();
		871	{
		872	fftw_real tmp353;
		873	fftw_real tmp354;
		874	fftw_real tmp368;
		875	fftw_real tmp369;
		876	ASSERT_ALIGNED_DOUBLE();
		877	tmp353 = tmp314 + tmp313;
		878	tmp354 = K707106781 * (tmp343 + tmp342);
		879	tmp355 = tmp353 + tmp354;
		880	tmp367 = tmp353 - tmp354;
		881	tmp368 = (K923879532 * tmp357) - (K382683432 * tmp356);
		882	tmp369 = (K382683432 * tmp359) + (K923879532 * tmp360);
		883	tmp370 = tmp368 - tmp369;
		884	tmp372 = tmp368 + tmp369;
		885	}
		886	{
		887	fftw_real tmp358;
		888	fftw_real tmp361;
		889	fftw_real tmp363;
		890	fftw_real tmp364;
		891	ASSERT_ALIGNED_DOUBLE();
		892	tmp358 = (K923879532 * tmp356) + (K382683432 * tmp357);
		893	tmp361 = (K923879532 * tmp359) - (K382683432 * tmp360);
		894	tmp362 = tmp358 + tmp361;
		895	tmp366 = tmp361 - tmp358;
		896	tmp363 = tmp339 + tmp340;
		897	tmp364 = K707106781 * (tmp318 + tmp321);
		898	tmp365 = tmp363 - tmp364;
		899	tmp371 = tmp363 + tmp364;
		900	}
		901	c_im(output[18 * ostride]) = tmp355 - tmp362;
		902	c_im(output[2 * ostride]) = tmp355 + tmp362;
		903	c_re(output[26 * ostride]) = tmp365 - tmp366;
		904	c_re(output[10 * ostride]) = tmp365 + tmp366;
		905	c_im(output[26 * ostride]) = tmp367 - tmp370;
		906	c_im(output[10 * ostride]) = tmp367 + tmp370;
		907	c_re(output[18 * ostride]) = tmp371 - tmp372;
		908	c_re(output[2 * ostride]) = tmp371 + tmp372;
		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 = (K195090322 * tmp228) + (K980785280 * tmp231);
		944	tmp250 = (K195090322 * tmp231) - (K980785280 * tmp228);
		945	tmp235 = tmp233 - tmp234;
		946	tmp238 = tmp236 - tmp237;
		947	tmp239 = (K195090322 * tmp235) - (K980785280 * tmp238);
		948	tmp251 = (K980785280 * tmp235) + (K195090322 * 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_im(output[23 * ostride]) = tmp225 - tmp240;
		958	c_im(output[7 * ostride]) = tmp225 + tmp240;
		959	tmp248 = tmp239 - tmp232;
		960	c_re(output[31 * ostride]) = tmp247 - tmp248;
		961	c_re(output[15 * ostride]) = tmp247 + tmp248;
		962	tmp252 = tmp250 - tmp251;
		963	c_im(output[31 * ostride]) = tmp249 - tmp252;
		964	c_im(output[15 * ostride]) = tmp249 + tmp252;
		965	tmp254 = tmp250 + tmp251;
		966	c_re(output[23 * ostride]) = tmp253 - tmp254;
		967	c_re(output[7 * 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 = (K831469612 * tmp258) + (K555570233 * tmp259);
		1004	tmp270 = (K831469612 * tmp259) - (K555570233 * tmp258);
		1005	tmp261 = tmp233 + tmp234;
		1006	tmp262 = tmp236 + tmp237;
		1007	tmp263 = (K831469612 * tmp261) - (K555570233 * tmp262);
		1008	tmp271 = (K555570233 * tmp261) + (K831469612 * 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_im(output[19 * ostride]) = tmp257 - tmp264;
		1018	c_im(output[3 * ostride]) = tmp257 + tmp264;
		1019	tmp268 = tmp263 - tmp260;
		1020	c_re(output[27 * ostride]) = tmp267 - tmp268;
		1021	c_re(output[11 * ostride]) = tmp267 + tmp268;
		1022	tmp272 = tmp270 - tmp271;
		1023	c_im(output[27 * ostride]) = tmp269 - tmp272;
		1024	c_im(output[11 * ostride]) = tmp269 + tmp272;
		1025	tmp274 = tmp270 + tmp271;
		1026	c_re(output[19 * ostride]) = tmp273 - tmp274;
		1027	c_re(output[3 * ostride]) = tmp273 + tmp274;
		1028	}
		1029	}
		1030	}
		1031
		1032	fftw_codelet_desc fftwi_no_twiddle_32_desc =
		1033	{
		1034	"fftwi_no_twiddle_32",
		1035	(void (*)()) fftwi_no_twiddle_32,
		1036	32,
		1037	FFTW_BACKWARD,
		1038	FFTW_NOTW,
		1039	716,
		1040	0,
		1041	(const int *) 0,
		1042	};

Subversion Repositories shark

(root)/shark/branches/xen/ports/fftw/fftw/fni_32.c - Rev 1618