41 0, 31650, 28469, 23705, 18050, 12266, 7041, 2873,
42 0, -1597, -2147, -1992, -1492, -933, -484, -188,
50 0, 31915, 29436, 25569, 20676, 15206, 9639, 4439,
51 0, -3390, -5579, -6549, -6414, -5392, -3773, -1874,
52 0, 1595, 2727, 3303, 3319, 2850, 2030, 1023,
53 0, -887, -1527, -1860, -1876, -1614, -1150, -579,
54 0, 501, 859, 1041, 1044, 892, 631, 315,
55 0, -266, -453, -543, -538, -455, -317, -156,
56 0, 130, 218, 258, 253, 212, 147, 72,
57 0, -59, -101, -122, -123, -106, -77, -40,
65 18022, 9912, 5451, 2998, 1649, 907, 499, 274, 151, 83
73 22938, 16057, 11240, 7868, 5508, 3856, 2699, 1889, 1322, 925
91 for (n = subframe_size - 1; n >= 0; n--) {
93 for (i = 0; i < 10; i++)
94 sum += filter_coeffs[i] * in[n - i - 1];
96 out[
n] = in[
n] + (sum >> 12);
111 const int16_t* residual, int16_t *residual_filt,
121 int corr_int_num, corr_int_den;
126 int16_t gain_num,gain_den;
127 int16_t sh_gain_num, sh_gain_den;
130 int16_t gain_long_num,gain_long_den;
131 int16_t sh_gain_long_num, sh_gain_long_den;
133 int16_t best_delay_int, best_delay_frac;
135 int16_t delayed_signal_offset;
136 int lt_filt_factor_a, lt_filt_factor_b;
138 int16_t * selected_signal;
139 const int16_t * selected_signal_const;
147 tmp |=
FFABS(residual[i]);
156 sig_scaled[i] = residual[i] >> shift;
159 sig_scaled[i] = residual[i] << -shift;
165 sig_scaled + RES_PREV_DATA_SIZE,
191 best_delay_int = pitch_delay_int - 1;
192 for (i = pitch_delay_int - 1; i <= pitch_delay_int + 1; i++) {
194 sig_scaled + RES_PREV_DATA_SIZE - i,
196 if (sum > corr_int_num) {
204 sig_scaled - best_delay_int + RES_PREV_DATA_SIZE,
214 &sig_scaled[RES_PREV_DATA_SIZE - best_delay_int],
216 ANALYZED_FRAC_DELAYS+1,
231 &delayed_signal[k][1],
233 corr_den[k][0] = sum + delayed_signal[k][0 ] * delayed_signal[k][0 ];
234 corr_den[k][1] = sum + delayed_signal[k][subframe_size] * delayed_signal[k][subframe_size];
236 tmp =
FFMAX3(tmp, corr_den[k][0], corr_den[k][1]);
239 sh_gain_den =
av_log2(tmp) - 14;
240 if (sh_gain_den >= 0) {
242 sh_gain_num =
FFMAX(sh_gain_den, sh_ener);
247 delayed_signal_offset = 1;
249 gain_den = corr_int_den >> sh_gain_den;
250 gain_num = corr_int_num >> sh_gain_num;
251 gain_num_square = gain_num * gain_num;
253 for (i = 0; i < 2; i++) {
254 int16_t gain_num_short, gain_den_short;
255 int gain_num_short_square;
259 sig_scaled + RES_PREV_DATA_SIZE,
261 gain_num_short =
FFMAX(sum >> sh_gain_num, 0);
268 gain_num_short_square = gain_num_short * gain_num_short;
269 gain_den_short = corr_den[k][i] >> sh_gain_den;
276 gain_num = gain_num_short;
277 gain_den = gain_den_short;
278 gain_num_square = gain_num_short_square;
279 delayed_signal_offset = i;
280 best_delay_frac = k + 1;
290 L64_temp0 = (int64_t)gain_num_square << ((sh_gain_num << 1) + 1);
291 L64_temp1 = ((int64_t)gain_den * ener) << (sh_gain_den + sh_ener);
292 if (L64_temp0 < L64_temp1)
300 memcpy(residual_filt, residual + RES_PREV_DATA_SIZE, subframe_size *
sizeof(int16_t));
305 if (best_delay_frac) {
308 &sig_scaled[RES_PREV_DATA_SIZE - best_delay_int + delayed_signal_offset],
316 sig_scaled + RES_PREV_DATA_SIZE,
321 sh_gain_long_num = 0;
326 sh_gain_long_num = tmp;
335 sh_gain_long_den = tmp;
340 L_temp0 = gain_num * gain_num;
343 L_temp1 = gain_long_num * gain_long_num;
346 tmp = ((sh_gain_long_num - sh_gain_num) << 1) - (sh_gain_long_den - sh_gain_den);
353 if (L_temp1 > L_temp0) {
355 selected_signal = residual_filt;
356 gain_num = gain_long_num;
357 gain_den = gain_long_den;
358 sh_gain_num = sh_gain_long_num;
359 sh_gain_den = sh_gain_long_den;
362 selected_signal = &delayed_signal[best_delay_frac-1][delayed_signal_offset];
366 for (i = 0; i < subframe_size; i++)
367 selected_signal[i] <<= shift;
369 for (i = 0; i < subframe_size; i++)
370 selected_signal[i] >>= -shift;
373 selected_signal_const = selected_signal;
376 selected_signal_const = residual + RES_PREV_DATA_SIZE - (best_delay_int + 1 - delayed_signal_offset);
378 tmp = sh_gain_num - sh_gain_den;
384 if (gain_num > gain_den)
389 lt_filt_factor_a = (gain_den << 15) / (gain_den + gain_num);
392 L64_temp0 = ((int64_t)gain_num) << (sh_gain_num - 1);
393 L64_temp1 = ((int64_t)gain_den) << sh_gain_den;
398 lt_filt_factor_b = 32767 - lt_filt_factor_a + 1;
401 selected_signal_const,
402 lt_filt_factor_a, lt_filt_factor_b,
403 1<<14, 15, subframe_size);
425 const int16_t *lp_gd, int16_t* speech,
450 if (
FFABS(rh1) > rh0 || !rh0)
454 for (i = 0; i < 20; i++)
455 gain_term +=
FFABS(lp_gn[i + 10]);
458 if (gain_term > 0x400) {
459 temp = 0x2000000 / gain_term;
460 for (i = 0; i < subframe_size; i++)
461 speech[i] = (speech[i] * temp + 0x4000) >> 15;
464 return -(rh1 << 15) / rh0;
477 int subframe_size, int16_t ht_prev_data)
484 if (refl_coeff > 0) {
493 ga = (fact << 15) / av_clip_int16(32768 -
FFABS(gt));
497 tmp = res_pst[subframe_size - 1];
499 for (i = subframe_size - 1; i >= 1; i--) {
500 tmp2 = (res_pst[i] << 15) + ((gt * res_pst[i-1]) << 1);
501 tmp2 = (tmp2 + 0x4000) >> 15;
503 tmp2 = (tmp2 * ga * 2 + fact) >> sh_fact;
506 tmp2 = (res_pst[0] << 15) + ((gt * ht_prev_data) << 1);
507 tmp2 = (tmp2 + 0x4000) >> 15;
508 tmp2 = (tmp2 * ga * 2 + fact) >> sh_fact;
515 const int16_t *lp_filter_coeffs,
int pitch_delay_int,
516 int16_t* residual, int16_t* res_filter_data,
517 int16_t* pos_filter_data, int16_t *speech,
int subframe_size)
526 memset(lp_gn, 0, 33 *
sizeof(int16_t));
529 for (i = 0; i < 10; i++)
533 for (i = 0; i < 10; i++)
537 memcpy(speech - 10, res_filter_data, 10 *
sizeof(int16_t));
540 memcpy(res_filter_data, speech + subframe_size - 10, 10 *
sizeof(int16_t));
545 residual, residual_filt_buf + 10,
552 tilt_comp_coeff =
get_tilt_comp(dsp, lp_gn, lp_gd, residual_filt_buf + 10, subframe_size);
556 residual_filt_buf + 10,
557 subframe_size, 10, 0, 0, 0x800);
558 memcpy(pos_filter_data, pos_filter_data + subframe_size, 10 *
sizeof(int16_t));
560 *ht_prev_data =
apply_tilt_comp(speech, pos_filter_data + 10, tilt_comp_coeff,
561 subframe_size, *ht_prev_data);
575 int subframe_size, int16_t gain_prev)
579 int exp_before, exp_after;
581 if(!gain_after && gain_before)
586 exp_before = 14 -
av_log2(gain_before);
587 gain_before =
bidir_sal(gain_before, exp_before);
589 exp_after = 14 -
av_log2(gain_after);
590 gain_after =
bidir_sal(gain_after, exp_after);
592 if (gain_before < gain_after) {
593 gain = (gain_before << 15) / gain_after;
594 gain =
bidir_sal(gain, exp_after - exp_before - 1);
596 gain = ((gain_before - gain_after) << 14) / gain_after + 0x4000;
597 gain =
bidir_sal(gain, exp_after - exp_before);
603 for (n = 0; n < subframe_size; n++) {
606 gain_prev = av_clip_int16(gain + gain_prev);
607 speech[
n] = av_clip_int16((speech[n] * gain_prev + 0x2000) >> 14);