39 #define AMR_USE_16BIT_TABLES
66 float prediction_error[4];
80 float hpf_31_mem[2], hpf_400_mem[2];
104 for (i = 0; i < 4; i++)
147 for (i = 0; i < 9; i++)
148 isf_q[i] =
dico1_isf[ind[0]][i] * (1.0
f / (1 << 15));
150 for (i = 0; i < 7; i++)
151 isf_q[i + 9] =
dico2_isf[ind[1]][i] * (1.0
f / (1 << 15));
153 for (i = 0; i < 5; i++)
156 for (i = 0; i < 4; i++)
159 for (i = 0; i < 7; i++)
174 for (i = 0; i < 9; i++)
175 isf_q[i] =
dico1_isf[ind[0]][i] * (1.0
f / (1 << 15));
177 for (i = 0; i < 7; i++)
178 isf_q[i + 9] =
dico2_isf[ind[1]][i] * (1.0
f / (1 << 15));
180 for (i = 0; i < 3; i++)
181 isf_q[i] +=
dico21_isf[ind[2]][i] * (1.0
f / (1 << 15));
183 for (i = 0; i < 3; i++)
184 isf_q[i + 3] +=
dico22_isf[ind[3]][i] * (1.0
f / (1 << 15));
186 for (i = 0; i < 3; i++)
187 isf_q[i + 6] +=
dico23_isf[ind[4]][i] * (1.0
f / (1 << 15));
189 for (i = 0; i < 3; i++)
190 isf_q[i + 9] +=
dico24_isf[ind[5]][i] * (1.0
f / (1 << 15));
192 for (i = 0; i < 4; i++)
193 isf_q[i + 12] +=
dico25_isf[ind[6]][i] * (1.0
f / (1 << 15));
211 isf_q[i] +=
isf_mean[i] * (1.0f / (1 << 15));
228 for (k = 0; k < 3; k++) {
231 isp_q[k][i] = (1.0 - c) * isp4_past[i] + c * isp_q[3][i];
247 uint8_t *base_lag_int,
int subframe)
249 if (subframe == 0 || subframe == 2) {
250 if (pitch_index < 376) {
251 *lag_int = (pitch_index + 137) >> 2;
252 *lag_frac = pitch_index - (*lag_int << 2) + 136;
253 }
else if (pitch_index < 440) {
254 *lag_int = (pitch_index + 257 - 376) >> 1;
255 *lag_frac = (pitch_index - (*lag_int << 1) + 256 - 376) << 1;
258 *lag_int = pitch_index - 280;
262 *base_lag_int = av_clip(*lag_int - 8 - (*lag_frac < 0),
268 *lag_int = (pitch_index + 1) >> 2;
269 *lag_frac = pitch_index - (*lag_int << 2);
270 *lag_int += *base_lag_int;
280 uint8_t *base_lag_int,
int subframe,
enum Mode mode)
282 if (subframe == 0 || (subframe == 2 && mode !=
MODE_6k60)) {
283 if (pitch_index < 116) {
284 *lag_int = (pitch_index + 69) >> 1;
285 *lag_frac = (pitch_index - (*lag_int << 1) + 68) << 1;
287 *lag_int = pitch_index - 24;
291 *base_lag_int = av_clip(*lag_int - 8 - (*lag_frac < 0),
294 *lag_int = (pitch_index + 1) >> 1;
295 *lag_frac = (pitch_index - (*lag_int << 1)) << 1;
296 *lag_int += *base_lag_int;
312 int pitch_lag_int, pitch_lag_frac;
325 pitch_lag_int += pitch_lag_frac > 0;
331 pitch_lag_frac + (pitch_lag_frac > 0 ? 0 : 4),
336 if (amr_subframe->
ltp) {
340 ctx->
pitch_vector[i] = 0.18 * exc[i - 1] + 0.64 * exc[i] +
342 memcpy(exc, ctx->
pitch_vector, AMRWB_SFR_SIZE *
sizeof(
float));
347 #define BIT_STR(x,lsb,len) (((x) >> (lsb)) & ((1 << (len)) - 1))
350 #define BIT_POS(x, p) (((x) >> (p)) & 1)
369 out[0] =
BIT_POS(code, m) ? -pos : pos;
377 out[0] =
BIT_POS(code, 2*m) ? -pos0 : pos0;
378 out[1] =
BIT_POS(code, 2*m) ? -pos1 : pos1;
379 out[1] = pos0 > pos1 ? -out[1] : out[1];
384 int half_2p =
BIT_POS(code, 2*m - 1) << (m - 1);
387 m - 1, off + half_2p);
393 int half_4p, subhalf_2p;
394 int b_offset = 1 << (m - 1);
396 switch (
BIT_STR(code, 4*m - 2, 2)) {
398 half_4p =
BIT_POS(code, 4*m - 3) << (m - 1);
399 subhalf_2p =
BIT_POS(code, 2*m - 3) << (m - 2);
402 m - 2, off + half_4p + subhalf_2p);
404 m - 1, off + half_4p);
410 m - 1, off + b_offset);
416 m - 1, off + b_offset);
422 m - 1, off + b_offset);
429 int half_3p =
BIT_POS(code, 5*m - 1) << (m - 1);
432 m - 1, off + half_3p);
439 int b_offset = 1 << (m - 1);
441 int half_more =
BIT_POS(code, 6*m - 5) << (m - 1);
442 int half_other = b_offset - half_more;
444 switch (
BIT_STR(code, 6*m - 4, 2)) {
447 m - 1, off + half_more);
449 m - 1, off + half_more);
453 m - 1, off + half_other);
455 m - 1, off + half_more);
459 m - 1, off + half_other);
461 m - 1, off + half_more);
467 m - 1, off + b_offset);
482 const uint16_t *pulse_lo,
const enum Mode mode)
487 int spacing = (mode ==
MODE_6k60) ? 2 : 4;
492 for (i = 0; i < 2; i++)
496 for (i = 0; i < 4; i++)
500 for (i = 0; i < 4; i++)
504 for (i = 0; i < 2; i++)
506 for (i = 2; i < 4; i++)
510 for (i = 0; i < 4; i++)
514 for (i = 0; i < 4; i++)
516 ((int) pulse_hi[i] << 14), 4, 1);
519 for (i = 0; i < 2; i++)
521 ((int) pulse_hi[i] << 10), 4, 1);
522 for (i = 2; i < 4; i++)
524 ((int) pulse_hi[i] << 14), 4, 1);
528 for (i = 0; i < 4; i++)
530 ((int) pulse_hi[i] << 11), 4, 1);
536 for (i = 0; i < 4; i++)
538 int pos = (
FFABS(sig_pos[i][j]) - 1) * spacing + i;
540 fixed_vector[pos] += sig_pos[i][j] < 0 ? -1.0 : 1.0;
553 float *fixed_gain_factor,
float *pitch_gain)
558 *pitch_gain = gains[0] * (1.0f / (1 << 14));
559 *fixed_gain_factor = gains[1] * (1.0f / (1 << 11));
576 fixed_vector[i] -= fixed_vector[i - 1] * ctx->
tilt_coef;
580 fixed_vector[i] += fixed_vector[i - ctx->
pitch_lag_int] * 0.85;
592 float *f_vector,
float f_gain)
599 return (p_ener - f_ener) / (p_ener + f_ener);
613 float *fixed_vector,
float *buf)
629 if (ir_filter_nr < 2)
634 for (i = 0; i < 6; i++)
650 if (ir_filter_nr < 2) {
682 acc += (isf[i] - isf_past[i]) * (isf[i] - isf_past[i]);
686 return FFMAX(0.0, 1.25 - acc * 0.8 * 512);
701 float voice_fac,
float stab_fac)
703 float sm_fac = 0.5 * (1 - voice_fac) * stab_fac;
709 if (fixed_gain < *prev_tr_gain) {
710 g0 =
FFMIN(*prev_tr_gain, fixed_gain + fixed_gain *
711 (6226 * (1.0
f / (1 << 15))));
713 g0 =
FFMAX(*prev_tr_gain, fixed_gain *
714 (27536 * (1.0
f / (1 << 15))));
718 return sm_fac * g0 + (1 - sm_fac) * fixed_gain;
730 float cpe = 0.125 * (1 + voice_fac);
731 float last = fixed_vector[0];
733 fixed_vector[0] -= cpe * fixed_vector[1];
736 float cur = fixed_vector[i];
738 fixed_vector[i] -= cpe * (last + fixed_vector[i + 1]);
742 fixed_vector[AMRWB_SFR_SIZE - 1] -= cpe * last;
756 float fixed_gain,
const float *fixed_vector,
776 energy, AMRWB_SFR_SIZE);
792 static void de_emphasis(
float *out,
float *in,
float m,
float mem[1])
796 out[0] = in[0] + m * mem[0];
799 out[i] = in[i] + out[i - 1] * m;
801 mem[0] = out[AMRWB_SFR_SIZE - 1];
816 int int_part = 0, frac_part;
819 for (j = 0; j < o_size / 5; j++) {
820 out[i] = in[int_part];
824 for (k = 1; k < 5; k++) {
844 uint16_t hb_idx, uint8_t vad)
856 return av_clipf((1.0 -
FFMAX(0.0, tilt)) * (1.25 - 0.25 * wsp), 0.1, 1.0);
869 const float *synth_exc,
float hb_gain)
879 energy * hb_gain * hb_gain,
891 for (i = 7; i <
LP_ORDER - 2; i++) {
892 float prod = (diff_isf[i] - mean) * (diff_isf[i - lag] - mean);
907 float diff_isf[
LP_ORDER - 2], diff_mean;
908 float *diff_hi = diff_isf -
LP_ORDER + 1;
913 isf[LP_ORDER_16k - 1] = isf[LP_ORDER - 1];
916 for (i = 0; i < LP_ORDER - 2; i++)
917 diff_isf[i] = isf[i + 1] - isf[i];
920 for (i = 2; i < LP_ORDER - 2; i++)
921 diff_mean += diff_isf[i] * (1.0
f / (LP_ORDER - 4));
925 for (i = 0; i < 3; i++) {
928 if (corr_lag[i] > corr_lag[i_max_corr])
933 for (i = LP_ORDER - 1; i < LP_ORDER_16k - 1; i++)
934 isf[i] = isf[i - 1] + isf[i - 1 - i_max_corr]
935 - isf[i - 2 - i_max_corr];
938 est = 7965 + (isf[2] - isf[3] - isf[4]) / 6.0;
939 scale = 0.5 * (
FFMIN(est, 7600) - isf[LP_ORDER - 2]) /
940 (isf[LP_ORDER_16k - 2] - isf[LP_ORDER - 2]);
942 for (i = LP_ORDER - 1; i < LP_ORDER_16k - 1; i++)
943 diff_hi[i] = scale * (isf[i] - isf[i - 1]);
946 for (i = LP_ORDER; i < LP_ORDER_16k - 1; i++)
947 if (diff_hi[i] + diff_hi[i - 1] < 5.0) {
948 if (diff_hi[i] > diff_hi[i - 1]) {
949 diff_hi[i - 1] = 5.0 - diff_hi[i];
951 diff_hi[i] = 5.0 - diff_hi[i - 1];
954 for (i = LP_ORDER - 1; i < LP_ORDER_16k - 1; i++)
955 isf[i] = isf[i - 1] + diff_hi[i] * (1.0
f / (1 << 15));
958 for (i = 0; i < LP_ORDER_16k - 1; i++)
976 for (i = 0; i <
size; i++) {
977 out[i] = lpc[i] * fac;
994 const float *exc,
const float *isf,
const float *isf_past)
1038 memcpy(data, mem, HB_FIR_SIZE *
sizeof(
float));
1044 out[i] += data[i + j] * fir_coef[j];
1047 memcpy(mem, data + AMRWB_SFR_SIZE_16k, HB_FIR_SIZE *
sizeof(
float));
1070 int *got_frame_ptr,
AVPacket *avpkt)
1074 const uint8_t *buf = avpkt->
data;
1075 int buf_size = avpkt->
size;
1076 int expected_fr_size, header_size;
1079 float fixed_gain_factor;
1080 float *synth_fixed_vector;
1081 float synth_fixed_gain;
1082 float voice_fac, stab_fac;
1105 if (buf_size < expected_fr_size) {
1107 "Frame too small (%d bytes). Truncated file?\n", buf_size);
1145 for (sub = 0; sub < 4; sub++)
1148 for (sub = 0; sub < 4; sub++) {
1173 ctx->
tilt_coef = voice_fac * 0.25 + 0.25;
1184 voice_fac, stab_fac);
1203 AMRWB_SFR_SIZE_16k);
1228 sub_buf[i] = (sub_buf[i] + hb_samples[i]) * (1.0f / (1 << 15));
1241 return expected_fr_size;