FFmpeg
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
mpegaudiodsp_mips_float.c
Go to the documentation of this file.
1 /*
2  * Copyright (c) 2012
3  * MIPS Technologies, Inc., California.
4  *
5  * Redistribution and use in source and binary forms, with or without
6  * modification, are permitted provided that the following conditions
7  * are met:
8  * 1. Redistributions of source code must retain the above copyright
9  * notice, this list of conditions and the following disclaimer.
10  * 2. Redistributions in binary form must reproduce the above copyright
11  * notice, this list of conditions and the following disclaimer in the
12  * documentation and/or other materials provided with the distribution.
13  * 3. Neither the name of the MIPS Technologies, Inc., nor the names of its
14  * contributors may be used to endorse or promote products derived from
15  * this software without specific prior written permission.
16  *
17  * THIS SOFTWARE IS PROVIDED BY THE MIPS TECHNOLOGIES, INC. ``AS IS'' AND
18  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
19  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
20  * ARE DISCLAIMED. IN NO EVENT SHALL THE MIPS TECHNOLOGIES, INC. BE LIABLE
21  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
22  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
23  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
24  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
25  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
26  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
27  * SUCH DAMAGE.
28  *
29  * Author: Bojan Zivkovic (bojan@mips.com)
30  *
31  * MPEG Audio decoder optimized for MIPS floating-point architecture
32  *
33  * This file is part of FFmpeg.
34  *
35  * FFmpeg is free software; you can redistribute it and/or
36  * modify it under the terms of the GNU Lesser General Public
37  * License as published by the Free Software Foundation; either
38  * version 2.1 of the License, or (at your option) any later version.
39  *
40  * FFmpeg is distributed in the hope that it will be useful,
41  * but WITHOUT ANY WARRANTY; without even the implied warranty of
42  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
43  * Lesser General Public License for more details.
44  *
45  * You should have received a copy of the GNU Lesser General Public
46  * License along with FFmpeg; if not, write to the Free Software
47  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
48  */
49 
50 /**
51  * @file
52  * Reference: libavcodec/mpegaudiodsp_template.c
53  * libavcodec/dct32.c
54  */
55 
56 #include <string.h>
57 
58 #include "libavutil/mips/asmdefs.h"
59 #include "libavcodec/mpegaudiodsp.h"
60 
61 static void ff_mpadsp_apply_window_mips_float(float *synth_buf, float *window,
62  int *dither_state, float *samples, int incr)
63 {
64  register const float *w, *w2, *p;
65  int j;
66  float *samples2;
67  float sum, sum2;
68  /* temporary variables */
69  int incr1 = incr << 2;
70  int t_sample;
71  float in1, in2, in3, in4, in5, in6, in7, in8;
72  float *p2;
73 
74  /* copy to avoid wrap */
75  memcpy(synth_buf + 512, synth_buf, 32 * sizeof(*synth_buf));
76 
77  /**
78  * instructions are scheduled to minimize pipeline stall.
79  * use of round_sample function from the original code is
80  * changed with appropriate assembly instructions.
81  */
82 
83  __asm__ volatile (
84  "lwc1 %[sum], 0(%[dither_state]) \t\n"
85  "sll %[t_sample], %[incr1], 5 \t\n"
86  "sub %[t_sample], %[t_sample], %[incr1] \n\t"
87  "li %[j], 4 \t\n"
88  "lwc1 %[in1], 0(%[window]) \t\n"
89  "lwc1 %[in2], 16*4(%[synth_buf]) \t\n"
90  "sw $zero, 0(%[dither_state]) \t\n"
91  "lwc1 %[in3], 64*4(%[window]) \t\n"
92  "lwc1 %[in4], 80*4(%[synth_buf]) \t\n"
93  PTR_ADDU "%[samples2],%[samples], %[t_sample] \t\n"
94  "madd.s %[sum], %[sum], %[in1], %[in2] \t\n"
95  "lwc1 %[in5], 128*4(%[window]) \t\n"
96  "lwc1 %[in6], 144*4(%[synth_buf]) \t\n"
97  "lwc1 %[in7], 192*4(%[window]) \t\n"
98  "madd.s %[sum], %[sum], %[in3], %[in4] \t\n"
99  "lwc1 %[in8], 208*4(%[synth_buf]) \t\n"
100  "lwc1 %[in1], 256*4(%[window]) \t\n"
101  "lwc1 %[in2], 272*4(%[synth_buf]) \t\n"
102  "madd.s %[sum], %[sum], %[in5], %[in6] \t\n"
103  "lwc1 %[in3], 320*4(%[window]) \t\n"
104  "lwc1 %[in4], 336*4(%[synth_buf]) \t\n"
105  "lwc1 %[in5], 384*4(%[window]) \t\n"
106  "madd.s %[sum], %[sum], %[in7], %[in8] \t\n"
107  "lwc1 %[in6], 400*4(%[synth_buf]) \t\n"
108  "lwc1 %[in7], 448*4(%[window]) \t\n"
109  "lwc1 %[in8], 464*4(%[synth_buf]) \t\n"
110  "madd.s %[sum], %[sum], %[in1], %[in2] \t\n"
111  "lwc1 %[in1], 32*4(%[window]) \t\n"
112  "lwc1 %[in2], 48*4(%[synth_buf]) \t\n"
113  "madd.s %[sum], %[sum], %[in3], %[in4] \t\n"
114  "lwc1 %[in3], 96*4(%[window]) \t\n"
115  "lwc1 %[in4], 112*4(%[synth_buf]) \t\n"
116  "madd.s %[sum], %[sum], %[in5], %[in6] \t\n"
117  "lwc1 %[in5], 160*4(%[window]) \t\n"
118  "lwc1 %[in6], 176*4(%[synth_buf]) \t\n"
119  "madd.s %[sum], %[sum], %[in7], %[in8] \t\n"
120  "lwc1 %[in7], 224*4(%[window]) \t\n"
121  "lwc1 %[in8], 240*4(%[synth_buf]) \t\n"
122  "nmsub.s %[sum], %[sum], %[in1], %[in2] \t\n"
123  "lwc1 %[in1], 288*4(%[window]) \t\n"
124  "lwc1 %[in2], 304*4(%[synth_buf]) \t\n"
125  "nmsub.s %[sum], %[sum], %[in3], %[in4] \t\n"
126  "lwc1 %[in3], 352*4(%[window]) \t\n"
127  "lwc1 %[in4], 368*4(%[synth_buf]) \t\n"
128  "nmsub.s %[sum], %[sum], %[in5], %[in6] \t\n"
129  "lwc1 %[in5], 416*4(%[window]) \t\n"
130  "lwc1 %[in6], 432*4(%[synth_buf]) \t\n"
131  "nmsub.s %[sum], %[sum], %[in7], %[in8] \t\n"
132  "lwc1 %[in7], 480*4(%[window]) \t\n"
133  "lwc1 %[in8], 496*4(%[synth_buf]) \t\n"
134  "nmsub.s %[sum], %[sum], %[in1], %[in2] \t\n"
135  PTR_ADDU "%[w], %[window], 4 \t\n"
136  "nmsub.s %[sum], %[sum], %[in3], %[in4] \t\n"
137  PTR_ADDU "%[w2], %[window], 124 \t\n"
138  PTR_ADDIU "%[p], %[synth_buf], 68 \t\n"
139  PTR_ADDIU "%[p2], %[synth_buf], 188 \t\n"
140  "nmsub.s %[sum], %[sum], %[in5], %[in6] \t\n"
141  "nmsub.s %[sum], %[sum], %[in7], %[in8] \t\n"
142  "swc1 %[sum], 0(%[samples]) \t\n"
143  PTR_ADDU "%[samples], %[samples], %[incr1] \t\n"
144 
145  /* calculate two samples at the same time to avoid one memory
146  access per two sample */
147 
148  "ff_mpadsp_apply_window_loop%=: \t\n"
149  "lwc1 %[in1], 0(%[w]) \t\n"
150  "lwc1 %[in2], 0(%[p]) \t\n"
151  "lwc1 %[in3], 0(%[w2]) \t\n"
152  "lwc1 %[in4], 64*4(%[w]) \t\n"
153  "lwc1 %[in5], 64*4(%[p]) \t\n"
154  "lwc1 %[in6], 64*4(%[w2]) \t\n"
155  "mul.s %[sum], %[in1], %[in2] \t\n"
156  "mul.s %[sum2], %[in2], %[in3] \t\n"
157  "lwc1 %[in1], 128*4(%[w]) \t\n"
158  "lwc1 %[in2], 128*4(%[p]) \t\n"
159  "madd.s %[sum], %[sum], %[in4], %[in5] \t\n"
160  "nmadd.s %[sum2], %[sum2], %[in5], %[in6] \t\n"
161  "lwc1 %[in3], 128*4(%[w2]) \t\n"
162  "lwc1 %[in4], 192*4(%[w]) \t\n"
163  "madd.s %[sum], %[sum], %[in1], %[in2] \t\n"
164  "lwc1 %[in5], 192*4(%[p]) \t\n"
165  "lwc1 %[in6], 192*4(%[w2]) \t\n"
166  "nmsub.s %[sum2], %[sum2], %[in2], %[in3] \t\n"
167  "lwc1 %[in1], 256*4(%[w]) \t\n"
168  "lwc1 %[in2], 256*4(%[p]) \t\n"
169  "madd.s %[sum], %[sum], %[in4], %[in5] \t\n"
170  "nmsub.s %[sum2], %[sum2], %[in5], %[in6] \t\n"
171  "lwc1 %[in3], 256*4(%[w2]) \t\n"
172  "lwc1 %[in4], 320*4(%[w]) \t\n"
173  "madd.s %[sum], %[sum], %[in1], %[in2] \t\n"
174  "lwc1 %[in5], 320*4(%[p]) \t\n"
175  "lwc1 %[in6], 320*4(%[w2]) \t\n"
176  "nmsub.s %[sum2], %[sum2], %[in2], %[in3] \t\n"
177  "lwc1 %[in1], 384*4(%[w]) \t\n"
178  "lwc1 %[in2], 384*4(%[p]) \t\n"
179  "madd.s %[sum], %[sum], %[in4], %[in5] \t\n"
180  "nmsub.s %[sum2], %[sum2], %[in5], %[in6] \t\n"
181  "lwc1 %[in3], 384*4(%[w2]) \t\n"
182  "lwc1 %[in4], 448*4(%[w]) \t\n"
183  "madd.s %[sum], %[sum], %[in1], %[in2] \t\n"
184  "lwc1 %[in5], 448*4(%[p]) \t\n"
185  "lwc1 %[in6], 448*4(%[w2]) \t\n"
186  "nmsub.s %[sum2], %[sum2], %[in2], %[in3] \t\n"
187  "madd.s %[sum], %[sum], %[in4], %[in5] \t\n"
188  "lwc1 %[in1], 32*4(%[w]) \t\n"
189  "lwc1 %[in2], 0(%[p2]) \t\n"
190  "nmsub.s %[sum2], %[sum2], %[in5], %[in6] \t\n"
191  "lwc1 %[in3], 32*4(%[w2]) \t\n"
192  "lwc1 %[in4], 96*4(%[w]) \t\n"
193  "lwc1 %[in5], 64*4(%[p2]) \t\n"
194  "nmsub.s %[sum], %[sum], %[in1], %[in2] \t\n"
195  "lwc1 %[in6], 96*4(%[w2]) \t\n"
196  "nmsub.s %[sum2], %[sum2], %[in2], %[in3] \t\n"
197  "lwc1 %[in1], 160*4(%[w]) \t\n"
198  "nmsub.s %[sum], %[sum], %[in4], %[in5] \t\n"
199  "lwc1 %[in2], 128*4(%[p2]) \t\n"
200  "nmsub.s %[sum2], %[sum2], %[in5], %[in6] \t\n"
201  "lwc1 %[in3], 160*4(%[w2]) \t\n"
202  "lwc1 %[in4], 224*4(%[w]) \t\n"
203  "lwc1 %[in5], 192*4(%[p2]) \t\n"
204  "nmsub.s %[sum], %[sum], %[in1], %[in2] \t\n"
205  "lwc1 %[in6], 224*4(%[w2]) \t\n"
206  "nmsub.s %[sum2], %[sum2], %[in2], %[in3] \t\n"
207  "lwc1 %[in1], 288*4(%[w]) \t\n"
208  "nmsub.s %[sum], %[sum], %[in4], %[in5] \t\n"
209  "lwc1 %[in2], 256*4(%[p2]) \t\n"
210  "nmsub.s %[sum2], %[sum2], %[in5], %[in6] \t\n"
211  "lwc1 %[in3], 288*4(%[w2]) \t\n"
212  "lwc1 %[in4], 352*4(%[w]) \t\n"
213  "lwc1 %[in5], 320*4(%[p2]) \t\n"
214  "nmsub.s %[sum], %[sum], %[in1], %[in2] \t\n"
215  "lwc1 %[in6], 352*4(%[w2]) \t\n"
216  "nmsub.s %[sum2], %[sum2], %[in2], %[in3] \t\n"
217  "lwc1 %[in1], 416*4(%[w]) \t\n"
218  "nmsub.s %[sum], %[sum], %[in4], %[in5] \t\n"
219  "lwc1 %[in2], 384*4(%[p2]) \t\n"
220  "nmsub.s %[sum2], %[sum2], %[in5], %[in6] \t\n"
221  "lwc1 %[in3], 416*4(%[w2]) \t\n"
222  "lwc1 %[in4], 480*4(%[w]) \t\n"
223  "lwc1 %[in5], 448*4(%[p2]) \t\n"
224  "nmsub.s %[sum], %[sum], %[in1], %[in2] \t\n"
225  "lwc1 %[in6], 480*4(%[w2]) \t\n"
226  "nmsub.s %[sum2], %[sum2], %[in2], %[in3] \t\n"
227  PTR_ADDIU "%[w], %[w], 4 \t\n"
228  "nmsub.s %[sum], %[sum], %[in4], %[in5] \t\n"
229  PTR_ADDIU "%[w2], %[w2], -4 \t\n"
230  "nmsub.s %[sum2], %[sum2], %[in5], %[in6] \t\n"
231  "addu %[j], %[j], 4 \t\n"
232  PTR_ADDIU "%[p], 4 \t\n"
233  "swc1 %[sum], 0(%[samples]) \t\n"
234  PTR_ADDIU "%[p2], -4 \t\n"
235  "swc1 %[sum2], 0(%[samples2]) \t\n"
236  PTR_ADDU "%[samples], %[samples], %[incr1] \t\n"
237  PTR_SUBU "%[samples2],%[samples2], %[incr1] \t\n"
238  "bne %[j], 64, ff_mpadsp_apply_window_loop%= \t\n"
239 
240  "lwc1 %[in1], 48*4(%[window]) \t\n"
241  "lwc1 %[in2], 32*4(%[synth_buf]) \t\n"
242  "lwc1 %[in3], 112*4(%[window]) \t\n"
243  "lwc1 %[in4], 96*4(%[synth_buf]) \t\n"
244  "lwc1 %[in5], 176*4(%[window]) \t\n"
245  "lwc1 %[in6], 160*4(%[synth_buf]) \t\n"
246  "mul.s %[sum], %[in1], %[in2] \t\n"
247  "lwc1 %[in7], 240*4(%[window]) \t\n"
248  "lwc1 %[in8], 224*4(%[synth_buf]) \t\n"
249  "lwc1 %[in1], 304*4(%[window]) \t\n"
250  "nmadd.s %[sum], %[sum], %[in3], %[in4] \t\n"
251  "lwc1 %[in2], 288*4(%[synth_buf]) \t\n"
252  "lwc1 %[in3], 368*4(%[window]) \t\n"
253  "lwc1 %[in4], 352*4(%[synth_buf]) \t\n"
254  "nmsub.s %[sum], %[sum], %[in5], %[in6] \t\n"
255  "nmsub.s %[sum], %[sum], %[in7], %[in8] \t\n"
256  "lwc1 %[in5], 432*4(%[window]) \t\n"
257  "lwc1 %[in6], 416*4(%[synth_buf]) \t\n"
258  "nmsub.s %[sum], %[sum], %[in1], %[in2] \t\n"
259  "lwc1 %[in7], 496*4(%[window]) \t\n"
260  "lwc1 %[in8], 480*4(%[synth_buf]) \t\n"
261  "nmsub.s %[sum], %[sum], %[in3], %[in4] \t\n"
262  "nmsub.s %[sum], %[sum], %[in5], %[in6] \t\n"
263  "nmsub.s %[sum], %[sum], %[in7], %[in8] \t\n"
264  "swc1 %[sum], 0(%[samples]) \t\n"
265 
266  : [sum] "=&f" (sum), [sum2] "=&f" (sum2),
267  [w2] "=&r" (w2), [w] "=&r" (w),
268  [p] "=&r" (p), [p2] "=&r" (p2), [j] "=&r" (j),
269  [samples] "+r" (samples), [samples2] "=&r" (samples2),
270  [in1] "=&f" (in1), [in2] "=&f" (in2),
271  [in3] "=&f" (in3), [in4] "=&f" (in4),
272  [in5] "=&f" (in5), [in6] "=&f" (in6),
273  [in7] "=&f" (in7), [in8] "=&f" (in8),
274  [t_sample] "=&r" (t_sample)
275  : [synth_buf] "r" (synth_buf), [window] "r" (window),
276  [dither_state] "r" (dither_state), [incr1] "r" (incr1)
277  : "memory"
278  );
279 }
280 
281 #if !HAVE_LOONGSON3
282 static void ff_dct32_mips_float(float *out, const float *tab)
283 {
284  float val0 , val1 , val2 , val3 , val4 , val5 , val6 , val7,
285  val8 , val9 , val10, val11, val12, val13, val14, val15,
286  val16, val17, val18, val19, val20, val21, val22, val23,
287  val24, val25, val26, val27, val28, val29, val30, val31;
288  float fTmp1, fTmp2, fTmp3, fTmp4, fTmp5, fTmp6, fTmp7, fTmp8,
289  fTmp9, fTmp10, fTmp11;
290 
291  /**
292  * instructions are scheduled to minimize pipeline stall.
293  */
294  __asm__ volatile (
295  "lwc1 %[fTmp1], 0*4(%[tab]) \n\t"
296  "lwc1 %[fTmp2], 31*4(%[tab]) \n\t"
297  "lwc1 %[fTmp3], 15*4(%[tab]) \n\t"
298  "lwc1 %[fTmp4], 16*4(%[tab]) \n\t"
299  "li.s %[fTmp7], 0.50241928618815570551 \n\t"
300  "add.s %[fTmp5], %[fTmp1], %[fTmp2] \n\t"
301  "sub.s %[fTmp8], %[fTmp1], %[fTmp2] \n\t"
302  "add.s %[fTmp6], %[fTmp3], %[fTmp4] \n\t"
303  "sub.s %[fTmp9], %[fTmp3], %[fTmp4] \n\t"
304  "li.s %[fTmp10], 0.50060299823519630134 \n\t"
305  "li.s %[fTmp11], 10.19000812354805681150 \n\t"
306  "mul.s %[fTmp8], %[fTmp8], %[fTmp10] \n\t"
307  "add.s %[val0], %[fTmp5], %[fTmp6] \n\t"
308  "sub.s %[val15], %[fTmp5], %[fTmp6] \n\t"
309  "lwc1 %[fTmp1], 7*4(%[tab]) \n\t"
310  "lwc1 %[fTmp2], 24*4(%[tab]) \n\t"
311  "madd.s %[val16], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
312  "nmsub.s %[val31], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
313  "mul.s %[val15], %[val15], %[fTmp7] \n\t"
314  "lwc1 %[fTmp3], 8*4(%[tab]) \n\t"
315  "lwc1 %[fTmp4], 23*4(%[tab]) \n\t"
316  "add.s %[fTmp5], %[fTmp1], %[fTmp2] \n\t"
317  "mul.s %[val31], %[val31], %[fTmp7] \n\t"
318  "sub.s %[fTmp8], %[fTmp1], %[fTmp2] \n\t"
319  "add.s %[fTmp6], %[fTmp3], %[fTmp4] \n\t"
320  "sub.s %[fTmp9], %[fTmp3], %[fTmp4] \n\t"
321  "li.s %[fTmp7], 5.10114861868916385802 \n\t"
322  "li.s %[fTmp10], 0.67480834145500574602 \n\t"
323  "li.s %[fTmp11], 0.74453627100229844977 \n\t"
324  "add.s %[val7], %[fTmp5], %[fTmp6] \n\t"
325  "sub.s %[val8], %[fTmp5], %[fTmp6] \n\t"
326  "mul.s %[fTmp8], %[fTmp8], %[fTmp10] \n\t"
327  "li.s %[fTmp1], 0.50979557910415916894 \n\t"
328  "sub.s %[fTmp2], %[val0], %[val7] \n\t"
329  "mul.s %[val8], %[val8], %[fTmp7] \n\t"
330  "madd.s %[val23], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
331  "nmsub.s %[val24], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
332  "add.s %[val0], %[val0], %[val7] \n\t"
333  "mul.s %[val7], %[fTmp1], %[fTmp2] \n\t"
334  "sub.s %[fTmp2], %[val15], %[val8] \n\t"
335  "add.s %[val8], %[val15], %[val8] \n\t"
336  "mul.s %[val24], %[val24], %[fTmp7] \n\t"
337  "sub.s %[fTmp3], %[val16], %[val23] \n\t"
338  "add.s %[val16], %[val16], %[val23] \n\t"
339  "mul.s %[val15], %[fTmp1], %[fTmp2] \n\t"
340  "sub.s %[fTmp4], %[val31], %[val24] \n\t"
341  "mul.s %[val23], %[fTmp1], %[fTmp3] \n\t"
342  "add.s %[val24], %[val31], %[val24] \n\t"
343  "mul.s %[val31], %[fTmp1], %[fTmp4] \n\t"
344 
345  : [fTmp1] "=&f" (fTmp1), [fTmp2] "=&f" (fTmp2), [fTmp3] "=&f" (fTmp3),
346  [fTmp4] "=&f" (fTmp4), [fTmp5] "=&f" (fTmp5), [fTmp6] "=&f" (fTmp6),
347  [fTmp7] "=&f" (fTmp7), [fTmp8] "=&f" (fTmp8), [fTmp9] "=&f" (fTmp9),
348  [fTmp10] "=&f" (fTmp10), [fTmp11] "=&f" (fTmp11),
349  [val0] "=f" (val0), [val7] "=f" (val7),
350  [val8] "=f" (val8), [val15] "=f" (val15),
351  [val16] "=f" (val16), [val23] "=f" (val23),
352  [val24] "=f" (val24), [val31] "=f" (val31)
353  : [tab] "r" (tab)
354  : "memory"
355  );
356 
357  __asm__ volatile (
358  "lwc1 %[fTmp1], 3*4(%[tab]) \n\t"
359  "lwc1 %[fTmp2], 28*4(%[tab]) \n\t"
360  "lwc1 %[fTmp3], 12*4(%[tab]) \n\t"
361  "lwc1 %[fTmp4], 19*4(%[tab]) \n\t"
362  "li.s %[fTmp7], 0.64682178335999012954 \n\t"
363  "add.s %[fTmp5], %[fTmp1], %[fTmp2] \n\t"
364  "sub.s %[fTmp8], %[fTmp1], %[fTmp2] \n\t"
365  "add.s %[fTmp6], %[fTmp3], %[fTmp4] \n\t"
366  "sub.s %[fTmp9], %[fTmp3], %[fTmp4] \n\t"
367  "li.s %[fTmp10], 0.53104259108978417447 \n\t"
368  "li.s %[fTmp11], 1.48416461631416627724 \n\t"
369  "mul.s %[fTmp8], %[fTmp8], %[fTmp10] \n\t"
370  "add.s %[val3], %[fTmp5], %[fTmp6] \n\t"
371  "sub.s %[val12], %[fTmp5], %[fTmp6] \n\t"
372  "lwc1 %[fTmp1], 4*4(%[tab]) \n\t"
373  "lwc1 %[fTmp2], 27*4(%[tab]) \n\t"
374  "madd.s %[val19], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
375  "nmsub.s %[val28], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
376  "mul.s %[val12], %[val12], %[fTmp7] \n\t"
377  "lwc1 %[fTmp3], 11*4(%[tab]) \n\t"
378  "lwc1 %[fTmp4], 20*4(%[tab]) \n\t"
379  "add.s %[fTmp5], %[fTmp1], %[fTmp2] \n\t"
380  "mul.s %[val28], %[val28], %[fTmp7] \n\t"
381  "sub.s %[fTmp8], %[fTmp1], %[fTmp2] \n\t"
382  "li.s %[fTmp7], 0.78815462345125022473 \n\t"
383  "add.s %[fTmp6], %[fTmp3], %[fTmp4] \n\t"
384  "sub.s %[fTmp9], %[fTmp3], %[fTmp4] \n\t"
385  "li.s %[fTmp10], 0.55310389603444452782 \n\t"
386  "li.s %[fTmp11], 1.16943993343288495515 \n\t"
387  "mul.s %[fTmp8], %[fTmp8], %[fTmp10] \n\t"
388  "add.s %[val4], %[fTmp5], %[fTmp6] \n\t"
389  "sub.s %[val11], %[fTmp5], %[fTmp6] \n\t"
390  "li.s %[fTmp1], 2.56291544774150617881 \n\t"
391  "madd.s %[val20], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
392  "nmsub.s %[val27], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
393  "mul.s %[val11], %[val11], %[fTmp7] \n\t"
394  "sub.s %[fTmp2], %[val3], %[val4] \n\t"
395  "add.s %[val3], %[val3], %[val4] \n\t"
396  "sub.s %[fTmp4], %[val19], %[val20] \n\t"
397  "mul.s %[val27], %[val27], %[fTmp7] \n\t"
398  "sub.s %[fTmp3], %[val12], %[val11] \n\t"
399  "mul.s %[val4], %[fTmp1], %[fTmp2] \n\t"
400  "add.s %[val11], %[val12], %[val11] \n\t"
401  "add.s %[val19], %[val19], %[val20] \n\t"
402  "mul.s %[val20], %[fTmp1], %[fTmp4] \n\t"
403  "mul.s %[val12], %[fTmp1], %[fTmp3] \n\t"
404  "sub.s %[fTmp2], %[val28], %[val27] \n\t"
405  "add.s %[val27], %[val28], %[val27] \n\t"
406  "mul.s %[val28], %[fTmp1], %[fTmp2] \n\t"
407 
408  : [fTmp1] "=&f" (fTmp1), [fTmp2] "=&f" (fTmp2), [fTmp3] "=&f" (fTmp3),
409  [fTmp4] "=&f" (fTmp4), [fTmp5] "=&f" (fTmp5), [fTmp6] "=&f" (fTmp6),
410  [fTmp7] "=&f" (fTmp7), [fTmp8] "=&f" (fTmp8), [fTmp9] "=&f" (fTmp9),
411  [fTmp10] "=&f" (fTmp10), [fTmp11] "=&f" (fTmp11),
412  [val3] "=f" (val3), [val4] "=f" (val4),
413  [val11] "=f" (val11), [val12] "=f" (val12),
414  [val19] "=f" (val19), [val20] "=f" (val20),
415  [val27] "=f" (val27), [val28] "=f" (val28)
416  : [tab] "r" (tab)
417  : "memory"
418  );
419 
420  __asm__ volatile (
421  "li.s %[fTmp1], 0.54119610014619698439 \n\t"
422  "sub.s %[fTmp2], %[val0], %[val3] \n\t"
423  "add.s %[val0], %[val0], %[val3] \n\t"
424  "sub.s %[fTmp3], %[val7], %[val4] \n\t"
425  "add.s %[val4], %[val7], %[val4] \n\t"
426  "sub.s %[fTmp4], %[val8], %[val11] \n\t"
427  "mul.s %[val3], %[fTmp1], %[fTmp2] \n\t"
428  "add.s %[val8], %[val8], %[val11] \n\t"
429  "mul.s %[val7], %[fTmp1], %[fTmp3] \n\t"
430  "sub.s %[fTmp2], %[val15], %[val12] \n\t"
431  "mul.s %[val11], %[fTmp1], %[fTmp4] \n\t"
432  "add.s %[val12], %[val15], %[val12] \n\t"
433  "mul.s %[val15], %[fTmp1], %[fTmp2] \n\t"
434 
435  : [val0] "+f" (val0), [val3] "+f" (val3),
436  [val4] "+f" (val4), [val7] "+f" (val7),
437  [val8] "+f" (val8), [val11] "+f" (val11),
438  [val12] "+f" (val12), [val15] "+f" (val15),
439  [fTmp1] "=f" (fTmp1), [fTmp2] "=&f" (fTmp2),
440  [fTmp3] "=&f" (fTmp3), [fTmp4] "=&f" (fTmp4)
441  :
442  );
443 
444  __asm__ volatile (
445  "sub.s %[fTmp2], %[val16], %[val19] \n\t"
446  "add.s %[val16], %[val16], %[val19] \n\t"
447  "sub.s %[fTmp3], %[val23], %[val20] \n\t"
448  "add.s %[val20], %[val23], %[val20] \n\t"
449  "sub.s %[fTmp4], %[val24], %[val27] \n\t"
450  "mul.s %[val19], %[fTmp1], %[fTmp2] \n\t"
451  "add.s %[val24], %[val24], %[val27] \n\t"
452  "mul.s %[val23], %[fTmp1], %[fTmp3] \n\t"
453  "sub.s %[fTmp2], %[val31], %[val28] \n\t"
454  "mul.s %[val27], %[fTmp1], %[fTmp4] \n\t"
455  "add.s %[val28], %[val31], %[val28] \n\t"
456  "mul.s %[val31], %[fTmp1], %[fTmp2] \n\t"
457 
458  : [fTmp2] "=&f" (fTmp2), [fTmp3] "=&f" (fTmp3), [fTmp4] "=&f" (fTmp4),
459  [val16] "+f" (val16), [val19] "+f" (val19), [val20] "+f" (val20),
460  [val23] "+f" (val23), [val24] "+f" (val24), [val27] "+f" (val27),
461  [val28] "+f" (val28), [val31] "+f" (val31)
462  : [fTmp1] "f" (fTmp1)
463  );
464 
465  __asm__ volatile (
466  "lwc1 %[fTmp1], 1*4(%[tab]) \n\t"
467  "lwc1 %[fTmp2], 30*4(%[tab]) \n\t"
468  "lwc1 %[fTmp3], 14*4(%[tab]) \n\t"
469  "lwc1 %[fTmp4], 17*4(%[tab]) \n\t"
470  "li.s %[fTmp7], 0.52249861493968888062 \n\t"
471  "add.s %[fTmp5], %[fTmp1], %[fTmp2] \n\t"
472  "sub.s %[fTmp8], %[fTmp1], %[fTmp2] \n\t"
473  "add.s %[fTmp6], %[fTmp3], %[fTmp4] \n\t"
474  "sub.s %[fTmp9], %[fTmp3], %[fTmp4] \n\t"
475  "li.s %[fTmp10], 0.50547095989754365998 \n\t"
476  "li.s %[fTmp11], 3.40760841846871878570 \n\t"
477  "mul.s %[fTmp8], %[fTmp8], %[fTmp10] \n\t"
478  "add.s %[val1], %[fTmp5], %[fTmp6] \n\t"
479  "sub.s %[val14], %[fTmp5], %[fTmp6] \n\t"
480  "lwc1 %[fTmp1], 6*4(%[tab]) \n\t"
481  "lwc1 %[fTmp2], 25*4(%[tab]) \n\t"
482  "madd.s %[val17], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
483  "nmsub.s %[val30], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
484  "mul.s %[val14], %[val14], %[fTmp7] \n\t"
485  "lwc1 %[fTmp3], 9*4(%[tab]) \n\t"
486  "lwc1 %[fTmp4], 22*4(%[tab]) \n\t"
487  "add.s %[fTmp5], %[fTmp1], %[fTmp2] \n\t"
488  "mul.s %[val30], %[val30], %[fTmp7] \n\t"
489  "sub.s %[fTmp8], %[fTmp1], %[fTmp2] \n\t"
490  "add.s %[fTmp6], %[fTmp3], %[fTmp4] \n\t"
491  "sub.s %[fTmp9], %[fTmp3], %[fTmp4] \n\t"
492  "li.s %[fTmp7], 1.72244709823833392782 \n\t"
493  "li.s %[fTmp10], 0.62250412303566481615 \n\t"
494  "li.s %[fTmp11], 0.83934964541552703873 \n\t"
495  "add.s %[val6], %[fTmp5], %[fTmp6] \n\t"
496  "sub.s %[val9], %[fTmp5], %[fTmp6] \n\t"
497  "mul.s %[fTmp8], %[fTmp8], %[fTmp10] \n\t"
498  "li.s %[fTmp1], 0.60134488693504528054 \n\t"
499  "sub.s %[fTmp2], %[val1], %[val6] \n\t"
500  "add.s %[val1], %[val1], %[val6] \n\t"
501  "mul.s %[val9], %[val9], %[fTmp7] \n\t"
502  "madd.s %[val22], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
503  "nmsub.s %[val25], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
504  "mul.s %[val6], %[fTmp1], %[fTmp2] \n\t"
505  "sub.s %[fTmp2], %[val14], %[val9] \n\t"
506  "add.s %[val9], %[val14], %[val9] \n\t"
507  "mul.s %[val25], %[val25], %[fTmp7] \n\t"
508  "sub.s %[fTmp3], %[val17], %[val22] \n\t"
509  "add.s %[val17], %[val17], %[val22] \n\t"
510  "mul.s %[val14], %[fTmp1], %[fTmp2] \n\t"
511  "sub.s %[fTmp2], %[val30], %[val25] \n\t"
512  "mul.s %[val22], %[fTmp1], %[fTmp3] \n\t"
513  "add.s %[val25], %[val30], %[val25] \n\t"
514  "mul.s %[val30], %[fTmp1], %[fTmp2] \n\t"
515 
516  : [fTmp1] "=&f" (fTmp1), [fTmp2] "=&f" (fTmp2), [fTmp3] "=&f" (fTmp3),
517  [fTmp4] "=&f" (fTmp4), [fTmp5] "=&f" (fTmp5), [fTmp6] "=&f" (fTmp6),
518  [fTmp7] "=&f" (fTmp7), [fTmp8] "=&f" (fTmp8), [fTmp9] "=&f" (fTmp9),
519  [fTmp10] "=&f" (fTmp10), [fTmp11] "=&f" (fTmp11),
520  [val1] "=f" (val1), [val6] "=f" (val6),
521  [val9] "=f" (val9), [val14] "=f" (val14),
522  [val17] "=f" (val17), [val22] "=f" (val22),
523  [val25] "=f" (val25), [val30] "=f" (val30)
524  : [tab] "r" (tab)
525  : "memory"
526  );
527 
528  __asm__ volatile (
529  "lwc1 %[fTmp1], 2*4(%[tab]) \n\t"
530  "lwc1 %[fTmp2], 29*4(%[tab]) \n\t"
531  "lwc1 %[fTmp3], 13*4(%[tab]) \n\t"
532  "lwc1 %[fTmp4], 18*4(%[tab]) \n\t"
533  "li.s %[fTmp7], 0.56694403481635770368 \n\t"
534  "add.s %[fTmp5], %[fTmp1], %[fTmp2] \n\t"
535  "sub.s %[fTmp8], %[fTmp1], %[fTmp2] \n\t"
536  "add.s %[fTmp6], %[fTmp3], %[fTmp4] \n\t"
537  "sub.s %[fTmp9], %[fTmp3], %[fTmp4] \n\t"
538  "li.s %[fTmp10], 0.51544730992262454697 \n\t"
539  "li.s %[fTmp11], 2.05778100995341155085 \n\t"
540  "mul.s %[fTmp8], %[fTmp8], %[fTmp10] \n\t"
541  "add.s %[val2], %[fTmp5], %[fTmp6] \n\t"
542  "sub.s %[val13], %[fTmp5], %[fTmp6] \n\t"
543  "lwc1 %[fTmp1], 5*4(%[tab]) \n\t"
544  "lwc1 %[fTmp2], 26*4(%[tab]) \n\t"
545  "madd.s %[val18], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
546  "nmsub.s %[val29], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
547  "mul.s %[val13], %[val13], %[fTmp7] \n\t"
548  "lwc1 %[fTmp3], 10*4(%[tab]) \n\t"
549  "lwc1 %[fTmp4], 21*4(%[tab]) \n\t"
550  "mul.s %[val29], %[val29], %[fTmp7] \n\t"
551  "add.s %[fTmp5], %[fTmp1], %[fTmp2] \n\t"
552  "sub.s %[fTmp8], %[fTmp1], %[fTmp2] \n\t"
553  "add.s %[fTmp6], %[fTmp3], %[fTmp4] \n\t"
554  "sub.s %[fTmp9], %[fTmp3], %[fTmp4] \n\t"
555  "li.s %[fTmp7], 1.06067768599034747134 \n\t"
556  "li.s %[fTmp10], 0.58293496820613387367 \n\t"
557  "li.s %[fTmp11], 0.97256823786196069369 \n\t"
558  "add.s %[val5], %[fTmp5], %[fTmp6] \n\t"
559  "sub.s %[val10], %[fTmp5], %[fTmp6] \n\t"
560  "mul.s %[fTmp8], %[fTmp8], %[fTmp10] \n\t"
561  "li.s %[fTmp1], 0.89997622313641570463 \n\t"
562  "sub.s %[fTmp2], %[val2], %[val5] \n\t"
563  "mul.s %[val10], %[val10], %[fTmp7] \n\t"
564  "madd.s %[val21], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
565  "nmsub.s %[val26], %[fTmp8], %[fTmp9], %[fTmp11] \n\t"
566  "add.s %[val2], %[val2], %[val5] \n\t"
567  "mul.s %[val5], %[fTmp1], %[fTmp2] \n\t"
568  "sub.s %[fTmp3], %[val13], %[val10] \n\t"
569  "add.s %[val10], %[val13], %[val10] \n\t"
570  "mul.s %[val26], %[val26], %[fTmp7] \n\t"
571  "sub.s %[fTmp4], %[val18], %[val21] \n\t"
572  "add.s %[val18], %[val18], %[val21] \n\t"
573  "mul.s %[val13], %[fTmp1], %[fTmp3] \n\t"
574  "sub.s %[fTmp2], %[val29], %[val26] \n\t"
575  "add.s %[val26], %[val29], %[val26] \n\t"
576  "mul.s %[val21], %[fTmp1], %[fTmp4] \n\t"
577  "mul.s %[val29], %[fTmp1], %[fTmp2] \n\t"
578 
579  : [fTmp1] "=&f" (fTmp1), [fTmp2] "=&f" (fTmp2), [fTmp3] "=&f" (fTmp3),
580  [fTmp4] "=&f" (fTmp4), [fTmp5] "=&f" (fTmp5), [fTmp6] "=&f" (fTmp6),
581  [fTmp7] "=&f" (fTmp7), [fTmp8] "=&f" (fTmp8), [fTmp9] "=&f" (fTmp9),
582  [fTmp10] "=&f" (fTmp10), [fTmp11] "=&f" (fTmp11),
583  [val2] "=f" (val2), [val5] "=f" (val5),
584  [val10] "=f" (val10), [val13] "=f" (val13),
585  [val18] "=f" (val18), [val21] "=f" (val21),
586  [val26] "=f" (val26), [val29] "=f" (val29)
587  : [tab] "r" (tab)
588  : "memory"
589  );
590 
591  __asm__ volatile (
592  "li.s %[fTmp1], 1.30656296487637652785 \n\t"
593  "sub.s %[fTmp2], %[val1], %[val2] \n\t"
594  "add.s %[val1], %[val1], %[val2] \n\t"
595  "sub.s %[fTmp3], %[val6], %[val5] \n\t"
596  "add.s %[val5], %[val6], %[val5] \n\t"
597  "sub.s %[fTmp4], %[val9], %[val10] \n\t"
598  "mul.s %[val2], %[fTmp1], %[fTmp2] \n\t"
599  "add.s %[val9], %[val9], %[val10] \n\t"
600  "mul.s %[val6], %[fTmp1], %[fTmp3] \n\t"
601  "sub.s %[fTmp2], %[val14], %[val13] \n\t"
602  "mul.s %[val10], %[fTmp1], %[fTmp4] \n\t"
603  "add.s %[val13], %[val14], %[val13] \n\t"
604  "mul.s %[val14], %[fTmp1], %[fTmp2] \n\t"
605 
606  : [fTmp1] "=f" (fTmp1), [fTmp2] "=&f" (fTmp2),
607  [fTmp3] "=&f" (fTmp3), [fTmp4] "=&f" (fTmp4),
608  [val1] "+f" (val1), [val2] "+f" (val2),
609  [val5] "+f" (val5), [val6] "+f" (val6),
610  [val9] "+f" (val9), [val10] "+f" (val10),
611  [val13] "+f" (val13), [val14] "+f" (val14)
612  :
613  );
614 
615  __asm__ volatile (
616  "sub.s %[fTmp2], %[val17], %[val18] \n\t"
617  "add.s %[val17], %[val17], %[val18] \n\t"
618  "sub.s %[fTmp3], %[val22], %[val21] \n\t"
619  "add.s %[val21], %[val22], %[val21] \n\t"
620  "sub.s %[fTmp4], %[val25], %[val26] \n\t"
621  "mul.s %[val18], %[fTmp1], %[fTmp2] \n\t"
622  "add.s %[val25], %[val25], %[val26] \n\t"
623  "mul.s %[val22], %[fTmp1], %[fTmp3] \n\t"
624  "sub.s %[fTmp2], %[val30], %[val29] \n\t"
625  "mul.s %[val26], %[fTmp1], %[fTmp4] \n\t"
626  "add.s %[val29], %[val30], %[val29] \n\t"
627  "mul.s %[val30], %[fTmp1], %[fTmp2] \n\t"
628 
629  : [fTmp2] "=&f" (fTmp2), [fTmp3] "=&f" (fTmp3), [fTmp4] "=&f" (fTmp4),
630  [val17] "+f" (val17), [val18] "+f" (val18), [val21] "+f" (val21),
631  [val22] "+f" (val22), [val25] "+f" (val25), [val26] "+f" (val26),
632  [val29] "+f" (val29), [val30] "+f" (val30)
633  : [fTmp1] "f" (fTmp1)
634  );
635 
636  __asm__ volatile (
637  "li.s %[fTmp1], 0.70710678118654752439 \n\t"
638  "sub.s %[fTmp2], %[val0], %[val1] \n\t"
639  "add.s %[val0], %[val0], %[val1] \n\t"
640  "sub.s %[fTmp3], %[val3], %[val2] \n\t"
641  "add.s %[val2], %[val3], %[val2] \n\t"
642  "sub.s %[fTmp4], %[val4], %[val5] \n\t"
643  "mul.s %[val1], %[fTmp1], %[fTmp2] \n\t"
644  "swc1 %[val0], 0(%[out]) \n\t"
645  "mul.s %[val3], %[fTmp3], %[fTmp1] \n\t"
646  "add.s %[val4], %[val4], %[val5] \n\t"
647  "mul.s %[val5], %[fTmp1], %[fTmp4] \n\t"
648  "swc1 %[val1], 16*4(%[out]) \n\t"
649  "sub.s %[fTmp2], %[val7], %[val6] \n\t"
650  "add.s %[val2], %[val2], %[val3] \n\t"
651  "swc1 %[val3], 24*4(%[out]) \n\t"
652  "add.s %[val6], %[val7], %[val6] \n\t"
653  "mul.s %[val7], %[fTmp1], %[fTmp2] \n\t"
654  "swc1 %[val2], 8*4(%[out]) \n\t"
655  "add.s %[val6], %[val6], %[val7] \n\t"
656  "swc1 %[val7], 28*4(%[out]) \n\t"
657  "add.s %[val4], %[val4], %[val6] \n\t"
658  "add.s %[val6], %[val6], %[val5] \n\t"
659  "add.s %[val5], %[val5], %[val7] \n\t"
660  "swc1 %[val4], 4*4(%[out]) \n\t"
661  "swc1 %[val5], 20*4(%[out]) \n\t"
662  "swc1 %[val6], 12*4(%[out]) \n\t"
663 
664  : [fTmp1] "=f" (fTmp1), [fTmp2] "=&f" (fTmp2),
665  [fTmp3] "=&f" (fTmp3), [fTmp4] "=&f" (fTmp4),
666  [val0] "+f" (val0), [val1] "+f" (val1),
667  [val2] "+f" (val2), [val3] "+f" (val3),
668  [val4] "+f" (val4), [val5] "+f" (val5),
669  [val6] "+f" (val6), [val7] "+f" (val7)
670  : [out] "r" (out)
671  );
672 
673  __asm__ volatile (
674  "sub.s %[fTmp2], %[val8], %[val9] \n\t"
675  "add.s %[val8], %[val8], %[val9] \n\t"
676  "sub.s %[fTmp3], %[val11], %[val10] \n\t"
677  "add.s %[val10], %[val11], %[val10] \n\t"
678  "sub.s %[fTmp4], %[val12], %[val13] \n\t"
679  "mul.s %[val9], %[fTmp1], %[fTmp2] \n\t"
680  "add.s %[val12], %[val12], %[val13] \n\t"
681  "mul.s %[val11], %[fTmp1], %[fTmp3] \n\t"
682  "sub.s %[fTmp2], %[val15], %[val14] \n\t"
683  "mul.s %[val13], %[fTmp1], %[fTmp4] \n\t"
684  "add.s %[val14], %[val15], %[val14] \n\t"
685  "add.s %[val10], %[val10], %[val11] \n\t"
686  "mul.s %[val15], %[fTmp1], %[fTmp2] \n\t"
687  "add.s %[val14], %[val14], %[val15] \n\t"
688  "add.s %[val12], %[val12], %[val14] \n\t"
689  "add.s %[val14], %[val14], %[val13] \n\t"
690  "add.s %[val13], %[val13], %[val15] \n\t"
691  "add.s %[val8], %[val8], %[val12] \n\t"
692  "add.s %[val12], %[val12], %[val10] \n\t"
693  "add.s %[val10], %[val10], %[val14] \n\t"
694  "add.s %[val14], %[val14], %[val9] \n\t"
695  "add.s %[val9], %[val9], %[val13] \n\t"
696  "add.s %[val13], %[val13], %[val11] \n\t"
697  "add.s %[val11], %[val11], %[val15] \n\t"
698  "swc1 %[val8], 2*4(%[out]) \n\t"
699  "swc1 %[val9], 18*4(%[out]) \n\t"
700  "swc1 %[val10], 10*4(%[out]) \n\t"
701  "swc1 %[val11], 26*4(%[out]) \n\t"
702  "swc1 %[val12], 6*4(%[out]) \n\t"
703  "swc1 %[val13], 22*4(%[out]) \n\t"
704  "swc1 %[val14], 14*4(%[out]) \n\t"
705  "swc1 %[val15], 30*4(%[out]) \n\t"
706 
707  : [fTmp2] "=&f" (fTmp2), [fTmp3] "=&f" (fTmp3), [fTmp4] "=&f" (fTmp4),
708  [val8] "+f" (val8), [val9] "+f" (val9), [val10] "+f" (val10),
709  [val11] "+f" (val11), [val12] "+f" (val12), [val13] "+f" (val13),
710  [val14] "+f" (val14), [val15] "+f" (val15)
711  : [fTmp1] "f" (fTmp1), [out] "r" (out)
712  );
713 
714  __asm__ volatile (
715  "sub.s %[fTmp2], %[val16], %[val17] \n\t"
716  "add.s %[val16], %[val16], %[val17] \n\t"
717  "sub.s %[fTmp3], %[val19], %[val18] \n\t"
718  "add.s %[val18], %[val19], %[val18] \n\t"
719  "sub.s %[fTmp4], %[val20], %[val21] \n\t"
720  "mul.s %[val17], %[fTmp1], %[fTmp2] \n\t"
721  "add.s %[val20], %[val20], %[val21] \n\t"
722  "mul.s %[val19], %[fTmp1], %[fTmp3] \n\t"
723  "sub.s %[fTmp2], %[val23], %[val22] \n\t"
724  "mul.s %[val21], %[fTmp1], %[fTmp4] \n\t"
725  "add.s %[val22], %[val23], %[val22] \n\t"
726  "add.s %[val18], %[val18], %[val19] \n\t"
727  "mul.s %[val23], %[fTmp1], %[fTmp2] \n\t"
728  "add.s %[val22], %[val22], %[val23] \n\t"
729  "add.s %[val20], %[val20], %[val22] \n\t"
730  "add.s %[val22], %[val22], %[val21] \n\t"
731  "add.s %[val21], %[val21], %[val23] \n\t"
732 
733  : [fTmp2] "=&f" (fTmp2), [fTmp3] "=&f" (fTmp3), [fTmp4] "=&f" (fTmp4),
734  [val16] "+f" (val16), [val17] "+f" (val17), [val18] "+f" (val18),
735  [val19] "+f" (val19), [val20] "+f" (val20), [val21] "+f" (val21),
736  [val22] "+f" (val22), [val23] "+f" (val23)
737  : [fTmp1] "f" (fTmp1)
738  );
739 
740  __asm__ volatile (
741  "sub.s %[fTmp2], %[val24], %[val25] \n\t"
742  "add.s %[val24], %[val24], %[val25] \n\t"
743  "sub.s %[fTmp3], %[val27], %[val26] \n\t"
744  "add.s %[val26], %[val27], %[val26] \n\t"
745  "sub.s %[fTmp4], %[val28], %[val29] \n\t"
746  "mul.s %[val25], %[fTmp1], %[fTmp2] \n\t"
747  "add.s %[val28], %[val28], %[val29] \n\t"
748  "mul.s %[val27], %[fTmp1], %[fTmp3] \n\t"
749  "sub.s %[fTmp2], %[val31], %[val30] \n\t"
750  "mul.s %[val29], %[fTmp1], %[fTmp4] \n\t"
751  "add.s %[val30], %[val31], %[val30] \n\t"
752  "add.s %[val26], %[val26], %[val27] \n\t"
753  "mul.s %[val31], %[fTmp1], %[fTmp2] \n\t"
754  "add.s %[val30], %[val30], %[val31] \n\t"
755  "add.s %[val28], %[val28], %[val30] \n\t"
756  "add.s %[val30], %[val30], %[val29] \n\t"
757  "add.s %[val29], %[val29], %[val31] \n\t"
758  "add.s %[val24], %[val24], %[val28] \n\t"
759  "add.s %[val28], %[val28], %[val26] \n\t"
760  "add.s %[val26], %[val26], %[val30] \n\t"
761  "add.s %[val30], %[val30], %[val25] \n\t"
762  "add.s %[val25], %[val25], %[val29] \n\t"
763  "add.s %[val29], %[val29], %[val27] \n\t"
764  "add.s %[val27], %[val27], %[val31] \n\t"
765 
766  : [fTmp2] "=&f" (fTmp2), [fTmp3] "=&f" (fTmp3), [fTmp4] "=&f" (fTmp4),
767  [val24] "+f" (val24), [val25] "+f" (val25), [val26] "+f" (val26),
768  [val27] "+f" (val27), [val28] "+f" (val28), [val29] "+f" (val29),
769  [val30] "+f" (val30), [val31] "+f" (val31)
770  : [fTmp1] "f" (fTmp1)
771  );
772 
773  out[ 1] = val16 + val24;
774  out[17] = val17 + val25;
775  out[ 9] = val18 + val26;
776  out[25] = val19 + val27;
777  out[ 5] = val20 + val28;
778  out[21] = val21 + val29;
779  out[13] = val22 + val30;
780  out[29] = val23 + val31;
781  out[ 3] = val24 + val20;
782  out[19] = val25 + val21;
783  out[11] = val26 + val22;
784  out[27] = val27 + val23;
785  out[ 7] = val28 + val18;
786  out[23] = val29 + val19;
787  out[15] = val30 + val17;
788  out[31] = val31;
789 }
790 #endif /* !HAVE_LOONGSON3 */
791 
792 static void imdct36_mips_float(float *out, float *buf, float *in, float *win)
793 {
794  float t0, t1, t2, t3, s0, s1, s2, s3;
795  float tmp[18];
796  /* temporary variables */
797  float in1, in2, in3, in4, in5, in6;
798  float out1, out2, out3, out4, out5;
799  float c1, c2, c3, c4, c5, c6, c7, c8, c9;
800 
801  /**
802  * all loops are unrolled totally, and instructions are scheduled to
803  * minimize pipeline stall. instructions of the first two loops are
804  * reorganized, in order to eliminate unnecessary readings and
805  * writings into array. values defined in macros and tables are
806  * eliminated - they are directly loaded in appropriate variables
807  */
808 
809  /* loop 1 and 2 */
810  __asm__ volatile (
811  "lwc1 %[in1], 17*4(%[in]) \t\n"
812  "lwc1 %[in2], 16*4(%[in]) \t\n"
813  "lwc1 %[in3], 15*4(%[in]) \t\n"
814  "lwc1 %[in4], 14*4(%[in]) \t\n"
815  "lwc1 %[in5], 13*4(%[in]) \t\n"
816  "lwc1 %[in6], 12*4(%[in]) \t\n"
817  "add.s %[out1], %[in1], %[in2] \t\n"
818  "add.s %[out2], %[in2], %[in3] \t\n"
819  "add.s %[out3], %[in3], %[in4] \t\n"
820  "add.s %[out4], %[in4], %[in5] \t\n"
821  "add.s %[out5], %[in5], %[in6] \t\n"
822  "lwc1 %[in1], 11*4(%[in]) \t\n"
823  "swc1 %[out2], 16*4(%[in]) \t\n"
824  "add.s %[out1], %[out1], %[out3] \t\n"
825  "swc1 %[out4], 14*4(%[in]) \t\n"
826  "add.s %[out3], %[out3], %[out5] \t\n"
827  "lwc1 %[in2], 10*4(%[in]) \t\n"
828  "lwc1 %[in3], 9*4(%[in]) \t\n"
829  "swc1 %[out1], 17*4(%[in]) \t\n"
830  "lwc1 %[in4], 8*4(%[in]) \t\n"
831  "swc1 %[out3], 15*4(%[in]) \t\n"
832  "add.s %[out1], %[in6], %[in1] \t\n"
833  "add.s %[out2], %[in1], %[in2] \t\n"
834  "add.s %[out3], %[in2], %[in3] \t\n"
835  "add.s %[out4], %[in3], %[in4] \t\n"
836  "lwc1 %[in5], 7*4(%[in]) \t\n"
837  "swc1 %[out1], 12*4(%[in]) \t\n"
838  "add.s %[out5], %[out5], %[out2] \t\n"
839  "swc1 %[out3], 10*4(%[in]) \t\n"
840  "add.s %[out2], %[out2], %[out4] \t\n"
841  "lwc1 %[in6], 6*4(%[in]) \t\n"
842  "lwc1 %[in1], 5*4(%[in]) \t\n"
843  "swc1 %[out5], 13*4(%[in]) \t\n"
844  "lwc1 %[in2], 4*4(%[in]) \t\n"
845  "swc1 %[out2], 11*4(%[in]) \t\n"
846  "add.s %[out5], %[in4], %[in5] \t\n"
847  "add.s %[out1], %[in5], %[in6] \t\n"
848  "add.s %[out2], %[in6], %[in1] \t\n"
849  "add.s %[out3], %[in1], %[in2] \t\n"
850  "lwc1 %[in3], 3*4(%[in]) \t\n"
851  "swc1 %[out5], 8*4(%[in]) \t\n"
852  "add.s %[out4], %[out4], %[out1] \t\n"
853  "swc1 %[out2], 6*4(%[in]) \t\n"
854  "add.s %[out1], %[out1], %[out3] \t\n"
855  "lwc1 %[in4], 2*4(%[in]) \t\n"
856  "lwc1 %[in5], 1*4(%[in]) \t\n"
857  "swc1 %[out4], 9*4(%[in]) \t\n"
858  "lwc1 %[in6], 0(%[in]) \t\n"
859  "swc1 %[out1], 7*4(%[in]) \t\n"
860  "add.s %[out4], %[in2], %[in3] \t\n"
861  "add.s %[out5], %[in3], %[in4] \t\n"
862  "add.s %[out1], %[in4], %[in5] \t\n"
863  "add.s %[out2], %[in5], %[in6] \t\n"
864  "swc1 %[out4], 4*4(%[in]) \t\n"
865  "add.s %[out3], %[out3], %[out5] \t\n"
866  "swc1 %[out1], 2*4(%[in]) \t\n"
867  "add.s %[out5], %[out5], %[out2] \t\n"
868  "swc1 %[out2], 1*4(%[in]) \t\n"
869  "swc1 %[out3], 5*4(%[in]) \t\n"
870  "swc1 %[out5], 3*4(%[in]) \t\n"
871 
872  : [in1] "=&f" (in1), [in2] "=&f" (in2),
873  [in3] "=&f" (in3), [in4] "=&f" (in4),
874  [in5] "=&f" (in5), [in6] "=&f" (in6),
875  [out1] "=&f" (out1), [out2] "=&f" (out2),
876  [out3] "=&f" (out3), [out4] "=&f" (out4),
877  [out5] "=&f" (out5)
878  : [in] "r" (in)
879  : "memory"
880  );
881 
882  /* loop 3 */
883  __asm__ volatile (
884  "li.s %[c1], 0.5 \t\n"
885  "lwc1 %[in1], 8*4(%[in]) \t\n"
886  "lwc1 %[in2], 16*4(%[in]) \t\n"
887  "lwc1 %[in3], 4*4(%[in]) \t\n"
888  "lwc1 %[in4], 0(%[in]) \t\n"
889  "lwc1 %[in5], 12*4(%[in]) \t\n"
890  "li.s %[c2], 0.93969262078590838405 \t\n"
891  "add.s %[t2], %[in1], %[in2] \t\n"
892  "add.s %[t0], %[in1], %[in3] \t\n"
893  "li.s %[c3], -0.76604444311897803520 \t\n"
894  "madd.s %[t3], %[in4], %[in5], %[c1] \t\n"
895  "sub.s %[t1], %[in4], %[in5] \t\n"
896  "sub.s %[t2], %[t2], %[in3] \t\n"
897  "mul.s %[t0], %[t0], %[c2] \t\n"
898  "li.s %[c4], -0.17364817766693034885 \t\n"
899  "li.s %[c5], -0.86602540378443864676 \t\n"
900  "li.s %[c6], 0.98480775301220805936 \t\n"
901  "nmsub.s %[out1], %[t1], %[t2], %[c1] \t\n"
902  "add.s %[out2], %[t1], %[t2] \t\n"
903  "add.s %[t2], %[in2], %[in3] \t\n"
904  "sub.s %[t1], %[in1], %[in2] \t\n"
905  "sub.s %[out3], %[t3], %[t0] \t\n"
906  "swc1 %[out1], 6*4(%[tmp]) \t\n"
907  "swc1 %[out2], 16*4(%[tmp]) \t\n"
908  "mul.s %[t2], %[t2], %[c3] \t\n"
909  "mul.s %[t1], %[t1], %[c4] \t\n"
910  "add.s %[out1], %[t3], %[t0] \t\n"
911  "lwc1 %[in1], 10*4(%[in]) \t\n"
912  "lwc1 %[in2], 14*4(%[in]) \t\n"
913  "sub.s %[out3], %[out3], %[t2] \t\n"
914  "add.s %[out2], %[t3], %[t2] \t\n"
915  "add.s %[out1], %[out1], %[t1] \t\n"
916  "lwc1 %[in3], 2*4(%[in]) \t\n"
917  "lwc1 %[in4], 6*4(%[in]) \t\n"
918  "swc1 %[out3], 10*4(%[tmp]) \t\n"
919  "sub.s %[out2], %[out2], %[t1] \t\n"
920  "swc1 %[out1], 2*4(%[tmp]) \t\n"
921  "add.s %[out1], %[in1], %[in2] \t\n"
922  "add.s %[t2], %[in1], %[in3] \t\n"
923  "sub.s %[t3], %[in1], %[in2] \t\n"
924  "swc1 %[out2], 14*4(%[tmp]) \t\n"
925  "li.s %[c7], -0.34202014332566873304 \t\n"
926  "sub.s %[out1], %[out1], %[in3] \t\n"
927  "mul.s %[t2], %[t2], %[c6] \t\n"
928  "mul.s %[t3], %[t3], %[c7] \t\n"
929  "li.s %[c8], 0.86602540378443864676 \t\n"
930  "mul.s %[t0], %[in4], %[c8] \t\n"
931  "mul.s %[out1], %[out1], %[c5] \t\n"
932  "add.s %[t1], %[in2], %[in3] \t\n"
933  "li.s %[c9], -0.64278760968653932632 \t\n"
934  "add.s %[out2], %[t2], %[t3] \t\n"
935  "lwc1 %[in1], 9*4(%[in]) \t\n"
936  "swc1 %[out1], 4*4(%[tmp]) \t\n"
937  "mul.s %[t1], %[t1], %[c9] \t\n"
938  "lwc1 %[in2], 17*4(%[in]) \t\n"
939  "add.s %[out2], %[out2], %[t0] \t\n"
940  "lwc1 %[in3], 5*4(%[in]) \t\n"
941  "lwc1 %[in4], 1*4(%[in]) \t\n"
942  "add.s %[out3], %[t2], %[t1] \t\n"
943  "sub.s %[out1], %[t3], %[t1] \t\n"
944  "swc1 %[out2], 0(%[tmp]) \t\n"
945  "lwc1 %[in5], 13*4(%[in]) \t\n"
946  "add.s %[t2], %[in1], %[in2] \t\n"
947  "sub.s %[out3], %[out3], %[t0] \t\n"
948  "sub.s %[out1], %[out1], %[t0] \t\n"
949  "add.s %[t0], %[in1], %[in3] \t\n"
950  "madd.s %[t3], %[in4], %[in5], %[c1] \t\n"
951  "sub.s %[t2], %[t2], %[in3] \t\n"
952  "swc1 %[out3], 12*4(%[tmp]) \t\n"
953  "swc1 %[out1], 8*4(%[tmp]) \t\n"
954  "sub.s %[t1], %[in4], %[in5] \t\n"
955  "mul.s %[t0], %[t0], %[c2] \t\n"
956  "nmsub.s %[out1], %[t1], %[t2], %[c1] \t\n"
957  "add.s %[out2], %[t1], %[t2] \t\n"
958  "add.s %[t2], %[in2], %[in3] \t\n"
959  "sub.s %[t1], %[in1], %[in2] \t\n"
960  "sub.s %[out3], %[t3], %[t0] \t\n"
961  "swc1 %[out1], 7*4(%[tmp]) \t\n"
962  "swc1 %[out2], 17*4(%[tmp]) \t\n"
963  "mul.s %[t2], %[t2], %[c3] \t\n"
964  "mul.s %[t1], %[t1], %[c4] \t\n"
965  "add.s %[out1], %[t3], %[t0] \t\n"
966  "lwc1 %[in1], 11*4(%[in]) \t\n"
967  "lwc1 %[in2], 15*4(%[in]) \t\n"
968  "sub.s %[out3], %[out3], %[t2] \t\n"
969  "add.s %[out2], %[t3], %[t2] \t\n"
970  "add.s %[out1], %[out1], %[t1] \t\n"
971  "lwc1 %[in3], 3*4(%[in]) \t\n"
972  "lwc1 %[in4], 7*4(%[in]) \t\n"
973  "swc1 %[out3], 11*4(%[tmp]) \t\n"
974  "sub.s %[out2], %[out2], %[t1] \t\n"
975  "swc1 %[out1], 3*4(%[tmp]) \t\n"
976  "add.s %[out3], %[in1], %[in2] \t\n"
977  "add.s %[t2], %[in1], %[in3] \t\n"
978  "sub.s %[t3], %[in1], %[in2] \t\n"
979  "swc1 %[out2], 15*4(%[tmp]) \t\n"
980  "mul.s %[t0], %[in4], %[c8] \t\n"
981  "sub.s %[out3], %[out3], %[in3] \t\n"
982  "mul.s %[t2], %[t2], %[c6] \t\n"
983  "mul.s %[t3], %[t3], %[c7] \t\n"
984  "add.s %[t1], %[in2], %[in3] \t\n"
985  "mul.s %[out3], %[out3], %[c5] \t\n"
986  "add.s %[out1], %[t2], %[t3] \t\n"
987  "mul.s %[t1], %[t1], %[c9] \t\n"
988  "swc1 %[out3], 5*4(%[tmp]) \t\n"
989  "add.s %[out1], %[out1], %[t0] \t\n"
990  "add.s %[out2], %[t2], %[t1] \t\n"
991  "sub.s %[out3], %[t3], %[t1] \t\n"
992  "swc1 %[out1], 1*4(%[tmp]) \t\n"
993  "sub.s %[out2], %[out2], %[t0] \t\n"
994  "sub.s %[out3], %[out3], %[t0] \t\n"
995  "swc1 %[out2], 13*4(%[tmp]) \t\n"
996  "swc1 %[out3], 9*4(%[tmp]) \t\n"
997 
998  : [t0] "=&f" (t0), [t1] "=&f" (t1),
999  [t2] "=&f" (t2), [t3] "=&f" (t3),
1000  [in1] "=&f" (in1), [in2] "=&f" (in2),
1001  [in3] "=&f" (in3), [in4] "=&f" (in4),
1002  [in5] "=&f" (in5),
1003  [out1] "=&f" (out1), [out2] "=&f" (out2),
1004  [out3] "=&f" (out3),
1005  [c1] "=&f" (c1), [c2] "=&f" (c2),
1006  [c3] "=&f" (c3), [c4] "=&f" (c4),
1007  [c5] "=&f" (c5), [c6] "=&f" (c6),
1008  [c7] "=&f" (c7), [c8] "=&f" (c8),
1009  [c9] "=&f" (c9)
1010  : [in] "r" (in), [tmp] "r" (tmp)
1011  : "memory"
1012  );
1013 
1014  /* loop 4 */
1015  __asm__ volatile (
1016  "lwc1 %[in1], 2*4(%[tmp]) \t\n"
1017  "lwc1 %[in2], 0(%[tmp]) \t\n"
1018  "lwc1 %[in3], 3*4(%[tmp]) \t\n"
1019  "lwc1 %[in4], 1*4(%[tmp]) \t\n"
1020  "li.s %[c1], 0.50190991877167369479 \t\n"
1021  "li.s %[c2], 5.73685662283492756461 \t\n"
1022  "add.s %[s0], %[in1], %[in2] \t\n"
1023  "sub.s %[s2], %[in1], %[in2] \t\n"
1024  "add.s %[s1], %[in3], %[in4] \t\n"
1025  "sub.s %[s3], %[in3], %[in4] \t\n"
1026  "lwc1 %[in1], 9*4(%[win]) \t\n"
1027  "lwc1 %[in2], 4*9*4(%[buf]) \t\n"
1028  "lwc1 %[in3], 8*4(%[win]) \t\n"
1029  "mul.s %[s1], %[s1], %[c1] \t\n"
1030  "mul.s %[s3], %[s3], %[c2] \t\n"
1031  "lwc1 %[in4], 4*8*4(%[buf]) \t\n"
1032  "lwc1 %[in5], 29*4(%[win]) \t\n"
1033  "lwc1 %[in6], 28*4(%[win]) \t\n"
1034  "add.s %[t0], %[s0], %[s1] \t\n"
1035  "sub.s %[t1], %[s0], %[s1] \t\n"
1036  "li.s %[c1], 0.51763809020504152469 \t\n"
1037  "li.s %[c2], 1.93185165257813657349 \t\n"
1038  "mul.s %[out3], %[in5], %[t0] \t\n"
1039  "madd.s %[out1], %[in2], %[in1], %[t1] \t\n"
1040  "madd.s %[out2], %[in4], %[in3], %[t1] \t\n"
1041  "mul.s %[out4], %[in6], %[t0] \t\n"
1042  "add.s %[t0], %[s2], %[s3] \t\n"
1043  "swc1 %[out3], 4*9*4(%[buf]) \t\n"
1044  "swc1 %[out1], 288*4(%[out]) \t\n"
1045  "swc1 %[out2], 256*4(%[out]) \t\n"
1046  "swc1 %[out4], 4*8*4(%[buf]) \t\n"
1047  "sub.s %[t1], %[s2], %[s3] \t\n"
1048  "lwc1 %[in1], 17*4(%[win]) \t\n"
1049  "lwc1 %[in2], 4*17*4(%[buf]) \t\n"
1050  "lwc1 %[in3], 0(%[win]) \t\n"
1051  "lwc1 %[in4], 0(%[buf]) \t\n"
1052  "lwc1 %[in5], 37*4(%[win]) \t\n"
1053  "lwc1 %[in6], 20*4(%[win]) \t\n"
1054  "madd.s %[out1], %[in2], %[in1], %[t1] \t\n"
1055  "lwc1 %[in1], 6*4(%[tmp]) \t\n"
1056  "madd.s %[out2], %[in4], %[in3], %[t1] \t\n"
1057  "mul.s %[out3], %[t0], %[in5] \t\n"
1058  "mul.s %[out4], %[t0], %[in6] \t\n"
1059  "swc1 %[out1], 544*4(%[out]) \t\n"
1060  "lwc1 %[in2], 4*4(%[tmp]) \t\n"
1061  "swc1 %[out2], 0(%[out]) \t\n"
1062  "swc1 %[out3], 4*17*4(%[buf]) \t\n"
1063  "swc1 %[out4], 0(%[buf]) \t\n"
1064  "lwc1 %[in3], 7*4(%[tmp]) \t\n"
1065  "add.s %[s0], %[in1], %[in2] \t\n"
1066  "sub.s %[s2], %[in1], %[in2] \t\n"
1067  "lwc1 %[in4], 5*4(%[tmp]) \t\n"
1068  "add.s %[s1], %[in3], %[in4] \t\n"
1069  "sub.s %[s3], %[in3], %[in4] \t\n"
1070  "lwc1 %[in1], 10*4(%[win]) \t\n"
1071  "lwc1 %[in2], 4*10*4(%[buf]) \t\n"
1072  "lwc1 %[in3], 7*4(%[win]) \t\n"
1073  "mul.s %[s1], %[s1], %[c1] \t\n"
1074  "mul.s %[s3], %[s3], %[c2] \t\n"
1075  "add.s %[t0], %[s0], %[s1] \t\n"
1076  "sub.s %[t1], %[s0], %[s1] \t\n"
1077  "lwc1 %[in4], 4*7*4(%[buf]) \t\n"
1078  "lwc1 %[in5], 30*4(%[win]) \t\n"
1079  "lwc1 %[in6], 27*4(%[win]) \t\n"
1080  "li.s %[c1], 0.55168895948124587824 \t\n"
1081  "madd.s %[out1], %[in2], %[in1], %[t1] \t\n"
1082  "madd.s %[out2], %[in4], %[in3], %[t1] \t\n"
1083  "mul.s %[out3], %[t0], %[in5] \t\n"
1084  "mul.s %[out4], %[t0], %[in6] \t\n"
1085  "add.s %[t0], %[s2], %[s3] \t\n"
1086  "swc1 %[out1], 320*4(%[out]) \t\n"
1087  "swc1 %[out2], 224*4(%[out]) \t\n"
1088  "swc1 %[out3], 4*10*4(%[buf]) \t\n"
1089  "swc1 %[out4], 4*7*4(%[buf]) \t\n"
1090  "sub.s %[t1], %[s2], %[s3] \t\n"
1091  "lwc1 %[in1], 16*4(%[win]) \t\n"
1092  "lwc1 %[in2], 4*16*4(%[buf]) \t\n"
1093  "lwc1 %[in3], 1*4(%[win]) \t\n"
1094  "lwc1 %[in4], 4*1*4(%[buf]) \t\n"
1095  "lwc1 %[in5], 36*4(%[win]) \t\n"
1096  "lwc1 %[in6], 21*4(%[win]) \t\n"
1097  "madd.s %[out1], %[in2], %[in1], %[t1] \t\n"
1098  "lwc1 %[in1], 10*4(%[tmp]) \t\n"
1099  "madd.s %[out2], %[in4], %[in3], %[t1] \t\n"
1100  "mul.s %[out3], %[in5], %[t0] \t\n"
1101  "mul.s %[out4], %[in6], %[t0] \t\n"
1102  "swc1 %[out1], 512*4(%[out]) \t\n"
1103  "lwc1 %[in2], 8*4(%[tmp]) \t\n"
1104  "swc1 %[out2], 32*4(%[out]) \t\n"
1105  "swc1 %[out3], 4*16*4(%[buf]) \t\n"
1106  "swc1 %[out4], 4*1*4(%[buf]) \t\n"
1107  "li.s %[c2], 1.18310079157624925896 \t\n"
1108  "add.s %[s0], %[in1], %[in2] \t\n"
1109  "sub.s %[s2], %[in1], %[in2] \t\n"
1110  "lwc1 %[in3], 11*4(%[tmp]) \t\n"
1111  "lwc1 %[in4], 9*4(%[tmp]) \t\n"
1112  "add.s %[s1], %[in3], %[in4] \t\n"
1113  "sub.s %[s3], %[in3], %[in4] \t\n"
1114  "lwc1 %[in1], 11*4(%[win]) \t\n"
1115  "lwc1 %[in2], 4*11*4(%[buf]) \t\n"
1116  "lwc1 %[in3], 6*4(%[win]) \t\n"
1117  "mul.s %[s1], %[s1], %[c1] \t\n"
1118  "mul.s %[s3], %[s3], %[c2] \t\n"
1119  "lwc1 %[in4], 4*6*4(%[buf]) \t\n"
1120  "lwc1 %[in5], 31*4(%[win]) \t\n"
1121  "lwc1 %[in6], 26*4(%[win]) \t\n"
1122  "add.s %[t0], %[s0], %[s1] \t\n"
1123  "sub.s %[t1], %[s0], %[s1] \t\n"
1124  "mul.s %[out3], %[t0], %[in5] \t\n"
1125  "mul.s %[out4], %[t0], %[in6] \t\n"
1126  "add.s %[t0], %[s2], %[s3] \t\n"
1127  "madd.s %[out1], %[in2], %[in1], %[t1] \t\n"
1128  "madd.s %[out2], %[in4], %[in3], %[t1] \t\n"
1129  "swc1 %[out3], 4*11*4(%[buf]) \t\n"
1130  "swc1 %[out4], 4*6*4(%[buf]) \t\n"
1131  "sub.s %[t1], %[s2], %[s3] \t\n"
1132  "swc1 %[out1], 352*4(%[out]) \t\n"
1133  "swc1 %[out2], 192*4(%[out]) \t\n"
1134  "lwc1 %[in1], 15*4(%[win]) \t\n"
1135  "lwc1 %[in2], 4*15*4(%[buf]) \t\n"
1136  "lwc1 %[in3], 2*4(%[win]) \t\n"
1137  "lwc1 %[in4], 4*2*4(%[buf]) \t\n"
1138  "lwc1 %[in5], 35*4(%[win]) \t\n"
1139  "lwc1 %[in6], 22*4(%[win]) \t\n"
1140  "madd.s %[out1], %[in2], %[in1], %[t1] \t\n"
1141  "lwc1 %[in1], 14*4(%[tmp]) \t\n"
1142  "madd.s %[out2], %[in4], %[in3], %[t1] \t\n"
1143  "mul.s %[out3], %[t0], %[in5] \t\n"
1144  "mul.s %[out4], %[t0], %[in6] \t\n"
1145  "swc1 %[out1], 480*4(%[out]) \t\n"
1146  "lwc1 %[in2], 12*4(%[tmp]) \t\n"
1147  "swc1 %[out2], 64*4(%[out]) \t\n"
1148  "swc1 %[out3], 4*15*4(%[buf]) \t\n"
1149  "swc1 %[out4], 4*2*4(%[buf]) \t\n"
1150  "lwc1 %[in3], 15*4(%[tmp]) \t\n"
1151  "add.s %[s0], %[in1], %[in2] \t\n"
1152  "sub.s %[s2], %[in1], %[in2] \t\n"
1153  "lwc1 %[in4], 13*4(%[tmp]) \t\n"
1154  "li.s %[c1], 0.61038729438072803416 \t\n"
1155  "li.s %[c2], 0.87172339781054900991 \t\n"
1156  "add.s %[s1], %[in3], %[in4] \t\n"
1157  "sub.s %[s3], %[in3], %[in4] \t\n"
1158  "lwc1 %[in1], 12*4(%[win]) \t\n"
1159  "lwc1 %[in2], 4*12*4(%[buf]) \t\n"
1160  "lwc1 %[in3], 5*4(%[win]) \t\n"
1161  "mul.s %[s1], %[s1], %[c1] \t\n"
1162  "mul.s %[s3], %[s3], %[c2] \t\n"
1163  "lwc1 %[in4], 4*5*4(%[buf]) \t\n"
1164  "lwc1 %[in5], 32*4(%[win]) \t\n"
1165  "lwc1 %[in6], 25*4(%[win]) \t\n"
1166  "add.s %[t0], %[s0], %[s1] \t\n"
1167  "sub.s %[t1], %[s0], %[s1] \t\n"
1168  "lwc1 %[s0], 16*4(%[tmp]) \t\n"
1169  "lwc1 %[s1], 17*4(%[tmp]) \t\n"
1170  "li.s %[c1], 0.70710678118654752439 \t\n"
1171  "mul.s %[out3], %[t0], %[in5] \t\n"
1172  "madd.s %[out1], %[in2], %[in1], %[t1] \t\n"
1173  "madd.s %[out2], %[in4], %[in3], %[t1] \t\n"
1174  "mul.s %[out4], %[t0], %[in6] \t\n"
1175  "add.s %[t0], %[s2], %[s3] \t\n"
1176  "swc1 %[out3], 4*12*4(%[buf]) \t\n"
1177  "swc1 %[out1], 384*4(%[out]) \t\n"
1178  "swc1 %[out2], 160*4(%[out]) \t\n"
1179  "swc1 %[out4], 4*5*4(%[buf]) \t\n"
1180  "sub.s %[t1], %[s2], %[s3] \t\n"
1181  "lwc1 %[in1], 14*4(%[win]) \t\n"
1182  "lwc1 %[in2], 4*14*4(%[buf]) \t\n"
1183  "lwc1 %[in3], 3*4(%[win]) \t\n"
1184  "lwc1 %[in4], 4*3*4(%[buf]) \t\n"
1185  "lwc1 %[in5], 34*4(%[win]) \t\n"
1186  "lwc1 %[in6], 23*4(%[win]) \t\n"
1187  "madd.s %[out1], %[in2], %[in1], %[t1] \t\n"
1188  "mul.s %[s1], %[s1], %[c1] \t\n"
1189  "madd.s %[out2], %[in4], %[in3], %[t1] \t\n"
1190  "mul.s %[out3], %[in5], %[t0] \t\n"
1191  "mul.s %[out4], %[in6], %[t0] \t\n"
1192  "swc1 %[out1], 448*4(%[out]) \t\n"
1193  "add.s %[t0], %[s0], %[s1] \t\n"
1194  "swc1 %[out2], 96*4(%[out]) \t\n"
1195  "swc1 %[out3], 4*14*4(%[buf]) \t\n"
1196  "swc1 %[out4], 4*3*4(%[buf]) \t\n"
1197  "sub.s %[t1], %[s0], %[s1] \t\n"
1198  "lwc1 %[in1], 13*4(%[win]) \t\n"
1199  "lwc1 %[in2], 4*13*4(%[buf]) \t\n"
1200  "lwc1 %[in3], 4*4(%[win]) \t\n"
1201  "lwc1 %[in4], 4*4*4(%[buf]) \t\n"
1202  "lwc1 %[in5], 33*4(%[win]) \t\n"
1203  "lwc1 %[in6], 24*4(%[win]) \t\n"
1204  "madd.s %[out1], %[in2], %[in1], %[t1] \t\n"
1205  "madd.s %[out2], %[in4], %[in3], %[t1] \t\n"
1206  "mul.s %[out3], %[t0], %[in5] \t\n"
1207  "mul.s %[out4], %[t0], %[in6] \t\n"
1208  "swc1 %[out1], 416*4(%[out]) \t\n"
1209  "swc1 %[out2], 128*4(%[out]) \t\n"
1210  "swc1 %[out3], 4*13*4(%[buf]) \t\n"
1211  "swc1 %[out4], 4*4*4(%[buf]) \t\n"
1212 
1213  : [c1] "=&f" (c1), [c2] "=&f" (c2),
1214  [in1] "=&f" (in1), [in2] "=&f" (in2),
1215  [in3] "=&f" (in3), [in4] "=&f" (in4),
1216  [in5] "=&f" (in5), [in6] "=&f" (in6),
1217  [out1] "=&f" (out1), [out2] "=&f" (out2),
1218  [out3] "=&f" (out3), [out4] "=&f" (out4),
1219  [t0] "=&f" (t0), [t1] "=&f" (t1),
1220  [t2] "=&f" (t2), [t3] "=&f" (t3),
1221  [s0] "=&f" (s0), [s1] "=&f" (s1),
1222  [s2] "=&f" (s2), [s3] "=&f" (s3)
1223  : [tmp] "r" (tmp), [win] "r" (win),
1224  [buf] "r" (buf), [out] "r" (out)
1225  : "memory"
1226  );
1227 }
1228 
1229 #if !HAVE_LOONGSON3
1230 static void ff_imdct36_blocks_mips_float(float *out, float *buf, float *in,
1231  int count, int switch_point, int block_type)
1232 {
1233  int j;
1234  for (j=0 ; j < count; j++) {
1235  /* apply window & overlap with previous buffer */
1236 
1237  /* select window */
1238  int win_idx = (switch_point && j < 2) ? 0 : block_type;
1239  float *win = ff_mdct_win_float[win_idx + (4 & -(j & 1))];
1240 
1241  imdct36_mips_float(out, buf, in, win);
1242 
1243  in += 18;
1244  buf += ((j&3) != 3 ? 1 : (72-3));
1245  out++;
1246  }
1247 }
1248 #endif /* !HAVE_LOONGSON3 */
1249 
1250 void ff_mpadsp_init_mipsfpu(MPADSPContext *s)
1251 {
1252  s->apply_window_float = ff_mpadsp_apply_window_mips_float;
1253 #if !HAVE_LOONGSON3
1254  s->imdct36_blocks_float = ff_imdct36_blocks_mips_float;
1255  s->dct32_float = ff_dct32_mips_float;
1256 #endif /* !HAVE_LOONGSON3 */
1257 }
s
const char * s
Definition: avisynth_c.h:631
MPADSPContext::imdct36_blocks_float
void(* imdct36_blocks_float)(float *out, float *buf, float *in, int count, int switch_point, int block_type)
Definition: mpegaudiodsp.h:33
asmdefs.h
MIPS assembly defines from sys/asm.h but rewritten for use with C inline assembly (rather than from w...
MPADSPContext::dct32_float
void(* dct32_float)(float *dst, const float *src)
Definition: mpegaudiodsp.h:30
ff_mdct_win_float
float ff_mdct_win_float[8][MDCT_BUF_SIZE]
t0
#define t0
Definition: regdef.h:28
c1
static const uint64_t c1
Definition: murmur3.c:49
ff_mpadsp_apply_window_mips_float
static void ff_mpadsp_apply_window_mips_float(float *synth_buf, float *window, int *dither_state, float *samples, int incr)
Definition: mpegaudiodsp_mips_float.c:61
s2
#define s2
Definition: regdef.h:39
imdct36_mips_float
static void imdct36_mips_float(float *out, float *buf, float *in, float *win)
Definition: mpegaudiodsp_mips_float.c:792
s0
#define s0
Definition: regdef.h:37
t1
#define t1
Definition: regdef.h:29
t3
#define t3
Definition: regdef.h:31
count
GLsizei count
Definition: opengl_enc.c:109
PTR_SUBU
#define PTR_SUBU
Definition: asmdefs.h:42
s3
#define s3
Definition: regdef.h:40
in
uint8_t pi<< 24) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0f/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S16, int16_t,(*(constint16_t *) pi >>8)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0f/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S32, int32_t,(*(constint32_t *) pi >>24)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0f/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_FLT, float, av_clip_uint8(lrintf(*(constfloat *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_FLT, float, av_clip_int16(lrintf(*(constfloat *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_FLT, float, av_clipl_int32(llrintf(*(constfloat *) pi *(1U<< 31)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_DBL, double, av_clip_uint8(lrint(*(constdouble *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_DBL, double, av_clip_int16(lrint(*(constdouble *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_DBL, double, av_clipl_int32(llrint(*(constdouble *) pi *(1U<< 31))))#defineSET_CONV_FUNC_GROUP(ofmt, ifmt) staticvoidset_generic_function(AudioConvert *ac){}voidff_audio_convert_free(AudioConvert **ac){if(!*ac) return;ff_dither_free(&(*ac) ->dc);av_freep(ac);}AudioConvert *ff_audio_convert_alloc(AVAudioResampleContext *avr, enumAVSampleFormatout_fmt, enumAVSampleFormatin_fmt, intchannels, intsample_rate, intapply_map){AudioConvert *ac;intin_planar, out_planar;ac=av_mallocz(sizeof(*ac));if(!ac) returnNULL;ac->avr=avr;ac->out_fmt=out_fmt;ac->in_fmt=in_fmt;ac->channels=channels;ac->apply_map=apply_map;if(avr->dither_method!=AV_RESAMPLE_DITHER_NONE &&av_get_packed_sample_fmt(out_fmt)==AV_SAMPLE_FMT_S16 &&av_get_bytes_per_sample(in_fmt)>2){ac->dc=ff_dither_alloc(avr, out_fmt, in_fmt, channels, sample_rate, apply_map);if(!ac->dc){av_free(ac);returnNULL;}returnac;}in_planar=ff_sample_fmt_is_planar(in_fmt, channels);out_planar=ff_sample_fmt_is_planar(out_fmt, channels);if(in_planar==out_planar){ac->func_type=CONV_FUNC_TYPE_FLAT;ac->planes=in_planar?ac->channels:1;}elseif(in_planar) ac->func_type=CONV_FUNC_TYPE_INTERLEAVE;elseac->func_type=CONV_FUNC_TYPE_DEINTERLEAVE;set_generic_function(ac);if(ARCH_AARCH64) ff_audio_convert_init_aarch64(ac);if(ARCH_ARM) ff_audio_convert_init_arm(ac);if(ARCH_X86) ff_audio_convert_init_x86(ac);returnac;}intff_audio_convert(AudioConvert *ac, AudioData *out, AudioData *in){intuse_generic=1;intlen=in->nb_samples;intp;if(ac->dc){av_log(ac->avr, AV_LOG_TRACE,"%dsamples-audio_convert:%sto%s(dithered)\n", len, av_get_sample_fmt_name(ac->in_fmt), av_get_sample_fmt_name(ac->out_fmt));returnff_convert_dither(ac-> in
Definition: audio_convert.c:194
buf
void * buf
Definition: avisynth_c.h:553
ff_dct32_mips_float
static void ff_dct32_mips_float(float *out, const float *tab)
Definition: mpegaudiodsp_mips_float.c:282
PTR_ADDIU
#define PTR_ADDIU
Definition: asmdefs.h:41
MPADSPContext::apply_window_float
void(* apply_window_float)(float *synth_buf, float *window, int *dither_state, float *samples, int incr)
Definition: mpegaudiodsp.h:26
ff_imdct36_blocks_mips_float
static void ff_imdct36_blocks_mips_float(float *out, float *buf, float *in, int count, int switch_point, int block_type)
Definition: mpegaudiodsp_mips_float.c:1230
s1
#define s1
Definition: regdef.h:38
ff_mpadsp_init_mipsfpu
void ff_mpadsp_init_mipsfpu(MPADSPContext *s)
Definition: mpegaudiodsp_mips_float.c:1250
c2
static const uint64_t c2
Definition: murmur3.c:50
tab
static const struct twinvq_data tab
Definition: twinvq_data.h:11135
out
uint8_t pi<< 24) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0f/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S16, int16_t,(*(constint16_t *) pi >>8)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0f/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S32, int32_t,(*(constint32_t *) pi >>24)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0f/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_FLT, float, av_clip_uint8(lrintf(*(constfloat *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_FLT, float, av_clip_int16(lrintf(*(constfloat *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_FLT, float, av_clipl_int32(llrintf(*(constfloat *) pi *(1U<< 31)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_DBL, double, av_clip_uint8(lrint(*(constdouble *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_DBL, double, av_clip_int16(lrint(*(constdouble *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_DBL, double, av_clipl_int32(llrint(*(constdouble *) pi *(1U<< 31))))#defineSET_CONV_FUNC_GROUP(ofmt, ifmt) staticvoidset_generic_function(AudioConvert *ac){}voidff_audio_convert_free(AudioConvert **ac){if(!*ac) return;ff_dither_free(&(*ac) ->dc);av_freep(ac);}AudioConvert *ff_audio_convert_alloc(AVAudioResampleContext *avr, enumAVSampleFormatout_fmt, enumAVSampleFormatin_fmt, intchannels, intsample_rate, intapply_map){AudioConvert *ac;intin_planar, out_planar;ac=av_mallocz(sizeof(*ac));if(!ac) returnNULL;ac->avr=avr;ac->out_fmt=out_fmt;ac->in_fmt=in_fmt;ac->channels=channels;ac->apply_map=apply_map;if(avr->dither_method!=AV_RESAMPLE_DITHER_NONE &&av_get_packed_sample_fmt(out_fmt)==AV_SAMPLE_FMT_S16 &&av_get_bytes_per_sample(in_fmt)>2){ac->dc=ff_dither_alloc(avr, out_fmt, in_fmt, channels, sample_rate, apply_map);if(!ac->dc){av_free(ac);returnNULL;}returnac;}in_planar=ff_sample_fmt_is_planar(in_fmt, channels);out_planar=ff_sample_fmt_is_planar(out_fmt, channels);if(in_planar==out_planar){ac->func_type=CONV_FUNC_TYPE_FLAT;ac->planes=in_planar?ac->channels:1;}elseif(in_planar) ac->func_type=CONV_FUNC_TYPE_INTERLEAVE;elseac->func_type=CONV_FUNC_TYPE_DEINTERLEAVE;set_generic_function(ac);if(ARCH_AARCH64) ff_audio_convert_init_aarch64(ac);if(ARCH_ARM) ff_audio_convert_init_arm(ac);if(ARCH_X86) ff_audio_convert_init_x86(ac);returnac;}intff_audio_convert(AudioConvert *ac, AudioData *out, AudioData *in){intuse_generic=1;intlen=in->nb_samples;intp;if(ac->dc){av_log(ac->avr, AV_LOG_TRACE,"%dsamples-audio_convert:%sto%s(dithered)\n", len, av_get_sample_fmt_name(ac->in_fmt), av_get_sample_fmt_name(ac->out_fmt));returnff_convert_dither(ac-> out
Definition: audio_convert.c:194
PTR_ADDU
#define PTR_ADDU
Definition: asmdefs.h:40
t2
#define t2
Definition: regdef.h:30
Generated on Wed Jun 10 2015 01:56:45 for FFmpeg by   doxygen 1.8.6