50 #define IMC_BLOCK_SIZE 64
51 #define IMC_FRAME_ID 0x21
106 float weights1[31], weights2[31];
111 #define VLC_TABLES_SIZE 9512
114 0, 640, 1156, 1732, 2308, 2852, 3396, 3924,
122 return 3.5 * atan((freq / 7500.0) * (freq / 7500.0)) + 13.0 * atan(freq * 0.00076);
127 double freqmin[32], freqmid[32], freqmax[32];
128 double scale = sampling_rate / (256.0 * 2.0 * 2.0);
129 double nyquist_freq = sampling_rate * 0.5;
130 double freq, bark, prev_bark = 0, tf, tb;
133 for (i = 0; i < 32; i++) {
138 tb = bark - prev_bark;
139 q->
weights1[i - 1] = pow(10.0, -1.0 * tb);
140 q->
weights2[i - 1] = pow(10.0, -2.7 * tb);
147 while (tf < nyquist_freq) {
159 if (tb <= bark - 0.5)
165 for (i = 0; i < 32; i++) {
167 for (j = 31; j > 0 && freq <= freqmid[j]; j--);
171 for (j = 0; j < 32 && freq >= freqmid[j]; j++);
190 for (j = 0; j < avctx->
channels; j++) {
193 for (i = 0; i <
BANDS; i++)
202 for (i = 0; i <
COEFFS; i++)
204 for (i = 0; i < COEFFS / 2; i++) {
205 q->
post_cos[i] = (1.0f / 32768) * cos(i / 256.0 * M_PI);
206 q->
post_sin[i] = (1.0f / 32768) * sin(i / 256.0 * M_PI);
208 r1 = sin((i * 4.0 + 1.0) / 1024.0 * M_PI);
209 r2 = cos((i * 4.0 + 1.0) / 1024.0 * M_PI);
213 q->
pre_coef2[i] = -(r1 - r2) * sqrt(2.0);
215 q->
pre_coef1[i] = -(r1 + r2) * sqrt(2.0);
222 for (i = 0; i < 30; i++)
226 for (i = 0; i < 4 ; i++) {
227 for (j = 0; j < 4; j++) {
229 huffman_vlc[i][j].
table_allocated = vlc_offsets[i * 4 + j + 1] - vlc_offsets[i * 4 + j];
259 float *flcoeffs2,
int *bandWidthT,
260 float *flcoeffs3,
float *flcoeffs5)
265 float snr_limit = 1.e-30;
269 for (i = 0; i <
BANDS; i++) {
270 flcoeffs5[i] = workT2[i] = 0.0;
272 workT1[i] = flcoeffs1[i] * flcoeffs1[i];
273 flcoeffs3[i] = 2.0 * flcoeffs2[i];
276 flcoeffs3[i] = -30000.0;
278 workT3[i] = bandWidthT[i] * workT1[i] * 0.01;
279 if (workT3[i] <= snr_limit)
283 for (i = 0; i <
BANDS; i++) {
284 for (cnt2 = i; cnt2 < q->
cyclTab[i]; cnt2++)
285 flcoeffs5[cnt2] = flcoeffs5[cnt2] + workT3[i];
286 workT2[cnt2 - 1] = workT2[cnt2 - 1] + workT3[i];
289 for (i = 1; i <
BANDS; i++) {
290 accum = (workT2[i - 1] + accum) * q->
weights1[i - 1];
291 flcoeffs5[i] += accum;
294 for (i = 0; i <
BANDS; i++)
297 for (i = 0; i <
BANDS; i++) {
298 for (cnt2 = i - 1; cnt2 > q->
cyclTab2[i]; cnt2--)
299 flcoeffs5[cnt2] += workT3[i];
300 workT2[cnt2+1] += workT3[i];
305 for (i = BANDS-2; i >= 0; i--) {
306 accum = (workT2[i+1] + accum) * q->
weights2[i];
307 flcoeffs5[i] += accum;
322 s = stream_format_code >> 1;
323 hufftab[0] = &huffman_vlc[s][0];
324 hufftab[1] = &huffman_vlc[s][1];
325 hufftab[2] = &huffman_vlc[s][2];
326 hufftab[3] = &huffman_vlc[s][3];
329 if (stream_format_code & 4)
333 for (i = start; i <
BANDS; i++) {
335 hufftab[cb_sel[i]]->
bits, 2);
336 if (levlCoeffs[i] == 17)
348 for (i = 1; i <
BANDS; i++)
353 float *flcoeffs1,
float *flcoeffs2)
359 flcoeffs1[0] = 20000.0 / pow (2, levlCoeffBuf[0] * 0.18945);
360 flcoeffs2[0] =
log2f(flcoeffs1[0]);
364 for (i = 1; i <
BANDS; i++) {
365 level = levlCoeffBuf[i];
372 else if (level <= 24)
378 tmp2 += 0.83048 *
level;
387 float *old_floor,
float *flcoeffs1,
395 for (i = 0; i <
BANDS; i++) {
397 if (levlCoeffBuf[i] < 16) {
398 flcoeffs1[i] =
imc_exp_tab2[levlCoeffBuf[i]] * old_floor[i];
399 flcoeffs2[i] = (levlCoeffBuf[i] - 7) * 0.83048 + flcoeffs2[i];
401 flcoeffs1[i] = old_floor[i];
407 float *flcoeffs1,
float *flcoeffs2)
413 flcoeffs1[pos] = 20000.0 / pow (2, levlCoeffBuf[0] * 0.18945);
414 flcoeffs2[pos] =
log2f(flcoeffs1[0]);
415 tmp = flcoeffs1[pos];
416 tmp2 = flcoeffs2[pos];
419 for (i = 0; i <
BANDS; i++) {
422 level = *levlCoeffBuf++;
423 flcoeffs1[i] = tmp *
powf(10.0, -level * 0.4375);
424 flcoeffs2[i] = tmp2 - 1.4533435415 *
level;
432 int stream_format_code,
int freebits,
int flag)
435 const float limit = -1.e20;
444 float lowest = 1.e10;
450 for (i = 0; i <
BANDS; i++)
453 for (i = 0; i < BANDS - 1; i++)
457 highest = highest * 0.25;
459 for (i = 0; i <
BANDS; i++) {
476 if (stream_format_code & 0x2) {
483 for (i = (stream_format_code & 0x2) ? 4 : 0; i < BANDS - 1; i++) {
492 summa = (summa * 0.5 - freebits) / iacc;
495 for (i = 0; i < BANDS / 2; i++) {
496 rres = summer - freebits;
497 if ((rres >= -8) && (rres <= 8))
503 for (j = (stream_format_code & 0x2) ? 4 : 0; j <
BANDS; j++) {
504 cwlen = av_clipf(((chctx->
flcoeffs4[j] * 0.5) - summa + 0.5), 0, 6);
515 if (freebits < summer)
522 summa = (float)(summer - freebits) / ((t1 + 1) * iacc) + summa;
525 for (i = (stream_format_code & 0x2) ? 4 : 0; i <
BANDS; i++) {
530 if (freebits > summer) {
531 for (i = 0; i <
BANDS; i++) {
532 workT[i] = (chctx->
bitsBandT[i] == 6) ? -1.e20
539 if (highest <= -1.e20)
545 for (i = 0; i <
BANDS; i++) {
546 if (workT[i] > highest) {
552 if (highest > -1.e20) {
553 workT[found_indx] -= 2.0;
555 workT[found_indx] = -1.e20;
557 for (j =
band_tab[found_indx]; j < band_tab[found_indx + 1] && (freebits > summer); j++) {
562 }
while (freebits > summer);
564 if (freebits < summer) {
565 for (i = 0; i <
BANDS; i++) {
569 if (stream_format_code & 0x2) {
575 while (freebits < summer) {
578 for (i = 0; i <
BANDS; i++) {
579 if (workT[i] < lowest) {
586 workT[low_indx] = lowest + 2.0;
589 workT[low_indx] = 1.e20;
591 for (j =
band_tab[low_indx]; j <
band_tab[low_indx+1] && (freebits < summer); j++) {
608 for (i = 0; i <
BANDS; i++) {
615 for (j = band_tab[i]; j < band_tab[i + 1]; j++) {
646 if (j < band_tab[i + 1]) {
667 for (i = 0; i <
BANDS; i++) {
668 workT[i] = (chctx->
bitsBandT[i] == 6) ? -1.e20
672 while (corrected < summer) {
673 if (highest <= -1.e20)
678 for (i = 0; i <
BANDS; i++) {
679 if (workT[i] > highest) {
685 if (highest > -1.e20) {
686 workT[found_indx] -= 2.0;
687 if (++(chctx->
bitsBandT[found_indx]) == 6)
688 workT[found_indx] = -1.e20;
690 for (j =
band_tab[found_indx]; j <
band_tab[found_indx+1] && (corrected < summer); j++) {
708 for (i = 0; i <
COEFFS / 2; i++) {
720 for (i = 0; i <
COEFFS / 2; i++) {
734 int stream_format_code)
737 int middle_value, cw_len, max_size;
738 const float *quantizer;
740 for (i = 0; i <
BANDS; i++) {
745 if (cw_len <= 0 || chctx->skipFlags[j])
748 max_size = 1 << cw_len;
749 middle_value = max_size >> 1;
775 int i, j, cw_len, cw;
777 for (i = 0; i <
BANDS; i++) {
786 av_dlog(
NULL,
"Band %i coeff %i cw_len %i\n", i, j, cw_len);
805 for (i = 0; i <
BANDS; i++) {
811 if ((((band_tab[i + 1] - band_tab[i]) * 1.5) > chctx->
sumLenArr[i]) && (chctx->
sumLenArr[i] > 0))
817 for (i = 0; i <
BANDS; i++) {
829 for (i = 0; i <
BANDS; i++) {
846 int stream_format_code;
847 int imc_hdr, i, j, ret;
850 int counter, bitscount;
856 if (imc_hdr & 0x18) {
863 if (stream_format_code & 0x04)
867 for (i = 0; i <
BANDS; i++)
875 if (stream_format_code & 0x1)
878 else if (stream_format_code & 0x1)
883 if (stream_format_code & 0x4)
893 if (stream_format_code & 0x1) {
894 for (i = 0; i <
BANDS; i++) {
901 for (i = 0; i <
BANDS; i++) {
910 for (i = 0; i < BANDS - 1; i++)
921 if (stream_format_code & 0x2) {
928 for (i = 1; i < 4; i++) {
929 if (stream_format_code & 0x1)
942 if (!(stream_format_code & 0x2))
954 if (stream_format_code & 0x1) {
955 for (i = 0; i <
BANDS; i++)
961 for (i = 0; i <
BANDS; i++) {
991 int *got_frame_ptr,
AVPacket *avpkt)
995 int buf_size = avpkt->
size;
1002 if (buf_size < IMC_BLOCK_SIZE * avctx->channels) {
1014 for (i = 0; i < avctx->
channels; i++) {