47 #define IMC_BLOCK_SIZE 64
48 #define IMC_FRAME_ID 0x21
98 #define VLC_TABLES_SIZE 9512
101 0, 640, 1156, 1732, 2308, 2852, 3396, 3924,
119 for(i = 0; i <
BANDS; i++)
124 for(i = 0; i <
COEFFS; i++)
126 for(i = 0; i < COEFFS/2; i++){
127 q->
post_cos[i] = (1.0f / 32768) * cos(i / 256.0 *
M_PI);
128 q->
post_sin[i] = (1.0f / 32768) * sin(i / 256.0 *
M_PI);
130 r1 = sin((i * 4.0 + 1.0) / 1024.0 *
M_PI);
131 r2 = cos((i * 4.0 + 1.0) / 1024.0 *
M_PI);
136 q->
pre_coef2[i] = -(r1 - r2) * sqrt(2.0);
140 q->
pre_coef1[i] = -(r1 + r2) * sqrt(2.0);
149 for(i = 0; i < 30; i++) {
154 for(i = 0; i < 4 ; i++) {
155 for(j = 0; j < 4; j++) {
157 huffman_vlc[i][j].
table_allocated = vlc_offsets[i * 4 + j + 1] - vlc_offsets[i * 4 + j];
180 float* flcoeffs3,
float* flcoeffs5)
185 float snr_limit = 1.e-30;
189 for(i = 0; i <
BANDS; i++) {
190 flcoeffs5[i] = workT2[i] = 0.0;
192 workT1[i] = flcoeffs1[i] * flcoeffs1[i];
193 flcoeffs3[i] = 2.0 * flcoeffs2[i];
196 flcoeffs3[i] = -30000.0;
198 workT3[i] = bandWidthT[i] * workT1[i] * 0.01;
199 if (workT3[i] <= snr_limit)
203 for(i = 0; i <
BANDS; i++) {
204 for(cnt2 = i; cnt2 <
cyclTab[i]; cnt2++)
205 flcoeffs5[cnt2] = flcoeffs5[cnt2] + workT3[i];
206 workT2[cnt2-1] = workT2[cnt2-1] + workT3[i];
209 for(i = 1; i <
BANDS; i++) {
211 flcoeffs5[i] += accum;
214 for(i = 0; i <
BANDS; i++)
217 for(i = 0; i <
BANDS; i++) {
218 for(cnt2 = i-1; cnt2 >
cyclTab2[i]; cnt2--)
219 flcoeffs5[cnt2] += workT3[i];
220 workT2[cnt2+1] += workT3[i];
225 for(i = BANDS-2; i >= 0; i--) {
227 flcoeffs5[i] += accum;
238 const uint8_t *cb_sel;
241 s = stream_format_code >> 1;
242 hufftab[0] = &huffman_vlc[s][0];
243 hufftab[1] = &huffman_vlc[s][1];
244 hufftab[2] = &huffman_vlc[s][2];
245 hufftab[3] = &huffman_vlc[s][3];
248 if(stream_format_code & 4)
252 for(i = start; i <
BANDS; i++){
254 if(levlCoeffs[i] == 17)
266 flcoeffs1[0] = 20000.0 / pow (2, levlCoeffBuf[0] * 0.18945);
267 flcoeffs2[0] = log(flcoeffs1[0])/log(2);
271 for(i = 1; i <
BANDS; i++) {
272 level = levlCoeffBuf[i];
279 else if (level <= 24)
285 tmp2 += 0.83048 *
level;
299 for(i = 0; i <
BANDS; i++) {
301 if(levlCoeffBuf[i] < 16) {
302 flcoeffs1[i] =
imc_exp_tab2[levlCoeffBuf[i]] * old_floor[i];
303 flcoeffs2[i] = (levlCoeffBuf[i]-7) * 0.83048 + flcoeffs2[i];
305 flcoeffs1[i] = old_floor[i];
315 const float limit = -1.e20;
324 float lowest = 1.e10;
330 for(i = 0; i <
BANDS; i++)
333 for(i = 0; i < BANDS-1; i++) {
338 highest = highest * 0.25;
340 for(i = 0; i <
BANDS; i++) {
357 if (stream_format_code & 0x2) {
364 for(i = (stream_format_code & 0x2)?4:0; i < BANDS-1; i++) {
373 summa = (summa * 0.5 - freebits) / iacc;
376 for(i = 0; i < BANDS/2; i++) {
377 rres = summer - freebits;
378 if((rres >= -8) && (rres <= 8))
break;
383 for(j = (stream_format_code & 0x2)?4:0; j <
BANDS; j++) {
384 cwlen = av_clipf(((q->
flcoeffs4[j] * 0.5) - summa + 0.5), 0, 6);
395 if (freebits < summer)
402 summa = (float)(summer - freebits) / ((t1 + 1) * iacc) + summa;
405 for(i = (stream_format_code & 0x2)?4:0; i <
BANDS; i++) {
410 if (freebits > summer) {
411 for(i = 0; i <
BANDS; i++) {
418 if (highest <= -1.e20)
424 for(i = 0; i <
BANDS; i++) {
425 if (workT[i] > highest) {
431 if (highest > -1.e20) {
432 workT[found_indx] -= 2.0;
434 workT[found_indx] = -1.e20;
436 for(j =
band_tab[found_indx]; j < band_tab[found_indx+1] && (freebits > summer); j++){
441 }
while (freebits > summer);
443 if (freebits < summer) {
444 for(i = 0; i <
BANDS; i++) {
447 if (stream_format_code & 0x2) {
453 while (freebits < summer){
456 for(i = 0; i <
BANDS; i++) {
457 if (workT[i] < lowest) {
463 workT[low_indx] = lowest + 2.0;
466 workT[low_indx] = 1.e20;
468 for(j =
band_tab[low_indx]; j <
band_tab[low_indx+1] && (freebits < summer); j++){
484 for(i = 0; i <
BANDS; i++) {
491 for(j = band_tab[i]; j < band_tab[i+1]; j++) {
540 for(i = 0; i <
BANDS; i++) {
544 while (corrected < summer) {
545 if(highest <= -1.e20)
550 for(i = 0; i <
BANDS; i++) {
551 if (workT[i] > highest) {
557 if (highest > -1.e20) {
558 workT[found_indx] -= 2.0;
560 workT[found_indx] = -1.e20;
562 for(j =
band_tab[found_indx]; j <
band_tab[found_indx+1] && (corrected < summer); j++) {
577 for(i=0; i <
COEFFS/2; i++){
589 for(i = 0; i <
COEFFS/2; i++){
600 int middle_value, cw_len, max_size;
601 const float* quantizer;
603 for(i = 0; i <
BANDS; i++) {
608 if (cw_len <= 0 || q->skipFlags[j])
611 max_size = 1 << cw_len;
612 middle_value = max_size >> 1;
637 int i, j, cw_len, cw;
639 for(i = 0; i <
BANDS; i++) {
662 int *got_frame_ptr,
AVPacket *avpkt)
664 const uint8_t *buf = avpkt->
data;
665 int buf_size = avpkt->
size;
669 int stream_format_code;
670 int imc_hdr, i, j, ret;
673 int counter, bitscount;
702 if(stream_format_code & 1){
703 av_log(avctx,
AV_LOG_ERROR,
"Stream code format %X is not supported\n", stream_format_code);
709 if (stream_format_code & 0x04)
722 if (stream_format_code & 0x4)
730 for (i=0 ; i<
BANDS ; i++) {
738 for(i = 0; i < BANDS-1; i++) {
747 if (stream_format_code & 0x2) {
754 for(i = 1; i < 4; i++){
771 for(i = 0; i <
BANDS; i++) {
777 if( (((band_tab[i+1] - band_tab[i]) * 1.5) > q->
sumLenArr[i]) && (q->
sumLenArr[i] > 0))
783 for(i = 0; i <
BANDS; i++) {
795 for(i = 0; i <
BANDS; i++) {
809 for(i = 0; i <
BANDS; i++) {