FFmpeg
aacsbr_mips.h
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  * Authors: Djordje Pesut (djordje@mips.com)
30  * Mirjana Vulin (mvulin@mips.com)
31  *
32  * This file is part of FFmpeg.
33  *
34  * FFmpeg is free software; you can redistribute it and/or
35  * modify it under the terms of the GNU Lesser General Public
36  * License as published by the Free Software Foundation; either
37  * version 2.1 of the License, or (at your option) any later version.
38  *
39  * FFmpeg is distributed in the hope that it will be useful,
40  * but WITHOUT ANY WARRANTY; without even the implied warranty of
41  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
42  * Lesser General Public License for more details.
43  *
44  * You should have received a copy of the GNU Lesser General Public
45  * License along with FFmpeg; if not, write to the Free Software
46  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
47  */
48 
49 /**
50  * @file
51  * Reference: libavcodec/aacsbr.c
52  */
53 
54 #ifndef AVCODEC_MIPS_AACSBR_MIPS_H
55 #define AVCODEC_MIPS_AACSBR_MIPS_H
56 
57 #include "libavcodec/aacdec.h"
58 #include "libavcodec/sbr.h"
59 #include "libavutil/mips/asmdefs.h"
60 
61 #if HAVE_INLINE_ASM
62 static void sbr_qmf_analysis_mips(AVFloatDSPContext *fdsp, AVTXContext *mdct, av_tx_fn mdct_fn,
63  SBRDSPContext *sbrdsp, const float *in, float *x,
64  float z[320], float W[2][32][32][2], int buf_idx)
65 {
66  int i;
67  float *w0;
68  float *w1;
69  int temp0, temp1, temp2, temp3, temp4, temp5, temp6, temp7;
70 
71  w0 = x;
72  w1 = x + 1024;
73  for(i = 0; i < 36; i++)
74  {
75  /* loop unrolled 8 times */
76  __asm__ volatile(
77  "lw %[temp0], 0(%[w1]) \n\t"
78  "lw %[temp1], 4(%[w1]) \n\t"
79  "lw %[temp2], 8(%[w1]) \n\t"
80  "lw %[temp3], 12(%[w1]) \n\t"
81  "lw %[temp4], 16(%[w1]) \n\t"
82  "lw %[temp5], 20(%[w1]) \n\t"
83  "lw %[temp6], 24(%[w1]) \n\t"
84  "lw %[temp7], 28(%[w1]) \n\t"
85  "sw %[temp0], 0(%[w0]) \n\t"
86  "sw %[temp1], 4(%[w0]) \n\t"
87  "sw %[temp2], 8(%[w0]) \n\t"
88  "sw %[temp3], 12(%[w0]) \n\t"
89  "sw %[temp4], 16(%[w0]) \n\t"
90  "sw %[temp5], 20(%[w0]) \n\t"
91  "sw %[temp6], 24(%[w0]) \n\t"
92  "sw %[temp7], 28(%[w0]) \n\t"
93  PTR_ADDIU " %[w0], %[w0], 32 \n\t"
94  PTR_ADDIU " %[w1], %[w1], 32 \n\t"
95 
96  : [w0]"+r"(w0), [w1]"+r"(w1),
97  [temp0]"=&r"(temp0), [temp1]"=&r"(temp1),
98  [temp2]"=&r"(temp2), [temp3]"=&r"(temp3),
99  [temp4]"=&r"(temp4), [temp5]"=&r"(temp5),
100  [temp6]"=&r"(temp6), [temp7]"=&r"(temp7)
101  :
102  : "memory"
103  );
104  }
105 
106  w0 = x + 288;
107  w1 = (float*)in;
108  for(i = 0; i < 128; i++)
109  {
110  /* loop unrolled 8 times */
111  __asm__ volatile(
112  "lw %[temp0], 0(%[w1]) \n\t"
113  "lw %[temp1], 4(%[w1]) \n\t"
114  "lw %[temp2], 8(%[w1]) \n\t"
115  "lw %[temp3], 12(%[w1]) \n\t"
116  "lw %[temp4], 16(%[w1]) \n\t"
117  "lw %[temp5], 20(%[w1]) \n\t"
118  "lw %[temp6], 24(%[w1]) \n\t"
119  "lw %[temp7], 28(%[w1]) \n\t"
120  "sw %[temp0], 0(%[w0]) \n\t"
121  "sw %[temp1], 4(%[w0]) \n\t"
122  "sw %[temp2], 8(%[w0]) \n\t"
123  "sw %[temp3], 12(%[w0]) \n\t"
124  "sw %[temp4], 16(%[w0]) \n\t"
125  "sw %[temp5], 20(%[w0]) \n\t"
126  "sw %[temp6], 24(%[w0]) \n\t"
127  "sw %[temp7], 28(%[w0]) \n\t"
128  PTR_ADDIU " %[w0], %[w0], 32 \n\t"
129  PTR_ADDIU " %[w1], %[w1], 32 \n\t"
130 
131  : [w0]"+r"(w0), [w1]"+r"(w1),
132  [temp0]"=&r"(temp0), [temp1]"=&r"(temp1),
133  [temp2]"=&r"(temp2), [temp3]"=&r"(temp3),
134  [temp4]"=&r"(temp4), [temp5]"=&r"(temp5),
135  [temp6]"=&r"(temp6), [temp7]"=&r"(temp7)
136  :
137  : "memory"
138  );
139  }
140 
141  for (i = 0; i < 32; i++) { // numTimeSlots*RATE = 16*2 as 960 sample frames
142  // are not supported
143  fdsp->vector_fmul_reverse(z, sbr_qmf_window_ds, x, 320);
144  sbrdsp->sum64x5(z);
145  sbrdsp->qmf_pre_shuffle(z);
146  mdct_fn(mdct, z, z+64, sizeof(float));
147  sbrdsp->qmf_post_shuffle(W[buf_idx][i], z);
148  x += 32;
149  }
150 }
151 
152 #if HAVE_MIPSFPU
153 #if !HAVE_MIPS32R6 && !HAVE_MIPS64R6
154 static void sbr_qmf_synthesis_mips(AVTXContext *mdct, av_tx_fn mdct_fn,
155  SBRDSPContext *sbrdsp, AVFloatDSPContext *fdsp,
156  float *out, float X[2][38][64],
157  float mdct_buf[2][64],
158  float *v0, int *v_off, const unsigned int div)
159 {
160  int i, n;
161  const float *sbr_qmf_window = div ? sbr_qmf_window_ds : sbr_qmf_window_us;
162  const int step = 128 >> div;
163  float *v;
164  float temp0, temp1, temp2, temp3, temp4, temp5, temp6, temp7, temp8, temp9, temp10, temp11, temp12, temp13;
165  float temp14, temp15, temp16, temp17, temp18, temp19;
166  float *vv0, *s0, *dst;
167  dst = out;
168 
169  for (i = 0; i < 32; i++) {
170  if (*v_off < step) {
171  int saved_samples = (1280 - 128) >> div;
172  memcpy(&v0[SBR_SYNTHESIS_BUF_SIZE - saved_samples], v0, saved_samples * sizeof(float));
173  *v_off = SBR_SYNTHESIS_BUF_SIZE - saved_samples - step;
174  } else {
175  *v_off -= step;
176  }
177  v = v0 + *v_off;
178  if (div) {
179  for (n = 0; n < 32; n++) {
180  X[0][i][ n] = -X[0][i][n];
181  X[0][i][32+n] = X[1][i][31-n];
182  }
183  mdct_fn(mdct, mdct_buf[0], X[0][i], sizeof(float));
184  sbrdsp->qmf_deint_neg(v, mdct_buf[0]);
185  } else {
186  sbrdsp->neg_odd_64(X[1][i]);
187  mdct_fn(mdct, mdct_buf[0], X[0][i], sizeof(float));
188  mdct_fn(mdct, mdct_buf[1], X[1][i], sizeof(float));
189  sbrdsp->qmf_deint_bfly(v, mdct_buf[1], mdct_buf[0]);
190  }
191 
192  if(div == 0)
193  {
194  float *v0_end;
195  vv0 = v;
196  v0_end = v + 60;
197  s0 = (float*)sbr_qmf_window;
198 
199  /* 10 calls of function vector_fmul_add merged into one loop
200  and loop unrolled 4 times */
201  __asm__ volatile(
202  ".set push \n\t"
203  ".set noreorder \n\t"
204  "lwc1 %[temp4], 0(%[v0]) \n\t"
205  "lwc1 %[temp5], 0(%[s0]) \n\t"
206  "lwc1 %[temp6], 4(%[v0]) \n\t"
207  "lwc1 %[temp7], 4(%[s0]) \n\t"
208  "lwc1 %[temp8], 8(%[v0]) \n\t"
209  "lwc1 %[temp9], 8(%[s0]) \n\t"
210  "lwc1 %[temp10], 12(%[v0]) \n\t"
211  "lwc1 %[temp11], 12(%[s0]) \n\t"
212  "lwc1 %[temp12], 768(%[v0]) \n\t"
213  "lwc1 %[temp13], 256(%[s0]) \n\t"
214  "lwc1 %[temp14], 772(%[v0]) \n\t"
215  "lwc1 %[temp15], 260(%[s0]) \n\t"
216  "lwc1 %[temp16], 776(%[v0]) \n\t"
217  "lwc1 %[temp17], 264(%[s0]) \n\t"
218  "lwc1 %[temp18], 780(%[v0]) \n\t"
219  "lwc1 %[temp19], 268(%[s0]) \n\t"
220  "1: \n\t"
221  "mul.s %[temp0], %[temp4], %[temp5] \n\t"
222  "lwc1 %[temp4], 1024(%[v0]) \n\t"
223  "mul.s %[temp1], %[temp6], %[temp7] \n\t"
224  "lwc1 %[temp5], 512(%[s0]) \n\t"
225  "mul.s %[temp2], %[temp8], %[temp9] \n\t"
226  "lwc1 %[temp6], 1028(%[v0]) \n\t"
227  "mul.s %[temp3], %[temp10], %[temp11] \n\t"
228  "lwc1 %[temp7], 516(%[s0]) \n\t"
229  "madd.s %[temp0], %[temp0], %[temp12], %[temp13] \n\t"
230  "lwc1 %[temp8], 1032(%[v0]) \n\t"
231  "madd.s %[temp1], %[temp1], %[temp14], %[temp15] \n\t"
232  "lwc1 %[temp9], 520(%[s0]) \n\t"
233  "madd.s %[temp2], %[temp2], %[temp16], %[temp17] \n\t"
234  "lwc1 %[temp10], 1036(%[v0]) \n\t"
235  "madd.s %[temp3], %[temp3], %[temp18], %[temp19] \n\t"
236  "lwc1 %[temp11], 524(%[s0]) \n\t"
237  "lwc1 %[temp12], 1792(%[v0]) \n\t"
238  "lwc1 %[temp13], 768(%[s0]) \n\t"
239  "lwc1 %[temp14], 1796(%[v0]) \n\t"
240  "lwc1 %[temp15], 772(%[s0]) \n\t"
241  "lwc1 %[temp16], 1800(%[v0]) \n\t"
242  "lwc1 %[temp17], 776(%[s0]) \n\t"
243  "lwc1 %[temp18], 1804(%[v0]) \n\t"
244  "lwc1 %[temp19], 780(%[s0]) \n\t"
245  "madd.s %[temp0], %[temp0], %[temp4], %[temp5] \n\t"
246  "lwc1 %[temp4], 2048(%[v0]) \n\t"
247  "madd.s %[temp1], %[temp1], %[temp6], %[temp7] \n\t"
248  "lwc1 %[temp5], 1024(%[s0]) \n\t"
249  "madd.s %[temp2], %[temp2], %[temp8], %[temp9] \n\t"
250  "lwc1 %[temp6], 2052(%[v0]) \n\t"
251  "madd.s %[temp3], %[temp3], %[temp10], %[temp11] \n\t"
252  "lwc1 %[temp7], 1028(%[s0]) \n\t"
253  "madd.s %[temp0], %[temp0], %[temp12], %[temp13] \n\t"
254  "lwc1 %[temp8], 2056(%[v0]) \n\t"
255  "madd.s %[temp1], %[temp1], %[temp14], %[temp15] \n\t"
256  "lwc1 %[temp9], 1032(%[s0]) \n\t"
257  "madd.s %[temp2], %[temp2], %[temp16], %[temp17] \n\t"
258  "lwc1 %[temp10], 2060(%[v0]) \n\t"
259  "madd.s %[temp3], %[temp3], %[temp18], %[temp19] \n\t"
260  "lwc1 %[temp11], 1036(%[s0]) \n\t"
261  "lwc1 %[temp12], 2816(%[v0]) \n\t"
262  "lwc1 %[temp13], 1280(%[s0]) \n\t"
263  "lwc1 %[temp14], 2820(%[v0]) \n\t"
264  "lwc1 %[temp15], 1284(%[s0]) \n\t"
265  "lwc1 %[temp16], 2824(%[v0]) \n\t"
266  "lwc1 %[temp17], 1288(%[s0]) \n\t"
267  "lwc1 %[temp18], 2828(%[v0]) \n\t"
268  "lwc1 %[temp19], 1292(%[s0]) \n\t"
269  "madd.s %[temp0], %[temp0], %[temp4], %[temp5] \n\t"
270  "lwc1 %[temp4], 3072(%[v0]) \n\t"
271  "madd.s %[temp1], %[temp1], %[temp6], %[temp7] \n\t"
272  "lwc1 %[temp5], 1536(%[s0]) \n\t"
273  "madd.s %[temp2], %[temp2], %[temp8], %[temp9] \n\t"
274  "lwc1 %[temp6], 3076(%[v0]) \n\t"
275  "madd.s %[temp3], %[temp3], %[temp10], %[temp11] \n\t"
276  "lwc1 %[temp7], 1540(%[s0]) \n\t"
277  "madd.s %[temp0], %[temp0], %[temp12], %[temp13] \n\t"
278  "lwc1 %[temp8], 3080(%[v0]) \n\t"
279  "madd.s %[temp1], %[temp1], %[temp14], %[temp15] \n\t"
280  "lwc1 %[temp9], 1544(%[s0]) \n\t"
281  "madd.s %[temp2], %[temp2], %[temp16], %[temp17] \n\t"
282  "lwc1 %[temp10], 3084(%[v0]) \n\t"
283  "madd.s %[temp3], %[temp3], %[temp18], %[temp19] \n\t"
284  "lwc1 %[temp11], 1548(%[s0]) \n\t"
285  "lwc1 %[temp12], 3840(%[v0]) \n\t"
286  "lwc1 %[temp13], 1792(%[s0]) \n\t"
287  "lwc1 %[temp14], 3844(%[v0]) \n\t"
288  "lwc1 %[temp15], 1796(%[s0]) \n\t"
289  "lwc1 %[temp16], 3848(%[v0]) \n\t"
290  "lwc1 %[temp17], 1800(%[s0]) \n\t"
291  "lwc1 %[temp18], 3852(%[v0]) \n\t"
292  "lwc1 %[temp19], 1804(%[s0]) \n\t"
293  "madd.s %[temp0], %[temp0], %[temp4], %[temp5] \n\t"
294  "lwc1 %[temp4], 4096(%[v0]) \n\t"
295  "madd.s %[temp1], %[temp1], %[temp6], %[temp7] \n\t"
296  "lwc1 %[temp5], 2048(%[s0]) \n\t"
297  "madd.s %[temp2], %[temp2], %[temp8], %[temp9] \n\t"
298  "lwc1 %[temp6], 4100(%[v0]) \n\t"
299  "madd.s %[temp3], %[temp3], %[temp10], %[temp11] \n\t"
300  "lwc1 %[temp7], 2052(%[s0]) \n\t"
301  "madd.s %[temp0], %[temp0], %[temp12], %[temp13] \n\t"
302  "lwc1 %[temp8], 4104(%[v0]) \n\t"
303  PTR_ADDIU "%[dst], %[dst], 16 \n\t"
304  "madd.s %[temp1], %[temp1], %[temp14], %[temp15] \n\t"
305  "lwc1 %[temp9], 2056(%[s0]) \n\t"
306  PTR_ADDIU " %[s0], %[s0], 16 \n\t"
307  "madd.s %[temp2], %[temp2], %[temp16], %[temp17] \n\t"
308  "lwc1 %[temp10], 4108(%[v0]) \n\t"
309  PTR_ADDIU " %[v0], %[v0], 16 \n\t"
310  "madd.s %[temp3], %[temp3], %[temp18], %[temp19] \n\t"
311  "lwc1 %[temp11], 2044(%[s0]) \n\t"
312  "lwc1 %[temp12], 4848(%[v0]) \n\t"
313  "lwc1 %[temp13], 2288(%[s0]) \n\t"
314  "lwc1 %[temp14], 4852(%[v0]) \n\t"
315  "lwc1 %[temp15], 2292(%[s0]) \n\t"
316  "lwc1 %[temp16], 4856(%[v0]) \n\t"
317  "lwc1 %[temp17], 2296(%[s0]) \n\t"
318  "lwc1 %[temp18], 4860(%[v0]) \n\t"
319  "lwc1 %[temp19], 2300(%[s0]) \n\t"
320  "madd.s %[temp0], %[temp0], %[temp4], %[temp5] \n\t"
321  "lwc1 %[temp4], 0(%[v0]) \n\t"
322  "madd.s %[temp1], %[temp1], %[temp6], %[temp7] \n\t"
323  "lwc1 %[temp5], 0(%[s0]) \n\t"
324  "madd.s %[temp2], %[temp2], %[temp8], %[temp9] \n\t"
325  "lwc1 %[temp6], 4(%[v0]) \n\t"
326  "madd.s %[temp3], %[temp3], %[temp10], %[temp11] \n\t"
327  "lwc1 %[temp7], 4(%[s0]) \n\t"
328  "madd.s %[temp0], %[temp0], %[temp12], %[temp13] \n\t"
329  "lwc1 %[temp8], 8(%[v0]) \n\t"
330  "madd.s %[temp1], %[temp1], %[temp14], %[temp15] \n\t"
331  "lwc1 %[temp9], 8(%[s0]) \n\t"
332  "madd.s %[temp2], %[temp2], %[temp16], %[temp17] \n\t"
333  "lwc1 %[temp10], 12(%[v0]) \n\t"
334  "madd.s %[temp3], %[temp3], %[temp18], %[temp19] \n\t"
335  "lwc1 %[temp11], 12(%[s0]) \n\t"
336  "lwc1 %[temp12], 768(%[v0]) \n\t"
337  "lwc1 %[temp13], 256(%[s0]) \n\t"
338  "lwc1 %[temp14], 772(%[v0]) \n\t"
339  "lwc1 %[temp15], 260(%[s0]) \n\t"
340  "lwc1 %[temp16], 776(%[v0]) \n\t"
341  "lwc1 %[temp17], 264(%[s0]) \n\t"
342  "lwc1 %[temp18], 780(%[v0]) \n\t"
343  "lwc1 %[temp19], 268(%[s0]) \n\t"
344  "swc1 %[temp0], -16(%[dst]) \n\t"
345  "swc1 %[temp1], -12(%[dst]) \n\t"
346  "swc1 %[temp2], -8(%[dst]) \n\t"
347  "bne %[v0], %[v0_end], 1b \n\t"
348  " swc1 %[temp3], -4(%[dst]) \n\t"
349  "mul.s %[temp0], %[temp4], %[temp5] \n\t"
350  "lwc1 %[temp4], 1024(%[v0]) \n\t"
351  "mul.s %[temp1], %[temp6], %[temp7] \n\t"
352  "lwc1 %[temp5], 512(%[s0]) \n\t"
353  "mul.s %[temp2], %[temp8], %[temp9] \n\t"
354  "lwc1 %[temp6], 1028(%[v0]) \n\t"
355  "mul.s %[temp3], %[temp10], %[temp11] \n\t"
356  "lwc1 %[temp7], 516(%[s0]) \n\t"
357  "madd.s %[temp0], %[temp0], %[temp12], %[temp13] \n\t"
358  "lwc1 %[temp8], 1032(%[v0]) \n\t"
359  "madd.s %[temp1], %[temp1], %[temp14], %[temp15] \n\t"
360  "lwc1 %[temp9], 520(%[s0]) \n\t"
361  "madd.s %[temp2], %[temp2], %[temp16], %[temp17] \n\t"
362  "lwc1 %[temp10], 1036(%[v0]) \n\t"
363  "madd.s %[temp3], %[temp3], %[temp18], %[temp19] \n\t"
364  "lwc1 %[temp11], 524(%[s0]) \n\t"
365  "lwc1 %[temp12], 1792(%[v0]) \n\t"
366  "lwc1 %[temp13], 768(%[s0]) \n\t"
367  "lwc1 %[temp14], 1796(%[v0]) \n\t"
368  "lwc1 %[temp15], 772(%[s0]) \n\t"
369  "lwc1 %[temp16], 1800(%[v0]) \n\t"
370  "lwc1 %[temp17], 776(%[s0]) \n\t"
371  "lwc1 %[temp18], 1804(%[v0]) \n\t"
372  "lwc1 %[temp19], 780(%[s0]) \n\t"
373  "madd.s %[temp0], %[temp0], %[temp4], %[temp5] \n\t"
374  "lwc1 %[temp4], 2048(%[v0]) \n\t"
375  "madd.s %[temp1], %[temp1], %[temp6], %[temp7] \n\t"
376  "lwc1 %[temp5], 1024(%[s0]) \n\t"
377  "madd.s %[temp2], %[temp2], %[temp8], %[temp9] \n\t"
378  "lwc1 %[temp6], 2052(%[v0]) \n\t"
379  "madd.s %[temp3], %[temp3], %[temp10], %[temp11] \n\t"
380  "lwc1 %[temp7], 1028(%[s0]) \n\t"
381  "madd.s %[temp0], %[temp0], %[temp12], %[temp13] \n\t"
382  "lwc1 %[temp8], 2056(%[v0]) \n\t"
383  "madd.s %[temp1], %[temp1], %[temp14], %[temp15] \n\t"
384  "lwc1 %[temp9], 1032(%[s0]) \n\t"
385  "madd.s %[temp2], %[temp2], %[temp16], %[temp17] \n\t"
386  "lwc1 %[temp10], 2060(%[v0]) \n\t"
387  "madd.s %[temp3], %[temp3], %[temp18], %[temp19] \n\t"
388  "lwc1 %[temp11], 1036(%[s0]) \n\t"
389  "lwc1 %[temp12], 2816(%[v0]) \n\t"
390  "lwc1 %[temp13], 1280(%[s0]) \n\t"
391  "lwc1 %[temp14], 2820(%[v0]) \n\t"
392  "lwc1 %[temp15], 1284(%[s0]) \n\t"
393  "lwc1 %[temp16], 2824(%[v0]) \n\t"
394  "lwc1 %[temp17], 1288(%[s0]) \n\t"
395  "lwc1 %[temp18], 2828(%[v0]) \n\t"
396  "lwc1 %[temp19], 1292(%[s0]) \n\t"
397  "madd.s %[temp0], %[temp0], %[temp4], %[temp5] \n\t"
398  "lwc1 %[temp4], 3072(%[v0]) \n\t"
399  "madd.s %[temp1], %[temp1], %[temp6], %[temp7] \n\t"
400  "lwc1 %[temp5], 1536(%[s0]) \n\t"
401  "madd.s %[temp2], %[temp2], %[temp8], %[temp9] \n\t"
402  "lwc1 %[temp6], 3076(%[v0]) \n\t"
403  "madd.s %[temp3], %[temp3], %[temp10], %[temp11] \n\t"
404  "lwc1 %[temp7], 1540(%[s0]) \n\t"
405  "madd.s %[temp0], %[temp0], %[temp12], %[temp13] \n\t"
406  "lwc1 %[temp8], 3080(%[v0]) \n\t"
407  "madd.s %[temp1], %[temp1], %[temp14], %[temp15] \n\t"
408  "lwc1 %[temp9], 1544(%[s0]) \n\t"
409  "madd.s %[temp2], %[temp2], %[temp16], %[temp17] \n\t"
410  "lwc1 %[temp10], 3084(%[v0]) \n\t"
411  "madd.s %[temp3], %[temp3], %[temp18], %[temp19] \n\t"
412  "lwc1 %[temp11], 1548(%[s0]) \n\t"
413  "lwc1 %[temp12], 3840(%[v0]) \n\t"
414  "lwc1 %[temp13], 1792(%[s0]) \n\t"
415  "lwc1 %[temp14], 3844(%[v0]) \n\t"
416  "lwc1 %[temp15], 1796(%[s0]) \n\t"
417  "lwc1 %[temp16], 3848(%[v0]) \n\t"
418  "lwc1 %[temp17], 1800(%[s0]) \n\t"
419  "lwc1 %[temp18], 3852(%[v0]) \n\t"
420  "lwc1 %[temp19], 1804(%[s0]) \n\t"
421  "madd.s %[temp0], %[temp0], %[temp4], %[temp5] \n\t"
422  "lwc1 %[temp4], 4096(%[v0]) \n\t"
423  "madd.s %[temp1], %[temp1], %[temp6], %[temp7] \n\t"
424  "lwc1 %[temp5], 2048(%[s0]) \n\t"
425  "madd.s %[temp2], %[temp2], %[temp8], %[temp9] \n\t"
426  "lwc1 %[temp6], 4100(%[v0]) \n\t"
427  "madd.s %[temp3], %[temp3], %[temp10], %[temp11] \n\t"
428  "lwc1 %[temp7], 2052(%[s0]) \n\t"
429  "madd.s %[temp0], %[temp0], %[temp12], %[temp13] \n\t"
430  "lwc1 %[temp8], 4104(%[v0]) \n\t"
431  "madd.s %[temp1], %[temp1], %[temp14], %[temp15] \n\t"
432  "lwc1 %[temp9], 2056(%[s0]) \n\t"
433  "madd.s %[temp2], %[temp2], %[temp16], %[temp17] \n\t"
434  "lwc1 %[temp10], 4108(%[v0]) \n\t"
435  "madd.s %[temp3], %[temp3], %[temp18], %[temp19] \n\t"
436  "lwc1 %[temp11], 2060(%[s0]) \n\t"
437  "lwc1 %[temp12], 4864(%[v0]) \n\t"
438  "lwc1 %[temp13], 2304(%[s0]) \n\t"
439  "lwc1 %[temp14], 4868(%[v0]) \n\t"
440  "lwc1 %[temp15], 2308(%[s0]) \n\t"
441  "madd.s %[temp0], %[temp0], %[temp4], %[temp5] \n\t"
442  "lwc1 %[temp16], 4872(%[v0]) \n\t"
443  "madd.s %[temp1], %[temp1], %[temp6], %[temp7] \n\t"
444  "lwc1 %[temp17], 2312(%[s0]) \n\t"
445  "madd.s %[temp2], %[temp2], %[temp8], %[temp9] \n\t"
446  "lwc1 %[temp18], 4876(%[v0]) \n\t"
447  "madd.s %[temp3], %[temp3], %[temp10], %[temp11] \n\t"
448  "lwc1 %[temp19], 2316(%[s0]) \n\t"
449  "madd.s %[temp0], %[temp0], %[temp12], %[temp13] \n\t"
450  PTR_ADDIU "%[dst], %[dst], 16 \n\t"
451  "madd.s %[temp1], %[temp1], %[temp14], %[temp15] \n\t"
452  "madd.s %[temp2], %[temp2], %[temp16], %[temp17] \n\t"
453  "madd.s %[temp3], %[temp3], %[temp18], %[temp19] \n\t"
454  "swc1 %[temp0], -16(%[dst]) \n\t"
455  "swc1 %[temp1], -12(%[dst]) \n\t"
456  "swc1 %[temp2], -8(%[dst]) \n\t"
457  "swc1 %[temp3], -4(%[dst]) \n\t"
458  ".set pop \n\t"
459 
460  : [dst]"+r"(dst), [v0]"+r"(vv0), [s0]"+r"(s0),
461  [temp0]"=&f"(temp0), [temp1]"=&f"(temp1), [temp2]"=&f"(temp2),
462  [temp3]"=&f"(temp3), [temp4]"=&f"(temp4), [temp5]"=&f"(temp5),
463  [temp6]"=&f"(temp6), [temp7]"=&f"(temp7), [temp8]"=&f"(temp8),
464  [temp9]"=&f"(temp9), [temp10]"=&f"(temp10), [temp11]"=&f"(temp11),
465  [temp12]"=&f"(temp12), [temp13]"=&f"(temp13), [temp14]"=&f"(temp14),
466  [temp15]"=&f"(temp15), [temp16]"=&f"(temp16), [temp17]"=&f"(temp17),
467  [temp18]"=&f"(temp18), [temp19]"=&f"(temp19)
468  : [v0_end]"r"(v0_end)
469  : "memory"
470  );
471  }
472  else
473  {
474  fdsp->vector_fmul (out, v , sbr_qmf_window , 64 >> div);
475  fdsp->vector_fmul_add(out, v + ( 192 >> div), sbr_qmf_window + ( 64 >> div), out , 64 >> div);
476  fdsp->vector_fmul_add(out, v + ( 256 >> div), sbr_qmf_window + (128 >> div), out , 64 >> div);
477  fdsp->vector_fmul_add(out, v + ( 448 >> div), sbr_qmf_window + (192 >> div), out , 64 >> div);
478  fdsp->vector_fmul_add(out, v + ( 512 >> div), sbr_qmf_window + (256 >> div), out , 64 >> div);
479  fdsp->vector_fmul_add(out, v + ( 704 >> div), sbr_qmf_window + (320 >> div), out , 64 >> div);
480  fdsp->vector_fmul_add(out, v + ( 768 >> div), sbr_qmf_window + (384 >> div), out , 64 >> div);
481  fdsp->vector_fmul_add(out, v + ( 960 >> div), sbr_qmf_window + (448 >> div), out , 64 >> div);
482  fdsp->vector_fmul_add(out, v + (1024 >> div), sbr_qmf_window + (512 >> div), out , 64 >> div);
483  fdsp->vector_fmul_add(out, v + (1216 >> div), sbr_qmf_window + (576 >> div), out , 64 >> div);
484  out += 64 >> div;
485  }
486  }
487 }
488 
489 #define sbr_qmf_analysis sbr_qmf_analysis_mips
490 #define sbr_qmf_synthesis sbr_qmf_synthesis_mips
491 
492 #endif /* !HAVE_MIPS32R6 && !HAVE_MIPS64R6 */
493 #endif /* HAVE_MIPSFPU */
494 #endif /* HAVE_INLINE_ASM */
495 
496 #endif /* AVCODEC_MIPS_AACSBR_MIPS_H */
SBRDSPContext
Definition: sbrdsp.h:27
out
FILE * out
Definition: movenc.c:54
AVTXContext
Definition: tx_priv.h:235
AVFloatDSPContext::vector_fmul_reverse
void(* vector_fmul_reverse)(float *dst, const float *src0, const float *src1, int len)
Calculate the entry wise product of two vectors of floats, and store the result in a vector of floats...
Definition: float_dsp.h:152
step
trying all byte sequences megabyte in length and selecting the best looking sequence will yield cases to try But a word about which is also called distortion Distortion can be quantified by almost any quality measurement one chooses the sum of squared differences is used but more complex methods that consider psychovisual effects can be used as well It makes no difference in this discussion First step
Definition: rate_distortion.txt:58
asmdefs.h
sbr.h
v0
#define v0
Definition: regdef.h:26
sbr_qmf_window_ds
static INTFLOAT sbr_qmf_window_ds[320]
Definition: aacsbrdata.h:45
av_tx_fn
void(* av_tx_fn)(AVTXContext *s, void *out, void *in, ptrdiff_t stride)
Function pointer to a function to perform the transform.
Definition: tx.h:151
SBRDSPContext::neg_odd_64
void(* neg_odd_64)(INTFLOAT *x)
Definition: sbrdsp.h:30
SBR_SYNTHESIS_BUF_SIZE
#define SBR_SYNTHESIS_BUF_SIZE
Definition: sbr.h:58
SBRDSPContext::qmf_deint_bfly
void(* qmf_deint_bfly)(INTFLOAT *v, const INTFLOAT *src0, const INTFLOAT *src1)
Definition: sbrdsp.h:34
SBRDSPContext::sum64x5
void(* sum64x5)(INTFLOAT *z)
Definition: sbrdsp.h:28
X
@ X
Definition: vf_addroi.c:27
SBRDSPContext::qmf_pre_shuffle
void(* qmf_pre_shuffle)(INTFLOAT *z)
Definition: sbrdsp.h:31
AVFloatDSPContext::vector_fmul
void(* vector_fmul)(float *dst, const float *src0, const float *src1, int len)
Calculate the entry wise product of two vectors of floats and store the result in a vector of floats.
Definition: float_dsp.h:36
AVFloatDSPContext
Definition: float_dsp.h:22
SBRDSPContext::qmf_deint_neg
void(* qmf_deint_neg)(INTFLOAT *v, const INTFLOAT *src)
Definition: sbrdsp.h:33
W
@ W
Definition: vf_addroi.c:27
i
#define i(width, name, range_min, range_max)
Definition: cbs_h2645.c:255
AVFloatDSPContext::vector_fmul_add
void(* vector_fmul_add)(float *dst, const float *src0, const float *src1, const float *src2, int len)
Calculate the entry wise product of two vectors of floats, add a third vector of floats and store the...
Definition: float_dsp.h:135
__asm__
__asm__(".macro parse_r var r\n\t" "\\var = -1\n\t" _IFC_REG(0) _IFC_REG(1) _IFC_REG(2) _IFC_REG(3) _IFC_REG(4) _IFC_REG(5) _IFC_REG(6) _IFC_REG(7) _IFC_REG(8) _IFC_REG(9) _IFC_REG(10) _IFC_REG(11) _IFC_REG(12) _IFC_REG(13) _IFC_REG(14) _IFC_REG(15) _IFC_REG(16) _IFC_REG(17) _IFC_REG(18) _IFC_REG(19) _IFC_REG(20) _IFC_REG(21) _IFC_REG(22) _IFC_REG(23) _IFC_REG(24) _IFC_REG(25) _IFC_REG(26) _IFC_REG(27) _IFC_REG(28) _IFC_REG(29) _IFC_REG(30) _IFC_REG(31) ".iflt \\var\n\t" ".error \"Unable to parse register name \\r\"\n\t" ".endif\n\t" ".endm")
SBRDSPContext::qmf_post_shuffle
void(* qmf_post_shuffle)(INTFLOAT W[32][2], const INTFLOAT *z)
Definition: sbrdsp.h:32
aacdec.h
sbr_qmf_window_us
static const INTFLOAT sbr_qmf_window_us[640]
Definition: aacsbrdata.h:48
PTR_ADDIU
#define PTR_ADDIU
Definition: asmdefs.h:50
s0
#define s0
Definition: regdef.h:37