WebSVN - shark - Blame - Rev 107 - /shark/trunk/ports/fftw/rfftw/frc

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

Subversion Repositories shark

(root)/shark/trunk/ports/fftw/rfftw/frc_64.c @ 1620 - Rev 107