48 #define IMC_BLOCK_SIZE 64
49 #define IMC_FRAME_ID 0x21
99 #define VLC_TABLES_SIZE 9512
102 0, 640, 1156, 1732, 2308, 2852, 3396, 3924,
120 for(i = 0; i <
BANDS; i++)
125 for(i = 0; i <
COEFFS; i++)
127 for(i = 0; i < COEFFS/2; i++){
128 q->
post_cos[i] = (1.0f / 32768) * cos(i / 256.0 *
M_PI);
129 q->
post_sin[i] = (1.0f / 32768) * sin(i / 256.0 *
M_PI);
131 r1 = sin((i * 4.0 + 1.0) / 1024.0 *
M_PI);
132 r2 = cos((i * 4.0 + 1.0) / 1024.0 *
M_PI);
137 q->
pre_coef2[i] = -(r1 - r2) * sqrt(2.0);
141 q->
pre_coef1[i] = -(r1 + r2) * sqrt(2.0);
150 for(i = 0; i < 30; i++) {
155 for(i = 0; i < 4 ; i++) {
156 for(j = 0; j < 4; j++) {
158 huffman_vlc[i][j].
table_allocated = vlc_offsets[i * 4 + j + 1] - vlc_offsets[i * 4 + j];
181 float* flcoeffs3,
float* flcoeffs5)
186 float snr_limit = 1.e-30;
190 for(i = 0; i <
BANDS; i++) {
191 flcoeffs5[i] = workT2[i] = 0.0;
193 workT1[i] = flcoeffs1[i] * flcoeffs1[i];
194 flcoeffs3[i] = 2.0 * flcoeffs2[i];
197 flcoeffs3[i] = -30000.0;
199 workT3[i] = bandWidthT[i] * workT1[i] * 0.01;
200 if (workT3[i] <= snr_limit)
204 for(i = 0; i <
BANDS; i++) {
205 for(cnt2 = i; cnt2 <
cyclTab[i]; cnt2++)
206 flcoeffs5[cnt2] = flcoeffs5[cnt2] + workT3[i];
207 workT2[cnt2-1] = workT2[cnt2-1] + workT3[i];
210 for(i = 1; i <
BANDS; i++) {
212 flcoeffs5[i] += accum;
215 for(i = 0; i <
BANDS; i++)
218 for(i = 0; i <
BANDS; i++) {
219 for(cnt2 = i-1; cnt2 >
cyclTab2[i]; cnt2--)
220 flcoeffs5[cnt2] += workT3[i];
221 workT2[cnt2+1] += workT3[i];
226 for(i = BANDS-2; i >= 0; i--) {
228 flcoeffs5[i] += accum;
239 const uint8_t *cb_sel;
242 s = stream_format_code >> 1;
243 hufftab[0] = &huffman_vlc[s][0];
244 hufftab[1] = &huffman_vlc[s][1];
245 hufftab[2] = &huffman_vlc[s][2];
246 hufftab[3] = &huffman_vlc[s][3];
249 if(stream_format_code & 4)
253 for(i = start; i <
BANDS; i++){
255 if(levlCoeffs[i] == 17)
267 flcoeffs1[0] = 20000.0 / pow (2, levlCoeffBuf[0] * 0.18945);
268 flcoeffs2[0] = log(flcoeffs1[0])/log(2);
272 for(i = 1; i <
BANDS; i++) {
273 level = levlCoeffBuf[i];
280 else if (level <= 24)
286 tmp2 += 0.83048 *
level;
300 for(i = 0; i <
BANDS; i++) {
302 if(levlCoeffBuf[i] < 16) {
303 flcoeffs1[i] =
imc_exp_tab2[levlCoeffBuf[i]] * old_floor[i];
304 flcoeffs2[i] = (levlCoeffBuf[i]-7) * 0.83048 + flcoeffs2[i];
306 flcoeffs1[i] = old_floor[i];
316 const float limit = -1.e20;
325 float lowest = 1.e10;
331 for(i = 0; i <
BANDS; i++)
334 for(i = 0; i < BANDS-1; i++) {
339 highest = highest * 0.25;
341 for(i = 0; i <
BANDS; i++) {
358 if (stream_format_code & 0x2) {
365 for(i = (stream_format_code & 0x2)?4:0; i < BANDS-1; i++) {
374 summa = (summa * 0.5 - freebits) / iacc;
377 for(i = 0; i < BANDS/2; i++) {
378 rres = summer - freebits;
379 if((rres >= -8) && (rres <= 8))
break;
384 for(j = (stream_format_code & 0x2)?4:0; j <
BANDS; j++) {
385 cwlen = av_clipf(((q->
flcoeffs4[j] * 0.5) - summa + 0.5), 0, 6);
396 if (freebits < summer)
403 summa = (float)(summer - freebits) / ((t1 + 1) * iacc) + summa;
406 for(i = (stream_format_code & 0x2)?4:0; i <
BANDS; i++) {
411 if (freebits > summer) {
412 for(i = 0; i <
BANDS; i++) {
419 if (highest <= -1.e20)
425 for(i = 0; i <
BANDS; i++) {
426 if (workT[i] > highest) {
432 if (highest > -1.e20) {
433 workT[found_indx] -= 2.0;
435 workT[found_indx] = -1.e20;
437 for(j =
band_tab[found_indx]; j < band_tab[found_indx+1] && (freebits > summer); j++){
442 }
while (freebits > summer);
444 if (freebits < summer) {
445 for(i = 0; i <
BANDS; i++) {
448 if (stream_format_code & 0x2) {
454 while (freebits < summer){
457 for(i = 0; i <
BANDS; i++) {
458 if (workT[i] < lowest) {
464 workT[low_indx] = lowest + 2.0;
467 workT[low_indx] = 1.e20;
469 for(j =
band_tab[low_indx]; j <
band_tab[low_indx+1] && (freebits < summer); j++){
485 for(i = 0; i <
BANDS; i++) {
492 for(j = band_tab[i]; j < band_tab[i+1]; j++) {
541 for(i = 0; i <
BANDS; i++) {
545 while (corrected < summer) {
546 if(highest <= -1.e20)
551 for(i = 0; i <
BANDS; i++) {
552 if (workT[i] > highest) {
558 if (highest > -1.e20) {
559 workT[found_indx] -= 2.0;
561 workT[found_indx] = -1.e20;
563 for(j =
band_tab[found_indx]; j <
band_tab[found_indx+1] && (corrected < summer); j++) {
578 for(i=0; i <
COEFFS/2; i++){
590 for(i = 0; i <
COEFFS/2; i++){
601 int middle_value, cw_len, max_size;
602 const float* quantizer;
604 for(i = 0; i <
BANDS; i++) {
609 if (cw_len <= 0 || q->skipFlags[j])
612 max_size = 1 << cw_len;
613 middle_value = max_size >> 1;
638 int i, j, cw_len, cw;
640 for(i = 0; i <
BANDS; i++) {
663 int *got_frame_ptr,
AVPacket *avpkt)
665 const uint8_t *buf = avpkt->
data;
666 int buf_size = avpkt->
size;
670 int stream_format_code;
671 int imc_hdr, i, j, ret;
674 int counter, bitscount;
703 if(stream_format_code & 1){
704 av_log(avctx,
AV_LOG_ERROR,
"Stream code format %X is not supported\n", stream_format_code);
710 if (stream_format_code & 0x04)
723 if (stream_format_code & 0x4)
731 for (i=0 ; i<
BANDS ; i++) {
739 for(i = 0; i < BANDS-1; i++) {
748 if (stream_format_code & 0x2) {
755 for(i = 1; i < 4; i++){
772 for(i = 0; i <
BANDS; i++) {
778 if( (((band_tab[i+1] - band_tab[i]) * 1.5) > q->
sumLenArr[i]) && (q->
sumLenArr[i] > 0))
784 for(i = 0; i <
BANDS; i++) {
796 for(i = 0; i <
BANDS; i++) {
810 for(i = 0; i <
BANDS; i++) {