65 float predictor_lspf[10];
66 float pitch_synthesis_filter_mem[303];
67 float pitch_pre_filter_mem[303];
68 float rnd_fir_filter_mem[180];
69 float formant_mem[170];
79 float postfilter_synth_mem[10];
96 for (i = 0; i < 10; i++)
119 float tmp_lspf, smooth, erasure_coeff;
120 const float *predictors;
130 for (i = 0; i < 10; i++) {
146 for (i = 0; i < 10; i++) {
148 lspf[i] = (i + 1) * (1 - erasure_coeff) / 11 +
149 erasure_coeff * predictors[i];
156 for (i = 1; i < 10; i++)
160 for (i = 9; i > 0; i--)
169 for (i = 0; i < 5; i++) {
176 if (lspf[9] <= .70 || lspf[9] >= .97)
178 for (i = 3; i < 10; i++)
179 if (fabs(lspf[i] - lspf[i - 2]) < .08)
182 if (lspf[9] <= .66 || lspf[9] >= .985)
184 for (i = 4; i < 10; i++)
185 if (fabs(lspf[i] - lspf[i - 4]) < .0931)
202 int i, subframes_count, g1[16];
207 case RATE_FULL: subframes_count = 16;
break;
208 case RATE_HALF: subframes_count = 4;
break;
209 default: subframes_count = 5;
211 for (i = 0; i < subframes_count; i++) {
214 g1[i] += av_clip((g1[i - 1] + g1[i - 2] + g1[i - 3]) / 3 - 6, 0, 32);
232 gain[6] = 0.4 * gain[3] + 0.6 * gain[4];
234 gain[4] = 0.8 * gain[2] + 0.2 * gain[3];
235 gain[3] = 0.2 * gain[1] + 0.8 * gain[2];
237 gain[1] = 0.6 * gain[0] + 0.4 * gain[1];
250 case 2 : g1[0] -= 1;
break;
251 case 3 : g1[0] -= 2;
break;
260 for (i = 1; i <= subframes_count; i++)
280 int i, diff, prev_diff = 0;
282 for (i = 1; i < 5; i++) {
283 diff = cbgain[i] - cbgain[i-1];
284 if (
FFABS(diff) > 10)
286 else if (
FFABS(diff - prev_diff) > 12)
318 uint16_t cbseed, cindex;
319 float *rnd, tmp_gain, fir_filter_value;
323 for (i = 0; i < 16; i++) {
326 for (j = 0; j < 10; j++)
331 for (i = 0; i < 4; i++) {
334 for (j = 0; j < 40; j++)
339 cbseed = (0x0003 & q->
frame.
lspv[4]) << 14 |
345 for (i = 0; i < 8; i++) {
347 for (k = 0; k < 20; k++) {
348 cbseed = 521 * cbseed + 259;
349 *rnd = (int16_t) cbseed;
352 fir_filter_value = 0.0;
353 for (j = 0; j < 10; j++)
355 (rnd[-j] + rnd[-20+j]);
358 *cdn_vector++ = tmp_gain * fir_filter_value;
367 for (i = 0; i < 8; i++) {
369 for (j = 0; j < 20; j++) {
370 cbseed = 521 * cbseed + 259;
371 *cdn_vector++ = tmp_gain * (int16_t) cbseed;
377 for (i = 0; i < 4; i++) {
379 for (j = 0; j < 40; j++)
384 memset(cdn_vector, 0, 160 *
sizeof(
float));
402 for (i = 0; i < 160; i += 40)
427 const float gain[4],
const uint8_t *lag,
428 const uint8_t pfrac[4])
431 float *v_lag, *v_out;
434 v_out = memory + 143;
436 for (i = 0; i < 4; i++) {
438 v_lag = memory + 143 + 40 * i - lag[i];
439 for (v_len = v_in + 40; v_in < v_len; v_in++) {
441 for (j = 0, *v_out = 0.; j < 4; j++)
446 *v_out = *v_in + gain[i] * *v_out;
452 memcpy(v_out, v_in, 40 *
sizeof(
float));
458 memmove(memory, memory + 160, 143 *
sizeof(
float));
472 const float *v_synthesis_filtered, *v_pre_filtered;
479 for (i = 0; i < 4; i++) {
485 float max_pitch_gain;
491 max_pitch_gain = 0.0;
494 max_pitch_gain = 1.0;
496 for (i = 0; i < 4; i++)
508 for (i = 0; i < 4; i++)
512 v_synthesis_filtered,
537 static void lspf2lpc(
const float *lspf,
float *lpc)
543 for (i = 0; i < 10; i++)
544 lsp[i] = cos(
M_PI * lspf[i]);
548 for (i = 0; i < 10; i++) {
549 lpc[i] *= bandwidth_expansion_coeff;
566 float *lpc,
const int subframe_num)
568 float interpolated_lspf[10];
572 weight = 0.25 * (subframe_num + 1);
580 weight, 1.0 - weight, 10);
621 if (bitrate > **buf) {
625 "Claimed bitrate and buffer size mismatch.\n");
629 }
else if (bitrate < **buf) {
631 "Buffer is too small for the claimed bitrate.\n");
637 "Bitrate byte is missing, guessing the bitrate from packet size.\n");
657 static const float pow_0_775[10] = {
658 0.775000, 0.600625, 0.465484, 0.360750, 0.279582,
659 0.216676, 0.167924, 0.130141, 0.100859, 0.078166
661 0.625000, 0.390625, 0.244141, 0.152588, 0.095367,
662 0.059605, 0.037253, 0.023283, 0.014552, 0.009095
664 float lpc_s[10], lpc_p[10], pole_out[170], zero_out[160];
667 for (n = 0; n < 10; n++) {
668 lpc_s[n] = lpc[n] * pow_0_625[n];
669 lpc_p[n] = lpc[n] * pow_0_775[n];
687 int *got_frame_ptr,
AVPacket *avpkt)
689 const uint8_t *buf = avpkt->
data;
690 int buf_size = avpkt->
size;
694 float quantized_lspf[10], lpc[10];
721 uint8_t *unpacked_data = (uint8_t *)&q->
frame;
727 for (; bitmaps < bitmaps_end; bitmaps++)
742 for (i = 0; i < 4; i++) {
773 for (i = 0; i < 4; i++) {