58 #define ENVELOPE_ADJUSTMENT_OFFSET 2 62 float X_low[32][40][2],
const float W[2][32][32][2],
66 int temp0, temp1, temp2, temp3, temp4, temp5, temp6, temp7;
67 float *p_x_low = &X_low[0][8][0];
68 float *p_w = (
float*)&W[buf_idx][0][0][0];
69 float *p_x1_low = &X_low[0][0][0];
70 float *p_w1 = (
float*)&W[1-buf_idx][24][0][0];
72 float *loop_end=p_x1_low + 2560;
77 "sw $0, 0(%[p_x1_low]) \n\t" 78 "sw $0, 4(%[p_x1_low]) \n\t" 79 "sw $0, 8(%[p_x1_low]) \n\t" 80 "sw $0, 12(%[p_x1_low]) \n\t" 81 "sw $0, 16(%[p_x1_low]) \n\t" 82 "sw $0, 20(%[p_x1_low]) \n\t" 83 "sw $0, 24(%[p_x1_low]) \n\t" 84 "sw $0, 28(%[p_x1_low]) \n\t" 85 PTR_ADDIU "%[p_x1_low],%[p_x1_low], 32 \n\t" 86 "bne %[p_x1_low], %[loop_end], 1b \n\t" 87 PTR_ADDIU "%[p_x1_low],%[p_x1_low], -10240 \n\t" 89 : [p_x1_low]
"+r"(p_x1_low)
90 : [loop_end]
"r"(loop_end)
94 for (k = 0; k < sbr->
kx[1]; k++) {
95 for (i = 0; i < 32; i+=4) {
98 "lw %[temp0], 0(%[p_w]) \n\t" 99 "lw %[temp1], 4(%[p_w]) \n\t" 100 "lw %[temp2], 256(%[p_w]) \n\t" 101 "lw %[temp3], 260(%[p_w]) \n\t" 102 "lw %[temp4], 512(%[p_w]) \n\t" 103 "lw %[temp5], 516(%[p_w]) \n\t" 104 "lw %[temp6], 768(%[p_w]) \n\t" 105 "lw %[temp7], 772(%[p_w]) \n\t" 106 "sw %[temp0], 0(%[p_x_low]) \n\t" 107 "sw %[temp1], 4(%[p_x_low]) \n\t" 108 "sw %[temp2], 8(%[p_x_low]) \n\t" 109 "sw %[temp3], 12(%[p_x_low]) \n\t" 110 "sw %[temp4], 16(%[p_x_low]) \n\t" 111 "sw %[temp5], 20(%[p_x_low]) \n\t" 112 "sw %[temp6], 24(%[p_x_low]) \n\t" 113 "sw %[temp7], 28(%[p_x_low]) \n\t" 114 PTR_ADDIU "%[p_x_low], %[p_x_low], 32 \n\t" 117 : [temp0]
"=&r"(temp0), [temp1]
"=&r"(temp1),
118 [temp2]
"=&r"(temp2), [temp3]
"=&r"(temp3),
119 [temp4]
"=&r"(temp4), [temp5]
"=&r"(temp5),
120 [temp6]
"=&r"(temp6), [temp7]
"=&r"(temp7),
121 [p_w]
"+r"(p_w), [p_x_low]
"+r"(p_x_low)
130 for (k = 0; k < sbr->
kx[0]; k++) {
131 for (i = 0; i < 2; i++) {
135 "lw %[temp0], 0(%[p_w1]) \n\t" 136 "lw %[temp1], 4(%[p_w1]) \n\t" 137 "lw %[temp2], 256(%[p_w1]) \n\t" 138 "lw %[temp3], 260(%[p_w1]) \n\t" 139 "lw %[temp4], 512(%[p_w1]) \n\t" 140 "lw %[temp5], 516(%[p_w1]) \n\t" 141 "lw %[temp6], 768(%[p_w1]) \n\t" 142 "lw %[temp7], 772(%[p_w1]) \n\t" 143 "sw %[temp0], 0(%[p_x1_low]) \n\t" 144 "sw %[temp1], 4(%[p_x1_low]) \n\t" 145 "sw %[temp2], 8(%[p_x1_low]) \n\t" 146 "sw %[temp3], 12(%[p_x1_low]) \n\t" 147 "sw %[temp4], 16(%[p_x1_low]) \n\t" 148 "sw %[temp5], 20(%[p_x1_low]) \n\t" 149 "sw %[temp6], 24(%[p_x1_low]) \n\t" 150 "sw %[temp7], 28(%[p_x1_low]) \n\t" 151 PTR_ADDIU "%[p_x1_low], %[p_x1_low], 32 \n\t" 154 : [temp0]
"=&r"(temp0), [temp1]
"=&r"(temp1),
155 [temp2]
"=&r"(temp2), [temp3]
"=&r"(temp3),
156 [temp4]
"=&r"(temp4), [temp5]
"=&r"(temp5),
157 [temp6]
"=&r"(temp6), [temp7]
"=&r"(temp7),
158 [p_w1]
"+r"(p_w1), [p_x1_low]
"+r"(p_x1_low)
170 const float Y0[38][64][2],
const float Y1[38][64][2],
171 const float X_low[32][40][2],
int ch)
175 int temp0, temp1, temp2, temp3;
176 const float *X_low1, *Y01, *Y11;
177 float *x1=&X[0][0][0];
184 "sw $0, 0(%[x1]) \n\t" 185 "sw $0, 4(%[x1]) \n\t" 186 "sw $0, 8(%[x1]) \n\t" 187 "sw $0, 12(%[x1]) \n\t" 188 "sw $0, 16(%[x1]) \n\t" 189 "sw $0, 20(%[x1]) \n\t" 190 "sw $0, 24(%[x1]) \n\t" 191 "sw $0, 28(%[x1]) \n\t" 193 "bne %[x1], %[j], 1b \n\t" 203 X_low1=&X_low[0][2][0];
205 for (k = 0; k < sbr->
kx[0]; k++) {
208 "move %[i], $zero \n\t" 210 "lw %[temp0], 0(%[X_low1]) \n\t" 211 "lw %[temp1], 4(%[X_low1]) \n\t" 212 "sw %[temp0], 0(%[x1]) \n\t" 213 "sw %[temp1], 9728(%[x1]) \n\t" 216 "addiu %[i], %[i], 1 \n\t" 217 "bne %[i], %[i_Temp], 2b \n\t" 219 : [x1]
"+r"(x1), [X_low1]
"+r"(X_low1), [i]
"=&r"(i),
220 [temp0]
"=&r"(temp0), [temp1]
"=&r"(temp1)
221 : [i_Temp]
"r"(i_Temp)
225 X_low1-=(i_Temp<<1)-80;
229 Y01=(
float*)&Y0[32][k][0];
231 for (; k < sbr->
kx[0] + sbr->
m[0]; k++) {
233 "move %[i], $zero \n\t" 235 "lw %[temp0], 0(%[Y01]) \n\t" 236 "lw %[temp1], 4(%[Y01]) \n\t" 237 "sw %[temp0], 0(%[x1]) \n\t" 238 "sw %[temp1], 9728(%[x1]) \n\t" 241 "addiu %[i], %[i], 1 \n\t" 242 "bne %[i], %[i_Temp], 3b \n\t" 244 : [x1]
"+r"(x1), [Y01]
"+r"(Y01), [i]
"=&r"(i),
245 [temp0]
"=&r"(temp0), [temp1]
"=&r"(temp1)
246 : [i_Temp]
"r"(i_Temp)
255 X_low1=&X_low[0][i_Temp+2][0];
258 for (k = 0; k < sbr->
kx[1]; k++) {
261 "move %[i], %[i_Temp] \n\t" 263 "lw %[temp0], 0(%[X_low1]) \n\t" 264 "lw %[temp1], 4(%[X_low1]) \n\t" 265 "sw %[temp0], 0(%[x1]) \n\t" 266 "sw %[temp1], 9728(%[x1]) \n\t" 269 "addiu %[i], %[i], 1 \n\t" 270 "bne %[i], %[temp3], 4b \n\t" 272 : [x1]
"+r"(x1), [X_low1]
"+r"(X_low1), [i]
"=&r"(i),
273 [temp0]
"=&r"(temp0), [temp1]
"=&r"(temp1),
275 : [i_Temp]
"r"(i_Temp), [temp3]
"r"(temp3)
278 x1 -= ((38-i_Temp)<<6)-1;
279 X_low1 -= ((38-i_Temp)<<1)- 80;
283 Y11=&Y1[i_Temp][k][0];
286 for (; k < sbr->
kx[1] + sbr->
m[1]; k++) {
289 "move %[i], %[i_Temp] \n\t" 291 "lw %[temp0], 0(%[Y11]) \n\t" 292 "lw %[temp1], 4(%[Y11]) \n\t" 293 "sw %[temp0], 0(%[x1]) \n\t" 294 "sw %[temp1], 9728(%[x1]) \n\t" 297 "addiu %[i], %[i], 1 \n\t" 298 "bne %[i], %[temp2], 5b \n\t" 300 : [x1]
"+r"(x1), [Y11]
"+r"(Y11), [i]
"=&r"(i),
301 [temp0]
"=&r"(temp0), [temp1]
"=&r"(temp1)
302 : [i_Temp]
"r"(i_Temp), [temp3]
"r"(temp3),
307 x1 -= ((32-i_Temp)<<6)-1;
308 Y11 -= ((32-i_Temp)<<7)-2;
314 #if !HAVE_MIPS32R6 && !HAVE_MIPS64R6 315 static void sbr_hf_assemble_mips(
float Y1[38][64][2],
316 const float X_high[64][40][2],
322 const int kx = sbr->
kx[1];
323 const int m_max = sbr->
m[1];
324 static const float h_smooth[5] = {
332 float (*g_temp)[48] = ch_data->
g_temp, (*q_temp)[48] = ch_data->
q_temp;
335 float *g_temp1, *q_temp1, *pok, *pok1;
336 float temp1, temp2, temp3, temp4;
340 for (i = 0; i < h_SL; i++) {
341 memcpy(g_temp[i + 2*ch_data->
t_env[0]], sbr->
gain[0], m_max *
sizeof(sbr->
gain[0][0]));
342 memcpy(q_temp[i + 2*ch_data->
t_env[0]], sbr->
q_m[0], m_max *
sizeof(sbr->
q_m[0][0]));
350 for (i = 2 * ch_data->
t_env[e]; i < 2 * ch_data->t_env[e + 1]; i++) {
351 g_temp1 = g_temp[h_SL + i];
353 q_temp1 = q_temp[h_SL + i];
357 for (j=0; j<(size>>2); j++) {
359 "lw %[temp1], 0(%[pok]) \n\t" 360 "lw %[temp2], 4(%[pok]) \n\t" 361 "lw %[temp3], 8(%[pok]) \n\t" 362 "lw %[temp4], 12(%[pok]) \n\t" 363 "sw %[temp1], 0(%[g_temp1]) \n\t" 364 "sw %[temp2], 4(%[g_temp1]) \n\t" 365 "sw %[temp3], 8(%[g_temp1]) \n\t" 366 "sw %[temp4], 12(%[g_temp1]) \n\t" 367 "lw %[temp1], 0(%[pok1]) \n\t" 368 "lw %[temp2], 4(%[pok1]) \n\t" 369 "lw %[temp3], 8(%[pok1]) \n\t" 370 "lw %[temp4], 12(%[pok1]) \n\t" 371 "sw %[temp1], 0(%[q_temp1]) \n\t" 372 "sw %[temp2], 4(%[q_temp1]) \n\t" 373 "sw %[temp3], 8(%[q_temp1]) \n\t" 374 "sw %[temp4], 12(%[q_temp1]) \n\t" 376 PTR_ADDIU "%[g_temp1], %[g_temp1], 16 \n\t" 378 PTR_ADDIU "%[q_temp1], %[q_temp1], 16 \n\t" 380 : [temp1]
"=&r"(temp1), [temp2]
"=&r"(temp2),
381 [temp3]
"=&r"(temp3), [temp4]
"=&r"(temp4),
382 [pok]
"+r"(pok), [g_temp1]
"+r"(g_temp1),
383 [pok1]
"+r"(pok1), [q_temp1]
"+r"(q_temp1)
389 for (j=0; j<(size&3); j++) {
391 "lw %[temp1], 0(%[pok]) \n\t" 392 "lw %[temp2], 0(%[pok1]) \n\t" 393 "sw %[temp1], 0(%[g_temp1]) \n\t" 394 "sw %[temp2], 0(%[q_temp1]) \n\t" 396 PTR_ADDIU "%[g_temp1], %[g_temp1], 4 \n\t" 398 PTR_ADDIU "%[q_temp1], %[q_temp1], 4 \n\t" 400 : [temp1]
"=&r"(temp1), [temp2]
"=&r"(temp2),
401 [temp3]
"=&r"(temp3), [temp4]
"=&r"(temp4),
402 [pok]
"+r"(pok), [g_temp1]
"+r"(g_temp1),
403 [pok1]
"+r"(pok1), [q_temp1]
"+r"(q_temp1)
412 for (i = 2 * ch_data->
t_env[e]; i < 2 * ch_data->t_env[e + 1]; i++) {
415 float *g_filt, *q_filt;
417 if (h_SL && e != e_a[0] && e != e_a[1]) {
421 for (m = 0; m < m_max; m++) {
422 const int idx1 = i + h_SL;
426 for (j = 0; j <= h_SL; j++) {
427 g_filt[m] += g_temp[idx1 - j][m] * h_smooth[j];
428 q_filt[m] += q_temp[idx1 - j][m] * h_smooth[j];
432 g_filt = g_temp[i + h_SL];
436 sbr->
dsp.
hf_g_filt(Y1[i] + kx, X_high + kx, g_filt, m_max,
439 if (e != e_a[0] && e != e_a[1]) {
444 int idx = indexsine&1;
445 int A = (1-((indexsine+(kx & 1))&2));
446 int B = (A^(-idx)) + idx;
447 float *
out = &Y1[i][kx][idx];
448 float *
in = sbr->
s_m[e];
449 float temp0, temp1, temp2, temp3, temp4, temp5;
450 float A_f = (float)A;
451 float B_f = (float)B;
453 for (m = 0; m+1 < m_max; m+=2) {
459 "lwc1 %[temp0], 0(%[in]) \n\t" 460 "lwc1 %[temp1], 4(%[in]) \n\t" 461 "madd.s %[temp4], %[temp2], %[temp0], %[A_f] \n\t" 462 "madd.s %[temp5], %[temp3], %[temp1], %[B_f] \n\t" 463 "swc1 %[temp4], 0(%[out]) \n\t" 464 "swc1 %[temp5], 8(%[out]) \n\t" 468 : [temp0]
"=&f" (temp0), [temp1]
"=&f"(temp1),
469 [temp4]
"=&f" (temp4), [temp5]
"=&f"(temp5),
470 [
in]
"+r"(
in), [out]
"+r"(out)
471 : [A_f]
"f"(A_f), [B_f]
"f"(B_f), [temp2]
"f"(temp2),
477 out[2*m ] += in[m ] *
A;
479 indexnoise = (indexnoise + m_max) & 0x1ff;
480 indexsine = (indexsine + 1) & 3;
488 float (*alpha0)[2],
float (*alpha1)[2],
489 const float X_low[32][40][2],
int k0)
492 float temp0, temp1, temp2, temp3, temp4, temp5, temp6, temp7,
c;
493 float *phi1, *alpha_1, *alpha_0, res1, res2, temp_real, temp_im;
497 for (k = 0; k < k0; k++) {
500 phi1 = &phi[0][0][0];
501 alpha_1 = &alpha1[k][0];
502 alpha_0 = &alpha0[k][0];
506 "lwc1 %[temp0], 40(%[phi1]) \n\t" 507 "lwc1 %[temp1], 16(%[phi1]) \n\t" 508 "lwc1 %[temp2], 24(%[phi1]) \n\t" 509 "lwc1 %[temp3], 28(%[phi1]) \n\t" 510 "mul.s %[dk], %[temp0], %[temp1] \n\t" 511 "lwc1 %[temp4], 0(%[phi1]) \n\t" 512 "mul.s %[res2], %[temp2], %[temp2] \n\t" 513 "lwc1 %[temp5], 4(%[phi1]) \n\t" 514 "madd.s %[res2], %[res2], %[temp3], %[temp3] \n\t" 515 "lwc1 %[temp6], 8(%[phi1]) \n\t" 516 "div.s %[res2], %[res2], %[c] \n\t" 517 "lwc1 %[temp0], 12(%[phi1]) \n\t" 518 "sub.s %[dk], %[dk], %[res2] \n\t" 520 : [temp0]
"=&f"(temp0), [temp1]
"=&f"(temp1), [temp2]
"=&f"(temp2),
521 [temp3]
"=&f"(temp3), [temp4]
"=&f"(temp4), [temp5]
"=&f"(temp5),
522 [temp6]
"=&f"(temp6), [res2]
"=&f"(res2), [dk]
"=&f"(dk)
523 : [phi1]
"r"(phi1), [
c]
"f"(
c)
532 "mul.s %[temp_real], %[temp4], %[temp2] \n\t" 533 "nmsub.s %[temp_real], %[temp_real], %[temp5], %[temp3] \n\t" 534 "nmsub.s %[temp_real], %[temp_real], %[temp6], %[temp1] \n\t" 535 "mul.s %[temp_im], %[temp4], %[temp3] \n\t" 536 "madd.s %[temp_im], %[temp_im], %[temp5], %[temp2] \n\t" 537 "nmsub.s %[temp_im], %[temp_im], %[temp0], %[temp1] \n\t" 538 "div.s %[temp_real], %[temp_real], %[dk] \n\t" 539 "div.s %[temp_im], %[temp_im], %[dk] \n\t" 540 "swc1 %[temp_real], 0(%[alpha_1]) \n\t" 541 "swc1 %[temp_im], 4(%[alpha_1]) \n\t" 543 : [temp_real]
"=&f" (temp_real), [temp_im]
"=&f"(temp_im)
544 : [phi1]
"r"(phi1), [temp0]
"f"(temp0), [temp1]
"f"(temp1),
545 [temp2]
"f"(temp2), [temp3]
"f"(temp3), [temp4]
"f"(temp4),
546 [temp5]
"f"(temp5), [temp6]
"f"(temp6),
547 [alpha_1]
"r"(alpha_1), [dk]
"f"(dk)
557 "lwc1 %[temp6], 0(%[alpha_1]) \n\t" 558 "lwc1 %[temp7], 4(%[alpha_1]) \n\t" 559 "mul.s %[temp_real], %[temp6], %[temp2] \n\t" 560 "add.s %[temp_real], %[temp_real], %[temp4] \n\t" 561 "madd.s %[temp_real], %[temp_real], %[temp7], %[temp3] \n\t" 562 "mul.s %[temp_im], %[temp7], %[temp2] \n\t" 563 "add.s %[temp_im], %[temp_im], %[temp5] \n\t" 564 "nmsub.s %[temp_im], %[temp_im], %[temp6], %[temp3] \n\t" 565 "div.s %[temp_real], %[temp_real], %[temp1] \n\t" 566 "div.s %[temp_im], %[temp_im], %[temp1] \n\t" 567 "neg.s %[temp_real], %[temp_real] \n\t" 568 "neg.s %[temp_im], %[temp_im] \n\t" 569 "swc1 %[temp_real], 0(%[alpha_0]) \n\t" 570 "swc1 %[temp_im], 4(%[alpha_0]) \n\t" 572 : [temp_real]
"=&f"(temp_real), [temp_im]
"=&f"(temp_im),
573 [temp6]
"=&f"(temp6), [temp7]
"=&f"(temp7),
574 [res1]
"=&f"(res1), [res2]
"=&f"(res2)
575 : [alpha_1]
"r"(alpha_1), [alpha_0]
"r"(alpha_0),
576 [temp0]
"f"(temp0), [temp1]
"f"(temp1), [temp2]
"f"(temp2),
577 [temp3]
"f"(temp3), [temp4]
"f"(temp4), [temp5]
"f"(temp5)
583 "lwc1 %[temp1], 0(%[alpha_1]) \n\t" 584 "lwc1 %[temp2], 4(%[alpha_1]) \n\t" 585 "lwc1 %[temp_real], 0(%[alpha_0]) \n\t" 586 "lwc1 %[temp_im], 4(%[alpha_0]) \n\t" 587 "mul.s %[res1], %[temp1], %[temp1] \n\t" 588 "madd.s %[res1], %[res1], %[temp2], %[temp2] \n\t" 589 "mul.s %[res2], %[temp_real], %[temp_real] \n\t" 590 "madd.s %[res2], %[res2], %[temp_im], %[temp_im] \n\t" 592 : [temp_real]
"=&f"(temp_real), [temp_im]
"=&f"(temp_im),
593 [temp1]
"=&f"(temp1), [temp2]
"=&f"(temp2),
594 [res1]
"=&f"(res1), [res2]
"=&f"(res2)
595 : [alpha_1]
"r"(alpha_1), [alpha_0]
"r"(alpha_0)
599 if (res1 >= 16.0f || res2 >= 16.0f) {
617 #if !HAVE_MIPS32R6 && !HAVE_MIPS64R6
MIPS assembly defines from sys/asm.h but rewritten for use with C inline assembly (rather than from w...
int(* sbr_lf_gen)(AACContext *ac, SpectralBandReplication *sbr, INTFLOAT X_low[32][40][2], const INTFLOAT W[2][32][32][2], int buf_idx)
unsigned bs_smoothing_mode
#define ENVELOPE_ADJUSTMENT_OFFSET
uint8_t pi<< 24) CONV_FUNC(AV_SAMPLE_FMT_S64, int64_t, AV_SAMPLE_FMT_U8,(uint64_t)((*(const uint8_t *) pi - 0x80U))<< 56) CONV_FUNC(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_U8,(*(const uint8_t *) pi - 0x80) *(1.0f/(1<< 7))) CONV_FUNC(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_U8,(*(const uint8_t *) pi - 0x80) *(1.0/(1<< 7))) CONV_FUNC(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S16,(*(const int16_t *) pi >>8)+0x80) CONV_FUNC(AV_SAMPLE_FMT_S64, int64_t, AV_SAMPLE_FMT_S16,(uint64_t)(*(const int16_t *) pi)<< 48) CONV_FUNC(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S16, *(const int16_t *) pi *(1.0f/(1<< 15))) CONV_FUNC(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S16, *(const int16_t *) pi *(1.0/(1<< 15))) CONV_FUNC(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S32,(*(const int32_t *) pi >>24)+0x80) CONV_FUNC(AV_SAMPLE_FMT_S64, int64_t, AV_SAMPLE_FMT_S32,(uint64_t)(*(const int32_t *) pi)<< 32) CONV_FUNC(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S32, *(const int32_t *) pi *(1.0f/(1U<< 31))) CONV_FUNC(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S32, *(const int32_t *) pi *(1.0/(1U<< 31))) CONV_FUNC(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S64,(*(const int64_t *) pi >>56)+0x80) CONV_FUNC(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S64, *(const int64_t *) pi *(1.0f/(INT64_C(1)<< 63))) CONV_FUNC(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S64, *(const int64_t *) pi *(1.0/(INT64_C(1)<< 63))) CONV_FUNC(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_FLT, av_clip_uint8(lrintf(*(const float *) pi *(1<< 7))+0x80)) CONV_FUNC(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_FLT, av_clip_int16(lrintf(*(const float *) pi *(1<< 15)))) CONV_FUNC(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_FLT, av_clipl_int32(llrintf(*(const float *) pi *(1U<< 31)))) CONV_FUNC(AV_SAMPLE_FMT_S64, int64_t, AV_SAMPLE_FMT_FLT, llrintf(*(const float *) pi *(INT64_C(1)<< 63))) CONV_FUNC(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_DBL, av_clip_uint8(lrint(*(const double *) pi *(1<< 7))+0x80)) CONV_FUNC(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_DBL, av_clip_int16(lrint(*(const double *) pi *(1<< 15)))) CONV_FUNC(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_DBL, av_clipl_int32(llrint(*(const double *) pi *(1U<< 31)))) CONV_FUNC(AV_SAMPLE_FMT_S64, int64_t, AV_SAMPLE_FMT_DBL, llrint(*(const double *) pi *(INT64_C(1)<< 63))) #define FMT_PAIR_FUNC(out, in) static conv_func_type *const fmt_pair_to_conv_functions[AV_SAMPLE_FMT_NB *AV_SAMPLE_FMT_NB]={ FMT_PAIR_FUNC(AV_SAMPLE_FMT_U8, AV_SAMPLE_FMT_U8), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S16, AV_SAMPLE_FMT_U8), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S32, AV_SAMPLE_FMT_U8), FMT_PAIR_FUNC(AV_SAMPLE_FMT_FLT, AV_SAMPLE_FMT_U8), FMT_PAIR_FUNC(AV_SAMPLE_FMT_DBL, AV_SAMPLE_FMT_U8), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S64, AV_SAMPLE_FMT_U8), FMT_PAIR_FUNC(AV_SAMPLE_FMT_U8, AV_SAMPLE_FMT_S16), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S16, AV_SAMPLE_FMT_S16), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S32, AV_SAMPLE_FMT_S16), FMT_PAIR_FUNC(AV_SAMPLE_FMT_FLT, AV_SAMPLE_FMT_S16), FMT_PAIR_FUNC(AV_SAMPLE_FMT_DBL, AV_SAMPLE_FMT_S16), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S64, AV_SAMPLE_FMT_S16), FMT_PAIR_FUNC(AV_SAMPLE_FMT_U8, AV_SAMPLE_FMT_S32), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S16, AV_SAMPLE_FMT_S32), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S32, AV_SAMPLE_FMT_S32), FMT_PAIR_FUNC(AV_SAMPLE_FMT_FLT, AV_SAMPLE_FMT_S32), FMT_PAIR_FUNC(AV_SAMPLE_FMT_DBL, AV_SAMPLE_FMT_S32), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S64, AV_SAMPLE_FMT_S32), FMT_PAIR_FUNC(AV_SAMPLE_FMT_U8, AV_SAMPLE_FMT_FLT), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S16, AV_SAMPLE_FMT_FLT), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S32, AV_SAMPLE_FMT_FLT), FMT_PAIR_FUNC(AV_SAMPLE_FMT_FLT, AV_SAMPLE_FMT_FLT), FMT_PAIR_FUNC(AV_SAMPLE_FMT_DBL, AV_SAMPLE_FMT_FLT), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S64, AV_SAMPLE_FMT_FLT), FMT_PAIR_FUNC(AV_SAMPLE_FMT_U8, AV_SAMPLE_FMT_DBL), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S16, AV_SAMPLE_FMT_DBL), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S32, AV_SAMPLE_FMT_DBL), FMT_PAIR_FUNC(AV_SAMPLE_FMT_FLT, AV_SAMPLE_FMT_DBL), FMT_PAIR_FUNC(AV_SAMPLE_FMT_DBL, AV_SAMPLE_FMT_DBL), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S64, AV_SAMPLE_FMT_DBL), FMT_PAIR_FUNC(AV_SAMPLE_FMT_U8, AV_SAMPLE_FMT_S64), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S16, AV_SAMPLE_FMT_S64), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S32, AV_SAMPLE_FMT_S64), FMT_PAIR_FUNC(AV_SAMPLE_FMT_FLT, AV_SAMPLE_FMT_S64), FMT_PAIR_FUNC(AV_SAMPLE_FMT_DBL, AV_SAMPLE_FMT_S64), FMT_PAIR_FUNC(AV_SAMPLE_FMT_S64, AV_SAMPLE_FMT_S64), };static void cpy1(uint8_t **dst, const uint8_t **src, int len){ memcpy(*dst, *src, len);} static void cpy2(uint8_t **dst, const uint8_t **src, int len){ memcpy(*dst, *src, 2 *len);} static void cpy4(uint8_t **dst, const uint8_t **src, int len){ memcpy(*dst, *src, 4 *len);} static void cpy8(uint8_t **dst, const uint8_t **src, int len){ memcpy(*dst, *src, 8 *len);} AudioConvert *swri_audio_convert_alloc(enum AVSampleFormat out_fmt, enum AVSampleFormat in_fmt, int channels, const int *ch_map, int flags) { AudioConvert *ctx;conv_func_type *f=fmt_pair_to_conv_functions[av_get_packed_sample_fmt(out_fmt)+AV_SAMPLE_FMT_NB *av_get_packed_sample_fmt(in_fmt)];if(!f) return NULL;ctx=av_mallocz(sizeof(*ctx));if(!ctx) return NULL;if(channels==1){ in_fmt=av_get_planar_sample_fmt(in_fmt);out_fmt=av_get_planar_sample_fmt(out_fmt);} ctx->channels=channels;ctx->conv_f=f;ctx->ch_map=ch_map;if(in_fmt==AV_SAMPLE_FMT_U8||in_fmt==AV_SAMPLE_FMT_U8P) memset(ctx->silence, 0x80, sizeof(ctx->silence));if(out_fmt==in_fmt &&!ch_map) { switch(av_get_bytes_per_sample(in_fmt)){ case 1:ctx->simd_f=cpy1;break;case 2:ctx->simd_f=cpy2;break;case 4:ctx->simd_f=cpy4;break;case 8:ctx->simd_f=cpy8;break;} } if(HAVE_X86ASM &&HAVE_MMX) swri_audio_convert_init_x86(ctx, out_fmt, in_fmt, channels);if(ARCH_ARM) swri_audio_convert_init_arm(ctx, out_fmt, in_fmt, channels);if(ARCH_AARCH64) swri_audio_convert_init_aarch64(ctx, out_fmt, in_fmt, channels);return ctx;} void swri_audio_convert_free(AudioConvert **ctx) { av_freep(ctx);} int swri_audio_convert(AudioConvert *ctx, AudioData *out, AudioData *in, int len) { int ch;int off=0;const int os=(out->planar ? 1 :out->ch_count) *out->bps;unsigned misaligned=0;av_assert0(ctx->channels==out->ch_count);if(ctx->in_simd_align_mask) { int planes=in->planar ? in->ch_count :1;unsigned m=0;for(ch=0;ch< planes;ch++) m|=(intptr_t) in->ch[ch];misaligned|=m &ctx->in_simd_align_mask;} if(ctx->out_simd_align_mask) { int planes=out->planar ? out->ch_count :1;unsigned m=0;for(ch=0;ch< planes;ch++) m|=(intptr_t) out->ch[ch];misaligned|=m &ctx->out_simd_align_mask;} if(ctx->simd_f &&!ctx->ch_map &&!misaligned){ off=len &~15;av_assert1(off >=0);av_assert1(off<=len);av_assert2(ctx->channels==SWR_CH_MAX||!in->ch[ctx->channels]);if(off >0){ if(out->planar==in->planar){ int planes=out->planar ? out->ch_count :1;for(ch=0;ch< planes;ch++){ ctx->simd_f(out-> ch ch
AAC_SIGNE kx[2]
kx', and kx respectively, kx is the first QMF subband where SBR is used.
AAC_SIGNE m[2]
M' and M respectively, M is the number of QMF subbands that use SBR.
uint8_t t_env_num_env_old
Envelope time border of the last envelope of the previous frame.
AAC Spectral Band Replication function declarations.
void(* sbr_hf_inverse_filter)(SBRDSPContext *dsp, INTFLOAT(*alpha0)[2], INTFLOAT(*alpha1)[2], const INTFLOAT X_low[32][40][2], int k0)
AAC definitions and structures.
int(* sbr_x_gen)(SpectralBandReplication *sbr, INTFLOAT X[2][38][64], const INTFLOAT Y0[38][64][2], const INTFLOAT Y1[38][64][2], const INTFLOAT X_low[32][40][2], int ch)
uint8_t pi<< 24) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_U8, uint8_t,(*(const uint8_t *) pi - 0x80) *(1.0f/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_U8, uint8_t,(*(const uint8_t *) pi - 0x80) *(1.0/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S16, int16_t,(*(const int16_t *) pi >> 8)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S16, int16_t, *(const int16_t *) pi *(1.0f/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S16, int16_t, *(const int16_t *) pi *(1.0/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S32, int32_t,(*(const int32_t *) pi >> 24)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S32, int32_t, *(const int32_t *) pi *(1.0f/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S32, int32_t, *(const int32_t *) pi *(1.0/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_FLT, float, av_clip_uint8(lrintf(*(const float *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_FLT, float, av_clip_int16(lrintf(*(const float *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_FLT, float, av_clipl_int32(llrintf(*(const float *) pi *(1U<< 31)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_DBL, double, av_clip_uint8(lrint(*(const double *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_DBL, double, av_clip_int16(lrint(*(const double *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_DBL, double, av_clipl_int32(llrint(*(const double *) pi *(1U<< 31)))) #define SET_CONV_FUNC_GROUP(ofmt, ifmt) static void set_generic_function(AudioConvert *ac) { } void ff_audio_convert_free(AudioConvert **ac) { if(! *ac) return;ff_dither_free(&(*ac) ->dc);av_freep(ac);} AudioConvert *ff_audio_convert_alloc(AVAudioResampleContext *avr, enum AVSampleFormat out_fmt, enum AVSampleFormat in_fmt, int channels, int sample_rate, int apply_map) { AudioConvert *ac;int in_planar, out_planar;ac=av_mallocz(sizeof(*ac));if(!ac) return NULL;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);return NULL;} return ac;} 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;} else if(in_planar) ac->func_type=CONV_FUNC_TYPE_INTERLEAVE;else ac->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);return ac;} int ff_audio_convert(AudioConvert *ac, AudioData *out, AudioData *in) { int use_generic=1;int len=in->nb_samples;int p;if(ac->dc) { av_log(ac->avr, AV_LOG_TRACE, "%d samples - audio_convert: %s to %s (dithered)\", len, av_get_sample_fmt_name(ac->in_fmt), av_get_sample_fmt_name(ac->out_fmt));return ff_convert_dither(ac-> in
AAC_FLOAT q_m[7][48]
Amplitude adjusted noise scalefactors.
void(* autocorrelate)(const INTFLOAT x[40][2], AAC_FLOAT phi[3][2][2])
uint8_t t_env[8]
Envelope time borders.
aacsbr functions pointers
AAC_FLOAT s_m[7][48]
Sinusoidal levels.
Spectral Band Replication per channel data.
void(* hf_apply_noise[4])(INTFLOAT(*Y)[2], const AAC_FLOAT *s_m, const AAC_FLOAT *q_filt, int noise, int kx, int m_max)
void(* sbr_hf_assemble)(INTFLOAT Y1[38][64][2], const INTFLOAT X_high[64][40][2], SpectralBandReplication *sbr, SBRData *ch_data, const int e_a[2])
void(* hf_g_filt)(INTFLOAT(*Y)[2], const INTFLOAT(*X_high)[40][2], const AAC_FLOAT *g_filt, int m_max, intptr_t ixh)
#define LOCAL_ALIGNED_16(t, v,...)
Spectral Band Replication.
void ff_aacsbr_func_ptr_init_mips(AACSBRContext *c)