48 #define IMC_BLOCK_SIZE 64
49 #define IMC_FRAME_ID 0x21
103 float weights1[31], weights2[31];
108 #define VLC_TABLES_SIZE 9512
111 0, 640, 1156, 1732, 2308, 2852, 3396, 3924,
119 return 3.5 * atan((freq / 7500.0) * (freq / 7500.0)) + 13.0 * atan(freq * 0.00076);
124 double freqmin[32], freqmid[32], freqmax[32];
125 double scale = sampling_rate / (256.0 * 2.0 * 2.0);
126 double nyquist_freq = sampling_rate * 0.5;
127 double freq, bark, prev_bark = 0, tf, tb;
130 for (i = 0; i < 32; i++) {
135 tb = bark - prev_bark;
136 q->
weights1[i - 1] = pow(10.0, -1.0 * tb);
137 q->
weights2[i - 1] = pow(10.0, -2.7 * tb);
144 while (tf < nyquist_freq) {
156 if (tb <= bark - 0.5)
162 for (i = 0; i < 32; i++) {
164 for (j = 31; j > 0 && freq <= freqmid[j]; j--);
168 for (j = 0; j < 32 && freq >= freqmid[j]; j++);
187 for (j = 0; j < avctx->
channels; j++) {
190 for (i = 0; i <
BANDS; i++)
199 for (i = 0; i <
COEFFS; i++)
201 for (i = 0; i < COEFFS / 2; i++) {
202 q->
post_cos[i] = (1.0f / 32768) * cos(i / 256.0 * M_PI);
203 q->
post_sin[i] = (1.0f / 32768) * sin(i / 256.0 * M_PI);
205 r1 = sin((i * 4.0 + 1.0) / 1024.0 * M_PI);
206 r2 = cos((i * 4.0 + 1.0) / 1024.0 * M_PI);
210 q->
pre_coef2[i] = -(r1 - r2) * sqrt(2.0);
212 q->
pre_coef1[i] = -(r1 + r2) * sqrt(2.0);
219 for (i = 0; i < 30; i++)
223 for (i = 0; i < 4 ; i++) {
224 for (j = 0; j < 4; j++) {
226 huffman_vlc[i][j].
table_allocated = vlc_offsets[i * 4 + j + 1] - vlc_offsets[i * 4 + j];
258 float *flcoeffs2,
int *bandWidthT,
259 float *flcoeffs3,
float *flcoeffs5)
264 float snr_limit = 1.e-30;
268 for (i = 0; i <
BANDS; i++) {
269 flcoeffs5[i] = workT2[i] = 0.0;
271 workT1[i] = flcoeffs1[i] * flcoeffs1[i];
272 flcoeffs3[i] = 2.0 * flcoeffs2[i];
275 flcoeffs3[i] = -30000.0;
277 workT3[i] = bandWidthT[i] * workT1[i] * 0.01;
278 if (workT3[i] <= snr_limit)
282 for (i = 0; i <
BANDS; i++) {
283 for (cnt2 = i; cnt2 < q->
cyclTab[i]; cnt2++)
284 flcoeffs5[cnt2] = flcoeffs5[cnt2] + workT3[i];
285 workT2[cnt2 - 1] = workT2[cnt2 - 1] + workT3[i];
288 for (i = 1; i <
BANDS; i++) {
289 accum = (workT2[i - 1] + accum) * q->
weights1[i - 1];
290 flcoeffs5[i] += accum;
293 for (i = 0; i <
BANDS; i++)
296 for (i = 0; i <
BANDS; i++) {
297 for (cnt2 = i - 1; cnt2 > q->
cyclTab2[i]; cnt2--)
298 flcoeffs5[cnt2] += workT3[i];
299 workT2[cnt2+1] += workT3[i];
304 for (i = BANDS-2; i >= 0; i--) {
305 accum = (workT2[i+1] + accum) * q->
weights2[i];
306 flcoeffs5[i] += accum;
321 s = stream_format_code >> 1;
322 hufftab[0] = &huffman_vlc[s][0];
323 hufftab[1] = &huffman_vlc[s][1];
324 hufftab[2] = &huffman_vlc[s][2];
325 hufftab[3] = &huffman_vlc[s][3];
328 if (stream_format_code & 4)
332 for (i = start; i <
BANDS; i++) {
334 hufftab[cb_sel[i]]->
bits, 2);
335 if (levlCoeffs[i] == 17)
341 float *flcoeffs1,
float *flcoeffs2)
347 flcoeffs1[0] = 20000.0 / pow (2, levlCoeffBuf[0] * 0.18945);
348 flcoeffs2[0] =
log2f(flcoeffs1[0]);
352 for (i = 1; i <
BANDS; i++) {
353 level = levlCoeffBuf[i];
360 else if (level <= 24)
366 tmp2 += 0.83048 *
level;
375 float *old_floor,
float *flcoeffs1,
383 for (i = 0; i <
BANDS; i++) {
385 if (levlCoeffBuf[i] < 16) {
386 flcoeffs1[i] =
imc_exp_tab2[levlCoeffBuf[i]] * old_floor[i];
387 flcoeffs2[i] = (levlCoeffBuf[i] - 7) * 0.83048 + flcoeffs2[i];
389 flcoeffs1[i] = old_floor[i];
398 int stream_format_code,
int freebits,
int flag)
401 const float limit = -1.e20;
410 float lowest = 1.e10;
416 for (i = 0; i <
BANDS; i++)
419 for (i = 0; i < BANDS - 1; i++)
423 highest = highest * 0.25;
425 for (i = 0; i <
BANDS; i++) {
442 if (stream_format_code & 0x2) {
449 for (i = (stream_format_code & 0x2) ? 4 : 0; i < BANDS - 1; i++) {
458 summa = (summa * 0.5 - freebits) / iacc;
461 for (i = 0; i < BANDS / 2; i++) {
462 rres = summer - freebits;
463 if ((rres >= -8) && (rres <= 8))
469 for (j = (stream_format_code & 0x2) ? 4 : 0; j <
BANDS; j++) {
470 cwlen = av_clipf(((chctx->
flcoeffs4[j] * 0.5) - summa + 0.5), 0, 6);
481 if (freebits < summer)
488 summa = (float)(summer - freebits) / ((t1 + 1) * iacc) + summa;
491 for (i = (stream_format_code & 0x2) ? 4 : 0; i <
BANDS; i++) {
496 if (freebits > summer) {
497 for (i = 0; i <
BANDS; i++) {
498 workT[i] = (chctx->
bitsBandT[i] == 6) ? -1.e20
505 if (highest <= -1.e20)
511 for (i = 0; i <
BANDS; i++) {
512 if (workT[i] > highest) {
518 if (highest > -1.e20) {
519 workT[found_indx] -= 2.0;
521 workT[found_indx] = -1.e20;
523 for (j =
band_tab[found_indx]; j < band_tab[found_indx + 1] && (freebits > summer); j++) {
528 }
while (freebits > summer);
530 if (freebits < summer) {
531 for (i = 0; i <
BANDS; i++) {
535 if (stream_format_code & 0x2) {
541 while (freebits < summer) {
544 for (i = 0; i <
BANDS; i++) {
545 if (workT[i] < lowest) {
552 workT[low_indx] = lowest + 2.0;
555 workT[low_indx] = 1.e20;
557 for (j =
band_tab[low_indx]; j <
band_tab[low_indx+1] && (freebits < summer); j++) {
574 for (i = 0; i <
BANDS; i++) {
581 for (j = band_tab[i]; j < band_tab[i + 1]; j++) {
612 if (j < band_tab[i + 1]) {
633 for (i = 0; i <
BANDS; i++) {
634 workT[i] = (chctx->
bitsBandT[i] == 6) ? -1.e20
638 while (corrected < summer) {
639 if (highest <= -1.e20)
644 for (i = 0; i <
BANDS; i++) {
645 if (workT[i] > highest) {
651 if (highest > -1.e20) {
652 workT[found_indx] -= 2.0;
653 if (++(chctx->
bitsBandT[found_indx]) == 6)
654 workT[found_indx] = -1.e20;
656 for (j =
band_tab[found_indx]; j <
band_tab[found_indx+1] && (corrected < summer); j++) {
674 for (i = 0; i <
COEFFS / 2; i++) {
686 for (i = 0; i <
COEFFS / 2; i++) {
700 int stream_format_code)
703 int middle_value, cw_len, max_size;
704 const float *quantizer;
706 for (i = 0; i <
BANDS; i++) {
711 if (cw_len <= 0 || chctx->skipFlags[j])
714 max_size = 1 << cw_len;
715 middle_value = max_size >> 1;
741 int i, j, cw_len, cw;
743 for (i = 0; i <
BANDS; i++) {
752 av_dlog(
NULL,
"Band %i coeff %i cw_len %i\n", i, j, cw_len);
768 int stream_format_code;
769 int imc_hdr, i, j, ret;
772 int counter, bitscount;
778 if (imc_hdr & 0x18) {
785 if (stream_format_code & 1) {
791 if (stream_format_code & 0x04)
795 for (i = 0; i <
BANDS; i++)
805 if (stream_format_code & 0x4)
815 for (i = 0; i <
BANDS; i++) {
823 for (i = 0; i < BANDS - 1; i++) {
832 if (stream_format_code & 0x2) {
839 for (i = 1; i < 4; i++) {
850 if (!(stream_format_code & 0x2))
862 for (i = 0; i <
BANDS; i++) {
868 if ((((band_tab[i + 1] - band_tab[i]) * 1.5) > chctx->
sumLenArr[i]) && (chctx->
sumLenArr[i] > 0))
874 for (i = 0; i <
BANDS; i++) {
886 for (i = 0; i <
BANDS; i++) {
900 for (i = 0; i <
BANDS; i++) {
930 int *got_frame_ptr,
AVPacket *avpkt)
933 int buf_size = avpkt->
size;
940 if (buf_size < IMC_BLOCK_SIZE * avctx->channels) {
952 for (i = 0; i < avctx->
channels; i++) {
int skipFlags[COEFFS]
skip coefficient decoding or not
struct IMCChannel IMCChannel
av_cold void ff_dsputil_init(DSPContext *c, AVCodecContext *avctx)
This structure describes decoded (raw) audio or video data.
int codewords[COEFFS]
raw codewords read from bitstream
static unsigned int get_bits(GetBitContext *s, int n)
Read 1-25 bits.
int skipFlagRaw[BANDS]
skip flags are stored in raw form or not
AVFrame * coded_frame
the picture in the bitstream
static const int vlc_offsets[17]
float mdct_sine_window[COEFFS]
MDCT tables.
void(* fft_permute)(struct FFTContext *s, FFTComplex *z)
Do the permutation needed BEFORE calling fft_calc().
static const uint8_t imc_huffman_lens[4][4][18]
int skipFlagCount[BANDS]
skipped coeffients per band
static const float imc_weights2[31]
av_dlog(ac->avr,"%d samples - audio_convert: %s to %s (%s)\n", len, av_get_sample_fmt_name(ac->in_fmt), av_get_sample_fmt_name(ac->out_fmt), use_generic?ac->func_descr_generic:ac->func_descr)
#define AV_CH_LAYOUT_STEREO
static void imc_read_level_coeffs(IMCContext *q, int stream_format_code, int *levlCoeffs)
static int decode(MimicContext *ctx, int quality, int num_coeffs, int is_iframe)
enum AVSampleFormat sample_fmt
audio sample format
int bandFlagsBuf[BANDS]
flags for each band
static av_cold int imc_decode_close(AVCodecContext *avctx)
static const int8_t cyclTab[32]
static int get_bits_count(const GetBitContext *s)
static const float imc_weights1[31]
bitstream reader API header.
static void imc_get_skip_coeff(IMCContext *q, IMCChannel *chctx)
#define LOCAL_ALIGNED_16(t, v,...)
static int init(AVCodecParserContext *s)
void(* bswap16_buf)(uint16_t *dst, const uint16_t *src, int len)
static const float *const imc_exp_tab2
void av_log_ask_for_sample(void *avc, const char *msg,...) av_printf_format(2
Log a generic warning message asking for a sample.
#define NULL_IF_CONFIG_SMALL(x)
Return NULL if CONFIG_SMALL is true, otherwise the argument without modification. ...
static void imc_imdct256(IMCContext *q, IMCChannel *chctx, int channels)
void av_log(void *avcl, int level, const char *fmt,...)
static void imc_calculate_coeffs(IMCContext *q, float *flcoeffs1, float *flcoeffs2, int *bandWidthT, float *flcoeffs3, float *flcoeffs5)
const char * name
Name of the codec implementation.
uint64_t channel_layout
Audio channel layout.
static const int8_t cyclTab2[32]
static int imc_decode_frame(AVCodecContext *avctx, void *data, int *got_frame_ptr, AVPacket *avpkt)
static void imc_decode_level_coefficients2(IMCContext *q, int *levlCoeffBuf, float *old_floor, float *flcoeffs1, float *flcoeffs2)
audio channel layout utility functions
static const uint16_t band_tab[33]
int ff_get_buffer(AVCodecContext *avctx, AVFrame *frame)
Get a buffer for a frame.
static int imc_get_coeffs(IMCContext *q, IMCChannel *chctx)
int bitsBandT[BANDS]
how many bits per codeword in band
static av_always_inline int get_vlc2(GetBitContext *s, VLC_TYPE(*table)[2], int bits, int max_depth)
Parse a vlc code.
#define INIT_VLC_USE_NEW_STATIC
float last_fft_im[COEFFS]
void ff_sine_window_init(float *window, int n)
Generate a sine window.
static const float xTab[14]
FFTComplex samples[COEFFS/2]
static int inverse_quant_coeff(IMCContext *q, IMCChannel *chctx, int stream_format_code)
int sample_rate
samples per second
main external API structure.
static void close(AVCodecParserContext *s)
void(* butterflies_float)(float *restrict v1, float *restrict v2, int len)
Calculate the sum and difference of two vectors of floats.
#define init_vlc(vlc, nb_bits, nb_codes,bits, bits_wrap, bits_size,codes, codes_wrap, codes_size,flags)
static double freq2bark(double freq)
static unsigned int get_bits1(GetBitContext *s)
void avcodec_get_frame_defaults(AVFrame *frame)
Set the fields of the given AVFrame to default values.
int bandWidthT[BANDS]
codewords per band
static int init_get_bits(GetBitContext *s, const uint8_t *buffer, int bit_size)
Initialize GetBitContext.
static VLC huffman_vlc[4][4]
static VLC_TYPE vlc_tables[VLC_TABLES_SIZE][2]
static const float imc_quantizer1[4][8]
static const uint16_t scale[4]
#define DECLARE_ALIGNED(n, t, v)
static av_cold void iac_generate_tabs(IMCContext *q, int sampling_rate)
common internal api header.
static int imc_decode_block(AVCodecContext *avctx, IMCContext *q, int ch)
#define AVERROR_INVALIDDATA
AVSampleFormat
Audio Sample Formats.
static int bit_allocation(IMCContext *q, IMCChannel *chctx, int stream_format_code, int freebits, int flag)
Perform bit allocation depending on bits available.
void(* fft_calc)(struct FFTContext *s, FFTComplex *z)
Do a complex FFT with the parameters defined in ff_fft_init().
#define AVERROR_PATCHWELCOME
static const uint8_t imc_cb_select[4][32]
int channels
number of audio channels
VLC_TYPE(* table)[2]
code, bits
static void imc_adjust_bit_allocation(IMCContext *q, IMCChannel *chctx, int summer)
Increase highest' band coefficient sizes as some bits won't be used.
static const float imc_quantizer2[2][56]
int sumLenArr[BANDS]
bits for all coeffs in band
static const uint8_t imc_huffman_sizes[4]
uint8_t ** extended_data
pointers to the data planes/channels.
#define AV_CH_LAYOUT_MONO
This structure stores compressed data.
static void imc_decode_level_coefficients(IMCContext *q, int *levlCoeffBuf, float *flcoeffs1, float *flcoeffs2)
int skipFlagBits[BANDS]
bits used to code skip flags
static av_cold int imc_decode_init(AVCodecContext *avctx)
int nb_samples
number of audio samples (per channel) described by this frame
static const float imc_exp_tab[32]
static const uint16_t imc_huffman_bits[4][4][18]
int CWlengthT[COEFFS]
how many bits in each codeword
if(!(ptr_align%ac->ptr_align)&&samples_align >=aligned_len)