41 #define AMR_USE_16BIT_TABLES
69 float prediction_error[4];
83 float hpf_31_mem[2], hpf_400_mem[2];
121 for (i = 0; i < 4; i++)
164 for (i = 0; i < 9; i++)
165 isf_q[i] =
dico1_isf[ind[0]][i] * (1.0f / (1 << 15));
167 for (i = 0; i < 7; i++)
168 isf_q[i + 9] =
dico2_isf[ind[1]][i] * (1.0f / (1 << 15));
170 for (i = 0; i < 5; i++)
173 for (i = 0; i < 4; i++)
176 for (i = 0; i < 7; i++)
191 for (i = 0; i < 9; i++)
192 isf_q[i] =
dico1_isf[ind[0]][i] * (1.0f / (1 << 15));
194 for (i = 0; i < 7; i++)
195 isf_q[i + 9] =
dico2_isf[ind[1]][i] * (1.0f / (1 << 15));
197 for (i = 0; i < 3; i++)
198 isf_q[i] +=
dico21_isf[ind[2]][i] * (1.0f / (1 << 15));
200 for (i = 0; i < 3; i++)
201 isf_q[i + 3] +=
dico22_isf[ind[3]][i] * (1.0f / (1 << 15));
203 for (i = 0; i < 3; i++)
204 isf_q[i + 6] +=
dico23_isf[ind[4]][i] * (1.0f / (1 << 15));
206 for (i = 0; i < 3; i++)
207 isf_q[i + 9] +=
dico24_isf[ind[5]][i] * (1.0f / (1 << 15));
209 for (i = 0; i < 4; i++)
210 isf_q[i + 12] +=
dico25_isf[ind[6]][i] * (1.0f / (1 << 15));
228 isf_q[i] +=
isf_mean[i] * (1.0f / (1 << 15));
245 for (k = 0; k < 3; k++) {
248 isp_q[k][i] = (1.0 - c) * isp4_past[i] + c * isp_q[3][i];
264 uint8_t *base_lag_int,
int subframe)
266 if (subframe == 0 || subframe == 2) {
267 if (pitch_index < 376) {
268 *lag_int = (pitch_index + 137) >> 2;
269 *lag_frac = pitch_index - (*lag_int << 2) + 136;
270 }
else if (pitch_index < 440) {
271 *lag_int = (pitch_index + 257 - 376) >> 1;
272 *lag_frac = (pitch_index - (*lag_int << 1) + 256 - 376) << 1;
275 *lag_int = pitch_index - 280;
279 *base_lag_int = av_clip(*lag_int - 8 - (*lag_frac < 0),
285 *lag_int = (pitch_index + 1) >> 2;
286 *lag_frac = pitch_index - (*lag_int << 2);
287 *lag_int += *base_lag_int;
299 if (subframe == 0 || (subframe == 2 && mode !=
MODE_6k60)) {
300 if (pitch_index < 116) {
301 *lag_int = (pitch_index + 69) >> 1;
302 *lag_frac = (pitch_index - (*lag_int << 1) + 68) << 1;
304 *lag_int = pitch_index - 24;
308 *base_lag_int = av_clip(*lag_int - 8 - (*lag_frac < 0),
311 *lag_int = (pitch_index + 1) >> 1;
312 *lag_frac = (pitch_index - (*lag_int << 1)) << 1;
313 *lag_int += *base_lag_int;
329 int pitch_lag_int, pitch_lag_frac;
342 pitch_lag_int += pitch_lag_frac > 0;
347 exc + 1 - pitch_lag_int,
349 pitch_lag_frac + (pitch_lag_frac > 0 ? 0 : 4),
354 if (amr_subframe->
ltp) {
358 ctx->
pitch_vector[i] = 0.18 * exc[i - 1] + 0.64 * exc[i] +
360 memcpy(exc, ctx->
pitch_vector, AMRWB_SFR_SIZE *
sizeof(
float));
365 #define BIT_STR(x,lsb,len) (((x) >> (lsb)) & ((1 << (len)) - 1))
368 #define BIT_POS(x, p) (((x) >> (p)) & 1)
387 out[0] =
BIT_POS(code, m) ? -pos : pos;
395 out[0] =
BIT_POS(code, 2*m) ? -pos0 : pos0;
396 out[1] =
BIT_POS(code, 2*m) ? -pos1 : pos1;
397 out[1] = pos0 > pos1 ? -out[1] : out[1];
402 int half_2p =
BIT_POS(code, 2*m - 1) << (m - 1);
405 m - 1, off + half_2p);
411 int half_4p, subhalf_2p;
412 int b_offset = 1 << (m - 1);
414 switch (
BIT_STR(code, 4*m - 2, 2)) {
416 half_4p =
BIT_POS(code, 4*m - 3) << (m - 1);
417 subhalf_2p =
BIT_POS(code, 2*m - 3) << (m - 2);
420 m - 2, off + half_4p + subhalf_2p);
422 m - 1, off + half_4p);
428 m - 1, off + b_offset);
434 m - 1, off + b_offset);
440 m - 1, off + b_offset);
447 int half_3p =
BIT_POS(code, 5*m - 1) << (m - 1);
450 m - 1, off + half_3p);
457 int b_offset = 1 << (m - 1);
459 int half_more =
BIT_POS(code, 6*m - 5) << (m - 1);
460 int half_other = b_offset - half_more;
462 switch (
BIT_STR(code, 6*m - 4, 2)) {
465 m - 1, off + half_more);
467 m - 1, off + half_more);
471 m - 1, off + half_other);
473 m - 1, off + half_more);
477 m - 1, off + half_other);
479 m - 1, off + half_more);
485 m - 1, off + b_offset);
500 const uint16_t *pulse_lo,
const enum Mode mode)
505 int spacing = (mode ==
MODE_6k60) ? 2 : 4;
510 for (i = 0; i < 2; i++)
514 for (i = 0; i < 4; i++)
518 for (i = 0; i < 4; i++)
522 for (i = 0; i < 2; i++)
524 for (i = 2; i < 4; i++)
528 for (i = 0; i < 4; i++)
532 for (i = 0; i < 4; i++)
534 ((int) pulse_hi[i] << 14), 4, 1);
537 for (i = 0; i < 2; i++)
539 ((int) pulse_hi[i] << 10), 4, 1);
540 for (i = 2; i < 4; i++)
542 ((int) pulse_hi[i] << 14), 4, 1);
546 for (i = 0; i < 4; i++)
548 ((int) pulse_hi[i] << 11), 4, 1);
554 for (i = 0; i < 4; i++)
556 int pos = (
FFABS(sig_pos[i][j]) - 1) * spacing + i;
558 fixed_vector[pos] += sig_pos[i][j] < 0 ? -1.0 : 1.0;
571 float *fixed_gain_factor,
float *pitch_gain)
576 *pitch_gain = gains[0] * (1.0f / (1 << 14));
577 *fixed_gain_factor = gains[1] * (1.0f / (1 << 11));
594 fixed_vector[i] -= fixed_vector[i - 1] * ctx->
tilt_coef;
598 fixed_vector[i] += fixed_vector[i - ctx->
pitch_lag_int] * 0.85;
611 float *f_vector,
float f_gain,
614 double p_ener = (double) ctx->
dot_productf(p_vector, p_vector,
617 double f_ener = (double) ctx->
dot_productf(f_vector, f_vector,
621 return (p_ener - f_ener) / (p_ener + f_ener);
635 float *fixed_vector,
float *buf)
651 if (ir_filter_nr < 2)
656 for (i = 0; i < 6; i++)
672 if (ir_filter_nr < 2) {
704 acc += (isf[i] - isf_past[i]) * (isf[i] - isf_past[i]);
708 return FFMAX(0.0, 1.25 - acc * 0.8 * 512);
723 float voice_fac,
float stab_fac)
725 float sm_fac = 0.5 * (1 - voice_fac) * stab_fac;
731 if (fixed_gain < *prev_tr_gain) {
732 g0 =
FFMIN(*prev_tr_gain, fixed_gain + fixed_gain *
733 (6226 * (1.0f / (1 << 15))));
735 g0 =
FFMAX(*prev_tr_gain, fixed_gain *
736 (27536 * (1.0f / (1 << 15))));
740 return sm_fac * g0 + (1 - sm_fac) * fixed_gain;
752 float cpe = 0.125 * (1 + voice_fac);
753 float last = fixed_vector[0];
755 fixed_vector[0] -= cpe * fixed_vector[1];
758 float cur = fixed_vector[i];
760 fixed_vector[i] -= cpe * (last + fixed_vector[i + 1]);
764 fixed_vector[AMRWB_SFR_SIZE - 1] -= cpe * last;
778 float fixed_gain,
const float *fixed_vector,
798 energy, AMRWB_SFR_SIZE);
818 out[0] = in[0] + m * mem[0];
821 out[i] = in[i] + out[i - 1] * m;
823 mem[0] = out[AMRWB_SFR_SIZE - 1];
839 int int_part = 0, frac_part;
842 for (j = 0; j < o_size / 5; j++) {
843 out[i] = in[int_part];
847 for (k = 1; k < 5; k++) {
880 return av_clipf((1.0 -
FFMAX(0.0, tilt)) * (1.25 - 0.25 * wsp), 0.1, 1.0);
893 const float *synth_exc,
float hb_gain)
903 energy * hb_gain * hb_gain,
915 for (i = 7; i <
LP_ORDER - 2; i++) {
916 float prod = (diff_isf[i] - mean) * (diff_isf[i - lag] - mean);
931 float diff_isf[
LP_ORDER - 2], diff_mean;
934 int i, j, i_max_corr;
936 isf[LP_ORDER_16k - 1] = isf[
LP_ORDER - 1];
940 diff_isf[i] = isf[i + 1] - isf[i];
943 for (i = 2; i < LP_ORDER - 2; i++)
944 diff_mean += diff_isf[i] * (1.0f / (LP_ORDER - 4));
948 for (i = 0; i < 3; i++) {
951 if (corr_lag[i] > corr_lag[i_max_corr])
956 for (i = LP_ORDER - 1; i < LP_ORDER_16k - 1; i++)
957 isf[i] = isf[i - 1] + isf[i - 1 - i_max_corr]
958 - isf[i - 2 - i_max_corr];
961 est = 7965 + (isf[2] - isf[3] - isf[4]) / 6.0;
962 scale = 0.5 * (
FFMIN(est, 7600) - isf[LP_ORDER - 2]) /
963 (isf[LP_ORDER_16k - 2] - isf[LP_ORDER - 2]);
965 for (i = LP_ORDER - 1, j = 0; i < LP_ORDER_16k - 1; i++, j++)
966 diff_isf[j] = scale * (isf[i] - isf[i - 1]);
969 for (i = 1; i < LP_ORDER_16k -
LP_ORDER; i++)
970 if (diff_isf[i] + diff_isf[i - 1] < 5.0) {
971 if (diff_isf[i] > diff_isf[i - 1]) {
972 diff_isf[i - 1] = 5.0 - diff_isf[i];
974 diff_isf[i] = 5.0 - diff_isf[i - 1];
977 for (i = LP_ORDER - 1, j = 0; i < LP_ORDER_16k - 1; i++, j++)
978 isf[i] = isf[i - 1] + diff_isf[j] * (1.0f / (1 << 15));
981 for (i = 0; i < LP_ORDER_16k - 1; i++)
999 for (i = 0; i <
size; i++) {
1000 out[i] = lpc[i] * fac;
1017 const float *exc,
const float *isf,
const float *isf_past)
1056 #ifndef hb_fir_filter
1063 memcpy(data, mem, HB_FIR_SIZE *
sizeof(
float));
1069 out[i] += data[i + j] * fir_coef[j];
1072 memcpy(mem, data + AMRWB_SFR_SIZE_16k, HB_FIR_SIZE *
sizeof(
float));
1096 int *got_frame_ptr,
AVPacket *avpkt)
1101 int buf_size = avpkt->
size;
1102 int expected_fr_size, header_size;
1105 float fixed_gain_factor;
1106 float *synth_fixed_vector;
1107 float synth_fixed_gain;
1108 float voice_fac, stab_fac;
1131 if (buf_size < expected_fr_size) {
1133 "Frame too small (%d bytes). Truncated file?\n", buf_size);
1171 for (sub = 0; sub < 4; sub++)
1174 for (sub = 0; sub < 4; sub++) {
1202 ctx->
tilt_coef = voice_fac * 0.25 + 0.25;
1213 voice_fac, stab_fac);
1257 sub_buf[i] = (sub_buf[i] + hb_samples[i]) * (1.0f / (1 << 15));
1270 return expected_fr_size;