44 #define EXPMAX ((19+EXPVLCBITS-1)/EXPVLCBITS)
46 #define HGAINVLCBITS 9
47 #define HGAINMAX ((13+HGAINVLCBITS-1)/HGAINVLCBITS)
143 t.v = ((u.v <<
LSP_POW_BITS) & ((1 << 23) - 1)) | (127 << 23);
154 wdel = M_PI / frame_len;
155 for(i=0;i<frame_len;i++)
182 float *out,
float *val_max_ptr,
186 float p, q, w, v, val_max;
205 *val_max_ptr = val_max;
217 if (i == 0 || i >= 8)
230 1.7782794100389e-04, 2.0535250264571e-04,
231 2.3713737056617e-04, 2.7384196342644e-04,
232 3.1622776601684e-04, 3.6517412725484e-04,
233 4.2169650342858e-04, 4.8696752516586e-04,
234 5.6234132519035e-04, 6.4938163157621e-04,
235 7.4989420933246e-04, 8.6596432336006e-04,
236 1.0000000000000e-03, 1.1547819846895e-03,
237 1.3335214321633e-03, 1.5399265260595e-03,
238 1.7782794100389e-03, 2.0535250264571e-03,
239 2.3713737056617e-03, 2.7384196342644e-03,
240 3.1622776601684e-03, 3.6517412725484e-03,
241 4.2169650342858e-03, 4.8696752516586e-03,
242 5.6234132519035e-03, 6.4938163157621e-03,
243 7.4989420933246e-03, 8.6596432336006e-03,
244 1.0000000000000e-02, 1.1547819846895e-02,
245 1.3335214321633e-02, 1.5399265260595e-02,
246 1.7782794100389e-02, 2.0535250264571e-02,
247 2.3713737056617e-02, 2.7384196342644e-02,
248 3.1622776601684e-02, 3.6517412725484e-02,
249 4.2169650342858e-02, 4.8696752516586e-02,
250 5.6234132519035e-02, 6.4938163157621e-02,
251 7.4989420933246e-02, 8.6596432336007e-02,
252 1.0000000000000e-01, 1.1547819846895e-01,
253 1.3335214321633e-01, 1.5399265260595e-01,
254 1.7782794100389e-01, 2.0535250264571e-01,
255 2.3713737056617e-01, 2.7384196342644e-01,
256 3.1622776601684e-01, 3.6517412725484e-01,
257 4.2169650342858e-01, 4.8696752516586e-01,
258 5.6234132519035e-01, 6.4938163157621e-01,
259 7.4989420933246e-01, 8.6596432336007e-01,
260 1.0000000000000e+00, 1.1547819846895e+00,
261 1.3335214321633e+00, 1.5399265260595e+00,
262 1.7782794100389e+00, 2.0535250264571e+00,
263 2.3713737056617e+00, 2.7384196342644e+00,
264 3.1622776601684e+00, 3.6517412725484e+00,
265 4.2169650342858e+00, 4.8696752516586e+00,
266 5.6234132519035e+00, 6.4938163157621e+00,
267 7.4989420933246e+00, 8.6596432336007e+00,
268 1.0000000000000e+01, 1.1547819846895e+01,
269 1.3335214321633e+01, 1.5399265260595e+01,
270 1.7782794100389e+01, 2.0535250264571e+01,
271 2.3713737056617e+01, 2.7384196342644e+01,
272 3.1622776601684e+01, 3.6517412725484e+01,
273 4.2169650342858e+01, 4.8696752516586e+01,
274 5.6234132519035e+01, 6.4938163157621e+01,
275 7.4989420933246e+01, 8.6596432336007e+01,
276 1.0000000000000e+02, 1.1547819846895e+02,
277 1.3335214321633e+02, 1.5399265260595e+02,
278 1.7782794100389e+02, 2.0535250264571e+02,
279 2.3713737056617e+02, 2.7384196342644e+02,
280 3.1622776601684e+02, 3.6517412725484e+02,
281 4.2169650342858e+02, 4.8696752516586e+02,
282 5.6234132519035e+02, 6.4938163157621e+02,
283 7.4989420933246e+02, 8.6596432336007e+02,
284 1.0000000000000e+03, 1.1547819846895e+03,
285 1.3335214321633e+03, 1.5399265260595e+03,
286 1.7782794100389e+03, 2.0535250264571e+03,
287 2.3713737056617e+03, 2.7384196342644e+03,
288 3.1622776601684e+03, 3.6517412725484e+03,
289 4.2169650342858e+03, 4.8696752516586e+03,
290 5.6234132519035e+03, 6.4938163157621e+03,
291 7.4989420933246e+03, 8.6596432336007e+03,
292 1.0000000000000e+04, 1.1547819846895e+04,
293 1.3335214321633e+04, 1.5399265260595e+04,
294 1.7782794100389e+04, 2.0535250264571e+04,
295 2.3713737056617e+04, 2.7384196342644e+04,
296 3.1622776601684e+04, 3.6517412725484e+04,
297 4.2169650342858e+04, 4.8696752516586e+04,
298 5.6234132519035e+04, 6.4938163157621e+04,
299 7.4989420933246e+04, 8.6596432336007e+04,
300 1.0000000000000e+05, 1.1547819846895e+05,
301 1.3335214321633e+05, 1.5399265260595e+05,
302 1.7782794100389e+05, 2.0535250264571e+05,
303 2.3713737056617e+05, 2.7384196342644e+05,
304 3.1622776601684e+05, 3.6517412725484e+05,
305 4.2169650342858e+05, 4.8696752516586e+05,
306 5.6234132519035e+05, 6.4938163157621e+05,
307 7.4989420933246e+05, 8.6596432336007e+05,
315 int last_exp, n, code;
318 uint32_t *q, *q_end, iv;
319 const float *ptab = pow_tab + 60;
320 const uint32_t *iptab = (
const uint32_t*)ptab;
329 iv = iptab[last_exp];
337 }
while ((n -= 4) > 0);
348 last_exp += code - 60;
355 iv = iptab[last_exp];
364 }
while ((n -= 4) > 0);
380 int block_len, bsize, n;
398 memcpy(out+n+block_len, in+n+block_len, n*
sizeof(
float));
416 memcpy(out, in, n*
sizeof(
float));
420 memset(out+n+block_len, 0, n*
sizeof(
float));
431 int n, v,
a, ch, bsize;
432 int coef_nb_bits, total_gain;
539 val = (int)0x80000000;
542 if (val == (
int)0x80000000) {
587 0, ptr, 0, nb_coefs[ch],
598 mdct_norm = 1.0 / (float)n4;
600 mdct_norm *= sqrt(n4);
608 float *coefs, *exponents,
mult, mult1,
noise;
609 int i, j, n, n1, last_high_band, esize;
615 mult = pow(10, total_gain * 0.05) / s->
max_exponent[ch];
617 coefs = s->
coefs[ch];
623 exponents[i<<bsize>>esize] * mult1;
639 for(i = 0;i < n; i++) {
640 v = exponents[i<<bsize>>esize];
643 exp_power[j] = e2 / n;
645 tprintf(s->
avctx,
"%d: power=%f (%d)\n", j, exp_power[j], n);
647 exponents += n<<bsize>>esize;
662 mult1 = sqrt(exp_power[j] / exp_power[last_high_band]);
667 for(i = 0;i < n; i++) {
671 exponents[i<<bsize>>esize] * mult1;
673 exponents += n<<bsize>>esize;
676 for(i = 0;i < n; i++) {
679 *coefs++ = ((*coefs1++) + noise) *
680 exponents[i<<bsize>>esize] *
mult;
682 exponents += n<<bsize>>esize;
688 mult1 = mult * exponents[((-1<<bsize))>>esize];
689 for(i = 0; i < n; i++) {
698 for(i = 0;i < n; i++) {
699 *coefs++ = coefs1[i] * exponents[i<<bsize>>esize] *
mult;
702 for(i = 0;i < n; i++)
779 memcpy(samples[ch] + samples_offset, s->
frame_out[ch],
786 dump_floats(s,
"samples", 6, samples[ch] + samples_offset, s->
frame_len);
794 int *got_frame_ptr,
AVPacket *avpkt)
797 int buf_size = avpkt->
size;
799 int nb_frames, bit_offset, i, pos,
len, ret;
804 tprintf(avctx,
"***decode_superframe:\n");
810 if (buf_size < avctx->block_align) {
812 "Input packet size too small (%d < %d)\n",
841 "Invalid last frame bit offset %d > buf size %d (%d)\n",
885 for(i=0;i<nb_frames;i++) {
895 len = buf_size - pos;
909 av_dlog(s->
avctx,
"%d %d %d %d outbytes:%td eaten:%d\n",
911 (int8_t *)samples - (int8_t *)data, avctx->
block_align);
const struct AVCodec * codec
static void wma_lsp_to_curve(WMACodecContext *s, float *out, float *val_max_ptr, int n, float *lsp)
NOTE: We use the same code as Vorbis here.
uint8_t last_superframe[MAX_CODED_SUPERFRAME_SIZE+FF_INPUT_BUFFER_PADDING_SIZE]
This structure describes decoded (raw) audio or video data.
static unsigned int get_bits(GetBitContext *s, int n)
Read 1-25 bits.
int next_block_len_bits
log2 of next block length
AVFrame * coded_frame
the picture in the bitstream
static void align_get_bits(GetBitContext *s)
static const float pow_tab[]
pow(10, i / 16.0) for i in -60..95
int ff_wma_run_level_decode(AVCodecContext *avctx, GetBitContext *gb, VLC *vlc, const float *level_table, const uint16_t *run_table, int version, WMACoef *ptr, int offset, int num_coefs, int block_len, int frame_len_bits, int coef_nb_bits)
Decode run level compressed coefficients.
int block_len
block length in samples
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)
float exponents[MAX_CHANNELS][BLOCK_MAX_SIZE]
static void wma_window(WMACodecContext *s, float *out)
Apply MDCT window and add into output.
int block_align
number of bytes per packet if constant and known or 0 Used by some WAV based audio codecs...
float lsp_pow_m_table2[(1<< LSP_POW_BITS)]
static int wma_decode_block(WMACodecContext *s)
float lsp_cos_table[BLOCK_MAX_SIZE]
int high_band_start[BLOCK_NB_SIZES]
index of first coef in high band
static int decode(MimicContext *ctx, int quality, int num_coeffs, int is_iframe)
enum AVSampleFormat sample_fmt
audio sample format
float WMACoef
type for decoded coefficients, int16_t would be enough for wma 1/2
const uint8_t ff_aac_scalefactor_bits[121]
uint8_t * extradata
some codecs need / can use extradata like Huffman tables.
int block_pos
current position in frame
static int decode_exp_vlc(WMACodecContext *s, int ch)
decode exponents coded with VLC codes
static int get_bits_count(const GetBitContext *s)
float lsp_pow_m_table1[(1<< LSP_POW_BITS)]
static int init(AVCodecParserContext *s)
void(* vector_fmul_add)(float *dst, const float *src0, const float *src1, const float *src2, int len)
int nb_block_sizes
number of block sizes
int ff_wma_total_gain_to_bits(int total_gain)
static int get_bits_left(GetBitContext *gb)
static float pow_m1_4(WMACodecContext *s, float x)
compute x^-0.25 with an exponent and mantissa table.
uint16_t exponent_bands[BLOCK_NB_SIZES][25]
uint8_t channel_coded[MAX_CHANNELS]
true if channel is coded
#define NULL_IF_CONFIG_SMALL(x)
Return NULL if CONFIG_SMALL is true, otherwise the argument without modification. ...
static void wma_lsp_to_curve_init(WMACodecContext *s, int frame_len)
void av_log(void *avcl, int level, const char *fmt,...)
const char * name
Name of the codec implementation.
static int wma_decode_superframe(AVCodecContext *avctx, void *data, int *got_frame_ptr, AVPacket *avpkt)
FFTSample output[BLOCK_MAX_SIZE *2]
const uint8_t ff_wma_hgain_huffbits[37]
static int wma_decode_init(AVCodecContext *avctx)
void(* imdct_calc)(struct FFTContext *s, FFTSample *output, const FFTSample *input)
int exponent_high_bands[BLOCK_NB_SIZES][HIGH_BAND_MAX_SIZE]
int ff_wma_end(AVCodecContext *avctx)
int high_band_values[MAX_CHANNELS][HIGH_BAND_MAX_SIZE]
float * windows[BLOCK_NB_SIZES]
#define MAX_CODED_SUPERFRAME_SIZE
int ff_get_buffer(AVCodecContext *avctx, AVFrame *frame)
Get a buffer for a frame.
const uint16_t ff_wma_hgain_huffcodes[37]
int version
1 = 0x160 (WMAV1), 2 = 0x161 (WMAV2)
static av_always_inline int get_vlc2(GetBitContext *s, VLC_TYPE(*table)[2], int bits, int max_depth)
Parse a vlc code.
int frame_len
frame length in samples
static int noise(AVBitStreamFilterContext *bsfc, AVCodecContext *avctx, const char *args, uint8_t **poutbuf, int *poutbuf_size, const uint8_t *buf, int buf_size, int keyframe)
static av_cold void flush(AVCodecContext *avctx)
int frame_len_bits
frame_len = 1 << frame_len_bits
static int wma_decode_frame(WMACodecContext *s, float **samples, int samples_offset)
#define HIGH_BAND_MAX_SIZE
int use_exp_vlc
exponent coding: 0 = lsp, 1 = vlc + delta
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)
float frame_out[MAX_CHANNELS][BLOCK_MAX_SIZE *2]
static int16_t mult(Float11 *f1, Float11 *f2)
int exponent_high_sizes[BLOCK_NB_SIZES]
static void decode_exp_lsp(WMACodecContext *s, int ch)
decode exponents coded with LSP coefficients (same idea as Vorbis)
static unsigned int get_bits1(GetBitContext *s)
void avcodec_get_frame_defaults(AVFrame *frame)
Set the fields of the given AVFrame to default values.
static void skip_bits(GetBitContext *s, int n)
#define FF_ARRAY_ELEMS(a)
int block_num
block number in current frame
int use_noise_coding
true if perceptual noise is added
static int init_get_bits(GetBitContext *s, const uint8_t *buffer, int bit_size)
Initialize GetBitContext.
#define FF_INPUT_BUFFER_PADDING_SIZE
int use_variable_block_len
uint8_t ms_stereo
true if mid/side stereo mode
FFTContext mdct_ctx[BLOCK_NB_SIZES]
const uint32_t ff_aac_scalefactor_code[121]
int exponents_bsize[MAX_CHANNELS]
log2 ratio frame/exp. length
void(* vector_fmul_reverse)(float *dst, const float *src0, const float *src1, int len)
float coefs[MAX_CHANNELS][BLOCK_MAX_SIZE]
int prev_block_len_bits
log2 of prev block length
int coefs_end[BLOCK_NB_SIZES]
max number of coded coefficients
struct WMACodecContext WMACodecContext
float lsp_pow_e_table[256]
const float ff_wma_lsp_codebook[NB_LSP_COEFS][16]
common internal api header.
#define AVERROR_INVALIDDATA
AVSampleFormat
Audio Sample Formats.
int channels
number of audio channels
VLC_TYPE(* table)[2]
code, bits
WMACoef coefs1[MAX_CHANNELS][BLOCK_MAX_SIZE]
static const struct twinvq_data tab
float max_exponent[MAX_CHANNELS]
int coefs_start
first coded coef
int block_len_bits
log2 of current block length
uint8_t ** extended_data
pointers to the data planes/channels.
This structure stores compressed data.
int nb_samples
number of audio samples (per channel) described by this frame
int high_band_coded[MAX_CHANNELS][HIGH_BAND_MAX_SIZE]
int ff_wma_init(AVCodecContext *avctx, int flags2)
float noise_table[NOISE_TAB_SIZE]
if(!(ptr_align%ac->ptr_align)&&samples_align >=aligned_len)