49 #define IMC_BLOCK_SIZE 64
50 #define IMC_FRAME_ID 0x21
105 float weights1[31], weights2[31];
110 #define VLC_TABLES_SIZE 9512
113 0, 640, 1156, 1732, 2308, 2852, 3396, 3924,
121 return 3.5 * atan((freq / 7500.0) * (freq / 7500.0)) + 13.0 * atan(freq * 0.00076);
126 double freqmin[32], freqmid[32], freqmax[32];
127 double scale = sampling_rate / (256.0 * 2.0 * 2.0);
128 double nyquist_freq = sampling_rate * 0.5;
129 double freq, bark, prev_bark = 0, tf, tb;
132 for (i = 0; i < 32; i++) {
137 tb = bark - prev_bark;
138 q->
weights1[i - 1] = pow(10.0, -1.0 * tb);
139 q->
weights2[i - 1] = pow(10.0, -2.7 * tb);
146 while (tf < nyquist_freq) {
158 if (tb <= bark - 0.5)
164 for (i = 0; i < 32; i++) {
166 for (j = 31; j > 0 && freq <= freqmid[j]; j--);
170 for (j = 0; j < 32 && freq >= freqmid[j]; j++);
189 for (j = 0; j < avctx->
channels; j++) {
192 for (i = 0; i <
BANDS; i++)
201 for (i = 0; i <
COEFFS; i++)
203 for (i = 0; i < COEFFS / 2; i++) {
204 q->
post_cos[i] = (1.0f / 32768) * cos(i / 256.0 * M_PI);
205 q->
post_sin[i] = (1.0f / 32768) * sin(i / 256.0 * M_PI);
207 r1 = sin((i * 4.0 + 1.0) / 1024.0 * M_PI);
208 r2 = cos((i * 4.0 + 1.0) / 1024.0 * M_PI);
212 q->
pre_coef2[i] = -(r1 - r2) * sqrt(2.0);
214 q->
pre_coef1[i] = -(r1 + r2) * sqrt(2.0);
221 for (i = 0; i < 30; i++)
225 for (i = 0; i < 4 ; i++) {
226 for (j = 0; j < 4; j++) {
228 huffman_vlc[i][j].
table_allocated = vlc_offsets[i * 4 + j + 1] - vlc_offsets[i * 4 + j];
261 float *flcoeffs2,
int *bandWidthT,
262 float *flcoeffs3,
float *flcoeffs5)
267 float snr_limit = 1.e-30;
271 for (i = 0; i <
BANDS; i++) {
272 flcoeffs5[i] = workT2[i] = 0.0;
274 workT1[i] = flcoeffs1[i] * flcoeffs1[i];
275 flcoeffs3[i] = 2.0 * flcoeffs2[i];
278 flcoeffs3[i] = -30000.0;
280 workT3[i] = bandWidthT[i] * workT1[i] * 0.01;
281 if (workT3[i] <= snr_limit)
285 for (i = 0; i <
BANDS; i++) {
286 for (cnt2 = i; cnt2 < q->
cyclTab[i]; cnt2++)
287 flcoeffs5[cnt2] = flcoeffs5[cnt2] + workT3[i];
288 workT2[cnt2 - 1] = workT2[cnt2 - 1] + workT3[i];
291 for (i = 1; i <
BANDS; i++) {
292 accum = (workT2[i - 1] + accum) * q->
weights1[i - 1];
293 flcoeffs5[i] += accum;
296 for (i = 0; i <
BANDS; i++)
299 for (i = 0; i <
BANDS; i++) {
300 for (cnt2 = i - 1; cnt2 > q->
cyclTab2[i]; cnt2--)
301 flcoeffs5[cnt2] += workT3[i];
302 workT2[cnt2+1] += workT3[i];
307 for (i = BANDS-2; i >= 0; i--) {
308 accum = (workT2[i+1] + accum) * q->
weights2[i];
309 flcoeffs5[i] += accum;
324 s = stream_format_code >> 1;
325 hufftab[0] = &huffman_vlc[s][0];
326 hufftab[1] = &huffman_vlc[s][1];
327 hufftab[2] = &huffman_vlc[s][2];
328 hufftab[3] = &huffman_vlc[s][3];
331 if (stream_format_code & 4)
335 for (i = start; i <
BANDS; i++) {
337 hufftab[cb_sel[i]]->
bits, 2);
338 if (levlCoeffs[i] == 17)
344 float *flcoeffs1,
float *flcoeffs2)
350 flcoeffs1[0] = 20000.0 / pow (2, levlCoeffBuf[0] * 0.18945);
351 flcoeffs2[0] =
log2f(flcoeffs1[0]);
355 for (i = 1; i <
BANDS; i++) {
356 level = levlCoeffBuf[i];
363 else if (level <= 24)
369 tmp2 += 0.83048 *
level;
378 float *old_floor,
float *flcoeffs1,
386 for (i = 0; i <
BANDS; i++) {
388 if (levlCoeffBuf[i] < 16) {
389 flcoeffs1[i] =
imc_exp_tab2[levlCoeffBuf[i]] * old_floor[i];
390 flcoeffs2[i] = (levlCoeffBuf[i] - 7) * 0.83048 + flcoeffs2[i];
392 flcoeffs1[i] = old_floor[i];
401 int stream_format_code,
int freebits,
int flag)
404 const float limit = -1.e20;
413 float lowest = 1.e10;
419 for (i = 0; i <
BANDS; i++)
422 for (i = 0; i < BANDS - 1; i++)
426 highest = highest * 0.25;
428 for (i = 0; i <
BANDS; i++) {
445 if (stream_format_code & 0x2) {
452 for (i = (stream_format_code & 0x2) ? 4 : 0; i < BANDS - 1; i++) {
461 summa = (summa * 0.5 - freebits) / iacc;
464 for (i = 0; i < BANDS / 2; i++) {
465 rres = summer - freebits;
466 if ((rres >= -8) && (rres <= 8))
472 for (j = (stream_format_code & 0x2) ? 4 : 0; j <
BANDS; j++) {
473 cwlen = av_clipf(((chctx->
flcoeffs4[j] * 0.5) - summa + 0.5), 0, 6);
484 if (freebits < summer)
491 summa = (float)(summer - freebits) / ((t1 + 1) * iacc) + summa;
494 for (i = (stream_format_code & 0x2) ? 4 : 0; i <
BANDS; i++) {
499 if (freebits > summer) {
500 for (i = 0; i <
BANDS; i++) {
501 workT[i] = (chctx->
bitsBandT[i] == 6) ? -1.e20
508 if (highest <= -1.e20)
514 for (i = 0; i <
BANDS; i++) {
515 if (workT[i] > highest) {
521 if (highest > -1.e20) {
522 workT[found_indx] -= 2.0;
524 workT[found_indx] = -1.e20;
526 for (j =
band_tab[found_indx]; j < band_tab[found_indx + 1] && (freebits > summer); j++) {
531 }
while (freebits > summer);
533 if (freebits < summer) {
534 for (i = 0; i <
BANDS; i++) {
538 if (stream_format_code & 0x2) {
544 while (freebits < summer) {
547 for (i = 0; i <
BANDS; i++) {
548 if (workT[i] < lowest) {
555 workT[low_indx] = lowest + 2.0;
558 workT[low_indx] = 1.e20;
560 for (j =
band_tab[low_indx]; j <
band_tab[low_indx+1] && (freebits < summer); j++) {
577 for (i = 0; i <
BANDS; i++) {
584 for (j = band_tab[i]; j < band_tab[i + 1]; j++) {
615 if (j < band_tab[i + 1]) {
636 for (i = 0; i <
BANDS; i++) {
637 workT[i] = (chctx->
bitsBandT[i] == 6) ? -1.e20
641 while (corrected < summer) {
642 if (highest <= -1.e20)
647 for (i = 0; i <
BANDS; i++) {
648 if (workT[i] > highest) {
654 if (highest > -1.e20) {
655 workT[found_indx] -= 2.0;
656 if (++(chctx->
bitsBandT[found_indx]) == 6)
657 workT[found_indx] = -1.e20;
659 for (j =
band_tab[found_indx]; j <
band_tab[found_indx+1] && (corrected < summer); j++) {
677 for (i = 0; i <
COEFFS / 2; i++) {
689 for (i = 0; i <
COEFFS / 2; i++) {
703 int stream_format_code)
706 int middle_value, cw_len, max_size;
707 const float *quantizer;
709 for (i = 0; i <
BANDS; i++) {
714 if (cw_len <= 0 || chctx->skipFlags[j])
717 max_size = 1 << cw_len;
718 middle_value = max_size >> 1;
744 int i, j, cw_len, cw;
746 for (i = 0; i <
BANDS; i++) {
755 av_dlog(
NULL,
"Band %i coeff %i cw_len %i\n", i, j, cw_len);
771 int stream_format_code;
772 int imc_hdr, i, j, ret;
775 int counter, bitscount;
781 if (imc_hdr & 0x18) {
788 if (stream_format_code & 1) {
794 if (stream_format_code & 0x04)
798 for (i = 0; i <
BANDS; i++)
808 if (stream_format_code & 0x4)
818 for (i = 0; i <
BANDS; i++) {
826 for (i = 0; i < BANDS - 1; i++) {
835 if (stream_format_code & 0x2) {
842 for (i = 1; i < 4; i++) {
853 if (!(stream_format_code & 0x2))
865 for (i = 0; i <
BANDS; i++) {
871 if ((((band_tab[i + 1] - band_tab[i]) * 1.5) > chctx->
sumLenArr[i]) && (chctx->
sumLenArr[i] > 0))
877 for (i = 0; i <
BANDS; i++) {
889 for (i = 0; i <
BANDS; i++) {
903 for (i = 0; i <
BANDS; i++) {
933 int *got_frame_ptr,
AVPacket *avpkt)
936 int buf_size = avpkt->
size;
943 if (buf_size < IMC_BLOCK_SIZE * avctx->channels) {
955 for (i = 0; i < avctx->
channels; i++) {