Mercurial > mplayer.hg
comparison libfaad2/sbr_dct.c @ 12527:4a370c80fe5c
update to the 2.0 release of faad, patch by adland
author | diego |
---|---|
date | Wed, 02 Jun 2004 22:59:04 +0000 |
parents | 3185f64f6350 |
children | d81145997036 |
comparison
equal
deleted
inserted
replaced
12526:e183ad37d24c | 12527:4a370c80fe5c |
---|---|
1 /* | 1 /* |
2 ** FAAD2 - Freeware Advanced Audio (AAC) Decoder including SBR decoding | 2 ** FAAD2 - Freeware Advanced Audio (AAC) Decoder including SBR decoding |
3 ** Copyright (C) 2003 M. Bakker, Ahead Software AG, http://www.nero.com | 3 ** Copyright (C) 2003-2004 M. Bakker, Ahead Software AG, http://www.nero.com |
4 ** | 4 ** |
5 ** This program is free software; you can redistribute it and/or modify | 5 ** This program is free software; you can redistribute it and/or modify |
6 ** it under the terms of the GNU General Public License as published by | 6 ** it under the terms of the GNU General Public License as published by |
7 ** the Free Software Foundation; either version 2 of the License, or | 7 ** the Free Software Foundation; either version 2 of the License, or |
8 ** (at your option) any later version. | 8 ** (at your option) any later version. |
9 ** | 9 ** |
10 ** This program is distributed in the hope that it will be useful, | 10 ** This program is distributed in the hope that it will be useful, |
11 ** but WITHOUT ANY WARRANTY; without even the implied warranty of | 11 ** but WITHOUT ANY WARRANTY; without even the implied warranty of |
12 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | 12 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
13 ** GNU General Public License for more details. | 13 ** GNU General Public License for more details. |
14 ** | 14 ** |
15 ** You should have received a copy of the GNU General Public License | 15 ** You should have received a copy of the GNU General Public License |
16 ** along with this program; if not, write to the Free Software | 16 ** along with this program; if not, write to the Free Software |
17 ** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. | 17 ** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. |
18 ** | 18 ** |
19 ** Any non-GPL usage of this software or parts of this software is strictly | 19 ** Any non-GPL usage of this software or parts of this software is strictly |
20 ** forbidden. | 20 ** forbidden. |
21 ** | 21 ** |
22 ** Commercial non-GPL licensing of this software is possible. | 22 ** Commercial non-GPL licensing of this software is possible. |
23 ** For more info contact Ahead Software through Mpeg4AAClicense@nero.com. | 23 ** For more info contact Ahead Software through Mpeg4AAClicense@nero.com. |
24 ** | 24 ** |
25 ** $Id: sbr_dct.c,v 1.4 2003/09/24 11:52:12 menno Exp $ | 25 ** $Id: sbr_dct.c,v 1.2 2003/10/03 22:22:27 alex Exp $ |
26 **/ | 26 **/ |
27 | 27 |
28 #include "common.h" | 28 #include "common.h" |
29 | 29 |
30 #ifdef SBR_DEC | 30 #ifdef SBR_DEC |
32 #ifdef _MSC_VER | 32 #ifdef _MSC_VER |
33 #pragma warning(disable:4305) | 33 #pragma warning(disable:4305) |
34 #pragma warning(disable:4244) | 34 #pragma warning(disable:4244) |
35 #endif | 35 #endif |
36 | 36 |
37 #define MUL_C_R(B,A) MUL_R_C(A,B) | |
38 | 37 |
39 #include "sbr_dct.h" | 38 #include "sbr_dct.h" |
40 | 39 |
41 #ifdef SBR_LOW_POWER | 40 #ifdef SBR_LOW_POWER |
42 | 41 |
69 real_t f241, f242, f243, f244, f245, f246, f247, f248, f249, f250; | 68 real_t f241, f242, f243, f244, f245, f246, f247, f248, f249, f250; |
70 real_t f251, f252, f253, f254, f255, f256, f257, f258, f259, f260; | 69 real_t f251, f252, f253, f254, f255, f256, f257, f258, f259, f260; |
71 real_t f261, f262, f263, f264, f265, f266, f267, f268, f269, f270; | 70 real_t f261, f262, f263, f264, f265, f266, f267, f268, f269, f270; |
72 real_t f271, f272; | 71 real_t f271, f272; |
73 | 72 |
74 f0 = MUL_C_R(COEF_CONST(0.7071067811865476), x[16]); | 73 f0 = MUL_C(COEF_CONST(0.7071067811865476), x[16]); |
75 f1 = x[0] - f0; | 74 f1 = x[0] - f0; |
76 f2 = x[0] + f0; | 75 f2 = x[0] + f0; |
77 f3 = x[8] + x[24]; | 76 f3 = x[8] + x[24]; |
78 f4 = MUL_C_R(COEF_CONST(1.3065629648763766), x[8]); | 77 f4 = MUL_C(COEF_CONST(1.3065629648763766), x[8]); |
79 f5 = MUL_C_R(COEF_CONST((-0.9238795325112866)), f3); | 78 f5 = MUL_C(COEF_CONST((-0.9238795325112866)), f3); |
80 f6 = MUL_C_R(COEF_CONST((-0.5411961001461967)), x[24]); | 79 f6 = MUL_C(COEF_CONST((-0.5411961001461967)), x[24]); |
81 f7 = f4 + f5; | 80 f7 = f4 + f5; |
82 f8 = f6 - f5; | 81 f8 = f6 - f5; |
83 f9 = f2 - f8; | 82 f9 = f2 - f8; |
84 f10 = f2 + f8; | 83 f10 = f2 + f8; |
85 f11 = f1 - f7; | 84 f11 = f1 - f7; |
86 f12 = f1 + f7; | 85 f12 = f1 + f7; |
87 f13 = x[4] + x[28]; | 86 f13 = x[4] + x[28]; |
88 f14 = MUL_C_R(COEF_CONST(1.1758756024193588), x[4]); | 87 f14 = MUL_C(COEF_CONST(1.1758756024193588), x[4]); |
89 f15 = MUL_C_R(COEF_CONST((-0.9807852804032304)), f13); | 88 f15 = MUL_C(COEF_CONST((-0.9807852804032304)), f13); |
90 f16 = MUL_C_R(COEF_CONST((-0.7856949583871021)), x[28]); | 89 f16 = MUL_C(COEF_CONST((-0.7856949583871021)), x[28]); |
91 f17 = f14 + f15; | 90 f17 = f14 + f15; |
92 f18 = f16 - f15; | 91 f18 = f16 - f15; |
93 f19 = x[12] + x[20]; | 92 f19 = x[12] + x[20]; |
94 f20 = MUL_C_R(COEF_CONST(1.3870398453221473), x[12]); | 93 f20 = MUL_C(COEF_CONST(1.3870398453221473), x[12]); |
95 f21 = MUL_C_R(COEF_CONST((-0.8314696123025455)), f19); | 94 f21 = MUL_C(COEF_CONST((-0.8314696123025455)), f19); |
96 f22 = MUL_C_R(COEF_CONST((-0.2758993792829436)), x[20]); | 95 f22 = MUL_C(COEF_CONST((-0.2758993792829436)), x[20]); |
97 f23 = f20 + f21; | 96 f23 = f20 + f21; |
98 f24 = f22 - f21; | 97 f24 = f22 - f21; |
99 f25 = f18 - f24; | 98 f25 = f18 - f24; |
100 f26 = f18 + f24; | 99 f26 = f18 + f24; |
101 f27 = MUL_C_R(COEF_CONST(0.7071067811865476), f25); | 100 f27 = MUL_C(COEF_CONST(0.7071067811865476), f25); |
102 f28 = f17 - f23; | 101 f28 = f17 - f23; |
103 f29 = f17 + f23; | 102 f29 = f17 + f23; |
104 f30 = MUL_C_R(COEF_CONST(0.7071067811865476), f29); | 103 f30 = MUL_C(COEF_CONST(0.7071067811865476), f29); |
105 f31 = f27 - f30; | 104 f31 = f27 - f30; |
106 f32 = f27 + f30; | 105 f32 = f27 + f30; |
107 f33 = f10 - f26; | 106 f33 = f10 - f26; |
108 f34 = f10 + f26; | 107 f34 = f10 + f26; |
109 f35 = f12 - f32; | 108 f35 = f12 - f32; |
111 f37 = f11 - f31; | 110 f37 = f11 - f31; |
112 f38 = f11 + f31; | 111 f38 = f11 + f31; |
113 f39 = f9 - f28; | 112 f39 = f9 - f28; |
114 f40 = f9 + f28; | 113 f40 = f9 + f28; |
115 f41 = x[2] + x[30]; | 114 f41 = x[2] + x[30]; |
116 f42 = MUL_C_R(COEF_CONST(1.0932018670017569), x[2]); | 115 f42 = MUL_C(COEF_CONST(1.0932018670017569), x[2]); |
117 f43 = MUL_C_R(COEF_CONST((-0.9951847266721969)), f41); | 116 f43 = MUL_C(COEF_CONST((-0.9951847266721969)), f41); |
118 f44 = MUL_C_R(COEF_CONST((-0.8971675863426368)), x[30]); | 117 f44 = MUL_C(COEF_CONST((-0.8971675863426368)), x[30]); |
119 f45 = f42 + f43; | 118 f45 = f42 + f43; |
120 f46 = f44 - f43; | 119 f46 = f44 - f43; |
121 f47 = x[6] + x[26]; | 120 f47 = x[6] + x[26]; |
122 f48 = MUL_C_R(COEF_CONST(1.2472250129866711), x[6]); | 121 f48 = MUL_C(COEF_CONST(1.2472250129866711), x[6]); |
123 f49 = MUL_C_R(COEF_CONST((-0.9569403357322089)), f47); | 122 f49 = MUL_C(COEF_CONST((-0.9569403357322089)), f47); |
124 f50 = MUL_C_R(COEF_CONST((-0.6666556584777469)), x[26]); | 123 f50 = MUL_C(COEF_CONST((-0.6666556584777469)), x[26]); |
125 f51 = f48 + f49; | 124 f51 = f48 + f49; |
126 f52 = f50 - f49; | 125 f52 = f50 - f49; |
127 f53 = x[10] + x[22]; | 126 f53 = x[10] + x[22]; |
128 f54 = MUL_C_R(COEF_CONST(1.3533180011743526), x[10]); | 127 f54 = MUL_C(COEF_CONST(1.3533180011743526), x[10]); |
129 f55 = MUL_C_R(COEF_CONST((-0.8819212643483551)), f53); | 128 f55 = MUL_C(COEF_CONST((-0.8819212643483551)), f53); |
130 f56 = MUL_C_R(COEF_CONST((-0.4105245275223575)), x[22]); | 129 f56 = MUL_C(COEF_CONST((-0.4105245275223575)), x[22]); |
131 f57 = f54 + f55; | 130 f57 = f54 + f55; |
132 f58 = f56 - f55; | 131 f58 = f56 - f55; |
133 f59 = x[14] + x[18]; | 132 f59 = x[14] + x[18]; |
134 f60 = MUL_C_R(COEF_CONST(1.4074037375263826), x[14]); | 133 f60 = MUL_C(COEF_CONST(1.4074037375263826), x[14]); |
135 f61 = MUL_C_R(COEF_CONST((-0.7730104533627369)), f59); | 134 f61 = MUL_C(COEF_CONST((-0.7730104533627369)), f59); |
136 f62 = MUL_C_R(COEF_CONST((-0.1386171691990913)), x[18]); | 135 f62 = MUL_C(COEF_CONST((-0.1386171691990913)), x[18]); |
137 f63 = f60 + f61; | 136 f63 = f60 + f61; |
138 f64 = f62 - f61; | 137 f64 = f62 - f61; |
139 f65 = f46 - f64; | 138 f65 = f46 - f64; |
140 f66 = f46 + f64; | 139 f66 = f46 + f64; |
141 f67 = f52 - f58; | 140 f67 = f52 - f58; |
142 f68 = f52 + f58; | 141 f68 = f52 + f58; |
143 f69 = f66 - f68; | 142 f69 = f66 - f68; |
144 f70 = f66 + f68; | 143 f70 = f66 + f68; |
145 f71 = MUL_C_R(COEF_CONST(0.7071067811865476), f69); | 144 f71 = MUL_C(COEF_CONST(0.7071067811865476), f69); |
146 f72 = f65 + f67; | 145 f72 = f65 + f67; |
147 f73 = MUL_C_R(COEF_CONST(1.3065629648763766), f65); | 146 f73 = MUL_C(COEF_CONST(1.3065629648763766), f65); |
148 f74 = MUL_C_R(COEF_CONST((-0.9238795325112866)), f72); | 147 f74 = MUL_C(COEF_CONST((-0.9238795325112866)), f72); |
149 f75 = MUL_C_R(COEF_CONST((-0.5411961001461967)), f67); | 148 f75 = MUL_C(COEF_CONST((-0.5411961001461967)), f67); |
150 f76 = f73 + f74; | 149 f76 = f73 + f74; |
151 f77 = f75 - f74; | 150 f77 = f75 - f74; |
152 f78 = f45 - f63; | 151 f78 = f45 - f63; |
153 f79 = f45 + f63; | 152 f79 = f45 + f63; |
154 f80 = f51 - f57; | 153 f80 = f51 - f57; |
155 f81 = f51 + f57; | 154 f81 = f51 + f57; |
156 f82 = f79 + f81; | 155 f82 = f79 + f81; |
157 f83 = MUL_C_R(COEF_CONST(1.3065629648763770), f79); | 156 f83 = MUL_C(COEF_CONST(1.3065629648763770), f79); |
158 f84 = MUL_C_R(COEF_CONST((-0.3826834323650904)), f82); | 157 f84 = MUL_C(COEF_CONST((-0.3826834323650904)), f82); |
159 f85 = MUL_C_R(COEF_CONST(0.5411961001461961), f81); | 158 f85 = MUL_C(COEF_CONST(0.5411961001461961), f81); |
160 f86 = f83 + f84; | 159 f86 = f83 + f84; |
161 f87 = f85 - f84; | 160 f87 = f85 - f84; |
162 f88 = f78 - f80; | 161 f88 = f78 - f80; |
163 f89 = f78 + f80; | 162 f89 = f78 + f80; |
164 f90 = MUL_C_R(COEF_CONST(0.7071067811865476), f89); | 163 f90 = MUL_C(COEF_CONST(0.7071067811865476), f89); |
165 f91 = f77 - f87; | 164 f91 = f77 - f87; |
166 f92 = f77 + f87; | 165 f92 = f77 + f87; |
167 f93 = f71 - f90; | 166 f93 = f71 - f90; |
168 f94 = f71 + f90; | 167 f94 = f71 + f90; |
169 f95 = f76 - f86; | 168 f95 = f76 - f86; |
183 f109 = f35 - f95; | 182 f109 = f35 - f95; |
184 f110 = f35 + f95; | 183 f110 = f35 + f95; |
185 f111 = f33 - f88; | 184 f111 = f33 - f88; |
186 f112 = f33 + f88; | 185 f112 = f33 + f88; |
187 f113 = x[1] + x[31]; | 186 f113 = x[1] + x[31]; |
188 f114 = MUL_C_R(COEF_CONST(1.0478631305325901), x[1]); | 187 f114 = MUL_C(COEF_CONST(1.0478631305325901), x[1]); |
189 f115 = MUL_C_R(COEF_CONST((-0.9987954562051724)), f113); | 188 f115 = MUL_C(COEF_CONST((-0.9987954562051724)), f113); |
190 f116 = MUL_C_R(COEF_CONST((-0.9497277818777548)), x[31]); | 189 f116 = MUL_C(COEF_CONST((-0.9497277818777548)), x[31]); |
191 f117 = f114 + f115; | 190 f117 = f114 + f115; |
192 f118 = f116 - f115; | 191 f118 = f116 - f115; |
193 f119 = x[5] + x[27]; | 192 f119 = x[5] + x[27]; |
194 f120 = MUL_C_R(COEF_CONST(1.2130114330978077), x[5]); | 193 f120 = MUL_C(COEF_CONST(1.2130114330978077), x[5]); |
195 f121 = MUL_C_R(COEF_CONST((-0.9700312531945440)), f119); | 194 f121 = MUL_C(COEF_CONST((-0.9700312531945440)), f119); |
196 f122 = MUL_C_R(COEF_CONST((-0.7270510732912803)), x[27]); | 195 f122 = MUL_C(COEF_CONST((-0.7270510732912803)), x[27]); |
197 f123 = f120 + f121; | 196 f123 = f120 + f121; |
198 f124 = f122 - f121; | 197 f124 = f122 - f121; |
199 f125 = x[9] + x[23]; | 198 f125 = x[9] + x[23]; |
200 f126 = MUL_C_R(COEF_CONST(1.3315443865537255), x[9]); | 199 f126 = MUL_C(COEF_CONST(1.3315443865537255), x[9]); |
201 f127 = MUL_C_R(COEF_CONST((-0.9039892931234433)), f125); | 200 f127 = MUL_C(COEF_CONST((-0.9039892931234433)), f125); |
202 f128 = MUL_C_R(COEF_CONST((-0.4764341996931612)), x[23]); | 201 f128 = MUL_C(COEF_CONST((-0.4764341996931612)), x[23]); |
203 f129 = f126 + f127; | 202 f129 = f126 + f127; |
204 f130 = f128 - f127; | 203 f130 = f128 - f127; |
205 f131 = x[13] + x[19]; | 204 f131 = x[13] + x[19]; |
206 f132 = MUL_C_R(COEF_CONST(1.3989068359730781), x[13]); | 205 f132 = MUL_C(COEF_CONST(1.3989068359730781), x[13]); |
207 f133 = MUL_C_R(COEF_CONST((-0.8032075314806453)), f131); | 206 f133 = MUL_C(COEF_CONST((-0.8032075314806453)), f131); |
208 f134 = MUL_C_R(COEF_CONST((-0.2075082269882124)), x[19]); | 207 f134 = MUL_C(COEF_CONST((-0.2075082269882124)), x[19]); |
209 f135 = f132 + f133; | 208 f135 = f132 + f133; |
210 f136 = f134 - f133; | 209 f136 = f134 - f133; |
211 f137 = x[17] + x[15]; | 210 f137 = x[17] + x[15]; |
212 f138 = MUL_C_R(COEF_CONST(1.4125100802019777), x[17]); | 211 f138 = MUL_C(COEF_CONST(1.4125100802019777), x[17]); |
213 f139 = MUL_C_R(COEF_CONST((-0.6715589548470187)), f137); | 212 f139 = MUL_C(COEF_CONST((-0.6715589548470187)), f137); |
214 f140 = MUL_C_R(COEF_CONST(0.0693921705079402), x[15]); | 213 f140 = MUL_C(COEF_CONST(0.0693921705079402), x[15]); |
215 f141 = f138 + f139; | 214 f141 = f138 + f139; |
216 f142 = f140 - f139; | 215 f142 = f140 - f139; |
217 f143 = x[21] + x[11]; | 216 f143 = x[21] + x[11]; |
218 f144 = MUL_C_R(COEF_CONST(1.3718313541934939), x[21]); | 217 f144 = MUL_C(COEF_CONST(1.3718313541934939), x[21]); |
219 f145 = MUL_C_R(COEF_CONST((-0.5141027441932219)), f143); | 218 f145 = MUL_C(COEF_CONST((-0.5141027441932219)), f143); |
220 f146 = MUL_C_R(COEF_CONST(0.3436258658070501), x[11]); | 219 f146 = MUL_C(COEF_CONST(0.3436258658070501), x[11]); |
221 f147 = f144 + f145; | 220 f147 = f144 + f145; |
222 f148 = f146 - f145; | 221 f148 = f146 - f145; |
223 f149 = x[25] + x[7]; | 222 f149 = x[25] + x[7]; |
224 f150 = MUL_C_R(COEF_CONST(1.2784339185752409), x[25]); | 223 f150 = MUL_C(COEF_CONST(1.2784339185752409), x[25]); |
225 f151 = MUL_C_R(COEF_CONST((-0.3368898533922200)), f149); | 224 f151 = MUL_C(COEF_CONST((-0.3368898533922200)), f149); |
226 f152 = MUL_C_R(COEF_CONST(0.6046542117908008), x[7]); | 225 f152 = MUL_C(COEF_CONST(0.6046542117908008), x[7]); |
227 f153 = f150 + f151; | 226 f153 = f150 + f151; |
228 f154 = f152 - f151; | 227 f154 = f152 - f151; |
229 f155 = x[29] + x[3]; | 228 f155 = x[29] + x[3]; |
230 f156 = MUL_C_R(COEF_CONST(1.1359069844201433), x[29]); | 229 f156 = MUL_C(COEF_CONST(1.1359069844201433), x[29]); |
231 f157 = MUL_C_R(COEF_CONST((-0.1467304744553624)), f155); | 230 f157 = MUL_C(COEF_CONST((-0.1467304744553624)), f155); |
232 f158 = MUL_C_R(COEF_CONST(0.8424460355094185), x[3]); | 231 f158 = MUL_C(COEF_CONST(0.8424460355094185), x[3]); |
233 f159 = f156 + f157; | 232 f159 = f156 + f157; |
234 f160 = f158 - f157; | 233 f160 = f158 - f157; |
235 f161 = f118 - f142; | 234 f161 = f118 - f142; |
236 f162 = f118 + f142; | 235 f162 = f118 + f142; |
237 f163 = f117 - f141; | 236 f163 = f117 - f141; |
247 f173 = f136 - f160; | 246 f173 = f136 - f160; |
248 f174 = f136 + f160; | 247 f174 = f136 + f160; |
249 f175 = f135 - f159; | 248 f175 = f135 - f159; |
250 f176 = f135 + f159; | 249 f176 = f135 + f159; |
251 f177 = f161 + f163; | 250 f177 = f161 + f163; |
252 f178 = MUL_C_R(COEF_CONST(1.1758756024193588), f161); | 251 f178 = MUL_C(COEF_CONST(1.1758756024193588), f161); |
253 f179 = MUL_C_R(COEF_CONST((-0.9807852804032304)), f177); | 252 f179 = MUL_C(COEF_CONST((-0.9807852804032304)), f177); |
254 f180 = MUL_C_R(COEF_CONST((-0.7856949583871021)), f163); | 253 f180 = MUL_C(COEF_CONST((-0.7856949583871021)), f163); |
255 f181 = f178 + f179; | 254 f181 = f178 + f179; |
256 f182 = f180 - f179; | 255 f182 = f180 - f179; |
257 f183 = f165 + f167; | 256 f183 = f165 + f167; |
258 f184 = MUL_C_R(COEF_CONST(1.3870398453221475), f165); | 257 f184 = MUL_C(COEF_CONST(1.3870398453221475), f165); |
259 f185 = MUL_C_R(COEF_CONST((-0.5555702330196022)), f183); | 258 f185 = MUL_C(COEF_CONST((-0.5555702330196022)), f183); |
260 f186 = MUL_C_R(COEF_CONST(0.2758993792829431), f167); | 259 f186 = MUL_C(COEF_CONST(0.2758993792829431), f167); |
261 f187 = f184 + f185; | 260 f187 = f184 + f185; |
262 f188 = f186 - f185; | 261 f188 = f186 - f185; |
263 f189 = f169 + f171; | 262 f189 = f169 + f171; |
264 f190 = MUL_C_R(COEF_CONST(0.7856949583871022), f169); | 263 f190 = MUL_C(COEF_CONST(0.7856949583871022), f169); |
265 f191 = MUL_C_R(COEF_CONST(0.1950903220161283), f189); | 264 f191 = MUL_C(COEF_CONST(0.1950903220161283), f189); |
266 f192 = MUL_C_R(COEF_CONST(1.1758756024193586), f171); | 265 f192 = MUL_C(COEF_CONST(1.1758756024193586), f171); |
267 f193 = f190 + f191; | 266 f193 = f190 + f191; |
268 f194 = f192 - f191; | 267 f194 = f192 - f191; |
269 f195 = f173 + f175; | 268 f195 = f173 + f175; |
270 f196 = MUL_C_R(COEF_CONST((-0.2758993792829430)), f173); | 269 f196 = MUL_C(COEF_CONST((-0.2758993792829430)), f173); |
271 f197 = MUL_C_R(COEF_CONST(0.8314696123025452), f195); | 270 f197 = MUL_C(COEF_CONST(0.8314696123025452), f195); |
272 f198 = MUL_C_R(COEF_CONST(1.3870398453221475), f175); | 271 f198 = MUL_C(COEF_CONST(1.3870398453221475), f175); |
273 f199 = f196 + f197; | 272 f199 = f196 + f197; |
274 f200 = f198 - f197; | 273 f200 = f198 - f197; |
275 f201 = f162 - f170; | 274 f201 = f162 - f170; |
276 f202 = f162 + f170; | 275 f202 = f162 + f170; |
277 f203 = f164 - f172; | 276 f203 = f164 - f172; |
287 f213 = f188 - f200; | 286 f213 = f188 - f200; |
288 f214 = f188 + f200; | 287 f214 = f188 + f200; |
289 f215 = f187 - f199; | 288 f215 = f187 - f199; |
290 f216 = f187 + f199; | 289 f216 = f187 + f199; |
291 f217 = f201 + f203; | 290 f217 = f201 + f203; |
292 f218 = MUL_C_R(COEF_CONST(1.3065629648763766), f201); | 291 f218 = MUL_C(COEF_CONST(1.3065629648763766), f201); |
293 f219 = MUL_C_R(COEF_CONST((-0.9238795325112866)), f217); | 292 f219 = MUL_C(COEF_CONST((-0.9238795325112866)), f217); |
294 f220 = MUL_C_R(COEF_CONST((-0.5411961001461967)), f203); | 293 f220 = MUL_C(COEF_CONST((-0.5411961001461967)), f203); |
295 f221 = f218 + f219; | 294 f221 = f218 + f219; |
296 f222 = f220 - f219; | 295 f222 = f220 - f219; |
297 f223 = f205 + f207; | 296 f223 = f205 + f207; |
298 f224 = MUL_C_R(COEF_CONST(0.5411961001461969), f205); | 297 f224 = MUL_C(COEF_CONST(0.5411961001461969), f205); |
299 f225 = MUL_C_R(COEF_CONST(0.3826834323650898), f223); | 298 f225 = MUL_C(COEF_CONST(0.3826834323650898), f223); |
300 f226 = MUL_C_R(COEF_CONST(1.3065629648763766), f207); | 299 f226 = MUL_C(COEF_CONST(1.3065629648763766), f207); |
301 f227 = f224 + f225; | 300 f227 = f224 + f225; |
302 f228 = f226 - f225; | 301 f228 = f226 - f225; |
303 f229 = f209 + f211; | 302 f229 = f209 + f211; |
304 f230 = MUL_C_R(COEF_CONST(1.3065629648763766), f209); | 303 f230 = MUL_C(COEF_CONST(1.3065629648763766), f209); |
305 f231 = MUL_C_R(COEF_CONST((-0.9238795325112866)), f229); | 304 f231 = MUL_C(COEF_CONST((-0.9238795325112866)), f229); |
306 f232 = MUL_C_R(COEF_CONST((-0.5411961001461967)), f211); | 305 f232 = MUL_C(COEF_CONST((-0.5411961001461967)), f211); |
307 f233 = f230 + f231; | 306 f233 = f230 + f231; |
308 f234 = f232 - f231; | 307 f234 = f232 - f231; |
309 f235 = f213 + f215; | 308 f235 = f213 + f215; |
310 f236 = MUL_C_R(COEF_CONST(0.5411961001461969), f213); | 309 f236 = MUL_C(COEF_CONST(0.5411961001461969), f213); |
311 f237 = MUL_C_R(COEF_CONST(0.3826834323650898), f235); | 310 f237 = MUL_C(COEF_CONST(0.3826834323650898), f235); |
312 f238 = MUL_C_R(COEF_CONST(1.3065629648763766), f215); | 311 f238 = MUL_C(COEF_CONST(1.3065629648763766), f215); |
313 f239 = f236 + f237; | 312 f239 = f236 + f237; |
314 f240 = f238 - f237; | 313 f240 = f238 - f237; |
315 f241 = f202 - f206; | 314 f241 = f202 - f206; |
316 f242 = f202 + f206; | 315 f242 = f202 + f206; |
317 f243 = f204 - f208; | 316 f243 = f204 - f208; |
328 f254 = f234 + f240; | 327 f254 = f234 + f240; |
329 f255 = f233 - f239; | 328 f255 = f233 - f239; |
330 f256 = f233 + f239; | 329 f256 = f233 + f239; |
331 f257 = f241 - f243; | 330 f257 = f241 - f243; |
332 f258 = f241 + f243; | 331 f258 = f241 + f243; |
333 f259 = MUL_C_R(COEF_CONST(0.7071067811865474), f257); | 332 f259 = MUL_C(COEF_CONST(0.7071067811865474), f257); |
334 f260 = MUL_C_R(COEF_CONST(0.7071067811865474), f258); | 333 f260 = MUL_C(COEF_CONST(0.7071067811865474), f258); |
335 f261 = f245 - f247; | 334 f261 = f245 - f247; |
336 f262 = f245 + f247; | 335 f262 = f245 + f247; |
337 f263 = MUL_C_R(COEF_CONST(0.7071067811865474), f261); | 336 f263 = MUL_C(COEF_CONST(0.7071067811865474), f261); |
338 f264 = MUL_C_R(COEF_CONST(0.7071067811865474), f262); | 337 f264 = MUL_C(COEF_CONST(0.7071067811865474), f262); |
339 f265 = f249 - f251; | 338 f265 = f249 - f251; |
340 f266 = f249 + f251; | 339 f266 = f249 + f251; |
341 f267 = MUL_C_R(COEF_CONST(0.7071067811865474), f265); | 340 f267 = MUL_C(COEF_CONST(0.7071067811865474), f265); |
342 f268 = MUL_C_R(COEF_CONST(0.7071067811865474), f266); | 341 f268 = MUL_C(COEF_CONST(0.7071067811865474), f266); |
343 f269 = f253 - f255; | 342 f269 = f253 - f255; |
344 f270 = f253 + f255; | 343 f270 = f253 + f255; |
345 f271 = MUL_C_R(COEF_CONST(0.7071067811865474), f269); | 344 f271 = MUL_C(COEF_CONST(0.7071067811865474), f269); |
346 f272 = MUL_C_R(COEF_CONST(0.7071067811865474), f270); | 345 f272 = MUL_C(COEF_CONST(0.7071067811865474), f270); |
347 y[31] = f98 - f242; | 346 y[31] = f98 - f242; |
348 y[0] = f98 + f242; | 347 y[0] = f98 + f242; |
349 y[30] = f100 - f250; | 348 y[30] = f100 - f250; |
350 y[1] = f100 + f250; | 349 y[1] = f100 + f250; |
351 y[29] = f102 - f254; | 350 y[29] = f102 - f254; |
386 real_t f21, f22, f23, f24, f25, f26, f27, f28, f29, f30; | 385 real_t f21, f22, f23, f24, f25, f26, f27, f28, f29, f30; |
387 real_t f31, f32, f33, f34, f35, f36, f37, f38, f39, f40; | 386 real_t f31, f32, f33, f34, f35, f36, f37, f38, f39, f40; |
388 real_t f41, f42, f43, f44, f45, f46, f47, f48, f49, f50; | 387 real_t f41, f42, f43, f44, f45, f46, f47, f48, f49, f50; |
389 real_t f51, f52, f53, f54, f55, f56, f57, f58, f59, f60; | 388 real_t f51, f52, f53, f54, f55, f56, f57, f58, f59, f60; |
390 real_t f61, f62, f65, f66, f67, f68, f71, f72, f73, f74; | 389 real_t f61, f62, f65, f66, f67, f68, f71, f72, f73, f74; |
391 real_t f75, f76, f77, f78, f79, f80, f81, f82, f85, f86; | 390 real_t f75, f76, f77, f78, f79, f80, f81, f82, f83, f85; |
392 real_t f87, f88, f91, f92, f93, f94, f95, f96, f97, f98; | 391 real_t f87, f88, f91, f92, f93, f94, f95, f96, f97, f98; |
393 real_t f99, f100, f101, f102, f103, f104, f105, f106, f107, f108; | 392 real_t f99, f100, f101, f102, f103, f104, f105, f106, f107, f108; |
394 real_t f109, f110, f111, f112, f113, f114, f115, f116, f117, f118; | 393 real_t f109, f110, f111, f112, f113, f114, f115, f116, f117, f118; |
395 real_t f119, f120, f121, f122, f123, f124, f125, f126, f127, f128; | 394 real_t f119, f120, f121, f122, f123, f124, f125, f126, f129, f130; |
396 real_t f129, f130, f133, f134, f135, f136, f139, f140, f141, f142; | 395 real_t f131, f132, f135, f136, f137, f138, f141, f142, f143, f144; |
397 real_t f145, f146, f147, f148, f151, f152, f153, f154, f155, f156; | 396 real_t f147, f148, f149, f150, f151, f152, f153, f154, f155, f156; |
398 real_t f157, f158, f159, f160, f161, f162, f163, f164, f165, f166; | 397 real_t f157, f158, f159, f160, f161, f162, f163, f164, f165, f166; |
399 real_t f167, f168, f169, f170, f171, f172, f173, f174, f175, f176; | 398 real_t f167, f168, f169, f170, f171, f172, f173, f174, f175, f176; |
400 real_t f177, f178, f179, f180, f181, f182, f183, f184, f185, f186; | 399 real_t f177, f178, f179, f180, f181, f182, f183, f184, f185, f186; |
401 real_t f187, f188, f189, f190, f191, f192, f193, f194, f195, f196; | 400 real_t f187, f188, f189, f190, f191, f192, f193, f194, f195, f196; |
402 real_t f197, f198, f199, f200, f201, f202, f203, f204, f205, f206; | 401 real_t f197, f198, f199, f200, f201, f202, f203, f204, f205, f206; |
403 real_t f207, f208, f209, f210, f211, f213, f214, f215, f216, f217; | 402 real_t f207, f208, f209, f210, f211, f212, f213, f214, f215, f216; |
404 real_t f218, f219, f220, f221, f222, f223, f224, f225, f226, f227; | 403 real_t f217, f218, f219, f220, f221, f222, f223, f224, f225, f226; |
405 real_t f228, f229, f230, f231, f232, f233, f234, f235, f236, f237; | 404 real_t f227, f228, f229, f230, f231, f232, f233, f234, f235, f236; |
406 real_t f238, f239, f240, f241, f242, f243, f244, f245, f246, f247; | 405 real_t f237, f238, f239, f240, f241, f242, f243, f244, f245, f246; |
407 real_t f248, f249, f250, f251, f252, f253, f254, f255, f256, f257; | 406 real_t f247, f248, f249, f250, f251, f252, f253, f254, f255, f256; |
408 real_t f258, f259, f260, f261, f262, f263, f264, f265, f266, f267; | 407 real_t f257, f258, f259, f260, f261, f262, f265, f266, f267, f268; |
409 real_t f268, f269, f270, f271, f272, f273, f274, f275, f276, f277; | 408 real_t f271, f272, f273, f274, f277, f278, f279, f280, f283, f284; |
410 real_t f279, f280, f295, f296, f297, f298, f299, f300, f301, f302; | 409 real_t f285, f286, f289, f290, f291, f292, f295, f296, f297, f298; |
411 real_t f303, f304, f305, f306, f307, f308, f309, f310, f311, f312; | 410 real_t f301, f302, f303, f304, f307, f308, f309, f310, f311, f312; |
412 real_t f313, f314, f315, f316, f317, f318, f319, f320, f321, f322; | 411 real_t f313, f314, f315, f316, f317, f318, f319, f320, f321, f322; |
413 real_t f323, f324, f325, f326, f327, f328, f329, f330, f331, f332; | 412 real_t f323, f324, f325, f326, f327, f328, f329, f330, f331, f332; |
414 real_t f333, f334, f335, f336, f337, f338, f339, f340, f341, f342; | 413 real_t f333, f334, f335, f336, f337, f338, f339, f340, f341, f342; |
415 real_t f343, f344, f345, f346, f347, f348, f349, f350, f351, f352; | 414 real_t f343, f344, f345, f346, f347, f348, f349, f350, f351, f352; |
416 real_t f353, f354, f355, f356, f357, f358, f359, f360, f361, f362; | 415 real_t f353, f354, f355, f356, f357, f358, f359, f360, f361, f362; |
418 real_t f373, f374, f375, f376, f377, f378, f379, f380, f381, f382; | 417 real_t f373, f374, f375, f376, f377, f378, f379, f380, f381, f382; |
419 real_t f383, f384, f385, f386, f387, f388, f389, f390, f391, f392; | 418 real_t f383, f384, f385, f386, f387, f388, f389, f390, f391, f392; |
420 real_t f393, f394, f395, f396, f397, f398, f399, f400, f401, f402; | 419 real_t f393, f394, f395, f396, f397, f398, f399, f400, f401, f402; |
421 real_t f403, f404, f405, f406, f407, f408, f409, f410, f411, f412; | 420 real_t f403, f404, f405, f406, f407, f408, f409, f410, f411, f412; |
422 real_t f413, f414, f415, f416, f417, f418, f419, f420, f421, f422; | 421 real_t f413, f414, f415, f416, f417, f418, f419, f420, f421, f422; |
423 real_t f423, f424, f425, f426, f427, f428, f429, f430, f431, f432; | 422 real_t f423, f424, f425, f426, f427, f428, f429, f430, f431, f433; |
424 real_t f433, f434, f435, f436, f437, f438, f439, f440, f441, f442; | 423 real_t f434, f435, f436, f437, f438, f439, f440, f441, f442, f443; |
425 real_t f443, f444, f445, f446, f447, f448, f449, f450, f451, f452; | 424 real_t f444, f445, f446, f447, f448, f449, f450, f451, f452, f453; |
426 real_t f453, f454, f455, f456, f457, f458, f459, f460, f461, f462; | 425 real_t f454, f455, f456, f457, f458, f459, f460, f461, f462, f463; |
427 real_t f463, f464, f465, f466, f467, f468, f469, f470, f471, f472; | 426 real_t f464, f465, f466, f467, f468, f469, f470, f471, f472, f473; |
428 real_t f473, f474, f475, f476, f477, f478, f479, f480, f481, f482; | 427 real_t f474, f475, f476, f477, f478, f479, f480, f481, f482, f483; |
429 real_t f483, f484, f485, f486, f487, f488, f489, f490, f491, f492; | 428 real_t f484, f485, f486, f487, f488, f489, f490, f491, f492, f493; |
430 real_t f493, f494, f495, f496, f497, f498, f499, f500, f501, f502; | 429 real_t f494, f495, f496, f497, f498, f499, f500, f501, f502, f503; |
431 real_t f503, f504, f505, f506, f507, f508, f509, f510, f511, f512; | 430 real_t f504, f505, f506, f507, f508, f509, f510, f511, f512, f513; |
432 real_t f513, f514, f515, f516, f517, f518, f519, f520, f521, f522; | 431 real_t f514, f515, f516, f517, f518, f519, f520, f521, f522, f523; |
433 real_t f523, f524, f525, f526, f527, f528, f529, f530, f531, f532; | 432 real_t f524, f525, f526, f527, f528, f529, f530, f531, f532, f533; |
434 real_t f533, f534, f535, f536, f537, f538, f539, f540, f541, f542; | 433 real_t f534, f535, f536, f537, f538, f539, f540, f541, f542, f543; |
435 real_t f543, f544, f545, f546, f547, f548, f549, f550, f551, f552; | 434 real_t f544, f546, f547, f548, f549, f550, f551, f552, f553, f554; |
436 real_t f553, f554, f557, f558, f559, f560, f563, f564, f565, f566; | 435 real_t f555, f556, f557, f558, f559, f560, f561, f562, f563, f564; |
437 real_t f569, f570, f571, f572, f575, f576, f577, f578, f581, f582; | 436 real_t f565, f566, f567, f568, f569, f570, f571, f572, f573, f574; |
438 real_t f583, f584, f587, f588, f589, f590, f593, f594, f595, f596; | 437 real_t f575, f576, f577, f578, f579, f580, f581, f582, f583, f584; |
439 real_t f599, f600, f601, f602, f605, f606, f607, f608, f611, f612; | 438 real_t f585, f586, f587, f588, f589, f590, f591, f592, f593, f594; |
440 real_t f613, f614, f617, f618, f619, f620, f623, f624, f625, f626; | 439 real_t f595, f596, f597, f598, f599, f600, f601, f602, f603, f604; |
441 real_t f629, f630, f631, f632, f635, f636, f637, f638, f641, f642; | 440 real_t f605, f606, f607, f608, f609, f610, f611, f612, f613, f614; |
442 real_t f643, f644; | 441 real_t f615, f616, f617, f618, f619, f620, f621, f622, f623, f624; |
443 static real_t t2[64]; | 442 real_t f625, f626, f627, f628; |
443 ALIGN static real_t t2[64]; | |
444 | 444 |
445 for (i0=0; i0<32; i0++) | 445 for (i0=0; i0<32; i0++) |
446 { | 446 { |
447 t2[2*i0+1] = x[i0] - x[-i0+63]; | 447 t2[2*i0+1] = x[i0] - x[-i0+63]; |
448 t2[2*i0] = x[i0] + x[-i0+63]; | 448 t2[2*i0] = x[i0] + x[-i0+63]; |
507 f59 = f51 + f57; | 507 f59 = f51 + f57; |
508 f60 = f53 - f55; | 508 f60 = f53 - f55; |
509 f61 = f53 + f55; | 509 f61 = f53 + f55; |
510 f62 = f59 - f61; | 510 f62 = f59 - f61; |
511 y[0] = f59 + f61; | 511 y[0] = f59 + f61; |
512 y[32] = MUL_C_R(COEF_CONST(0.7071067811865476), f62); | 512 y[32] = MUL_C(COEF_CONST(0.7071067811865476), f62); |
513 f65 = f58 + f60; | 513 f65 = f58 + f60; |
514 f66 = MUL_C_R(COEF_CONST(1.3065629648763766), f58); | 514 f66 = MUL_C(COEF_CONST(1.3065629648763766), f58); |
515 f67 = MUL_C_R(COEF_CONST((-0.9238795325112866)), f65); | 515 f67 = MUL_C(COEF_CONST((-0.9238795325112866)), f65); |
516 f68 = MUL_C_R(COEF_CONST((-0.5411961001461967)), f60); | 516 f68 = MUL_C(COEF_CONST((-0.5411961001461967)), f60); |
517 y[48] = f66 + f67; | 517 y[48] = f66 + f67; |
518 y[16] = f68 - f67; | 518 y[16] = f68 - f67; |
519 f71 = f52 - f54; | 519 f71 = f50 + f56; |
520 f72 = f52 + f54; | 520 f72 = MUL_C(COEF_CONST(1.1758756024193588), f50); |
521 f73 = MUL_C_R(COEF_CONST(0.7071067811865476), f72); | 521 f73 = MUL_C(COEF_CONST((-0.9807852804032304)), f71); |
522 f74 = MUL_C_R(COEF_CONST(0.7071067811865476), f71); | 522 f74 = MUL_C(COEF_CONST((-0.7856949583871021)), f56); |
523 f75 = f50 - f73; | 523 f75 = f72 + f73; |
524 f76 = f50 + f73; | 524 f76 = f74 - f73; |
525 f77 = f56 - f74; | 525 f77 = f52 + f54; |
526 f78 = f56 + f74; | 526 f78 = MUL_C(COEF_CONST(1.3870398453221473), f52); |
527 f79 = f78 + f76; | 527 f79 = MUL_C(COEF_CONST((-0.8314696123025455)), f77); |
528 f80 = MUL_C_R(COEF_CONST((-0.7856949583871021)), f78); | 528 f80 = MUL_C(COEF_CONST((-0.2758993792829436)), f54); |
529 f81 = MUL_C_R(COEF_CONST(0.9807852804032304), f79); | 529 f81 = f78 + f79; |
530 f82 = MUL_C_R(COEF_CONST(1.1758756024193588), f76); | 530 f82 = f80 - f79; |
531 y[8] = f80 + f81; | 531 f83 = f76 - f82; |
532 y[56] = f82 - f81; | 532 y[8] = f76 + f82; |
533 f85 = f77 + f75; | 533 f85 = MUL_C(COEF_CONST(0.7071067811865476), f83); |
534 f86 = MUL_C_R(COEF_CONST(0.2758993792829431), f77); | 534 y[56] = f75 - f81; |
535 f87 = MUL_C_R(COEF_CONST(0.5555702330196022), f85); | 535 f87 = f75 + f81; |
536 f88 = MUL_C_R(COEF_CONST(1.3870398453221475), f75); | 536 f88 = MUL_C(COEF_CONST(0.7071067811865476), f87); |
537 y[40] = f86 + f87; | 537 y[40] = f85 - f88; |
538 y[24] = f88 - f87; | 538 y[24] = f85 + f88; |
539 f91 = f40 - f42; | 539 f91 = f36 - f38; |
540 f92 = f40 + f42; | 540 f92 = f36 + f38; |
541 f93 = MUL_C_R(COEF_CONST(0.7071067811865476), f92); | 541 f93 = f40 - f42; |
542 f94 = MUL_C_R(COEF_CONST(0.7071067811865476), f91); | 542 f94 = f40 + f42; |
543 f95 = f38 - f44; | 543 f95 = f44 - f46; |
544 f96 = f38 + f44; | 544 f96 = f44 + f46; |
545 f97 = MUL_C_R(COEF_CONST(0.7071067811865476), f96); | 545 f97 = MUL_C(COEF_CONST(0.7071067811865476), f94); |
546 f98 = MUL_C_R(COEF_CONST(0.7071067811865476), f95); | 546 f98 = f34 - f97; |
547 f99 = f34 - f93; | 547 f99 = f34 + f97; |
548 f100 = f34 + f93; | 548 f100 = f92 + f96; |
549 f101 = f48 - f94; | 549 f101 = MUL_C(COEF_CONST(1.3065629648763766), f92); |
550 f102 = f48 + f94; | 550 f102 = MUL_C(COEF_CONST((-0.9238795325112866)), f100); |
551 f103 = f36 - f97; | 551 f103 = MUL_C(COEF_CONST((-0.5411961001461967)), f96); |
552 f104 = f36 + f97; | 552 f104 = f101 + f102; |
553 f105 = f46 - f98; | 553 f105 = f103 - f102; |
554 f106 = f46 + f98; | 554 f106 = f99 - f105; |
555 f107 = f106 + f104; | 555 f107 = f99 + f105; |
556 f108 = MUL_C_R(COEF_CONST((-0.5411961001461969)), f106); | 556 f108 = f98 - f104; |
557 f109 = MUL_C_R(COEF_CONST(0.9238795325112867), f107); | 557 f109 = f98 + f104; |
558 f110 = MUL_C_R(COEF_CONST(1.3065629648763766), f104); | 558 f110 = f91 + f95; |
559 f111 = f108 + f109; | 559 f111 = MUL_C(COEF_CONST(1.3065629648763770), f91); |
560 f112 = f110 - f109; | 560 f112 = MUL_C(COEF_CONST((-0.3826834323650904)), f110); |
561 f113 = f105 + f103; | 561 f113 = MUL_C(COEF_CONST(0.5411961001461961), f95); |
562 f114 = MUL_C_R(COEF_CONST(1.3065629648763770), f105); | 562 f114 = f111 + f112; |
563 f115 = MUL_C_R(COEF_CONST((-0.3826834323650904)), f113); | 563 f115 = f113 - f112; |
564 f116 = MUL_C_R(COEF_CONST(0.5411961001461961), f103); | 564 f116 = MUL_C(COEF_CONST(0.7071067811865476), f93); |
565 f117 = f114 + f115; | 565 f117 = f116 - f48; |
566 f118 = f116 - f115; | 566 f118 = f116 + f48; |
567 f119 = f100 - f111; | 567 f119 = f115 - f118; |
568 f120 = f100 + f111; | 568 f120 = f115 + f118; |
569 f121 = f102 - f112; | 569 f121 = f114 - f117; |
570 f122 = f102 + f112; | 570 f122 = f114 + f117; |
571 f123 = f99 - f117; | 571 f123 = f120 + f107; |
572 f124 = f99 + f117; | 572 f124 = MUL_C(COEF_CONST((-0.8971675863426361)), f120); |
573 f125 = f101 - f118; | 573 f125 = MUL_C(COEF_CONST(0.9951847266721968), f123); |
574 f126 = f101 + f118; | 574 f126 = MUL_C(COEF_CONST(1.0932018670017576), f107); |
575 f127 = f122 + f120; | 575 y[4] = f124 + f125; |
576 f128 = MUL_C_R(COEF_CONST((-0.8971675863426361)), f122); | 576 y[60] = f126 - f125; |
577 f129 = MUL_C_R(COEF_CONST(0.9951847266721968), f127); | 577 f129 = f122 + f109; |
578 f130 = MUL_C_R(COEF_CONST(1.0932018670017576), f120); | 578 f130 = MUL_C(COEF_CONST((-0.6666556584777466)), f122); |
579 y[4] = f128 + f129; | 579 f131 = MUL_C(COEF_CONST(0.9569403357322089), f129); |
580 y[60] = f130 - f129; | 580 f132 = MUL_C(COEF_CONST(1.2472250129866713), f109); |
581 f133 = f126 + f124; | 581 y[12] = f130 + f131; |
582 f134 = MUL_C_R(COEF_CONST((-0.4105245275223571)), f126); | 582 y[52] = f132 - f131; |
583 f135 = MUL_C_R(COEF_CONST(0.8819212643483549), f133); | 583 f135 = f121 + f108; |
584 f136 = MUL_C_R(COEF_CONST(1.3533180011743529), f124); | 584 f136 = MUL_C(COEF_CONST((-0.4105245275223571)), f121); |
585 y[20] = f134 + f135; | 585 f137 = MUL_C(COEF_CONST(0.8819212643483549), f135); |
586 y[44] = f136 - f135; | 586 f138 = MUL_C(COEF_CONST(1.3533180011743529), f108); |
587 f139 = f121 + f119; | 587 y[20] = f136 + f137; |
588 f140 = MUL_C_R(COEF_CONST(0.1386171691990915), f121); | 588 y[44] = f138 - f137; |
589 f141 = MUL_C_R(COEF_CONST(0.6343932841636455), f139); | 589 f141 = f119 + f106; |
590 f142 = MUL_C_R(COEF_CONST(1.4074037375263826), f119); | 590 f142 = MUL_C(COEF_CONST((-0.1386171691990915)), f119); |
591 y[36] = f140 + f141; | 591 f143 = MUL_C(COEF_CONST(0.7730104533627370), f141); |
592 y[28] = f142 - f141; | 592 f144 = MUL_C(COEF_CONST(1.4074037375263826), f106); |
593 f145 = f125 + f123; | 593 y[28] = f142 + f143; |
594 f146 = MUL_C_R(COEF_CONST(0.6666556584777466), f125); | 594 y[36] = f144 - f143; |
595 f147 = MUL_C_R(COEF_CONST(0.2902846772544623), f145); | 595 f147 = f16 - f18; |
596 f148 = MUL_C_R(COEF_CONST(1.2472250129866711), f123); | 596 f148 = f16 + f18; |
597 y[52] = f146 + f147; | 597 f149 = MUL_C(COEF_CONST(0.7071067811865476), f148); |
598 y[12] = f148 - f147; | 598 f150 = MUL_C(COEF_CONST(0.7071067811865476), f147); |
599 f151 = f2 + f32; | 599 f151 = f10 - f24; |
600 f152 = MUL_C_R(COEF_CONST(1.0478631305325901), f2); | 600 f152 = f10 + f24; |
601 f153 = MUL_C_R(COEF_CONST((-0.9987954562051724)), f151); | 601 f153 = MUL_C(COEF_CONST(0.7071067811865476), f152); |
602 f154 = MUL_C_R(COEF_CONST((-0.9497277818777548)), f32); | 602 f154 = MUL_C(COEF_CONST(0.7071067811865476), f151); |
603 f155 = f152 + f153; | 603 f155 = f14 - f20; |
604 f156 = f154 - f153; | 604 f156 = f14 + f20; |
605 f157 = f4 + f30; | 605 f157 = MUL_C(COEF_CONST(0.7071067811865476), f156); |
606 f158 = MUL_C_R(COEF_CONST(1.1359069844201428), f4); | 606 f158 = MUL_C(COEF_CONST(0.7071067811865476), f155); |
607 f159 = MUL_C_R(COEF_CONST((-0.9891765099647809)), f157); | 607 f159 = f12 - f22; |
608 f160 = MUL_C_R(COEF_CONST((-0.8424460355094190)), f30); | 608 f160 = f12 + f22; |
609 f161 = f158 + f159; | 609 f161 = MUL_C(COEF_CONST(0.7071067811865476), f160); |
610 f162 = f160 - f159; | 610 f162 = MUL_C(COEF_CONST(0.7071067811865476), f159); |
611 f163 = f6 + f28; | 611 f163 = f2 - f149; |
612 f164 = MUL_C_R(COEF_CONST(1.2130114330978077), f6); | 612 f164 = f2 + f149; |
613 f165 = MUL_C_R(COEF_CONST((-0.9700312531945440)), f163); | 613 f165 = f32 - f150; |
614 f166 = MUL_C_R(COEF_CONST((-0.7270510732912803)), f28); | 614 f166 = f32 + f150; |
615 f167 = f164 + f165; | 615 f167 = f8 - f153; |
616 f168 = f166 - f165; | 616 f168 = f8 + f153; |
617 f169 = f8 + f26; | 617 f169 = f26 - f154; |
618 f170 = MUL_C_R(COEF_CONST(1.2784339185752405), f8); | 618 f170 = f26 + f154; |
619 f171 = MUL_C_R(COEF_CONST((-0.9415440651830209)), f169); | 619 f171 = f4 - f157; |
620 f172 = MUL_C_R(COEF_CONST((-0.6046542117908014)), f26); | 620 f172 = f4 + f157; |
621 f173 = f170 + f171; | 621 f173 = f30 - f158; |
622 f174 = f172 - f171; | 622 f174 = f30 + f158; |
623 f175 = f10 + f24; | 623 f175 = f6 - f161; |
624 f176 = MUL_C_R(COEF_CONST(1.3315443865537255), f10); | 624 f176 = f6 + f161; |
625 f177 = MUL_C_R(COEF_CONST((-0.9039892931234433)), f175); | 625 f177 = f28 - f162; |
626 f178 = MUL_C_R(COEF_CONST((-0.4764341996931612)), f24); | 626 f178 = f28 + f162; |
627 f179 = f176 + f177; | 627 f179 = f170 + f168; |
628 f180 = f178 - f177; | 628 f180 = MUL_C(COEF_CONST((-0.5411961001461969)), f170); |
629 f181 = f12 + f22; | 629 f181 = MUL_C(COEF_CONST(0.9238795325112867), f179); |
630 f182 = MUL_C_R(COEF_CONST(1.3718313541934939), f12); | 630 f182 = MUL_C(COEF_CONST(1.3065629648763766), f168); |
631 f183 = MUL_C_R(COEF_CONST((-0.8577286100002722)), f181); | 631 f183 = f180 + f181; |
632 f184 = MUL_C_R(COEF_CONST((-0.3436258658070507)), f22); | 632 f184 = f182 - f181; |
633 f185 = f182 + f183; | 633 f185 = f169 + f167; |
634 f186 = f184 - f183; | 634 f186 = MUL_C(COEF_CONST(1.3065629648763770), f169); |
635 f187 = f14 + f20; | 635 f187 = MUL_C(COEF_CONST((-0.3826834323650904)), f185); |
636 f188 = MUL_C_R(COEF_CONST(1.3989068359730781), f14); | 636 f188 = MUL_C(COEF_CONST(0.5411961001461961), f167); |
637 f189 = MUL_C_R(COEF_CONST((-0.8032075314806453)), f187); | 637 f189 = f186 + f187; |
638 f190 = MUL_C_R(COEF_CONST((-0.2075082269882124)), f20); | 638 f190 = f188 - f187; |
639 f191 = f188 + f189; | 639 f191 = f178 + f176; |
640 f192 = f190 - f189; | 640 f192 = MUL_C(COEF_CONST((-0.5411961001461969)), f178); |
641 f193 = f16 + f18; | 641 f193 = MUL_C(COEF_CONST(0.9238795325112867), f191); |
642 f194 = MUL_C_R(COEF_CONST(1.4125100802019774), f16); | 642 f194 = MUL_C(COEF_CONST(1.3065629648763766), f176); |
643 f195 = MUL_C_R(COEF_CONST((-0.7409511253549591)), f193); | 643 f195 = f192 + f193; |
644 f196 = MUL_C_R(COEF_CONST((-0.0693921705079408)), f18); | 644 f196 = f194 - f193; |
645 f197 = f194 + f195; | 645 f197 = f177 + f175; |
646 f198 = f196 - f195; | 646 f198 = MUL_C(COEF_CONST(1.3065629648763770), f177); |
647 f199 = f156 - f198; | 647 f199 = MUL_C(COEF_CONST((-0.3826834323650904)), f197); |
648 f200 = f156 + f198; | 648 f200 = MUL_C(COEF_CONST(0.5411961001461961), f175); |
649 f201 = f162 - f192; | 649 f201 = f198 + f199; |
650 f202 = f162 + f192; | 650 f202 = f200 - f199; |
651 f203 = f168 - f186; | 651 f203 = f164 - f183; |
652 f204 = f168 + f186; | 652 f204 = f164 + f183; |
653 f205 = f174 - f180; | 653 f205 = f166 - f184; |
654 f206 = f174 + f180; | 654 f206 = f166 + f184; |
655 f207 = f200 - f206; | 655 f207 = f163 - f189; |
656 f208 = f200 + f206; | 656 f208 = f163 + f189; |
657 f209 = f202 - f204; | 657 f209 = f165 - f190; |
658 f210 = f202 + f204; | 658 f210 = f165 + f190; |
659 f211 = f208 - f210; | 659 f211 = f172 - f195; |
660 y[2] = f208 + f210; | 660 f212 = f172 + f195; |
661 f213 = MUL_C_R(COEF_CONST(0.7071067811865476), f211); | 661 f213 = f174 - f196; |
662 f214 = f207 + f209; | 662 f214 = f174 + f196; |
663 f215 = MUL_C_R(COEF_CONST(1.3065629648763766), f207); | 663 f215 = f171 - f201; |
664 f216 = MUL_C_R(COEF_CONST((-0.9238795325112866)), f214); | 664 f216 = f171 + f201; |
665 f217 = MUL_C_R(COEF_CONST((-0.5411961001461967)), f209); | 665 f217 = f173 - f202; |
666 f218 = f215 + f216; | 666 f218 = f173 + f202; |
667 f219 = f217 - f216; | 667 f219 = f214 + f212; |
668 f220 = f201 - f203; | 668 f220 = MUL_C(COEF_CONST((-0.7856949583871021)), f214); |
669 f221 = f201 + f203; | 669 f221 = MUL_C(COEF_CONST(0.9807852804032304), f219); |
670 f222 = MUL_C_R(COEF_CONST(0.7071067811865476), f221); | 670 f222 = MUL_C(COEF_CONST(1.1758756024193588), f212); |
671 f223 = MUL_C_R(COEF_CONST(0.7071067811865476), f220); | 671 f223 = f220 + f221; |
672 f224 = f199 - f222; | 672 f224 = f222 - f221; |
673 f225 = f199 + f222; | 673 f225 = f218 + f216; |
674 f226 = f205 - f223; | 674 f226 = MUL_C(COEF_CONST(0.2758993792829431), f218); |
675 f227 = f205 + f223; | 675 f227 = MUL_C(COEF_CONST(0.5555702330196022), f225); |
676 f228 = f227 + f225; | 676 f228 = MUL_C(COEF_CONST(1.3870398453221475), f216); |
677 f229 = MUL_C_R(COEF_CONST((-0.7856949583871021)), f227); | 677 f229 = f226 + f227; |
678 f230 = MUL_C_R(COEF_CONST(0.9807852804032304), f228); | 678 f230 = f228 - f227; |
679 f231 = MUL_C_R(COEF_CONST(1.1758756024193588), f225); | 679 f231 = f213 + f211; |
680 f232 = f229 + f230; | 680 f232 = MUL_C(COEF_CONST(1.1758756024193591), f213); |
681 f233 = f231 - f230; | 681 f233 = MUL_C(COEF_CONST((-0.1950903220161287)), f231); |
682 f234 = f226 + f224; | 682 f234 = MUL_C(COEF_CONST(0.7856949583871016), f211); |
683 f235 = MUL_C_R(COEF_CONST(0.2758993792829431), f226); | 683 f235 = f232 + f233; |
684 f236 = MUL_C_R(COEF_CONST(0.5555702330196022), f234); | 684 f236 = f234 - f233; |
685 f237 = MUL_C_R(COEF_CONST(1.3870398453221475), f224); | 685 f237 = f217 + f215; |
686 f238 = f235 + f236; | 686 f238 = MUL_C(COEF_CONST(1.3870398453221473), f217); |
687 f239 = f237 - f236; | 687 f239 = MUL_C(COEF_CONST((-0.8314696123025455)), f237); |
688 f240 = f155 - f197; | 688 f240 = MUL_C(COEF_CONST((-0.2758993792829436)), f215); |
689 f241 = f155 + f197; | 689 f241 = f238 + f239; |
690 f242 = f161 - f191; | 690 f242 = f240 - f239; |
691 f243 = f161 + f191; | 691 f243 = f204 - f223; |
692 f244 = f167 - f185; | 692 f244 = f204 + f223; |
693 f245 = f167 + f185; | 693 f245 = f206 - f224; |
694 f246 = f173 - f179; | 694 f246 = f206 + f224; |
695 f247 = f173 + f179; | 695 f247 = f208 - f229; |
696 f248 = f245 - f243; | 696 f248 = f208 + f229; |
697 f249 = f245 + f243; | 697 f249 = f210 - f230; |
698 f250 = MUL_C_R(COEF_CONST(0.7071067811865476), f249); | 698 f250 = f210 + f230; |
699 f251 = f247 - f250; | 699 f251 = f203 - f235; |
700 f252 = f247 + f250; | 700 f252 = f203 + f235; |
701 f253 = MUL_C_R(COEF_CONST(0.7071067811865476), f248); | 701 f253 = f205 - f236; |
702 f254 = f253 - f241; | 702 f254 = f205 + f236; |
703 f255 = f253 + f241; | 703 f255 = f207 - f241; |
704 f256 = f255 + f252; | 704 f256 = f207 + f241; |
705 f257 = MUL_C_R(COEF_CONST((-0.7856949583871021)), f255); | 705 f257 = f209 - f242; |
706 f258 = MUL_C_R(COEF_CONST(0.9807852804032304), f256); | 706 f258 = f209 + f242; |
707 f259 = MUL_C_R(COEF_CONST(1.1758756024193588), f252); | 707 f259 = f246 + f244; |
708 f260 = f257 + f258; | 708 f260 = MUL_C(COEF_CONST((-0.9497277818777543)), f246); |
709 f261 = f259 - f258; | 709 f261 = MUL_C(COEF_CONST(0.9987954562051724), f259); |
710 f262 = f254 + f251; | 710 f262 = MUL_C(COEF_CONST(1.0478631305325905), f244); |
711 f263 = MUL_C_R(COEF_CONST((-0.2758993792829430)), f254); | 711 y[2] = f260 + f261; |
712 f264 = MUL_C_R(COEF_CONST(0.8314696123025452), f262); | 712 y[62] = f262 - f261; |
713 f265 = MUL_C_R(COEF_CONST(1.3870398453221475), f251); | 713 f265 = f250 + f248; |
714 f266 = f263 + f264; | 714 f266 = MUL_C(COEF_CONST((-0.7270510732912801)), f250); |
715 f267 = f265 - f264; | 715 f267 = MUL_C(COEF_CONST(0.9700312531945440), f265); |
716 f268 = f240 - f246; | 716 f268 = MUL_C(COEF_CONST(1.2130114330978079), f248); |
717 f269 = f240 + f246; | 717 y[10] = f266 + f267; |
718 f270 = f242 - f244; | 718 y[54] = f268 - f267; |
719 f271 = f242 + f244; | 719 f271 = f254 + f252; |
720 f272 = f269 + f271; | 720 f272 = MUL_C(COEF_CONST((-0.4764341996931611)), f254); |
721 f273 = MUL_C_R(COEF_CONST(1.3065629648763770), f269); | 721 f273 = MUL_C(COEF_CONST(0.9039892931234433), f271); |
722 f274 = MUL_C_R(COEF_CONST((-0.3826834323650904)), f272); | 722 f274 = MUL_C(COEF_CONST(1.3315443865537255), f252); |
723 f275 = MUL_C_R(COEF_CONST(0.5411961001461961), f271); | 723 y[18] = f272 + f273; |
724 f276 = f273 + f274; | 724 y[46] = f274 - f273; |
725 f277 = f275 - f274; | 725 f277 = f258 + f256; |
726 y[62] = f268 - f270; | 726 f278 = MUL_C(COEF_CONST((-0.2075082269882114)), f258); |
727 f279 = f268 + f270; | 727 f279 = MUL_C(COEF_CONST(0.8032075314806448), f277); |
728 f280 = MUL_C_R(COEF_CONST(0.7071067811865476), f279); | 728 f280 = MUL_C(COEF_CONST(1.3989068359730783), f256); |
729 y[10] = f232 - f260; | 729 y[26] = f278 + f279; |
730 y[6] = f232 + f260; | 730 y[38] = f280 - f279; |
731 y[18] = f219 - f277; | 731 f283 = f245 + f243; |
732 y[14] = f219 + f277; | 732 f284 = MUL_C(COEF_CONST(0.0693921705079408), f245); |
733 y[26] = f239 + f266; | 733 f285 = MUL_C(COEF_CONST(0.6715589548470183), f283); |
734 y[22] = f239 - f266; | 734 f286 = MUL_C(COEF_CONST(1.4125100802019774), f243); |
735 y[34] = f213 - f280; | 735 y[34] = f284 + f285; |
736 y[30] = f213 + f280; | 736 y[30] = f286 - f285; |
737 y[42] = f238 - f267; | 737 f289 = f249 + f247; |
738 y[38] = f238 + f267; | 738 f290 = MUL_C(COEF_CONST(0.3436258658070505), f249); |
739 y[50] = f218 - f276; | 739 f291 = MUL_C(COEF_CONST(0.5141027441932217), f289); |
740 y[46] = f218 + f276; | 740 f292 = MUL_C(COEF_CONST(1.3718313541934939), f247); |
741 y[58] = f233 + f261; | 741 y[42] = f290 + f291; |
742 y[54] = f233 - f261; | 742 y[22] = f292 - f291; |
743 f295 = t2[3] - t2[5]; | 743 f295 = f253 + f251; |
744 f296 = t2[3] + t2[5]; | 744 f296 = MUL_C(COEF_CONST(0.6046542117908007), f253); |
745 f297 = t2[7] - t2[9]; | 745 f297 = MUL_C(COEF_CONST(0.3368898533922201), f295); |
746 f298 = t2[7] + t2[9]; | 746 f298 = MUL_C(COEF_CONST(1.2784339185752409), f251); |
747 f299 = t2[11] - t2[13]; | 747 y[50] = f296 + f297; |
748 f300 = t2[11] + t2[13]; | 748 y[14] = f298 - f297; |
749 f301 = t2[15] - t2[17]; | 749 f301 = f257 + f255; |
750 f302 = t2[15] + t2[17]; | 750 f302 = MUL_C(COEF_CONST(0.8424460355094192), f257); |
751 f303 = t2[19] - t2[21]; | 751 f303 = MUL_C(COEF_CONST(0.1467304744553618), f301); |
752 f304 = t2[19] + t2[21]; | 752 f304 = MUL_C(COEF_CONST(1.1359069844201428), f255); |
753 f305 = t2[23] - t2[25]; | 753 y[58] = f302 + f303; |
754 f306 = t2[23] + t2[25]; | 754 y[6] = f304 - f303; |
755 f307 = t2[27] - t2[29]; | 755 f307 = t2[1] + t2[63]; |
756 f308 = t2[27] + t2[29]; | 756 f308 = MUL_C(COEF_CONST(1.0242400472191162), t2[1]); |
757 f309 = t2[31] - t2[33]; | 757 f309 = MUL_C(COEF_CONST((-0.9996988186962043)), f307); |
758 f310 = t2[31] + t2[33]; | 758 f310 = MUL_C(COEF_CONST((-0.9751575901732922)), t2[63]); |
759 f311 = t2[35] - t2[37]; | 759 f311 = f308 + f309; |
760 f312 = t2[35] + t2[37]; | 760 f312 = f310 - f309; |
761 f313 = t2[39] - t2[41]; | 761 f313 = t2[3] + t2[61]; |
762 f314 = t2[39] + t2[41]; | 762 f314 = MUL_C(COEF_CONST(1.0708550202783571),t2[3]); |
763 f315 = t2[43] - t2[45]; | 763 f315 = MUL_C(COEF_CONST((-0.9972904566786902)), f313); |
764 f316 = t2[43] + t2[45]; | 764 f316 = MUL_C(COEF_CONST((-0.9237258930790232)), t2[61]); |
765 f317 = t2[47] - t2[49]; | 765 f317 = f314 + f315; |
766 f318 = t2[47] + t2[49]; | 766 f318 = f316 - f315; |
767 f319 = t2[51] - t2[53]; | 767 f319 = t2[5] + t2[59]; |
768 f320 = t2[51] + t2[53]; | 768 f320 = MUL_C(COEF_CONST(1.1148902097979256), t2[5]); |
769 f321 = t2[55] - t2[57]; | 769 f321 = MUL_C(COEF_CONST((-0.9924795345987101)), f319); |
770 f322 = t2[55] + t2[57]; | 770 f322 = MUL_C(COEF_CONST((-0.8700688593994945)), t2[59]); |
771 f323 = t2[59] - t2[61]; | 771 f323 = f320 + f321; |
772 f324 = t2[59] + t2[61]; | 772 f324 = f322 - f321; |
773 f325 = MUL_C_R(COEF_CONST(0.7071067811865476), f310); | 773 f325 = t2[7] + t2[57]; |
774 f326 = t2[1] - f325; | 774 f326 = MUL_C(COEF_CONST(1.1562395311492426), t2[7]); |
775 f327 = t2[1] + f325; | 775 f327 = MUL_C(COEF_CONST((-0.9852776423889412)), f325); |
776 f328 = f302 + f318; | 776 f328 = MUL_C(COEF_CONST((-0.8143157536286398)), t2[57]); |
777 f329 = MUL_C_R(COEF_CONST(1.3065629648763766), f302); | 777 f329 = f326 + f327; |
778 f330 = MUL_C_R(COEF_CONST((-0.9238795325112866)), f328); | 778 f330 = f328 - f327; |
779 f331 = MUL_C_R(COEF_CONST((-0.5411961001461967)), f318); | 779 f331 = t2[9] + t2[55]; |
780 f332 = f329 + f330; | 780 f332 = MUL_C(COEF_CONST(1.1948033701953984), t2[9]); |
781 f333 = f331 - f330; | 781 f333 = MUL_C(COEF_CONST((-0.9757021300385286)), f331); |
782 f334 = f327 - f333; | 782 f334 = MUL_C(COEF_CONST((-0.7566008898816589)), t2[55]); |
783 f335 = f327 + f333; | 783 f335 = f332 + f333; |
784 f336 = f326 - f332; | 784 f336 = f334 - f333; |
785 f337 = f326 + f332; | 785 f337 = t2[11] + t2[53]; |
786 f338 = f306 - f314; | 786 f338 = MUL_C(COEF_CONST(1.2304888232703382), t2[11]); |
787 f339 = f306 + f314; | 787 f339 = MUL_C(COEF_CONST((-0.9637760657954400)), f337); |
788 f340 = MUL_C_R(COEF_CONST(0.7071067811865476), f339); | 788 f340 = MUL_C(COEF_CONST((-0.6970633083205418)), t2[53]); |
789 f341 = f298 - f340; | 789 f341 = f338 + f339; |
790 f342 = f298 + f340; | 790 f342 = f340 - f339; |
791 f343 = MUL_C_R(COEF_CONST(0.7071067811865476), f338); | 791 f343 = t2[13] + t2[51]; |
792 f344 = f343 - f322; | 792 f344 = MUL_C(COEF_CONST(1.2632099209919279), t2[13]); |
793 f345 = f343 + f322; | 793 f345 = MUL_C(COEF_CONST((-0.9495281805930368)), f343); |
794 f346 = f345 + f342; | 794 f346 = MUL_C(COEF_CONST((-0.6358464401941457)), t2[51]); |
795 f347 = MUL_C_R(COEF_CONST((-0.7856949583871021)), f345); | 795 f347 = f344 + f345; |
796 f348 = MUL_C_R(COEF_CONST(0.9807852804032304), f346); | 796 f348 = f346 - f345; |
797 f349 = MUL_C_R(COEF_CONST(1.1758756024193588), f342); | 797 f349 = t2[15] + t2[49]; |
798 f350 = f347 + f348; | 798 f350 = MUL_C(COEF_CONST(1.2928878353697266), t2[15]); |
799 f351 = f349 - f348; | 799 f351 = MUL_C(COEF_CONST((-0.9329927988347391)), f349); |
800 f352 = f344 + f341; | 800 f352 = MUL_C(COEF_CONST((-0.5730977622997515)), t2[49]); |
801 f353 = MUL_C_R(COEF_CONST((-0.2758993792829430)), f344); | 801 f353 = f350 + f351; |
802 f354 = MUL_C_R(COEF_CONST(0.8314696123025452), f352); | 802 f354 = f352 - f351; |
803 f355 = MUL_C_R(COEF_CONST(1.3870398453221475), f341); | 803 f355 = t2[17] + t2[47]; |
804 f356 = f353 + f354; | 804 f356 = MUL_C(COEF_CONST(1.3194510697085207), t2[17]); |
805 f357 = f355 - f354; | 805 f357 = MUL_C(COEF_CONST((-0.9142097557035306)), f355); |
806 f358 = f335 - f350; | 806 f358 = MUL_C(COEF_CONST((-0.5089684416985405)), t2[47]); |
807 f359 = f335 + f350; | 807 f359 = f356 + f357; |
808 f360 = f337 - f356; | 808 f360 = f358 - f357; |
809 f361 = f337 + f356; | 809 f361 = t2[19] + t2[45]; |
810 f362 = f336 - f357; | 810 f362 = MUL_C(COEF_CONST(1.3428356308501219), t2[19]); |
811 f363 = f336 + f357; | 811 f363 = MUL_C(COEF_CONST((-0.8932243011955153)), f361); |
812 f364 = f334 - f351; | 812 f364 = MUL_C(COEF_CONST((-0.4436129715409087)), t2[45]); |
813 f365 = f334 + f351; | 813 f365 = f362 + f363; |
814 f366 = MUL_C_R(COEF_CONST(5.1011486186891641), f296); | 814 f366 = f364 - f363; |
815 f367 = MUL_C_R(COEF_CONST(1.7224470982383342), f300); | 815 f367 = t2[21] + t2[43]; |
816 f368 = MUL_C_R(COEF_CONST(1.0606776859903475), f304); | 816 f368 = MUL_C(COEF_CONST(1.3629851833384954), t2[21]); |
817 f369 = MUL_C_R(COEF_CONST(0.7881546234512502), f308); | 817 f369 = MUL_C(COEF_CONST((-0.8700869911087115)), f367); |
818 f370 = MUL_C_R(COEF_CONST(0.5024192861881557), f324); | 818 f370 = MUL_C(COEF_CONST((-0.3771887988789276)), t2[43]); |
819 f371 = MUL_C_R(COEF_CONST(0.5224986149396889), f320); | 819 f371 = f368 + f369; |
820 f372 = MUL_C_R(COEF_CONST(0.5669440348163577), f316); | 820 f372 = f370 - f369; |
821 f373 = MUL_C_R(COEF_CONST(0.6468217833599901), f312); | 821 f373 = t2[23] + t2[41]; |
822 f374 = f366 - f370; | 822 f374 = MUL_C(COEF_CONST(1.3798511851368040), t2[23]); |
823 f375 = f366 + f370; | 823 f375 = MUL_C(COEF_CONST((-0.8448535652497072)), f373); |
824 f376 = f367 - f371; | 824 f376 = MUL_C(COEF_CONST((-0.3098559453626103)), t2[41]); |
825 f377 = f367 + f371; | 825 f377 = f374 + f375; |
826 f378 = f368 - f372; | 826 f378 = f376 - f375; |
827 f379 = f368 + f372; | 827 f379 = t2[25] + t2[39]; |
828 f380 = f369 - f373; | 828 f380 = MUL_C(COEF_CONST(1.3933930045694289), t2[25]); |
829 f381 = f369 + f373; | 829 f381 = MUL_C(COEF_CONST((-0.8175848131515840)), f379); |
830 f382 = MUL_C_R(COEF_CONST(0.5097955791041592), f375); | 830 f382 = MUL_C(COEF_CONST((-0.2417766217337392)), t2[39]); |
831 f383 = MUL_C_R(COEF_CONST(0.6013448869350453), f377); | 831 f383 = f380 + f381; |
832 f384 = MUL_C_R(COEF_CONST(0.8999762231364156), f379); | 832 f384 = f382 - f381; |
833 f385 = MUL_C_R(COEF_CONST(2.5629154477415055), f381); | 833 f385 = t2[27] + t2[37]; |
834 f386 = f382 + f385; | 834 f386 = MUL_C(COEF_CONST(1.4035780182072330), t2[27]); |
835 f387 = f382 - f385; | 835 f387 = MUL_C(COEF_CONST((-0.7883464276266061)), f385); |
836 f388 = f383 + f384; | 836 f388 = MUL_C(COEF_CONST((-0.1731148370459791)), t2[37]); |
837 f389 = f384 - f383; | 837 f389 = f386 + f387; |
838 f390 = f387 - f389; | 838 f390 = f388 - f387; |
839 f391 = f387 + f389; | 839 f391 = t2[29] + t2[35]; |
840 f392 = MUL_C_R(COEF_CONST(0.7071067811865476), f390); | 840 f392 = MUL_C(COEF_CONST(1.4103816894602614), t2[29]); |
841 f393 = f386 - f388; | 841 f393 = MUL_C(COEF_CONST((-0.7572088465064846)), f391); |
842 f394 = MUL_C_R(COEF_CONST(1.3065629648763766), f386); | 842 f394 = MUL_C(COEF_CONST((-0.1040360035527078)), t2[35]); |
843 f395 = MUL_C_R(COEF_CONST((-0.9238795325112866)), f393); | 843 f395 = f392 + f393; |
844 f396 = MUL_C_R(COEF_CONST((-0.5411961001461967)), f388); | 844 f396 = f394 - f393; |
845 f397 = f394 + f395; | 845 f397 = t2[31] + t2[33]; |
846 f398 = f395 + f396; | 846 f398 = MUL_C(COEF_CONST(1.4137876276885337), t2[31]); |
847 f399 = f391 - f398; | 847 f399 = MUL_C(COEF_CONST((-0.7242470829514670)), f397); |
848 f400 = f392 - f398; | 848 f400 = MUL_C(COEF_CONST((-0.0347065382144002)), t2[33]); |
849 f401 = f392 + f397; | 849 f401 = f398 + f399; |
850 f402 = f380 - f374; | 850 f402 = f400 - f399; |
851 f403 = f374 + f380; | 851 f403 = f312 - f402; |
852 f404 = f378 - f376; | 852 f404 = f312 + f402; |
853 f405 = f376 + f378; | 853 f405 = f318 - f396; |
854 f406 = f403 + f405; | 854 f406 = f318 + f396; |
855 f407 = MUL_C_R(COEF_CONST(1.3065629648763770), f403); | 855 f407 = f324 - f390; |
856 f408 = MUL_C_R(COEF_CONST((-0.3826834323650904)), f406); | 856 f408 = f324 + f390; |
857 f409 = MUL_C_R(COEF_CONST(0.5411961001461961), f405); | 857 f409 = f330 - f384; |
858 f410 = f407 + f408; | 858 f410 = f330 + f384; |
859 f411 = f408 - f409; | 859 f411 = f336 - f378; |
860 f412 = f402 - f404; | 860 f412 = f336 + f378; |
861 f413 = f402 + f404; | 861 f413 = f342 - f372; |
862 f414 = MUL_C_R(COEF_CONST(0.7071067811865476), f413); | 862 f414 = f342 + f372; |
863 f415 = f411 + f397; | 863 f415 = f348 - f366; |
864 f416 = f401 + f411; | 864 f416 = f348 + f366; |
865 f417 = f414 + f401; | 865 f417 = f354 - f360; |
866 f418 = f400 + f414; | 866 f418 = f354 + f360; |
867 f419 = f400 - f410; | 867 f419 = f404 - f418; |
868 f420 = f399 - f410; | 868 f420 = f404 + f418; |
869 f421 = f412 + f399; | 869 f421 = f406 - f416; |
870 f422 = f359 - f397; | 870 f422 = f406 + f416; |
871 f423 = f359 + f397; | 871 f423 = f408 - f414; |
872 f424 = f361 + f415; | 872 f424 = f408 + f414; |
873 f425 = f361 - f415; | 873 f425 = f410 - f412; |
874 f426 = f363 - f416; | 874 f426 = f410 + f412; |
875 f427 = f363 + f416; | 875 f427 = f420 - f426; |
876 f428 = f365 + f417; | 876 f428 = f420 + f426; |
877 f429 = f365 - f417; | 877 f429 = f422 - f424; |
878 f430 = f364 - f418; | 878 f430 = f422 + f424; |
879 f431 = f364 + f418; | 879 f431 = f428 - f430; |
880 f432 = f362 + f419; | 880 y[1] = f428 + f430; |
881 f433 = f362 - f419; | 881 f433 = MUL_C(COEF_CONST(0.7071067811865476), f431); |
882 f434 = f360 - f420; | 882 f434 = f427 + f429; |
883 f435 = f360 + f420; | 883 f435 = MUL_C(COEF_CONST(1.3065629648763766), f427); |
884 f436 = f358 + f421; | 884 f436 = MUL_C(COEF_CONST((-0.9238795325112866)), f434); |
885 f437 = f358 - f421; | 885 f437 = MUL_C(COEF_CONST((-0.5411961001461967)), f429); |
886 f438 = MUL_C_R(COEF_CONST(5.1011486186891641), f295); | 886 f438 = f435 + f436; |
887 f439 = MUL_C_R(COEF_CONST(1.7224470982383342), f299); | 887 f439 = f437 - f436; |
888 f440 = MUL_C_R(COEF_CONST(1.0606776859903475), f303); | 888 f440 = f419 + f425; |
889 f441 = MUL_C_R(COEF_CONST(0.7881546234512502), f307); | 889 f441 = MUL_C(COEF_CONST(1.1758756024193588), f419); |
890 f442 = MUL_C_R(COEF_CONST(0.5024192861881557), f323); | 890 f442 = MUL_C(COEF_CONST((-0.9807852804032304)), f440); |
891 f443 = MUL_C_R(COEF_CONST(0.5224986149396889), f319); | 891 f443 = MUL_C(COEF_CONST((-0.7856949583871021)), f425); |
892 f444 = MUL_C_R(COEF_CONST(0.5669440348163577), f315); | 892 f444 = f441 + f442; |
893 f445 = MUL_C_R(COEF_CONST(0.6468217833599901), f311); | 893 f445 = f443 - f442; |
894 f446 = f438 + f442; | 894 f446 = f421 + f423; |
895 f447 = f438 - f442; | 895 f447 = MUL_C(COEF_CONST(1.3870398453221473), f421); |
896 f448 = f439 + f443; | 896 f448 = MUL_C(COEF_CONST((-0.8314696123025455)), f446); |
897 f449 = f443 - f439; | 897 f449 = MUL_C(COEF_CONST((-0.2758993792829436)), f423); |
898 f450 = f440 + f444; | 898 f450 = f447 + f448; |
899 f451 = f440 - f444; | 899 f451 = f449 - f448; |
900 f452 = f441 + f445; | 900 f452 = f445 - f451; |
901 f453 = f445 - f441; | 901 f453 = f445 + f451; |
902 f454 = MUL_C_R(COEF_CONST(0.5097955791041592), f447); | 902 f454 = MUL_C(COEF_CONST(0.7071067811865476), f452); |
903 f455 = MUL_C_R(COEF_CONST(0.6013448869350453), f449); | 903 f455 = f444 - f450; |
904 f456 = MUL_C_R(COEF_CONST(0.8999762231364156), f451); | 904 f456 = f444 + f450; |
905 f457 = MUL_C_R(COEF_CONST(2.5629154477415055), f453); | 905 f457 = MUL_C(COEF_CONST(0.7071067811865476), f456); |
906 f458 = f454 + f457; | 906 f458 = f454 - f457; |
907 f459 = f454 - f457; | 907 f459 = f454 + f457; |
908 f460 = f455 + f456; | 908 f460 = f405 - f407; |
909 f461 = f456 - f455; | 909 f461 = f405 + f407; |
910 f462 = f459 - f461; | 910 f462 = f409 - f411; |
911 f463 = f459 + f461; | 911 f463 = f409 + f411; |
912 f464 = MUL_C_R(COEF_CONST(0.7071067811865476), f462); | 912 f464 = f413 - f415; |
913 f465 = f458 - f460; | 913 f465 = f413 + f415; |
914 f466 = MUL_C_R(COEF_CONST(1.3065629648763766), f458); | 914 f466 = MUL_C(COEF_CONST(0.7071067811865476), f463); |
915 f467 = MUL_C_R(COEF_CONST((-0.9238795325112866)), f465); | 915 f467 = f403 - f466; |
916 f468 = MUL_C_R(COEF_CONST((-0.5411961001461967)), f460); | 916 f468 = f403 + f466; |
917 f469 = f466 + f467; | 917 f469 = f461 + f465; |
918 f470 = f467 + f468; | 918 f470 = MUL_C(COEF_CONST(1.3065629648763766), f461); |
919 f471 = f463 - f470; | 919 f471 = MUL_C(COEF_CONST((-0.9238795325112866)), f469); |
920 f472 = f464 - f470; | 920 f472 = MUL_C(COEF_CONST((-0.5411961001461967)), f465); |
921 f473 = f464 + f469; | 921 f473 = f470 + f471; |
922 f474 = f446 + f452; | 922 f474 = f472 - f471; |
923 f475 = f452 - f446; | 923 f475 = f468 - f474; |
924 f476 = f448 + f450; | 924 f476 = f468 + f474; |
925 f477 = f448 - f450; | 925 f477 = f467 - f473; |
926 f478 = f475 + f477; | 926 f478 = f467 + f473; |
927 f479 = MUL_C_R(COEF_CONST(1.3065629648763770), f475); | 927 f479 = f460 + f464; |
928 f480 = MUL_C_R(COEF_CONST((-0.3826834323650904)), f478); | 928 f480 = MUL_C(COEF_CONST(1.3065629648763770), f460); |
929 f481 = MUL_C_R(COEF_CONST(0.5411961001461961), f477); | 929 f481 = MUL_C(COEF_CONST((-0.3826834323650904)), f479); |
930 f482 = f479 + f480; | 930 f482 = MUL_C(COEF_CONST(0.5411961001461961), f464); |
931 f483 = f481 - f480; | 931 f483 = f480 + f481; |
932 f484 = f474 + f476; | 932 f484 = f482 - f481; |
933 f485 = f476 - f474; | 933 f485 = MUL_C(COEF_CONST(0.7071067811865476), f462); |
934 f486 = MUL_C_R(COEF_CONST(0.7071067811865476), f485); | 934 f486 = f485 - f417; |
935 f487 = f483 + f469; | 935 f487 = f485 + f417; |
936 f488 = f473 + f483; | 936 f488 = f484 - f487; |
937 f489 = f486 + f473; | 937 f489 = f484 + f487; |
938 f490 = f472 + f486; | 938 f490 = f483 - f486; |
939 f491 = f482 + f472; | 939 f491 = f483 + f486; |
940 f492 = f471 + f482; | 940 f492 = f489 + f476; |
941 f493 = f471 - f484; | 941 f493 = MUL_C(COEF_CONST((-0.8971675863426361)), f489); |
942 f494 = MUL_C_R(COEF_CONST(0.7071067811865476), f309); | 942 f494 = MUL_C(COEF_CONST(0.9951847266721968), f492); |
943 f495 = t2[63] - f494; | 943 f495 = MUL_C(COEF_CONST(1.0932018670017576), f476); |
944 f496 = t2[63] + f494; | 944 f496 = f493 + f494; |
945 f497 = f317 + f301; | 945 f497 = f495 - f494; |
946 f498 = MUL_C_R(COEF_CONST(1.3065629648763766), f317); | 946 f498 = f491 + f478; |
947 f499 = MUL_C_R(COEF_CONST((-0.9238795325112866)), f497); | 947 f499 = MUL_C(COEF_CONST((-0.6666556584777466)), f491); |
948 f500 = MUL_C_R(COEF_CONST((-0.5411961001461967)), f301); | 948 f500 = MUL_C(COEF_CONST(0.9569403357322089), f498); |
949 f501 = f498 + f499; | 949 f501 = MUL_C(COEF_CONST(1.2472250129866713), f478); |
950 f502 = f500 - f499; | 950 f502 = f499 + f500; |
951 f503 = f496 - f502; | 951 f503 = f501 - f500; |
952 f504 = f496 + f502; | 952 f504 = f490 + f477; |
953 f505 = f495 - f501; | 953 f505 = MUL_C(COEF_CONST((-0.4105245275223571)), f490); |
954 f506 = f495 + f501; | 954 f506 = MUL_C(COEF_CONST(0.8819212643483549), f504); |
955 f507 = MUL_C_R(COEF_CONST(0.5097955791041592), f321); | 955 f507 = MUL_C(COEF_CONST(1.3533180011743529), f477); |
956 f508 = MUL_C_R(COEF_CONST(0.6013448869350453), f313); | 956 f508 = f505 + f506; |
957 f509 = MUL_C_R(COEF_CONST(0.8999762231364156), f305); | 957 f509 = f507 - f506; |
958 f510 = MUL_C_R(COEF_CONST(2.5629154477415055), f297); | 958 f510 = f488 + f475; |
959 f511 = f507 - f510; | 959 f511 = MUL_C(COEF_CONST((-0.1386171691990915)), f488); |
960 f512 = f507 + f510; | 960 f512 = MUL_C(COEF_CONST(0.7730104533627370), f510); |
961 f513 = f508 - f509; | 961 f513 = MUL_C(COEF_CONST(1.4074037375263826), f475); |
962 f514 = f508 + f509; | 962 f514 = f511 + f512; |
963 f515 = f512 - f514; | 963 f515 = f513 - f512; |
964 f516 = f512 + f514; | 964 f516 = f311 + f401; |
965 f517 = MUL_C_R(COEF_CONST(0.7071067811865476), f515); | 965 f517 = f311 - f401; |
966 f518 = f511 + f513; | 966 f518 = f317 + f395; |
967 f519 = MUL_C_R(COEF_CONST(1.3065629648763766), f511); | 967 f519 = f395 - f317; |
968 f520 = MUL_C_R(COEF_CONST((-0.9238795325112866)), f518); | 968 f520 = f323 + f389; |
969 f521 = MUL_C_R(COEF_CONST((-0.5411961001461967)), f513); | 969 f521 = f323 - f389; |
970 f522 = f519 + f520; | 970 f522 = f329 + f383; |
971 f523 = f521 - f520; | 971 f523 = f383 - f329; |
972 f524 = f516 + f523; | 972 f524 = f335 + f377; |
973 f525 = f523 + f517; | 973 f525 = f335 - f377; |
974 f526 = f517 + f522; | 974 f526 = f341 + f371; |
975 f527 = f504 - f524; | 975 f527 = f371 - f341; |
976 f528 = f504 + f524; | 976 f528 = f347 + f365; |
977 f529 = f506 - f525; | 977 f529 = f347 - f365; |
978 f530 = f506 + f525; | 978 f530 = f353 + f359; |
979 f531 = f505 - f526; | 979 f531 = f359 - f353; |
980 f532 = f505 + f526; | 980 f532 = f517 - f531; |
981 f533 = f503 - f522; | 981 f533 = f517 + f531; |
982 f534 = f503 + f522; | 982 f534 = f519 - f529; |
983 f535 = f493 + f528; | 983 f535 = f519 + f529; |
984 f536 = f528 - f493; | 984 f536 = f521 - f527; |
985 f537 = f492 + f530; | 985 f537 = f521 + f527; |
986 f538 = f492 - f530; | 986 f538 = f523 - f525; |
987 f539 = f491 + f532; | 987 f539 = f523 + f525; |
988 f540 = f532 - f491; | 988 f540 = f533 - f539; |
989 f541 = f490 + f534; | 989 f541 = f533 + f539; |
990 f542 = f490 - f534; | 990 f542 = f535 - f537; |
991 f543 = f489 + f533; | 991 f543 = f535 + f537; |
992 f544 = f533 - f489; | 992 f544 = f541 - f543; |
993 f545 = f488 + f531; | 993 y[63] = f541 + f543; |
994 f546 = f488 - f531; | 994 f546 = MUL_C(COEF_CONST(0.7071067811865476), f544); |
995 f547 = f487 + f529; | 995 f547 = f540 + f542; |
996 f548 = f529 - f487; | 996 f548 = MUL_C(COEF_CONST(1.3065629648763766), f540); |
997 f549 = f469 + f527; | 997 f549 = MUL_C(COEF_CONST((-0.9238795325112866)), f547); |
998 f550 = f469 - f527; | 998 f550 = MUL_C(COEF_CONST((-0.5411961001461967)), f542); |
999 f551 = f536 + f423; | 999 f551 = f548 + f549; |
1000 f552 = MUL_C_R(COEF_CONST((-0.9751575901732920)), f536); | 1000 f552 = f550 - f549; |
1001 f553 = MUL_C_R(COEF_CONST(0.9996988186962043), f551); | 1001 f553 = f532 + f538; |
1002 f554 = MUL_C_R(COEF_CONST(1.0242400472191164), f423); | 1002 f554 = MUL_C(COEF_CONST(1.1758756024193588), f532); |
1003 y[1] = f552 + f553; | 1003 f555 = MUL_C(COEF_CONST((-0.9807852804032304)), f553); |
1004 y[63] = f554 - f553; | 1004 f556 = MUL_C(COEF_CONST((-0.7856949583871021)), f538); |
1005 f557 = f538 + f425; | 1005 f557 = f554 + f555; |
1006 f558 = MUL_C_R(COEF_CONST((-0.9237258930790228)), f538); | 1006 f558 = f556 - f555; |
1007 f559 = MUL_C_R(COEF_CONST(0.9972904566786902), f557); | 1007 f559 = f534 + f536; |
1008 f560 = MUL_C_R(COEF_CONST(1.0708550202783576), f425); | 1008 f560 = MUL_C(COEF_CONST(1.3870398453221473), f534); |
1009 y[3] = f558 + f559; | 1009 f561 = MUL_C(COEF_CONST((-0.8314696123025455)), f559); |
1010 y[61] = f560 - f559; | 1010 f562 = MUL_C(COEF_CONST((-0.2758993792829436)), f536); |
1011 f563 = f540 + f427; | 1011 f563 = f560 + f561; |
1012 f564 = MUL_C_R(COEF_CONST((-0.8700688593994936)), f540); | 1012 f564 = f562 - f561; |
1013 f565 = MUL_C_R(COEF_CONST(0.9924795345987100), f563); | 1013 f565 = f558 - f564; |
1014 f566 = MUL_C_R(COEF_CONST(1.1148902097979263), f427); | 1014 f566 = f558 + f564; |
1015 y[5] = f564 + f565; | 1015 f567 = MUL_C(COEF_CONST(0.7071067811865476), f565); |
1016 y[59] = f566 - f565; | 1016 f568 = f557 - f563; |
1017 f569 = f542 + f429; | 1017 f569 = f557 + f563; |
1018 f570 = MUL_C_R(COEF_CONST((-0.8143157536286398)), f542); | 1018 f570 = MUL_C(COEF_CONST(0.7071067811865476), f569); |
1019 f571 = MUL_C_R(COEF_CONST(0.9852776423889412), f569); | 1019 f571 = f567 - f570; |
1020 f572 = MUL_C_R(COEF_CONST(1.1562395311492426), f429); | 1020 f572 = f567 + f570; |
1021 y[7] = f570 + f571; | 1021 f573 = MUL_C(COEF_CONST(0.5024192861881557), f516); |
1022 y[57] = f572 - f571; | 1022 f574 = MUL_C(COEF_CONST(0.5224986149396889), f518); |
1023 f575 = f544 + f431; | 1023 f575 = MUL_C(COEF_CONST(0.5669440348163577), f520); |
1024 f576 = MUL_C_R(COEF_CONST((-0.7566008898816587)), f544); | 1024 f576 = MUL_C(COEF_CONST(0.6468217833599901), f522); |
1025 f577 = MUL_C_R(COEF_CONST(0.9757021300385286), f575); | 1025 f577 = MUL_C(COEF_CONST(0.7881546234512502), f524); |
1026 f578 = MUL_C_R(COEF_CONST(1.1948033701953984), f431); | 1026 f578 = MUL_C(COEF_CONST(1.0606776859903471), f526); |
1027 y[9] = f576 + f577; | 1027 f579 = MUL_C(COEF_CONST(1.7224470982383342), f528); |
1028 y[55] = f578 - f577; | 1028 f580 = MUL_C(COEF_CONST(5.1011486186891553), f530); |
1029 f581 = f546 + f433; | 1029 f581 = f573 + f580; |
1030 f582 = MUL_C_R(COEF_CONST((-0.6970633083205414)), f546); | 1030 f582 = f573 - f580; |
1031 f583 = MUL_C_R(COEF_CONST(0.9637760657954398), f581); | 1031 f583 = f574 + f579; |
1032 f584 = MUL_C_R(COEF_CONST(1.2304888232703384), f433); | 1032 f584 = f579 - f574; |
1033 y[11] = f582 + f583; | 1033 f585 = f575 + f578; |
1034 y[53] = f584 - f583; | 1034 f586 = f575 - f578; |
1035 f587 = f548 + f435; | 1035 f587 = f576 + f577; |
1036 f588 = MUL_C_R(COEF_CONST((-0.6358464401941451)), f548); | 1036 f588 = f577 - f576; |
1037 f589 = MUL_C_R(COEF_CONST(0.9495281805930367), f587); | 1037 f589 = f582 - f588; |
1038 f590 = MUL_C_R(COEF_CONST(1.2632099209919283), f435); | 1038 f590 = f582 + f588; |
1039 y[13] = f588 + f589; | 1039 f591 = f584 - f586; |
1040 y[51] = f590 - f589; | 1040 f592 = f584 + f586; |
1041 f593 = f550 + f437; | 1041 f593 = f590 - f592; |
1042 f594 = MUL_C_R(COEF_CONST((-0.5730977622997506)), f550); | 1042 f594 = f590 + f592; |
1043 f595 = MUL_C_R(COEF_CONST(0.9329927988347389), f593); | 1043 f595 = MUL_C(COEF_CONST(0.7071067811865476), f593); |
1044 f596 = MUL_C_R(COEF_CONST(1.2928878353697271), f437); | 1044 f596 = f589 + f591; |
1045 y[15] = f594 + f595; | 1045 f597 = MUL_C(COEF_CONST(1.3065629648763766), f589); |
1046 y[49] = f596 - f595; | 1046 f598 = MUL_C(COEF_CONST((-0.9238795325112866)), f596); |
1047 f599 = f549 + f436; | 1047 f599 = MUL_C(COEF_CONST((-0.5411961001461967)), f591); |
1048 f600 = MUL_C_R(COEF_CONST((-0.5089684416985408)), f549); | 1048 f600 = f597 + f598; |
1049 f601 = MUL_C_R(COEF_CONST(0.9142097557035307), f599); | 1049 f601 = f599 - f598; |
1050 f602 = MUL_C_R(COEF_CONST(1.3194510697085207), f436); | 1050 f602 = f583 + f585; |
1051 y[17] = f600 + f601; | 1051 f603 = f585 - f583; |
1052 y[47] = f602 - f601; | 1052 f604 = MUL_C(COEF_CONST(0.7071067811865476), f603); |
1053 f605 = f434 - f547; | 1053 f605 = MUL_C(COEF_CONST(0.7071067811865476), f602); |
1054 f606 = MUL_C_R(COEF_CONST((-0.4436129715409087)), f547); | 1054 f606 = f581 - f604; |
1055 f607 = MUL_C_R(COEF_CONST(0.8932243011955153), f605); | 1055 f607 = f581 + f604; |
1056 f608 = MUL_C_R(COEF_CONST(1.3428356308501219), f434); | 1056 f608 = f605 - f587; |
1057 y[19] = f607 - f606; | 1057 f609 = f587 + f605; |
1058 y[45] = f608 - f607; | 1058 f610 = f607 - f609; |
1059 f611 = f545 + f432; | 1059 f611 = MUL_C(COEF_CONST((-0.7856949583871021)), f609); |
1060 f612 = MUL_C_R(COEF_CONST((-0.3771887988789273)), f545); | 1060 f612 = MUL_C(COEF_CONST(0.9807852804032304), f610); |
1061 f613 = MUL_C_R(COEF_CONST(0.8700869911087114), f611); | 1061 f613 = MUL_C(COEF_CONST(1.1758756024193588), f607); |
1062 f614 = MUL_C_R(COEF_CONST(1.3629851833384954), f432); | 1062 f614 = f612 - f611; |
1063 y[21] = f612 + f613; | 1063 f615 = f613 - f612; |
1064 y[43] = f614 - f613; | 1064 f616 = f608 + f606; |
1065 f617 = f430 - f543; | 1065 f617 = MUL_C(COEF_CONST(0.2758993792829431), f608); |
1066 f618 = MUL_C_R(COEF_CONST((-0.3098559453626097)), f543); | 1066 f618 = MUL_C(COEF_CONST(0.5555702330196022), f616); |
1067 f619 = MUL_C_R(COEF_CONST(0.8448535652497070), f617); | 1067 f619 = MUL_C(COEF_CONST(1.3870398453221475), f606); |
1068 f620 = MUL_C_R(COEF_CONST(1.3798511851368043), f430); | 1068 f620 = f617 + f618; |
1069 y[23] = f619 - f618; | 1069 f621 = f619 - f618; |
1070 y[41] = f620 - f619; | 1070 f622 = f594 + f614; |
1071 f623 = f541 + f428; | 1071 f623 = f614 + f601; |
1072 f624 = MUL_C_R(COEF_CONST((-0.2417766217337384)), f541); | 1072 f624 = f601 + f621; |
1073 f625 = MUL_C_R(COEF_CONST(0.8175848131515837), f623); | 1073 f625 = f621 + f595; |
1074 f626 = MUL_C_R(COEF_CONST(1.3933930045694289), f428); | 1074 f626 = f595 + f620; |
1075 y[25] = f624 + f625; | 1075 f627 = f620 + f600; |
1076 y[39] = f626 - f625; | 1076 f628 = f600 + f615; |
1077 f629 = f426 - f539; | 1077 y[5] = f496 - f615; |
1078 f630 = MUL_C_R(COEF_CONST((-0.1731148370459794)), f539); | 1078 y[3] = f496 + f615; |
1079 f631 = MUL_C_R(COEF_CONST(0.7883464276266062), f629); | 1079 y[9] = f453 - f568; |
1080 f632 = MUL_C_R(COEF_CONST(1.4035780182072330), f426); | 1080 y[7] = f453 + f568; |
1081 y[27] = f631 - f630; | 1081 y[13] = f502 - f628; |
1082 y[37] = f632 - f631; | 1082 y[11] = f502 + f628; |
1083 f635 = f537 + f424; | 1083 y[17] = f439 - f551; |
1084 f636 = MUL_C_R(COEF_CONST((-0.1040360035527077)), f537); | 1084 y[15] = f439 + f551; |
1085 f637 = MUL_C_R(COEF_CONST(0.7572088465064845), f635); | 1085 y[21] = f508 - f627; |
1086 f638 = MUL_C_R(COEF_CONST(1.4103816894602612), f424); | 1086 y[19] = f508 + f627; |
1087 y[29] = f636 + f637; | 1087 y[25] = f459 - f571; |
1088 y[35] = f638 - f637; | 1088 y[23] = f459 + f571; |
1089 f641 = f422 - f535; | 1089 y[29] = f514 - f626; |
1090 f642 = MUL_C_R(COEF_CONST((-0.0347065382144000)), f535); | 1090 y[27] = f514 + f626; |
1091 f643 = MUL_C_R(COEF_CONST(0.7242470829514669), f641); | 1091 y[33] = f433 - f546; |
1092 f644 = MUL_C_R(COEF_CONST(1.4137876276885337), f422); | 1092 y[31] = f433 + f546; |
1093 y[31] = f643 - f642; | 1093 y[37] = f515 - f625; |
1094 y[33] = f644 - f643; | 1094 y[35] = f515 + f625; |
1095 y[41] = f458 - f572; | |
1096 y[39] = f458 + f572; | |
1097 y[45] = f509 - f624; | |
1098 y[43] = f509 + f624; | |
1099 y[49] = f438 - f552; | |
1100 y[47] = f438 + f552; | |
1101 y[53] = f503 - f623; | |
1102 y[51] = f503 + f623; | |
1103 y[57] = f455 - f566; | |
1104 y[55] = f455 + f566; | |
1105 y[61] = f497 - f622; | |
1106 y[59] = f497 + f622; | |
1107 } | |
1108 | |
1109 void DST2_64_unscaled(real_t *y, real_t *x) | |
1110 { | |
1111 int16_t i0; | |
1112 real_t f2, f3, f4, f5, f6, f7; | |
1113 real_t f8, f9, f10, f11, f12, f13; | |
1114 real_t f14, f15, f16, f17, f18, f19; | |
1115 real_t f20, f21, f22, f23, f24, f25; | |
1116 real_t f26, f27, f28, f29, f30, f31; | |
1117 real_t f32, f33, f34, f35, f36, f37; | |
1118 real_t f38, f39, f40, f41, f42, f43; | |
1119 real_t f44, f45, f46, f47, f48, f49; | |
1120 real_t f50, f51, f52, f53, f54, f55; | |
1121 real_t f56, f57, f58, f59, f60, f61; | |
1122 real_t f62, f63, f64, f65, f66, f67; | |
1123 real_t f68, f69, f70, f71, f72, f73; | |
1124 real_t f74, f75, f76, f77, f78, f79; | |
1125 real_t f80, f81, f82, f83, f84, f85; | |
1126 real_t f86, f87, f88, f89, f90, f91; | |
1127 real_t f92, f93, f94, f95, f96, f97; | |
1128 real_t f98, f99, f100, f101, f102, f103; | |
1129 real_t f104, f105, f106, f107, f108, f109; | |
1130 real_t f110, f111, f112, f113, f114, f115; | |
1131 real_t f116, f117, f118, f119, f120, f121; | |
1132 real_t f122, f123, f124, f125, f126, f127; | |
1133 real_t f128, f129, f130, f131, f132, f133; | |
1134 real_t f134, f135, f136, f137, f138, f139; | |
1135 real_t f140, f141, f142, f143, f144, f145; | |
1136 real_t f146, f147, f148, f149, f150, f151; | |
1137 real_t f152, f153, f154, f155, f156, f157; | |
1138 real_t f158, f159, f160, f161, f162, f163; | |
1139 real_t f164, f165, f166, f167, f168, f169; | |
1140 real_t f170, f171, f172, f173, f174, f175; | |
1141 real_t f176, f177, f178, f179, f180, f181; | |
1142 real_t f182, f183, f184, f185, f186, f187; | |
1143 real_t f188, f189, f190, f191, f192, f193; | |
1144 real_t f194, f195, f196, f197, f198, f199; | |
1145 real_t f200, f201, f202, f203, f204, f205; | |
1146 real_t f206, f207, f208, f209, f210, f211; | |
1147 real_t f212, f213, f214, f215, f216, f217; | |
1148 real_t f218, f219, f220, f221, f222, f223; | |
1149 real_t f224, f225, f226, f227, f228, f229; | |
1150 real_t f230, f231, f232, f233, f234, f235; | |
1151 real_t f236, f237, f238, f239, f240, f241; | |
1152 real_t f242, f243, f244, f245, f246, f247; | |
1153 real_t f248, f249, f250, f251, f252, f253; | |
1154 real_t f254, f255, f256, f257, f258, f259; | |
1155 real_t f260, f261, f264, f265, f266, f267; | |
1156 real_t f270, f271, f272, f273, f276, f277; | |
1157 real_t f278, f279, f282, f283, f284, f285; | |
1158 real_t f288, f289, f290, f291, f294, f295; | |
1159 real_t f296, f297, f300, f301, f302, f303; | |
1160 real_t f306, f307, f308, f309, f312, f313; | |
1161 real_t f314, f315, f318, f319, f320, f321; | |
1162 real_t f324, f325, f326, f327, f330, f331; | |
1163 real_t f332, f333, f336, f337, f338, f339; | |
1164 real_t f342, f343, f344, f345, f348, f349; | |
1165 real_t f350, f351, f354, f355, f356, f357; | |
1166 real_t f358, f359, f360, f361, f362, f363; | |
1167 real_t f364, f365, f366, f367, f368, f369; | |
1168 real_t f370, f371, f372, f373, f374, f375; | |
1169 real_t f376, f377, f378, f379, f380, f381; | |
1170 real_t f382, f383, f384, f385, f386, f387; | |
1171 real_t f388, f389, f390, f391, f392, f393; | |
1172 real_t f394, f395, f396, f397, f398, f399; | |
1173 real_t f400, f401, f402, f403, f404, f405; | |
1174 real_t f406, f407, f408, f409, f410, f411; | |
1175 real_t f412, f413, f414, f415, f416, f417; | |
1176 real_t f418, f419, f420, f421, f422, f423; | |
1177 real_t f424, f425, f426, f427, f428, f429; | |
1178 real_t f430, f431, f432, f433, f434, f435; | |
1179 real_t f436, f437, f438, f439, f440, f441; | |
1180 real_t f442, f443, f444, f445, f446, f447; | |
1181 real_t f448, f449, f450, f451, f452, f453; | |
1182 real_t f454, f455, f456, f457, f458, f459; | |
1183 real_t f460, f461, f462, f463, f464, f465; | |
1184 real_t f466, f467, f468, f469, f470, f471; | |
1185 real_t f472, f473, f474, f475, f476, f477; | |
1186 real_t f478, f479, f480, f481, f482, f483; | |
1187 real_t f484, f485, f486, f487, f488, f489; | |
1188 real_t f490, f491, f492, f493, f494, f495; | |
1189 real_t f496, f497, f498, f499, f500, f501; | |
1190 real_t f504, f505, f506, f507, f510, f511; | |
1191 real_t f512, f513, f516, f517, f518, f519; | |
1192 real_t f522, f523, f524, f525, f528, f529; | |
1193 real_t f530, f531, f534, f535, f536, f537; | |
1194 real_t f540, f541, f542, f543, f546, f547; | |
1195 real_t f548, f549, f550, f551, f552, f553; | |
1196 real_t f554, f555, f556, f557, f558, f559; | |
1197 real_t f560, f561, f562, f563, f564, f565; | |
1198 real_t f566, f567, f568, f569, f570, f571; | |
1199 real_t f572, f573, f574, f577, f578, f579; | |
1200 real_t f580, f583, f584, f585, f586, f587; | |
1201 real_t f588, f589, f590, f591, f592, f593; | |
1202 real_t f594, f595, f596, f597, f598, f603; | |
1203 real_t f604, f605, f606, f607, f608, f609; | |
1204 real_t f610, f611, f612, f613, f614, f615; | |
1205 real_t f616, f617, f618, f619, f620, f621; | |
1206 real_t f622, f623, f624, f625, f626, f627; | |
1207 real_t f628, f629, f630, f631, f632, f633; | |
1208 real_t f634, f635, f636, f637, f638, f639; | |
1209 real_t f640, f641, f642, f643, f644, f645; | |
1210 real_t f646, f647, f648, f649, f650; | |
1211 ALIGN static real_t t2[64]; | |
1212 | |
1213 for (i0=0; i0<32; i0++) | |
1214 { | |
1215 t2[2*i0+1] = x[i0] - x[-i0+63]; | |
1216 t2[2*i0] = x[i0] + x[-i0+63]; | |
1217 } | |
1218 f2 = t2[2] + t2[4]; | |
1219 f3 = t2[4] - t2[2]; | |
1220 f4 = t2[6] + t2[8]; | |
1221 f5 = t2[8] - t2[6]; | |
1222 f6 = t2[10] + t2[12]; | |
1223 f7 = t2[12] - t2[10]; | |
1224 f8 = t2[14] + t2[16]; | |
1225 f9 = t2[16] - t2[14]; | |
1226 f10 = t2[18] + t2[20]; | |
1227 f11 = t2[20] - t2[18]; | |
1228 f12 = t2[22] + t2[24]; | |
1229 f13 = t2[24] - t2[22]; | |
1230 f14 = t2[26] + t2[28]; | |
1231 f15 = t2[28] - t2[26]; | |
1232 f16 = t2[30] + t2[32]; | |
1233 f17 = t2[32] - t2[30]; | |
1234 f18 = t2[34] + t2[36]; | |
1235 f19 = t2[36] - t2[34]; | |
1236 f20 = t2[38] + t2[40]; | |
1237 f21 = t2[40] - t2[38]; | |
1238 f22 = t2[42] + t2[44]; | |
1239 f23 = t2[44] - t2[42]; | |
1240 f24 = t2[46] + t2[48]; | |
1241 f25 = t2[48] - t2[46]; | |
1242 f26 = t2[50] + t2[52]; | |
1243 f27 = t2[52] - t2[50]; | |
1244 f28 = t2[54] + t2[56]; | |
1245 f29 = t2[56] - t2[54]; | |
1246 f30 = t2[58] + t2[60]; | |
1247 f31 = t2[60] - t2[58]; | |
1248 f32 = MUL_C(COEF_CONST(0.7071067811865476), f17); | |
1249 f33 = t2[0] - f32; | |
1250 f34 = t2[0] + f32; | |
1251 f35 = f9 + f25; | |
1252 f36 = MUL_C(COEF_CONST(1.3065629648763766), f9); | |
1253 f37 = MUL_C(COEF_CONST((-0.9238795325112866)), f35); | |
1254 f38 = MUL_C(COEF_CONST((-0.5411961001461967)), f25); | |
1255 f39 = f36 + f37; | |
1256 f40 = f38 - f37; | |
1257 f41 = f34 - f40; | |
1258 f42 = f34 + f40; | |
1259 f43 = f33 - f39; | |
1260 f44 = f33 + f39; | |
1261 f45 = MUL_C(COEF_CONST(2.5629154477415064), f5); | |
1262 f46 = MUL_C(COEF_CONST(0.8999762231364158), f13); | |
1263 f47 = MUL_C(COEF_CONST(0.5097955791041592), f29); | |
1264 f48 = MUL_C(COEF_CONST(0.6013448869350453), f21); | |
1265 f49 = f45 - f47; | |
1266 f50 = f45 + f47; | |
1267 f51 = f46 - f48; | |
1268 f52 = f46 + f48; | |
1269 f53 = f50 + f52; | |
1270 f54 = MUL_C(COEF_CONST(1.3065629648763770), f50); | |
1271 f55 = MUL_C(COEF_CONST((-0.3826834323650904)), f53); | |
1272 f56 = MUL_C(COEF_CONST(0.5411961001461961), f52); | |
1273 f57 = f54 + f55; | |
1274 f58 = f56 - f55; | |
1275 f59 = f51 - f49; | |
1276 f60 = f49 + f51; | |
1277 f61 = MUL_C(COEF_CONST(0.7071067811865476), f60); | |
1278 f62 = f58 - f61; | |
1279 f63 = f57 - f61; | |
1280 f64 = f59 + f57; | |
1281 f65 = f42 - f58; | |
1282 f66 = f42 + f58; | |
1283 f67 = f44 + f62; | |
1284 f68 = f44 - f62; | |
1285 f69 = f43 - f63; | |
1286 f70 = f43 + f63; | |
1287 f71 = f41 + f64; | |
1288 f72 = f41 - f64; | |
1289 f73 = f7 - f11; | |
1290 f74 = f7 + f11; | |
1291 f75 = f15 - f19; | |
1292 f76 = f15 + f19; | |
1293 f77 = f23 - f27; | |
1294 f78 = f23 + f27; | |
1295 f79 = MUL_C(COEF_CONST(0.7071067811865476), f76); | |
1296 f80 = f3 - f79; | |
1297 f81 = f3 + f79; | |
1298 f82 = f74 + f78; | |
1299 f83 = MUL_C(COEF_CONST(1.3065629648763766), f74); | |
1300 f84 = MUL_C(COEF_CONST((-0.9238795325112866)), f82); | |
1301 f85 = MUL_C(COEF_CONST((-0.5411961001461967)), f78); | |
1302 f86 = f83 + f84; | |
1303 f87 = f85 - f84; | |
1304 f88 = f81 - f87; | |
1305 f89 = f81 + f87; | |
1306 f90 = f80 - f86; | |
1307 f91 = f80 + f86; | |
1308 f92 = MUL_C(COEF_CONST(0.7071067811865476), f75); | |
1309 f93 = f31 - f92; | |
1310 f94 = f31 + f92; | |
1311 f95 = f77 + f73; | |
1312 f96 = MUL_C(COEF_CONST(1.3065629648763766), f77); | |
1313 f97 = MUL_C(COEF_CONST((-0.9238795325112866)), f95); | |
1314 f98 = MUL_C(COEF_CONST((-0.5411961001461967)), f73); | |
1315 f99 = f96 + f97; | |
1316 f100 = f98 - f97; | |
1317 f101 = f94 - f100; | |
1318 f102 = f94 + f100; | |
1319 f103 = f93 - f99; | |
1320 f104 = f93 + f99; | |
1321 f105 = f102 + f89; | |
1322 f106 = MUL_C(COEF_CONST((-0.8971675863426361)), f102); | |
1323 f107 = MUL_C(COEF_CONST(0.9951847266721968), f105); | |
1324 f108 = MUL_C(COEF_CONST(1.0932018670017576), f89); | |
1325 f109 = f106 + f107; | |
1326 f110 = f108 - f107; | |
1327 f111 = f91 - f104; | |
1328 f112 = MUL_C(COEF_CONST((-0.6666556584777466)), f104); | |
1329 f113 = MUL_C(COEF_CONST(0.9569403357322089), f111); | |
1330 f114 = MUL_C(COEF_CONST(1.2472250129866713), f91); | |
1331 f115 = f113 - f112; | |
1332 f116 = f114 - f113; | |
1333 f117 = f103 + f90; | |
1334 f118 = MUL_C(COEF_CONST((-0.4105245275223571)), f103); | |
1335 f119 = MUL_C(COEF_CONST(0.8819212643483549), f117); | |
1336 f120 = MUL_C(COEF_CONST(1.3533180011743529), f90); | |
1337 f121 = f118 + f119; | |
1338 f122 = f120 - f119; | |
1339 f123 = f88 - f101; | |
1340 f124 = MUL_C(COEF_CONST((-0.1386171691990915)), f101); | |
1341 f125 = MUL_C(COEF_CONST(0.7730104533627370), f123); | |
1342 f126 = MUL_C(COEF_CONST(1.4074037375263826), f88); | |
1343 f127 = f125 - f124; | |
1344 f128 = f126 - f125; | |
1345 f129 = f66 - f109; | |
1346 f130 = f66 + f109; | |
1347 f131 = f68 - f115; | |
1348 f132 = f68 + f115; | |
1349 f133 = f70 - f121; | |
1350 f134 = f70 + f121; | |
1351 f135 = f72 - f127; | |
1352 f136 = f72 + f127; | |
1353 f137 = f71 - f128; | |
1354 f138 = f71 + f128; | |
1355 f139 = f69 - f122; | |
1356 f140 = f69 + f122; | |
1357 f141 = f67 - f116; | |
1358 f142 = f67 + f116; | |
1359 f143 = f65 - f110; | |
1360 f144 = f65 + f110; | |
1361 f145 = f26 + f30; | |
1362 f146 = f22 + f26; | |
1363 f147 = f18 + f22; | |
1364 f148 = f14 + f18; | |
1365 f149 = f10 + f14; | |
1366 f150 = f6 + f10; | |
1367 f151 = f2 + f6; | |
1368 f152 = MUL_C(COEF_CONST(0.7071067811865476), f148); | |
1369 f153 = f152 - f30; | |
1370 f154 = f30 + f152; | |
1371 f155 = f146 + f150; | |
1372 f156 = MUL_C(COEF_CONST(1.3065629648763766), f146); | |
1373 f157 = MUL_C(COEF_CONST((-0.9238795325112866)), f155); | |
1374 f158 = MUL_C(COEF_CONST((-0.5411961001461967)), f150); | |
1375 f159 = f156 + f157; | |
1376 f160 = f157 - f158; | |
1377 f161 = f154 + f160; | |
1378 f162 = f160 - f154; | |
1379 f163 = f153 + f159; | |
1380 f164 = f153 - f159; | |
1381 f165 = f147 + f145; | |
1382 f166 = f149 + f147; | |
1383 f167 = f151 + f149; | |
1384 f168 = MUL_C(COEF_CONST(0.7071067811865476), f166); | |
1385 f169 = f168 - f145; | |
1386 f170 = f145 + f168; | |
1387 f171 = f165 + f167; | |
1388 f172 = MUL_C(COEF_CONST(1.3065629648763766), f165); | |
1389 f173 = MUL_C(COEF_CONST((-0.9238795325112866)), f171); | |
1390 f174 = MUL_C(COEF_CONST((-0.5411961001461967)), f167); | |
1391 f175 = f172 + f173; | |
1392 f176 = f173 - f174; | |
1393 f177 = f170 + f176; | |
1394 f178 = f176 - f170; | |
1395 f179 = f169 + f175; | |
1396 f180 = f169 - f175; | |
1397 f181 = MUL_C(COEF_CONST(0.5097955791041592), f178); | |
1398 f182 = MUL_C(COEF_CONST(0.6013448869350453), f180); | |
1399 f183 = MUL_C(COEF_CONST(0.8999762231364156), f179); | |
1400 f184 = MUL_C(COEF_CONST(2.5629154477415055), f177); | |
1401 f185 = f162 - f181; | |
1402 f186 = f162 + f181; | |
1403 f187 = f164 - f182; | |
1404 f188 = f164 + f182; | |
1405 f189 = f163 - f183; | |
1406 f190 = f163 + f183; | |
1407 f191 = f184 - f161; | |
1408 f192 = f161 + f184; | |
1409 f193 = MUL_C(COEF_CONST(0.5024192861881557), f186); | |
1410 f194 = MUL_C(COEF_CONST(0.5224986149396889), f188); | |
1411 f195 = MUL_C(COEF_CONST(0.5669440348163577), f190); | |
1412 f196 = MUL_C(COEF_CONST(0.6468217833599901), f192); | |
1413 f197 = MUL_C(COEF_CONST(0.7881546234512502), f191); | |
1414 f198 = MUL_C(COEF_CONST(1.0606776859903471), f189); | |
1415 f199 = MUL_C(COEF_CONST(1.7224470982383342), f187); | |
1416 f200 = MUL_C(COEF_CONST(5.1011486186891553), f185); | |
1417 f201 = MUL_C(COEF_CONST(0.7071067811865476), f16); | |
1418 f202 = f201 - t2[62]; | |
1419 f203 = t2[62] + f201; | |
1420 f204 = f24 + f8; | |
1421 f205 = MUL_C(COEF_CONST(1.3065629648763766), f24); | |
1422 f206 = MUL_C(COEF_CONST((-0.9238795325112866)), f204); | |
1423 f207 = MUL_C(COEF_CONST((-0.5411961001461967)), f8); | |
1424 f208 = f205 + f206; | |
1425 f209 = f206 - f207; | |
1426 f210 = f203 + f209; | |
1427 f211 = f209 - f203; | |
1428 f212 = f202 + f208; | |
1429 f213 = f202 - f208; | |
1430 f214 = f20 + f28; | |
1431 f215 = f12 + f20; | |
1432 f216 = f4 + f12; | |
1433 f217 = MUL_C(COEF_CONST(0.7071067811865476), f215); | |
1434 f218 = f217 - f28; | |
1435 f219 = f28 + f217; | |
1436 f220 = f214 + f216; | |
1437 f221 = MUL_C(COEF_CONST(1.3065629648763766), f214); | |
1438 f222 = MUL_C(COEF_CONST((-0.9238795325112866)), f220); | |
1439 f223 = MUL_C(COEF_CONST((-0.5411961001461967)), f216); | |
1440 f224 = f221 + f222; | |
1441 f225 = f222 - f223; | |
1442 f226 = f219 + f225; | |
1443 f227 = f225 - f219; | |
1444 f228 = f218 + f224; | |
1445 f229 = f218 - f224; | |
1446 f230 = MUL_C(COEF_CONST(0.5097955791041592), f227); | |
1447 f231 = MUL_C(COEF_CONST(0.6013448869350453), f229); | |
1448 f232 = MUL_C(COEF_CONST(0.8999762231364156), f228); | |
1449 f233 = MUL_C(COEF_CONST(2.5629154477415055), f226); | |
1450 f234 = f211 - f230; | |
1451 f235 = f211 + f230; | |
1452 f236 = f213 - f231; | |
1453 f237 = f213 + f231; | |
1454 f238 = f212 - f232; | |
1455 f239 = f212 + f232; | |
1456 f240 = f233 - f210; | |
1457 f241 = f210 + f233; | |
1458 f242 = f193 - f235; | |
1459 f243 = f193 + f235; | |
1460 f244 = f237 - f194; | |
1461 f245 = f194 + f237; | |
1462 f246 = f195 - f239; | |
1463 f247 = f195 + f239; | |
1464 f248 = f196 - f241; | |
1465 f249 = f196 + f241; | |
1466 f250 = f197 - f240; | |
1467 f251 = f197 + f240; | |
1468 f252 = f238 - f198; | |
1469 f253 = f198 + f238; | |
1470 f254 = f199 - f236; | |
1471 f255 = f199 + f236; | |
1472 f256 = f234 - f200; | |
1473 f257 = f200 + f234; | |
1474 f258 = f243 + f130; | |
1475 f259 = MUL_C(COEF_CONST((-0.9751575901732920)), f243); | |
1476 f260 = MUL_C(COEF_CONST(0.9996988186962043), f258); | |
1477 f261 = MUL_C(COEF_CONST(1.0242400472191164), f130); | |
1478 y[62] = f259 + f260; | |
1479 y[0] = f261 - f260; | |
1480 f264 = f132 - f245; | |
1481 f265 = MUL_C(COEF_CONST((-0.9237258930790228)), f245); | |
1482 f266 = MUL_C(COEF_CONST(0.9972904566786902), f264); | |
1483 f267 = MUL_C(COEF_CONST(1.0708550202783576), f132); | |
1484 y[60] = f266 - f265; | |
1485 y[2] = f267 - f266; | |
1486 f270 = f247 + f134; | |
1487 f271 = MUL_C(COEF_CONST((-0.8700688593994936)), f247); | |
1488 f272 = MUL_C(COEF_CONST(0.9924795345987100), f270); | |
1489 f273 = MUL_C(COEF_CONST(1.1148902097979263), f134); | |
1490 y[58] = f271 + f272; | |
1491 y[4] = f273 - f272; | |
1492 f276 = f249 + f136; | |
1493 f277 = MUL_C(COEF_CONST((-0.8143157536286398)), f249); | |
1494 f278 = MUL_C(COEF_CONST(0.9852776423889412), f276); | |
1495 f279 = MUL_C(COEF_CONST(1.1562395311492426), f136); | |
1496 y[56] = f277 + f278; | |
1497 y[6] = f279 - f278; | |
1498 f282 = f251 + f138; | |
1499 f283 = MUL_C(COEF_CONST((-0.7566008898816587)), f251); | |
1500 f284 = MUL_C(COEF_CONST(0.9757021300385286), f282); | |
1501 f285 = MUL_C(COEF_CONST(1.1948033701953984), f138); | |
1502 y[54] = f283 + f284; | |
1503 y[8] = f285 - f284; | |
1504 f288 = f140 - f253; | |
1505 f289 = MUL_C(COEF_CONST((-0.6970633083205414)), f253); | |
1506 f290 = MUL_C(COEF_CONST(0.9637760657954398), f288); | |
1507 f291 = MUL_C(COEF_CONST(1.2304888232703384), f140); | |
1508 y[52] = f290 - f289; | |
1509 y[10] = f291 - f290; | |
1510 f294 = f255 + f142; | |
1511 f295 = MUL_C(COEF_CONST((-0.6358464401941451)), f255); | |
1512 f296 = MUL_C(COEF_CONST(0.9495281805930367), f294); | |
1513 f297 = MUL_C(COEF_CONST(1.2632099209919283), f142); | |
1514 y[50] = f295 + f296; | |
1515 y[12] = f297 - f296; | |
1516 f300 = f144 - f257; | |
1517 f301 = MUL_C(COEF_CONST((-0.5730977622997506)), f257); | |
1518 f302 = MUL_C(COEF_CONST(0.9329927988347389), f300); | |
1519 f303 = MUL_C(COEF_CONST(1.2928878353697271), f144); | |
1520 y[48] = f302 - f301; | |
1521 y[14] = f303 - f302; | |
1522 f306 = f256 + f143; | |
1523 f307 = MUL_C(COEF_CONST((-0.5089684416985408)), f256); | |
1524 f308 = MUL_C(COEF_CONST(0.9142097557035307), f306); | |
1525 f309 = MUL_C(COEF_CONST(1.3194510697085207), f143); | |
1526 y[46] = f307 + f308; | |
1527 y[16] = f309 - f308; | |
1528 f312 = f254 + f141; | |
1529 f313 = MUL_C(COEF_CONST((-0.4436129715409087)), f254); | |
1530 f314 = MUL_C(COEF_CONST(0.8932243011955153), f312); | |
1531 f315 = MUL_C(COEF_CONST(1.3428356308501219), f141); | |
1532 y[44] = f313 + f314; | |
1533 y[18] = f315 - f314; | |
1534 f318 = f252 + f139; | |
1535 f319 = MUL_C(COEF_CONST((-0.3771887988789273)), f252); | |
1536 f320 = MUL_C(COEF_CONST(0.8700869911087114), f318); | |
1537 f321 = MUL_C(COEF_CONST(1.3629851833384954), f139); | |
1538 y[42] = f319 + f320; | |
1539 y[20] = f321 - f320; | |
1540 f324 = f250 + f137; | |
1541 f325 = MUL_C(COEF_CONST((-0.3098559453626097)), f250); | |
1542 f326 = MUL_C(COEF_CONST(0.8448535652497070), f324); | |
1543 f327 = MUL_C(COEF_CONST(1.3798511851368043), f137); | |
1544 y[40] = f325 + f326; | |
1545 y[22] = f327 - f326; | |
1546 f330 = f248 + f135; | |
1547 f331 = MUL_C(COEF_CONST((-0.2417766217337384)), f248); | |
1548 f332 = MUL_C(COEF_CONST(0.8175848131515837), f330); | |
1549 f333 = MUL_C(COEF_CONST(1.3933930045694289), f135); | |
1550 y[38] = f331 + f332; | |
1551 y[24] = f333 - f332; | |
1552 f336 = f246 + f133; | |
1553 f337 = MUL_C(COEF_CONST((-0.1731148370459794)), f246); | |
1554 f338 = MUL_C(COEF_CONST(0.7883464276266062), f336); | |
1555 f339 = MUL_C(COEF_CONST(1.4035780182072330), f133); | |
1556 y[36] = f337 + f338; | |
1557 y[26] = f339 - f338; | |
1558 f342 = f244 + f131; | |
1559 f343 = MUL_C(COEF_CONST((-0.1040360035527077)), f244); | |
1560 f344 = MUL_C(COEF_CONST(0.7572088465064845), f342); | |
1561 f345 = MUL_C(COEF_CONST(1.4103816894602612), f131); | |
1562 y[34] = f343 + f344; | |
1563 y[28] = f345 - f344; | |
1564 f348 = f242 + f129; | |
1565 f349 = MUL_C(COEF_CONST((-0.0347065382144000)), f242); | |
1566 f350 = MUL_C(COEF_CONST(0.7242470829514669), f348); | |
1567 f351 = MUL_C(COEF_CONST(1.4137876276885337), f129); | |
1568 y[32] = f349 + f350; | |
1569 y[30] = f351 - f350; | |
1570 f354 = t2[1] - t2[63]; | |
1571 f355 = t2[1] + t2[63]; | |
1572 f356 = t2[3] - t2[61]; | |
1573 f357 = t2[3] + t2[61]; | |
1574 f358 = t2[5] - t2[59]; | |
1575 f359 = t2[5] + t2[59]; | |
1576 f360 = t2[7] - t2[57]; | |
1577 f361 = t2[7] + t2[57]; | |
1578 f362 = t2[9] - t2[55]; | |
1579 f363 = t2[9] + t2[55]; | |
1580 f364 = t2[11] - t2[53]; | |
1581 f365 = t2[11] + t2[53]; | |
1582 f366 = t2[13] - t2[51]; | |
1583 f367 = t2[13] + t2[51]; | |
1584 f368 = t2[15] - t2[49]; | |
1585 f369 = t2[15] + t2[49]; | |
1586 f370 = t2[17] - t2[47]; | |
1587 f371 = t2[17] + t2[47]; | |
1588 f372 = t2[19] - t2[45]; | |
1589 f373 = t2[19] + t2[45]; | |
1590 f374 = t2[21] - t2[43]; | |
1591 f375 = t2[21] + t2[43]; | |
1592 f376 = t2[23] - t2[41]; | |
1593 f377 = t2[23] + t2[41]; | |
1594 f378 = t2[25] - t2[39]; | |
1595 f379 = t2[25] + t2[39]; | |
1596 f380 = t2[27] - t2[37]; | |
1597 f381 = t2[27] + t2[37]; | |
1598 f382 = t2[29] - t2[35]; | |
1599 f383 = t2[29] + t2[35]; | |
1600 f384 = t2[31] - t2[33]; | |
1601 f385 = t2[31] + t2[33]; | |
1602 f386 = f369 + f371; | |
1603 f387 = f371 - f369; | |
1604 f388 = MUL_C(COEF_CONST(0.7071067811865476), f387); | |
1605 f389 = MUL_C(COEF_CONST(0.7071067811865476), f386); | |
1606 f390 = f363 + f377; | |
1607 f391 = f363 - f377; | |
1608 f392 = MUL_C(COEF_CONST(0.7071067811865476), f391); | |
1609 f393 = MUL_C(COEF_CONST(0.7071067811865476), f390); | |
1610 f394 = f367 + f373; | |
1611 f395 = f367 - f373; | |
1612 f396 = MUL_C(COEF_CONST(0.7071067811865476), f395); | |
1613 f397 = MUL_C(COEF_CONST(0.7071067811865476), f394); | |
1614 f398 = f365 + f375; | |
1615 f399 = f375 - f365; | |
1616 f400 = MUL_C(COEF_CONST(0.7071067811865476), f399); | |
1617 f401 = MUL_C(COEF_CONST(0.7071067811865476), f398); | |
1618 f402 = f355 - f388; | |
1619 f403 = f355 + f388; | |
1620 f404 = f389 - f385; | |
1621 f405 = f385 + f389; | |
1622 f406 = f361 + f392; | |
1623 f407 = f392 - f361; | |
1624 f408 = f379 - f393; | |
1625 f409 = f379 + f393; | |
1626 f410 = f357 + f396; | |
1627 f411 = f396 - f357; | |
1628 f412 = f383 - f397; | |
1629 f413 = f383 + f397; | |
1630 f414 = f359 - f400; | |
1631 f415 = f359 + f400; | |
1632 f416 = f401 - f381; | |
1633 f417 = f381 + f401; | |
1634 f418 = f409 + f407; | |
1635 f419 = MUL_C(COEF_CONST((-0.5411961001461969)), f409); | |
1636 f420 = MUL_C(COEF_CONST(0.9238795325112867), f418); | |
1637 f421 = MUL_C(COEF_CONST(1.3065629648763766), f407); | |
1638 f422 = f419 + f420; | |
1639 f423 = f421 - f420; | |
1640 f424 = f408 - f406; | |
1641 f425 = MUL_C(COEF_CONST(1.3065629648763770), f408); | |
1642 f426 = MUL_C(COEF_CONST((-0.3826834323650904)), f424); | |
1643 f427 = MUL_C(COEF_CONST(0.5411961001461961), f406); | |
1644 f428 = f425 + f426; | |
1645 f429 = f426 + f427; | |
1646 f430 = f415 - f417; | |
1647 f431 = MUL_C(COEF_CONST((-0.5411961001461969)), f417); | |
1648 f432 = MUL_C(COEF_CONST(0.9238795325112867), f430); | |
1649 f433 = MUL_C(COEF_CONST(1.3065629648763766), f415); | |
1650 f434 = f432 - f431; | |
1651 f435 = f433 - f432; | |
1652 f436 = f416 + f414; | |
1653 f437 = MUL_C(COEF_CONST(1.3065629648763770), f416); | |
1654 f438 = MUL_C(COEF_CONST((-0.3826834323650904)), f436); | |
1655 f439 = MUL_C(COEF_CONST(0.5411961001461961), f414); | |
1656 f440 = f437 + f438; | |
1657 f441 = f439 - f438; | |
1658 f442 = f403 - f422; | |
1659 f443 = f403 + f422; | |
1660 f444 = f405 + f423; | |
1661 f445 = f423 - f405; | |
1662 f446 = f402 - f428; | |
1663 f447 = f402 + f428; | |
1664 f448 = f404 + f429; | |
1665 f449 = f404 - f429; | |
1666 f450 = f411 - f434; | |
1667 f451 = f411 + f434; | |
1668 f452 = f413 - f435; | |
1669 f453 = f413 + f435; | |
1670 f454 = f410 + f440; | |
1671 f455 = f440 - f410; | |
1672 f456 = f412 - f441; | |
1673 f457 = f412 + f441; | |
1674 f458 = f453 + f451; | |
1675 f459 = MUL_C(COEF_CONST((-0.7856949583871021)), f453); | |
1676 f460 = MUL_C(COEF_CONST(0.9807852804032304), f458); | |
1677 f461 = MUL_C(COEF_CONST(1.1758756024193588), f451); | |
1678 f462 = f459 + f460; | |
1679 f463 = f461 - f460; | |
1680 f464 = f457 + f455; | |
1681 f465 = MUL_C(COEF_CONST(0.2758993792829431), f457); | |
1682 f466 = MUL_C(COEF_CONST(0.5555702330196022), f464); | |
1683 f467 = MUL_C(COEF_CONST(1.3870398453221475), f455); | |
1684 f468 = f465 + f466; | |
1685 f469 = f467 - f466; | |
1686 f470 = f452 + f450; | |
1687 f471 = MUL_C(COEF_CONST(1.1758756024193591), f452); | |
1688 f472 = MUL_C(COEF_CONST((-0.1950903220161287)), f470); | |
1689 f473 = MUL_C(COEF_CONST(0.7856949583871016), f450); | |
1690 f474 = f471 + f472; | |
1691 f475 = f473 - f472; | |
1692 f476 = f456 - f454; | |
1693 f477 = MUL_C(COEF_CONST(1.3870398453221473), f456); | |
1694 f478 = MUL_C(COEF_CONST((-0.8314696123025455)), f476); | |
1695 f479 = MUL_C(COEF_CONST((-0.2758993792829436)), f454); | |
1696 f480 = f477 + f478; | |
1697 f481 = f478 + f479; | |
1698 f482 = f443 - f462; | |
1699 f483 = f443 + f462; | |
1700 f484 = f445 - f463; | |
1701 f485 = f445 + f463; | |
1702 f486 = f447 - f468; | |
1703 f487 = f447 + f468; | |
1704 f488 = f449 - f469; | |
1705 f489 = f449 + f469; | |
1706 f490 = f442 - f474; | |
1707 f491 = f442 + f474; | |
1708 f492 = f444 + f475; | |
1709 f493 = f475 - f444; | |
1710 f494 = f446 - f480; | |
1711 f495 = f446 + f480; | |
1712 f496 = f448 + f481; | |
1713 f497 = f448 - f481; | |
1714 f498 = f485 + f483; | |
1715 f499 = MUL_C(COEF_CONST((-0.9497277818777543)), f485); | |
1716 f500 = MUL_C(COEF_CONST(0.9987954562051724), f498); | |
1717 f501 = MUL_C(COEF_CONST(1.0478631305325905), f483); | |
1718 y[61] = f499 + f500; | |
1719 y[1] = f501 - f500; | |
1720 f504 = f489 + f487; | |
1721 f505 = MUL_C(COEF_CONST((-0.7270510732912801)), f489); | |
1722 f506 = MUL_C(COEF_CONST(0.9700312531945440), f504); | |
1723 f507 = MUL_C(COEF_CONST(1.2130114330978079), f487); | |
1724 y[53] = f505 + f506; | |
1725 y[9] = f507 - f506; | |
1726 f510 = f493 + f491; | |
1727 f511 = MUL_C(COEF_CONST((-0.4764341996931611)), f493); | |
1728 f512 = MUL_C(COEF_CONST(0.9039892931234433), f510); | |
1729 f513 = MUL_C(COEF_CONST(1.3315443865537255), f491); | |
1730 y[45] = f511 + f512; | |
1731 y[17] = f513 - f512; | |
1732 f516 = f497 + f495; | |
1733 f517 = MUL_C(COEF_CONST((-0.2075082269882114)), f497); | |
1734 f518 = MUL_C(COEF_CONST(0.8032075314806448), f516); | |
1735 f519 = MUL_C(COEF_CONST(1.3989068359730783), f495); | |
1736 y[37] = f517 + f518; | |
1737 y[25] = f519 - f518; | |
1738 f522 = f484 + f482; | |
1739 f523 = MUL_C(COEF_CONST(0.0693921705079408), f484); | |
1740 f524 = MUL_C(COEF_CONST(0.6715589548470183), f522); | |
1741 f525 = MUL_C(COEF_CONST(1.4125100802019774), f482); | |
1742 y[29] = f523 + f524; | |
1743 y[33] = f525 - f524; | |
1744 f528 = f488 + f486; | |
1745 f529 = MUL_C(COEF_CONST(0.3436258658070505), f488); | |
1746 f530 = MUL_C(COEF_CONST(0.5141027441932217), f528); | |
1747 f531 = MUL_C(COEF_CONST(1.3718313541934939), f486); | |
1748 y[21] = f529 + f530; | |
1749 y[41] = f531 - f530; | |
1750 f534 = f490 - f492; | |
1751 f535 = MUL_C(COEF_CONST(0.6046542117908007), f492); | |
1752 f536 = MUL_C(COEF_CONST(0.3368898533922201), f534); | |
1753 f537 = MUL_C(COEF_CONST(1.2784339185752409), f490); | |
1754 y[13] = f536 - f535; | |
1755 y[49] = f537 - f536; | |
1756 f540 = f496 + f494; | |
1757 f541 = MUL_C(COEF_CONST(0.8424460355094192), f496); | |
1758 f542 = MUL_C(COEF_CONST(0.1467304744553618), f540); | |
1759 f543 = MUL_C(COEF_CONST(1.1359069844201428), f494); | |
1760 y[5] = f541 + f542; | |
1761 y[57] = f543 - f542; | |
1762 f546 = f354 + f384; | |
1763 f547 = f354 - f384; | |
1764 f548 = f356 + f382; | |
1765 f549 = f382 - f356; | |
1766 f550 = f358 + f380; | |
1767 f551 = f358 - f380; | |
1768 f552 = f360 + f378; | |
1769 f553 = f378 - f360; | |
1770 f554 = f362 + f376; | |
1771 f555 = f362 - f376; | |
1772 f556 = f364 + f374; | |
1773 f557 = f374 - f364; | |
1774 f558 = f366 + f372; | |
1775 f559 = f366 - f372; | |
1776 f560 = f368 + f370; | |
1777 f561 = f370 - f368; | |
1778 f562 = f547 - f561; | |
1779 f563 = f547 + f561; | |
1780 f564 = f549 - f559; | |
1781 f565 = f549 + f559; | |
1782 f566 = f551 - f557; | |
1783 f567 = f551 + f557; | |
1784 f568 = f553 - f555; | |
1785 f569 = f553 + f555; | |
1786 f570 = f563 - f569; | |
1787 f571 = f563 + f569; | |
1788 f572 = f565 - f567; | |
1789 f573 = f565 + f567; | |
1790 f574 = f571 - f573; | |
1791 y[63] = f571 + f573; | |
1792 y[31] = MUL_C(COEF_CONST(0.7071067811865476), f574); | |
1793 f577 = f570 + f572; | |
1794 f578 = MUL_C(COEF_CONST(1.3065629648763766), f570); | |
1795 f579 = MUL_C(COEF_CONST((-0.9238795325112866)), f577); | |
1796 f580 = MUL_C(COEF_CONST((-0.5411961001461967)), f572); | |
1797 y[15] = f578 + f579; | |
1798 y[47] = f580 - f579; | |
1799 f583 = f564 + f562; | |
1800 f584 = f566 + f564; | |
1801 f585 = f568 + f566; | |
1802 f586 = MUL_C(COEF_CONST(0.7071067811865476), f584); | |
1803 f587 = f562 - f586; | |
1804 f588 = f562 + f586; | |
1805 f589 = f583 + f585; | |
1806 f590 = MUL_C(COEF_CONST(1.3065629648763766), f583); | |
1807 f591 = MUL_C(COEF_CONST((-0.9238795325112866)), f589); | |
1808 f592 = MUL_C(COEF_CONST((-0.5411961001461967)), f585); | |
1809 f593 = f590 + f591; | |
1810 f594 = f592 - f591; | |
1811 f595 = f588 - f594; | |
1812 f596 = f588 + f594; | |
1813 f597 = f587 - f593; | |
1814 f598 = f587 + f593; | |
1815 y[55] = MUL_C(COEF_CONST(0.5097955791041592), f596); | |
1816 y[39] = MUL_C(COEF_CONST(0.6013448869350453), f598); | |
1817 y[23] = MUL_C(COEF_CONST(0.8999762231364156), f597); | |
1818 y[7] = MUL_C(COEF_CONST(2.5629154477415055), f595); | |
1819 f603 = MUL_C(COEF_CONST(0.5024192861881557), f546); | |
1820 f604 = MUL_C(COEF_CONST(0.5224986149396889), f548); | |
1821 f605 = MUL_C(COEF_CONST(0.5669440348163577), f550); | |
1822 f606 = MUL_C(COEF_CONST(0.6468217833599901), f552); | |
1823 f607 = MUL_C(COEF_CONST(0.7881546234512502), f554); | |
1824 f608 = MUL_C(COEF_CONST(1.0606776859903471), f556); | |
1825 f609 = MUL_C(COEF_CONST(1.7224470982383342), f558); | |
1826 f610 = MUL_C(COEF_CONST(5.1011486186891553), f560); | |
1827 f611 = f603 + f610; | |
1828 f612 = f603 - f610; | |
1829 f613 = f604 + f609; | |
1830 f614 = f609 - f604; | |
1831 f615 = f605 + f608; | |
1832 f616 = f605 - f608; | |
1833 f617 = f606 + f607; | |
1834 f618 = f607 - f606; | |
1835 f619 = f612 - f618; | |
1836 f620 = f612 + f618; | |
1837 f621 = f614 - f616; | |
1838 f622 = f614 + f616; | |
1839 f623 = f620 - f622; | |
1840 f624 = f620 + f622; | |
1841 f625 = MUL_C(COEF_CONST(0.7071067811865476), f623); | |
1842 f626 = f619 + f621; | |
1843 f627 = MUL_C(COEF_CONST(1.3065629648763766), f619); | |
1844 f628 = MUL_C(COEF_CONST((-0.9238795325112866)), f626); | |
1845 f629 = MUL_C(COEF_CONST((-0.5411961001461967)), f621); | |
1846 f630 = f627 + f628; | |
1847 f631 = f629 - f628; | |
1848 f632 = f611 - f613; | |
1849 f633 = f615 - f613; | |
1850 f634 = f615 - f617; | |
1851 f635 = MUL_C(COEF_CONST(0.7071067811865476), f633); | |
1852 f636 = f611 - f635; | |
1853 f637 = f611 + f635; | |
1854 f638 = f632 + f634; | |
1855 f639 = MUL_C(COEF_CONST(1.3065629648763766), f632); | |
1856 f640 = MUL_C(COEF_CONST((-0.9238795325112866)), f638); | |
1857 f641 = MUL_C(COEF_CONST((-0.5411961001461967)), f634); | |
1858 f642 = f639 + f640; | |
1859 f643 = f641 - f640; | |
1860 f644 = f637 - f643; | |
1861 f645 = f637 + f643; | |
1862 f646 = f636 - f642; | |
1863 f647 = f636 + f642; | |
1864 f648 = MUL_C(COEF_CONST(0.5097955791041592), f645); | |
1865 f649 = MUL_C(COEF_CONST(0.6013448869350453), f647); | |
1866 f650 = MUL_C(COEF_CONST(0.8999762231364156), f646); | |
1867 y[3] = MUL_C(COEF_CONST(2.5629154477415055), f644); | |
1868 y[59] = f624 + f648; | |
1869 y[51] = f648 + f631; | |
1870 y[43] = f631 + f649; | |
1871 y[35] = f649 + f625; | |
1872 y[27] = f625 + f650; | |
1873 y[19] = f650 + f630; | |
1874 y[11] = f630 + y[3]; | |
1095 } | 1875 } |
1096 | 1876 |
1097 #else | 1877 #else |
1098 | 1878 |
1099 void DCT4_64(real_t *y, real_t *x) | 1879 void DCT4_64(real_t *y, real_t *x) |
1100 { | 1880 { |
1101 int16_t i0; | 1881 int16_t i0; |
1102 real_t f2; | 1882 ALIGN static real_t t2[64]; |
1103 real_t f3; | |
1104 real_t f4; | |
1105 real_t f5; | |
1106 real_t f6; | |
1107 real_t f7; | |
1108 real_t f8; | |
1109 real_t f9; | |
1110 real_t f10; | |
1111 real_t f11; | |
1112 real_t f12; | |
1113 real_t f13; | |
1114 real_t f14; | |
1115 real_t f15; | |
1116 real_t f16; | |
1117 real_t f17; | |
1118 real_t f18; | |
1119 real_t f19; | |
1120 real_t f20; | |
1121 real_t f21; | |
1122 real_t f22; | |
1123 real_t f23; | |
1124 real_t f24; | |
1125 real_t f25; | |
1126 real_t f26; | |
1127 real_t f27; | |
1128 real_t f28; | |
1129 real_t f29; | |
1130 real_t f30; | |
1131 real_t f31; | |
1132 real_t f32; | |
1133 real_t f33; | |
1134 real_t f34; | |
1135 real_t f35; | |
1136 real_t f36; | |
1137 real_t f37; | |
1138 real_t f38; | |
1139 real_t f39; | |
1140 real_t f40; | |
1141 real_t f41; | |
1142 real_t f42; | |
1143 real_t f43; | |
1144 real_t f44; | |
1145 real_t f45; | |
1146 real_t f46; | |
1147 real_t f47; | |
1148 real_t f48; | |
1149 real_t f49; | |
1150 real_t f50; | |
1151 real_t f51; | |
1152 real_t f52; | |
1153 real_t f53; | |
1154 real_t f54; | |
1155 real_t f55; | |
1156 real_t f56; | |
1157 real_t f57; | |
1158 real_t f58; | |
1159 real_t f59; | |
1160 real_t f60; | |
1161 real_t f61; | |
1162 real_t f62; | |
1163 real_t f63; | |
1164 real_t f64; | |
1165 real_t f65; | |
1166 real_t f66; | |
1167 real_t f67; | |
1168 real_t f68; | |
1169 real_t f69; | |
1170 real_t f70; | |
1171 real_t f71; | |
1172 real_t f72; | |
1173 real_t f73; | |
1174 real_t f74; | |
1175 real_t f75; | |
1176 real_t f76; | |
1177 real_t f77; | |
1178 real_t f78; | |
1179 real_t f79; | |
1180 real_t f80; | |
1181 real_t f81; | |
1182 real_t f82; | |
1183 real_t f83; | |
1184 real_t f84; | |
1185 real_t f85; | |
1186 real_t f86; | |
1187 real_t f87; | |
1188 real_t f88; | |
1189 real_t f89; | |
1190 real_t f90; | |
1191 real_t f91; | |
1192 real_t f92; | |
1193 real_t f93; | |
1194 real_t f94; | |
1195 real_t f95; | |
1196 real_t f96; | |
1197 real_t f97; | |
1198 real_t f98; | |
1199 real_t f99; | |
1200 real_t f100; | |
1201 real_t f101; | |
1202 real_t f102; | |
1203 real_t f103; | |
1204 real_t f104; | |
1205 real_t f105; | |
1206 real_t f106; | |
1207 real_t f107; | |
1208 real_t f108; | |
1209 real_t f109; | |
1210 real_t f110; | |
1211 real_t f111; | |
1212 real_t f112; | |
1213 real_t f113; | |
1214 real_t f114; | |
1215 real_t f115; | |
1216 real_t f116; | |
1217 real_t f117; | |
1218 real_t f118; | |
1219 real_t f119; | |
1220 real_t f120; | |
1221 real_t f121; | |
1222 real_t f122; | |
1223 real_t f123; | |
1224 real_t f124; | |
1225 real_t f125; | |
1226 real_t f126; | |
1227 real_t f127; | |
1228 real_t f128; | |
1229 real_t f129; | |
1230 real_t f130; | |
1231 real_t f131; | |
1232 real_t f132; | |
1233 real_t f133; | |
1234 real_t f134; | |
1235 real_t f135; | |
1236 real_t f136; | |
1237 real_t f137; | |
1238 real_t f138; | |
1239 real_t f139; | |
1240 real_t f140; | |
1241 real_t f141; | |
1242 real_t f142; | |
1243 real_t f143; | |
1244 real_t f144; | |
1245 real_t f145; | |
1246 real_t f146; | |
1247 real_t f147; | |
1248 real_t f148; | |
1249 real_t f149; | |
1250 real_t f150; | |
1251 real_t f151; | |
1252 real_t f152; | |
1253 real_t f153; | |
1254 real_t f154; | |
1255 real_t f155; | |
1256 real_t f156; | |
1257 real_t f157; | |
1258 real_t f158; | |
1259 real_t f159; | |
1260 real_t f160; | |
1261 real_t f161; | |
1262 real_t f162; | |
1263 real_t f163; | |
1264 real_t f164; | |
1265 real_t f165; | |
1266 real_t f166; | |
1267 real_t f167; | |
1268 real_t f168; | |
1269 real_t f169; | |
1270 real_t f170; | |
1271 real_t f171; | |
1272 real_t f172; | |
1273 real_t f173; | |
1274 real_t f174; | |
1275 real_t f175; | |
1276 real_t f176; | |
1277 real_t f177; | |
1278 real_t f178; | |
1279 real_t f179; | |
1280 real_t f180; | |
1281 real_t f181; | |
1282 real_t f182; | |
1283 real_t f183; | |
1284 real_t f184; | |
1285 real_t f185; | |
1286 real_t f186; | |
1287 real_t f187; | |
1288 real_t f188; | |
1289 real_t f189; | |
1290 real_t f190; | |
1291 real_t f191; | |
1292 real_t f192; | |
1293 real_t f193; | |
1294 real_t f194; | |
1295 real_t f195; | |
1296 real_t f196; | |
1297 real_t f197; | |
1298 real_t f198; | |
1299 real_t f199; | |
1300 real_t f200; | |
1301 real_t f201; | |
1302 real_t f202; | |
1303 real_t f203; | |
1304 real_t f204; | |
1305 real_t f205; | |
1306 real_t f206; | |
1307 real_t f207; | |
1308 real_t f208; | |
1309 real_t f209; | |
1310 real_t f210; | |
1311 real_t f211; | |
1312 real_t f212; | |
1313 real_t f213; | |
1314 real_t f214; | |
1315 real_t f215; | |
1316 real_t f216; | |
1317 real_t f217; | |
1318 real_t f218; | |
1319 real_t f219; | |
1320 real_t f220; | |
1321 real_t f221; | |
1322 real_t f222; | |
1323 real_t f223; | |
1324 real_t f224; | |
1325 real_t f225; | |
1326 real_t f226; | |
1327 real_t f227; | |
1328 real_t f228; | |
1329 real_t f229; | |
1330 real_t f230; | |
1331 real_t f231; | |
1332 real_t f232; | |
1333 real_t f233; | |
1334 real_t f234; | |
1335 real_t f235; | |
1336 real_t f236; | |
1337 real_t f237; | |
1338 real_t f238; | |
1339 real_t f239; | |
1340 real_t f240; | |
1341 real_t f241; | |
1342 real_t f242; | |
1343 real_t f243; | |
1344 real_t f244; | |
1345 real_t f245; | |
1346 real_t f246; | |
1347 real_t f247; | |
1348 real_t f248; | |
1349 real_t f249; | |
1350 real_t f250; | |
1351 real_t f251; | |
1352 real_t f252; | |
1353 real_t f253; | |
1354 real_t f254; | |
1355 real_t f255; | |
1356 real_t f256; | |
1357 real_t f257; | |
1358 real_t f258; | |
1359 real_t f259; | |
1360 real_t f260; | |
1361 real_t f261; | |
1362 real_t f262; | |
1363 real_t f263; | |
1364 real_t f264; | |
1365 real_t f265; | |
1366 real_t f266; | |
1367 real_t f267; | |
1368 real_t f268; | |
1369 real_t f269; | |
1370 real_t f270; | |
1371 real_t f271; | |
1372 real_t f272; | |
1373 real_t f273; | |
1374 real_t f274; | |
1375 real_t f275; | |
1376 real_t f276; | |
1377 real_t f277; | |
1378 real_t f278; | |
1379 real_t f279; | |
1380 real_t f280; | |
1381 real_t f281; | |
1382 real_t f282; | |
1383 real_t f283; | |
1384 real_t f284; | |
1385 real_t f285; | |
1386 real_t f286; | |
1387 real_t f287; | |
1388 real_t f288; | |
1389 real_t f289; | |
1390 real_t f290; | |
1391 real_t f291; | |
1392 real_t f292; | |
1393 real_t f293; | |
1394 real_t f294; | |
1395 real_t f295; | |
1396 real_t f296; | |
1397 real_t f297; | |
1398 real_t f298; | |
1399 real_t f299; | |
1400 real_t f300; | |
1401 real_t f301; | |
1402 real_t f302; | |
1403 real_t f303; | |
1404 real_t f304; | |
1405 real_t f305; | |
1406 real_t f306; | |
1407 real_t f307; | |
1408 real_t f308; | |
1409 real_t f309; | |
1410 real_t f310; | |
1411 real_t f311; | |
1412 real_t f312; | |
1413 real_t f313; | |
1414 real_t f314; | |
1415 real_t f315; | |
1416 real_t f316; | |
1417 real_t f317; | |
1418 real_t f318; | |
1419 real_t f319; | |
1420 real_t f320; | |
1421 real_t f321; | |
1422 real_t f322; | |
1423 real_t f323; | |
1424 real_t f324; | |
1425 real_t f325; | |
1426 real_t f326; | |
1427 real_t f327; | |
1428 real_t f328; | |
1429 real_t f329; | |
1430 real_t f330; | |
1431 real_t f331; | |
1432 real_t f332; | |
1433 real_t f333; | |
1434 real_t f334; | |
1435 real_t f335; | |
1436 real_t f336; | |
1437 real_t f337; | |
1438 real_t f338; | |
1439 real_t f339; | |
1440 real_t f340; | |
1441 real_t f341; | |
1442 real_t f342; | |
1443 real_t f343; | |
1444 real_t f344; | |
1445 real_t f345; | |
1446 real_t f346; | |
1447 real_t f347; | |
1448 real_t f348; | |
1449 real_t f349; | |
1450 real_t f350; | |
1451 real_t f351; | |
1452 real_t f352; | |
1453 real_t f353; | |
1454 real_t f354; | |
1455 real_t f355; | |
1456 real_t f356; | |
1457 real_t f357; | |
1458 real_t f358; | |
1459 real_t f359; | |
1460 real_t f360; | |
1461 real_t f361; | |
1462 real_t f362; | |
1463 real_t f363; | |
1464 real_t f364; | |
1465 real_t f365; | |
1466 real_t f366; | |
1467 real_t f367; | |
1468 real_t f368; | |
1469 real_t f369; | |
1470 real_t f370; | |
1471 real_t f371; | |
1472 real_t f372; | |
1473 real_t f373; | |
1474 real_t f374; | |
1475 real_t f375; | |
1476 real_t f376; | |
1477 real_t f377; | |
1478 real_t f378; | |
1479 real_t f379; | |
1480 real_t f380; | |
1481 real_t f381; | |
1482 real_t f382; | |
1483 real_t f383; | |
1484 real_t f384; | |
1485 real_t f385; | |
1486 real_t f386; | |
1487 real_t f387; | |
1488 real_t f388; | |
1489 real_t f389; | |
1490 real_t f390; | |
1491 real_t f391; | |
1492 real_t f392; | |
1493 real_t f393; | |
1494 real_t f394; | |
1495 real_t f395; | |
1496 real_t f396; | |
1497 real_t f397; | |
1498 real_t f398; | |
1499 real_t f399; | |
1500 real_t f400; | |
1501 real_t f401; | |
1502 real_t f402; | |
1503 real_t f403; | |
1504 real_t f404; | |
1505 real_t f405; | |
1506 real_t f406; | |
1507 real_t f407; | |
1508 real_t f408; | |
1509 real_t f409; | |
1510 real_t f410; | |
1511 real_t f411; | |
1512 real_t f412; | |
1513 real_t f413; | |
1514 real_t f414; | |
1515 real_t f415; | |
1516 real_t f416; | |
1517 real_t f417; | |
1518 real_t f418; | |
1519 real_t f419; | |
1520 real_t f420; | |
1521 real_t f421; | |
1522 real_t f422; | |
1523 real_t f423; | |
1524 real_t f424; | |
1525 real_t f425; | |
1526 real_t f426; | |
1527 real_t f427; | |
1528 real_t f428; | |
1529 real_t f429; | |
1530 real_t f430; | |
1531 real_t f431; | |
1532 real_t f432; | |
1533 real_t f433; | |
1534 real_t f434; | |
1535 real_t f435; | |
1536 real_t f436; | |
1537 real_t f437; | |
1538 real_t f438; | |
1539 real_t f439; | |
1540 real_t f440; | |
1541 real_t f441; | |
1542 real_t f442; | |
1543 real_t f443; | |
1544 real_t f444; | |
1545 real_t f445; | |
1546 real_t f446; | |
1547 real_t f447; | |
1548 real_t f448; | |
1549 real_t f449; | |
1550 real_t f450; | |
1551 real_t f451; | |
1552 real_t f452; | |
1553 real_t f453; | |
1554 real_t f454; | |
1555 real_t f455; | |
1556 real_t f456; | |
1557 real_t f457; | |
1558 real_t f458; | |
1559 real_t f459; | |
1560 real_t f460; | |
1561 real_t f461; | |
1562 real_t f462; | |
1563 real_t f463; | |
1564 real_t f464; | |
1565 real_t f465; | |
1566 real_t f466; | |
1567 real_t f467; | |
1568 real_t f468; | |
1569 real_t f469; | |
1570 real_t f470; | |
1571 real_t f471; | |
1572 real_t f472; | |
1573 real_t f473; | |
1574 real_t f474; | |
1575 real_t f475; | |
1576 real_t f476; | |
1577 real_t f477; | |
1578 real_t f478; | |
1579 real_t f479; | |
1580 real_t f480; | |
1581 real_t f481; | |
1582 real_t f482; | |
1583 real_t f483; | |
1584 real_t f484; | |
1585 real_t f485; | |
1586 real_t f486; | |
1587 real_t f487; | |
1588 real_t f488; | |
1589 real_t f489; | |
1590 real_t f490; | |
1591 real_t f491; | |
1592 real_t f492; | |
1593 real_t f493; | |
1594 real_t f494; | |
1595 real_t f495; | |
1596 real_t f496; | |
1597 real_t f497; | |
1598 real_t f498; | |
1599 real_t f499; | |
1600 real_t f500; | |
1601 real_t f501; | |
1602 real_t f502; | |
1603 real_t f503; | |
1604 real_t f504; | |
1605 real_t f505; | |
1606 real_t f506; | |
1607 real_t f507; | |
1608 real_t f508; | |
1609 real_t f509; | |
1610 real_t f510; | |
1611 real_t f511; | |
1612 real_t f512; | |
1613 real_t f513; | |
1614 real_t f514; | |
1615 real_t f515; | |
1616 real_t f516; | |
1617 real_t f517; | |
1618 real_t f518; | |
1619 real_t f519; | |
1620 real_t f520; | |
1621 real_t f521; | |
1622 real_t f522; | |
1623 real_t f523; | |
1624 real_t f524; | |
1625 real_t f525; | |
1626 real_t f526; | |
1627 real_t f527; | |
1628 real_t f528; | |
1629 real_t f529; | |
1630 real_t f530; | |
1631 real_t f531; | |
1632 real_t f532; | |
1633 real_t f533; | |
1634 real_t f534; | |
1635 real_t f535; | |
1636 real_t f536; | |
1637 real_t f537; | |
1638 real_t f538; | |
1639 real_t f539; | |
1640 real_t f540; | |
1641 real_t f541; | |
1642 real_t f542; | |
1643 real_t f543; | |
1644 real_t f544; | |
1645 real_t f545; | |
1646 real_t f546; | |
1647 real_t f547; | |
1648 real_t f548; | |
1649 real_t f549; | |
1650 real_t f550; | |
1651 real_t f551; | |
1652 real_t f552; | |
1653 real_t f553; | |
1654 real_t f554; | |
1655 real_t f555; | |
1656 real_t f556; | |
1657 real_t f557; | |
1658 real_t f558; | |
1659 real_t f559; | |
1660 real_t f560; | |
1661 real_t f561; | |
1662 real_t f562; | |
1663 real_t f563; | |
1664 real_t f564; | |
1665 real_t f565; | |
1666 real_t f566; | |
1667 real_t f567; | |
1668 real_t f568; | |
1669 real_t f569; | |
1670 real_t f570; | |
1671 real_t f571; | |
1672 real_t f572; | |
1673 real_t f573; | |
1674 real_t f574; | |
1675 real_t f575; | |
1676 real_t f576; | |
1677 real_t f577; | |
1678 real_t f578; | |
1679 real_t f579; | |
1680 real_t f580; | |
1681 real_t f581; | |
1682 real_t f582; | |
1683 real_t f583; | |
1684 real_t f584; | |
1685 real_t f585; | |
1686 real_t f586; | |
1687 real_t f587; | |
1688 real_t f588; | |
1689 real_t f589; | |
1690 real_t f590; | |
1691 real_t f591; | |
1692 real_t f592; | |
1693 real_t f593; | |
1694 real_t f594; | |
1695 real_t f595; | |
1696 real_t f596; | |
1697 real_t f597; | |
1698 real_t f598; | |
1699 real_t f599; | |
1700 real_t f600; | |
1701 real_t f601; | |
1702 real_t f602; | |
1703 real_t f603; | |
1704 real_t f604; | |
1705 real_t f605; | |
1706 real_t f606; | |
1707 real_t f607; | |
1708 real_t f608; | |
1709 real_t f609; | |
1710 real_t f610; | |
1711 real_t f611; | |
1712 real_t f612; | |
1713 real_t f613; | |
1714 real_t f614; | |
1715 real_t f615; | |
1716 real_t f618; | |
1717 real_t f619; | |
1718 real_t f620; | |
1719 real_t f621; | |
1720 real_t f624; | |
1721 real_t f625; | |
1722 real_t f626; | |
1723 real_t f627; | |
1724 real_t f630; | |
1725 real_t f631; | |
1726 real_t f632; | |
1727 real_t f633; | |
1728 real_t f636; | |
1729 real_t f637; | |
1730 real_t f638; | |
1731 real_t f639; | |
1732 real_t f642; | |
1733 real_t f643; | |
1734 real_t f644; | |
1735 real_t f645; | |
1736 real_t f648; | |
1737 real_t f649; | |
1738 real_t f650; | |
1739 real_t f651; | |
1740 real_t f654; | |
1741 real_t f655; | |
1742 real_t f656; | |
1743 real_t f657; | |
1744 real_t f660; | |
1745 real_t f661; | |
1746 real_t f662; | |
1747 real_t f663; | |
1748 real_t f666; | |
1749 real_t f667; | |
1750 real_t f668; | |
1751 real_t f669; | |
1752 real_t f672; | |
1753 real_t f673; | |
1754 real_t f674; | |
1755 real_t f675; | |
1756 real_t f678; | |
1757 real_t f679; | |
1758 real_t f680; | |
1759 real_t f681; | |
1760 real_t f684; | |
1761 real_t f685; | |
1762 real_t f686; | |
1763 real_t f687; | |
1764 real_t f690; | |
1765 real_t f691; | |
1766 real_t f692; | |
1767 real_t f693; | |
1768 real_t f696; | |
1769 real_t f697; | |
1770 real_t f698; | |
1771 real_t f699; | |
1772 real_t f702; | |
1773 real_t f703; | |
1774 real_t f704; | |
1775 real_t f705; | |
1776 real_t f708; | |
1777 real_t f709; | |
1778 real_t f710; | |
1779 real_t f711; | |
1780 real_t f714; | |
1781 real_t f715; | |
1782 real_t f716; | |
1783 real_t f717; | |
1784 real_t f720; | |
1785 real_t f721; | |
1786 real_t f722; | |
1787 real_t f723; | |
1788 real_t f726; | |
1789 real_t f727; | |
1790 real_t f728; | |
1791 real_t f729; | |
1792 real_t f732; | |
1793 real_t f733; | |
1794 real_t f734; | |
1795 real_t f735; | |
1796 real_t f738; | |
1797 real_t f739; | |
1798 real_t f740; | |
1799 real_t f741; | |
1800 real_t f744; | |
1801 real_t f745; | |
1802 real_t f746; | |
1803 real_t f747; | |
1804 real_t f750; | |
1805 real_t f751; | |
1806 real_t f752; | |
1807 real_t f753; | |
1808 real_t f756; | |
1809 real_t f757; | |
1810 real_t f758; | |
1811 real_t f759; | |
1812 real_t f762; | |
1813 real_t f763; | |
1814 real_t f764; | |
1815 real_t f765; | |
1816 real_t f768; | |
1817 real_t f769; | |
1818 real_t f770; | |
1819 real_t f771; | |
1820 real_t f774; | |
1821 real_t f775; | |
1822 real_t f776; | |
1823 real_t f777; | |
1824 real_t f780; | |
1825 real_t f781; | |
1826 real_t f782; | |
1827 real_t f783; | |
1828 real_t f786; | |
1829 real_t f787; | |
1830 real_t f788; | |
1831 real_t f789; | |
1832 real_t f792; | |
1833 real_t f793; | |
1834 real_t f794; | |
1835 real_t f795; | |
1836 real_t f798; | |
1837 real_t f799; | |
1838 real_t f800; | |
1839 real_t f801; | |
1840 static real_t t2[64]; | |
1841 | 1883 |
1842 t2[0] = x[0]; | 1884 t2[0] = x[0]; |
1843 for (i0=0; i0<31; i0++) | 1885 for (i0=0; i0<31; i0++) |
1844 { | 1886 { |
1845 t2[2*i0+1] = x[2*i0+1] - x[2*i0+2]; | 1887 t2[2*i0+1] = x[2*i0+1] - x[2*i0+2]; |
1846 t2[2*i0+2] = x[2*i0+1] + x[2*i0+2]; | 1888 t2[2*i0+2] = x[2*i0+1] + x[2*i0+2]; |
1847 } | 1889 } |
1848 t2[63] = x[63]; | 1890 t2[63] = x[63]; |
1891 | |
1892 DCT4_64_kernel(y, t2); | |
1893 } | |
1894 | |
1895 void DCT4_64_kernel(real_t *y, real_t *t2) | |
1896 { | |
1897 real_t f2, f3, f4, f5, f6, f7, f8; | |
1898 real_t f9, f10, f11, f12, f13, f14, f15; | |
1899 real_t f16, f17, f18, f19, f20, f21, f22; | |
1900 real_t f23, f24, f25, f26, f27, f28, f29; | |
1901 real_t f30, f31, f32, f33, f34, f35, f36; | |
1902 real_t f37, f38, f39, f40, f41, f42, f43; | |
1903 real_t f44, f45, f46, f47, f48, f49, f50; | |
1904 real_t f51, f52, f53, f54, f55, f56, f57; | |
1905 real_t f58, f59, f60, f61, f62, f63, f64; | |
1906 real_t f65, f66, f67, f68, f69, f70, f71; | |
1907 real_t f72, f73, f74, f75, f76, f77, f78; | |
1908 real_t f79, f80, f81, f82, f83, f84, f85; | |
1909 real_t f86, f87, f88, f89, f90, f91, f92; | |
1910 real_t f93, f94, f95, f96, f97, f98, f99; | |
1911 real_t f100, f101, f102, f103, f104, f105, f106; | |
1912 real_t f107, f108, f109, f110, f111, f112, f113; | |
1913 real_t f114, f115, f116, f117, f118, f119, f120; | |
1914 real_t f121, f122, f123, f124, f125, f126, f127; | |
1915 real_t f128, f129, f130, f131, f132, f133, f134; | |
1916 real_t f135, f136, f137, f138, f139, f140, f141; | |
1917 real_t f142, f143, f144, f145, f146, f147, f148; | |
1918 real_t f149, f150, f151, f152, f153, f154, f155; | |
1919 real_t f156, f157, f158, f159, f160, f161, f162; | |
1920 real_t f163, f164, f165, f166, f167, f168, f169; | |
1921 real_t f170, f171, f172, f173, f174, f175, f176; | |
1922 real_t f177, f178, f179, f180, f181, f182, f183; | |
1923 real_t f184, f185, f186, f187, f188, f189, f190; | |
1924 real_t f191, f192, f193, f194, f195, f196, f197; | |
1925 real_t f198, f199, f200, f201, f202, f203, f204; | |
1926 real_t f205, f206, f207, f208, f209, f210, f211; | |
1927 real_t f212, f213, f214, f215, f216, f217, f218; | |
1928 real_t f219, f220, f221, f222, f223, f224, f225; | |
1929 real_t f226, f227, f228, f229, f230, f231, f232; | |
1930 real_t f233, f234, f235, f236, f237, f238, f239; | |
1931 real_t f240, f241, f242, f243, f244, f245, f246; | |
1932 real_t f247, f248, f249, f250, f251, f252, f253; | |
1933 real_t f254, f255, f256, f257, f258, f259, f260; | |
1934 real_t f261, f262, f263, f264, f265, f266, f267; | |
1935 real_t f268, f269, f270, f271, f272, f273, f274; | |
1936 real_t f275, f276, f277, f278, f279, f280, f281; | |
1937 real_t f282, f283, f284, f285, f286, f287, f288; | |
1938 real_t f289, f290, f291, f292, f293, f294, f295; | |
1939 real_t f296, f297, f298, f299, f300, f301, f302; | |
1940 real_t f303, f304, f305, f306, f307, f308, f309; | |
1941 real_t f310, f311, f312, f313, f314, f315, f316; | |
1942 real_t f317, f318, f319, f320, f321, f322, f323; | |
1943 real_t f324, f325, f326, f327, f328, f329, f330; | |
1944 real_t f331, f332, f333, f334, f335, f336, f337; | |
1945 real_t f338, f339, f340, f341, f342, f343, f344; | |
1946 real_t f345, f346, f347, f348, f349, f350, f351; | |
1947 real_t f352, f353, f354, f355, f356, f357, f358; | |
1948 real_t f359, f360, f361, f362, f363, f364, f365; | |
1949 real_t f366, f367, f368, f369, f370, f371, f372; | |
1950 real_t f373, f374, f375, f376, f377, f378, f379; | |
1951 real_t f380, f381, f382, f383, f384, f385, f386; | |
1952 real_t f387, f388, f389, f390, f391, f392, f393; | |
1953 real_t f394, f395, f396, f397, f398, f399, f400; | |
1954 real_t f401, f402, f403, f404, f405, f406, f407; | |
1955 real_t f408, f409, f410, f411, f412, f413, f414; | |
1956 real_t f415, f416, f417, f418, f419, f420, f421; | |
1957 real_t f422, f423, f424, f425, f426, f427, f428; | |
1958 real_t f429, f430, f431, f432, f433, f434, f435; | |
1959 real_t f436, f437, f438, f439, f440, f441, f442; | |
1960 real_t f443, f444, f445, f446, f447, f448, f449; | |
1961 real_t f450, f451, f452, f453, f454, f455, f456; | |
1962 real_t f457, f458, f459, f460, f461, f462, f463; | |
1963 real_t f464, f465, f466, f467, f468, f469, f470; | |
1964 real_t f471, f472, f473, f474, f475, f476, f477; | |
1965 real_t f478, f479, f480, f481, f482, f483, f484; | |
1966 real_t f485, f486, f487, f488, f489, f490, f491; | |
1967 real_t f492, f493, f494, f495, f496, f497, f498; | |
1968 real_t f499, f500, f501, f502, f503, f504, f505; | |
1969 real_t f506, f507, f508, f509, f510, f511, f512; | |
1970 real_t f513, f514, f515, f516, f517, f518, f519; | |
1971 real_t f520, f521, f522, f523, f524, f525, f526; | |
1972 real_t f527, f528, f529, f530, f531, f532, f533; | |
1973 real_t f534, f535, f536, f537, f538, f539, f540; | |
1974 real_t f541, f542, f543, f544, f545, f546, f547; | |
1975 real_t f548, f549, f550, f551, f552, f553, f554; | |
1976 real_t f555, f556, f557, f558, f559, f560, f561; | |
1977 real_t f562, f563, f564, f565, f566, f567, f568; | |
1978 real_t f569, f570, f571, f572, f573, f574, f575; | |
1979 real_t f576, f577, f578, f579, f580, f581, f582; | |
1980 real_t f583, f584, f585, f586, f587, f588, f589; | |
1981 real_t f590, f591, f592, f593, f594, f595, f596; | |
1982 real_t f597, f598, f599, f600, f601, f602, f603; | |
1983 real_t f604, f605, f606, f607, f608, f609, f610; | |
1984 real_t f611, f612, f613, f614, f615, f618, f619; | |
1985 real_t f620, f621, f624, f625, f626, f627, f630; | |
1986 real_t f631, f632, f633, f636, f637, f638, f639; | |
1987 real_t f642, f643, f644, f645, f648, f649, f650; | |
1988 real_t f651, f654, f655, f656, f657, f660, f661; | |
1989 real_t f662, f663, f666, f667, f668, f669, f672; | |
1990 real_t f673, f674, f675, f678, f679, f680, f681; | |
1991 real_t f684, f685, f686, f687, f690, f691, f692; | |
1992 real_t f693, f696, f697, f698, f699, f702, f703; | |
1993 real_t f704, f705, f708, f709, f710, f711, f714; | |
1994 real_t f715, f716, f717, f720, f721, f722, f723; | |
1995 real_t f726, f727, f728, f729, f732, f733, f734; | |
1996 real_t f735, f738, f739, f740, f741, f744, f745; | |
1997 real_t f746, f747, f750, f751, f752, f753, f756; | |
1998 real_t f757, f758, f759, f762, f763, f764, f765; | |
1999 real_t f768, f769, f770, f771, f774, f775, f776; | |
2000 real_t f777, f780, f781, f782, f783, f786, f787; | |
2001 real_t f788, f789, f792, f793, f794, f795, f798; | |
2002 real_t f799, f800, f801; | |
2003 | |
1849 f2 = 0.7071067811865476 * t2[32]; | 2004 f2 = 0.7071067811865476 * t2[32]; |
1850 f3 = x[0] - f2; | 2005 f3 = t2[0] - f2; |
1851 f4 = x[0] + f2; | 2006 f4 = t2[0] + f2; |
1852 f5 = t2[16] + t2[48]; | 2007 f5 = t2[16] + t2[48]; |
1853 f6 = 1.3065629648763766 * t2[16]; | 2008 f6 = 1.3065629648763766 * t2[16]; |
1854 f7 = (-0.9238795325112866) * f5; | 2009 f7 = (-0.9238795325112866) * f5; |
1855 f8 = (-0.5411961001461967) * t2[48]; | 2010 f8 = (-0.5411961001461967) * t2[48]; |
1856 f9 = f6 + f7; | 2011 f9 = f6 + f7; |
2310 f463 = f447 - f449; | 2465 f463 = f447 - f449; |
2311 f464 = f447 + f449; | 2466 f464 = f447 + f449; |
2312 f465 = 0.7071067811865474 * f463; | 2467 f465 = 0.7071067811865474 * f463; |
2313 f466 = 0.7071067811865474 * f464; | 2468 f466 = 0.7071067811865474 * f464; |
2314 f467 = 0.7071067811865476 * t2[31]; | 2469 f467 = 0.7071067811865476 * t2[31]; |
2315 f468 = x[63] - f467; | 2470 f468 = t2[63] - f467; |
2316 f469 = x[63] + f467; | 2471 f469 = t2[63] + f467; |
2317 f470 = t2[47] + t2[15]; | 2472 f470 = t2[47] + t2[15]; |
2318 f471 = 1.3065629648763766 * t2[47]; | 2473 f471 = 1.3065629648763766 * t2[47]; |
2319 f472 = (-0.9238795325112866) * f470; | 2474 f472 = (-0.9238795325112866) * f470; |
2320 f473 = (-0.5411961001461967) * t2[15]; | 2475 f473 = (-0.5411961001461967) * t2[15]; |
2321 f474 = f471 + f472; | 2476 f474 = f471 + f472; |