45 #include "libavutil/channel_layout.h"
46 #include "libavutil/lfg.h"
58 #define MONO 0x1000001
59 #define STEREO 0x1000002
60 #define JOINT_STEREO 0x1000003
61 #define MC_COOK 0x2000000 // multichannel Cook, not supported
63 #define SUBBAND_SIZE 20
64 #define MAX_SUBPACKETS 5
89 float mono_previous_buffer1[1024];
90 float mono_previous_buffer2[1024];
100 typedef struct cook {
105 void (*scalar_dequant)(
struct cook *q,
int index,
int quant_index,
106 int *subband_coef_index,
int *subband_coef_sign,
109 void (*decouple)(
struct cook *q,
113 float *decode_buffer,
114 float *mlt_buffer1,
float *mlt_buffer2);
116 void (*imlt_window)(
struct cook *q,
float *buffer1,
117 cook_gains *gains_ptr,
float *previous_buffer);
120 int gain_index,
int gain_index_next);
122 void (*saturate_output)(
struct cook *q,
float *
out);
139 VLC envelope_quant_index[13];
144 float gain_table[23];
150 float decode_buffer_1[1024];
151 float decode_buffer_2[1024];
152 float decode_buffer_0[1060];
168 for (i = -63; i < 64; i++) {
179 for (i = 0; i < 23; i++)
190 for (i = 0; i < 13; i++) {
196 for (i = 0; i < 7; i++) {
226 for (j = 0; j < mlt_size; j++)
243 for (i = 0; i < 5; i++)
249 #define DECODE_BYTES_PAD1(bytes) (3 - ((bytes) + 3) % 4)
250 #define DECODE_BYTES_PAD2(bytes) ((bytes) % 4 + DECODE_BYTES_PAD1(2 * (bytes)))
274 static const uint32_t
tab[4] = {
281 uint32_t *obuf = (uint32_t *) out;
288 off = (intptr_t) inbuffer & 3;
289 buf = (
const uint32_t *) (inbuffer - off);
292 for (i = 0; i < bytes / 4; i++)
293 obuf[i] = c ^ buf[i];
312 for (i = 0; i < 13; i++)
314 for (i = 0; i < 7; i++)
346 gaininfo[i++] = gain;
359 int *quant_index_table)
363 quant_index_table[0] =
get_bits(&q->
gb, 6) - 6;
379 quant_index_table[i] = quant_index_table[i - 1] + j - 12;
380 if (quant_index_table[i] > 63 || quant_index_table[i] < -63) {
382 "Invalid quantizer %d at position %d, outside [-63, 63] range\n",
383 quant_index_table[i], i);
400 int *category,
int *category_index)
402 int exp_idx, bias, tmpbias1, tmpbias2, bits_left, num_bits,
index,
v, i, j;
403 int exp_index2[102] = { 0 };
404 int exp_index1[102] = { 0 };
406 int tmp_categorize_array[128 * 2] = { 0 };
419 for (i = 32; i > 0; i = i / 2) {
423 exp_idx = av_clip((i - quant_index_table[index] + bias) / 2, 0, 7);
427 if (num_bits >= bits_left - 32)
434 exp_idx = av_clip((bias - quant_index_table[i]) / 2, 0, 7);
436 exp_index1[i] = exp_idx;
437 exp_index2[i] = exp_idx;
439 tmpbias1 = tmpbias2 = num_bits;
442 if (tmpbias1 + tmpbias2 > 2 * bits_left) {
446 if (exp_index1[i] < 7) {
447 v = (-2 * exp_index1[i]) - quant_index_table[i] + bias;
456 tmp_categorize_array[tmp_categorize_array1_idx++] =
index;
464 if (exp_index2[i] > 0) {
465 v = (-2 * exp_index2[i]) - quant_index_table[i] + bias;
474 tmp_categorize_array[--tmp_categorize_array2_idx] =
index;
482 category[i] = exp_index2[i];
485 category_index[i] = tmp_categorize_array[tmp_categorize_array2_idx++];
502 int idx = category_index[i];
519 int *subband_coef_index,
int *subband_coef_sign,
526 if (subband_coef_index[i]) {
528 if (subband_coef_sign[i])
548 int *subband_coef_index,
int *subband_coef_sign)
551 int vlc, vd, tmp, result;
555 for (i = 0; i <
vpr_tab[category]; i++) {
561 for (j = vd - 1; j >= 0; j--) {
563 subband_coef_index[vd * i + j] = vlc - tmp * (
kmax_tab[category] + 1);
566 for (j = 0; j < vd; j++) {
567 if (subband_coef_index[i * vd + j]) {
572 subband_coef_sign[i * vd + j] = 0;
575 subband_coef_sign[i * vd + j] = 0;
592 int *quant_index_table,
float *mlt_buffer)
604 index = category[
band];
605 if (category[band] < 7) {
606 if (
unpack_SQVH(q, p, category[band], subband_coef_index, subband_coef_sign)) {
609 category[band + j] = 7;
613 memset(subband_coef_index, 0,
sizeof(subband_coef_index));
614 memset(subband_coef_sign, 0,
sizeof(subband_coef_sign));
617 subband_coef_index, subband_coef_sign,
629 int category_index[128] = { 0 };
630 int category[128] = { 0 };
631 int quant_index_table[102];
637 categorize(q, p, quant_index_table, category, category_index);
658 int gain_index,
int gain_index_next)
662 fc1 =
pow2tab[gain_index + 63];
664 if (gain_index == gain_index_next) {
668 fc2 = q->
gain_table[11 + (gain_index_next - gain_index)];
685 cook_gains *gains_ptr,
float *previous_buffer)
697 inbuffer[i] = inbuffer[i] * fc * q->
mlt_window[i] -
713 cook_gains *gains_ptr,
float *previous_buffer)
722 q->
imlt_window(q, buffer1, gains_ptr, previous_buffer);
725 for (i = 0; i < 8; i++)
726 if (gains_ptr->
now[i] || gains_ptr->
now[i + 1])
728 gains_ptr->
now[i], gains_ptr->
now[i + 1]);
731 memcpy(previous_buffer, buffer0,
748 int length = end - start + 1;
754 for (i = 0; i <
length; i++)
759 for (i = 0; i <
length; i++) {
765 decouple_tab[start + i] =
v;
785 float *decode_buffer,
786 float *mlt_buffer1,
float *mlt_buffer2)
791 mlt_buffer1[SUBBAND_SIZE * subband + j] = f1 * decode_buffer[tmp_idx];
792 mlt_buffer2[SUBBAND_SIZE * subband + j] = f2 * decode_buffer[tmp_idx];
804 float *mlt_buffer_left,
float *mlt_buffer_right)
811 const float *cplscale;
816 memset(mlt_buffer_left, 0, 1024 *
sizeof(*mlt_buffer_left));
817 memset(mlt_buffer_right, 0, 1024 *
sizeof(*mlt_buffer_right));
825 mlt_buffer_left[i * 20 + j] = decode_buffer[i * 40 + j];
826 mlt_buffer_right[i * 20 + j] = decode_buffer[i * 40 + 20 + j];
835 idx -= decouple_tab[cpl_tmp];
837 f1 = cplscale[decouple_tab[cpl_tmp] + 1];
839 q->
decouple(q, p, i, f1, f2, decode_buffer,
840 mlt_buffer_left, mlt_buffer_right);
896 cook_gains *gains_ptr,
float *previous_buffer,
899 imlt_gain(q, decode_buffer, gains_ptr, previous_buffer);
914 const uint8_t *inbuffer,
float **outbuffer)
916 int sub_packet_size = p->
size;
938 outbuffer ? outbuffer[p->
ch_idx] : NULL);
944 outbuffer ? outbuffer[p->
ch_idx + 1] : NULL);
948 outbuffer ? outbuffer[p->
ch_idx + 1] : NULL);
956 int *got_frame_ptr,
AVPacket *avpkt)
960 int buf_size = avpkt->
size;
962 float **samples = NULL;
967 if (buf_size < avctx->block_align)
986 "frame subpacket size total > avctx->block_align!\n");
997 "subpacket[%i] size %i js %i %i block_align %i\n",
1025 #define PRINT(a, b) av_dlog(q->avctx, " %s = %d\n", a, b);
1057 unsigned int channel_mask = 0;
1058 int samples_per_frame = 0;
1063 if (extradata_size <= 0) {
1080 while (edata_ptr < edata_ptr_end) {
1083 if (extradata_size >= 8) {
1085 samples_per_frame = bytestream_get_be16(&edata_ptr);
1087 extradata_size -= 8;
1089 if (extradata_size >= 8) {
1090 bytestream_get_be32(&edata_ptr);
1098 extradata_size -= 8;
1151 if (extradata_size >= 4)
1275 dump_cook_context(q);
static void mlt_compensate_output(COOKContext *q, float *decode_buffer, cook_gains *gains_ptr, float *previous_buffer, float *out)
Final part of subpacket decoding: Apply modulated lapped transform, gain compensation, clip and convert to integer.
static av_cold void init_cplscales_table(COOKContext *q)
static const int cplband[51]
av_cold void ff_dsputil_init(DSPContext *c, AVCodecContext *avctx)
#define AVERROR_PATCHWELCOME
This structure describes decoded (raw) audio or video data.
void(* scalar_dequant)(struct cook *q, int index, int quant_index, int *subband_coef_index, int *subband_coef_sign, float *mlt_p)
static const uint16_t envelope_quant_index_huffcodes[13][24]
static unsigned int get_bits(GetBitContext *s, int n)
Read 1-25 bits.
float decode_buffer_1[1024]
static av_cold int init(AVCodecContext *avctx)
static const int kmax_tab[7]
static const int expbits_tab[8]
static void categorize(COOKContext *q, COOKSubpacket *p, const int *quant_index_table, int *category, int *category_index)
Calculate the category and category_index vector.
static const float *const cplscales[5]
uint8_t ** extended_data
pointers to the data planes/channels.
#define DECLARE_ALIGNED(n, t, v)
void av_log(void *avcl, int level, const char *fmt,...) av_printf_format(3
Send the specified message to the log if the level is less than or equal to the current av_log_level...
static av_cold void init_pow2table(void)
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
VLC envelope_quant_index[13]
int block_align
number of bytes per packet if constant and known or 0 Used by some WAV based audio codecs...
uint8_t * extradata
some codecs need / can use extradata like Huffman tables.
static const uint8_t *const ccpl_huffbits[5]
static const int vhsize_tab[7]
static const float quant_centroid_tab[7][14]
void(* imlt_window)(struct cook *q, float *buffer1, cook_gains *gains_ptr, float *previous_buffer)
static void imlt_gain(COOKContext *q, float *inbuffer, cook_gains *gains_ptr, float *previous_buffer)
The modulated lapped transform, this takes transform coefficients and transforms them into timedomain...
static av_cold void init_gain_table(COOKContext *q)
void void avpriv_request_sample(void *avc, const char *msg,...) av_printf_format(2
Log a generic warning message about a missing feature.
enum AVSampleFormat sample_fmt
audio sample format
uint8_t * decoded_bytes_buffer
static int decouple_info(COOKContext *q, COOKSubpacket *p, int *decouple_tab)
function for getting the jointstereo coupling information
float mono_previous_buffer1[1024]
static void decode_vectors(COOKContext *q, COOKSubpacket *p, int *category, int *quant_index_table, float *mlt_buffer)
Fill the mlt_buffer with mlt coefficients.
static unsigned int av_lfg_get(AVLFG *c)
Get the next random unsigned 32-bit number using an ALFG.
static const uint8_t offset[511][2]
static void expand_category(COOKContext *q, int *category, int *category_index)
Expand the category vector.
static av_cold int end(AVCodecContext *avctx)
const char * name
Name of the codec implementation.
#define CODEC_CAP_DR1
Codec uses get_buffer() for allocating buffers and supports custom allocators.
static void interpolate(float *out, float v1, float v2, int size)
static int get_bits_count(const GetBitContext *s)
bitstream reader API header.
const float * cplscales[5]
#define FFSWAP(type, a, b)
static int decode_subpacket(COOKContext *q, COOKSubpacket *p, const uint8_t *inbuffer, float **outbuffer)
Cook subpacket decoding.
static int cook_decode_frame(AVCodecContext *avctx, void *data, int *got_frame_ptr, AVPacket *avpkt)
static void decode_bytes_and_gain(COOKContext *q, COOKSubpacket *p, const uint8_t *inbuffer, cook_gains *gains_ptr)
First part of subpacket decoding: decode raw stream bytes and read gain info.
#define DECODE_BYTES_PAD1(bytes)
#define AV_LOG_ERROR
Something went wrong and cannot losslessly be recovered.
static const int vd_tab[7]
void av_free(void *ptr)
Free a memory block which has been allocated with av_malloc(z)() or av_realloc(). ...
static const float dither_tab[9]
#define NULL_IF_CONFIG_SMALL(x)
Return NULL if CONFIG_SMALL is true, otherwise the argument without modification. ...
#define AV_LOG_DEBUG
Stuff which is only useful for libav* developers.
int av_get_channel_layout_nb_channels(uint64_t channel_layout)
Return the number of channels in the channel layout.
static const uint16_t *const ccpl_huffcodes[5]
float mono_previous_buffer2[1024]
static int decode_envelope(COOKContext *q, COOKSubpacket *p, int *quant_index_table)
Create the quant index table needed for the envelope.
Libavcodec external API header.
uint64_t channel_layout
Audio channel layout.
static void saturate_output_float(COOKContext *q, float *out)
Saturate the output signal and interleave.
void(* imdct_calc)(struct FFTContext *s, FFTSample *output, const FFTSample *input)
static const int vhvlcsize_tab[7]
#define FF_INPUT_BUFFER_PADDING_SIZE
Required number of additionally allocated bytes at the end of the input bitstream for decoding...
#define FF_ARRAY_ELEMS(a)
static int unpack_SQVH(COOKContext *q, COOKSubpacket *p, int category, int *subband_coef_index, int *subband_coef_sign)
Unpack the subband_coef_index and subband_coef_sign vectors.
int bit_rate
the average bitrate
static av_cold int init_cook_mlt(COOKContext *q)
static int mono_decode(COOKContext *q, COOKSubpacket *p, float *mlt_buffer)
void * av_malloc(size_t size) av_malloc_attrib 1(1)
Allocate a block of size bytes with alignment suitable for all memory accesses (including vectors if ...
static av_cold int cook_decode_init(AVCodecContext *avctx)
Cook initialization.
static av_always_inline int get_vlc2(GetBitContext *s, VLC_TYPE(*table)[2], int bits, int max_depth)
Parse a vlc code.
static const uint16_t *const cvh_huffcodes[7]
void ff_sine_window_init(float *window, int n)
Generate a sine window.
AVSampleFormat
Audio Sample Formats.
static void interpolate_float(COOKContext *q, float *buffer, int gain_index, int gain_index_next)
the actual requantization of the timedomain samples
void av_lfg_init(AVLFG *c, unsigned int seed)
void(* interpolate)(struct cook *q, float *buffer, int gain_index, int gain_index_next)
typedef void(RENAME(mix_any_func_type))
static av_cold int init_cook_vlc_tables(COOKContext *q)
int sample_rate
samples per second
void(* vector_clipf)(float *dst, const float *src, float min, float max, int len)
main external API structure.
static void close(AVCodecParserContext *s)
float mono_mdct_output[2048]
int ff_get_buffer(AVCodecContext *avctx, AVFrame *frame, int flags)
Get a buffer for a frame.
#define init_vlc(vlc, nb_bits, nb_codes,bits, bits_wrap, bits_size,codes, codes_wrap, codes_size,flags)
static unsigned int get_bits1(GetBitContext *s)
void(* decouple)(struct cook *q, COOKSubpacket *p, int subband, float f1, float f2, float *decode_buffer, float *mlt_buffer1, float *mlt_buffer2)
static int joint_decode(COOKContext *q, COOKSubpacket *p, float *mlt_buffer_left, float *mlt_buffer_right)
function for decoding joint stereo data
static av_cold int cook_decode_close(AVCodecContext *avctx)
static float pow2tab[127]
static int init_get_bits(GetBitContext *s, const uint8_t *buffer, int bit_size)
Initialize GetBitContext.
float decode_buffer_0[1060]
COOKSubpacket subpacket[MAX_SUBPACKETS]
float decode_buffer_2[1024]
static float rootpow2tab[127]
static const uint8_t envelope_quant_index_huffbits[13][24]
static const uint8_t *const cvh_huffbits[7]
static int decode_bytes(const uint8_t *inbuffer, uint8_t *out, int bytes)
Cook indata decoding, every 32 bits are XORed with 0x37c511f2.
static void scalar_dequant_float(COOKContext *q, int index, int quant_index, int *subband_coef_index, int *subband_coef_sign, float *mlt_p)
The real requantization of the mltcoefs.
common internal api header.
#define AVERROR_INVALIDDATA
static const int invradix_tab[7]
uint8_t pi<< 24) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0f/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S16, int16_t,(*(constint16_t *) pi >>8)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0f/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S32, int32_t,(*(constint32_t *) pi >>24)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0f/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_FLT, float, av_clip_uint8(lrintf(*(constfloat *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_FLT, float, av_clip_int16(lrintf(*(constfloat *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_FLT, float, av_clipl_int32(llrintf(*(constfloat *) pi *(1U<< 31)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_DBL, double, av_clip_uint8(lrint(*(constdouble *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_DBL, double, av_clip_int16(lrint(*(constdouble *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_DBL, double, av_clipl_int32(llrint(*(constdouble *) pi *(1U<< 31))))#defineSET_CONV_FUNC_GROUP(ofmt, ifmt) staticvoidset_generic_function(AudioConvert *ac){}voidff_audio_convert_free(AudioConvert **ac){return;ff_dither_free(&(*ac) ->dc);av_freep(ac);}AudioConvert *ff_audio_convert_alloc(AVAudioResampleContext *avr, enumAVSampleFormatout_fmt, enumAVSampleFormatin_fmt, intchannels, intsample_rate, intapply_map){AudioConvert *ac;intin_planar, out_planar;ac=av_mallocz(sizeof(*ac));returnNULL;ac->avr=avr;ac->out_fmt=out_fmt;ac->in_fmt=in_fmt;ac->channels=channels;ac->apply_map=apply_map;if(avr->dither_method!=AV_RESAMPLE_DITHER_NONE &&av_get_packed_sample_fmt(out_fmt)==AV_SAMPLE_FMT_S16 &&av_get_bytes_per_sample(in_fmt)>2){ac->dc=ff_dither_alloc(avr, out_fmt, in_fmt, channels, sample_rate, apply_map);if(!ac->dc){av_free(ac);returnNULL;}returnac;}in_planar=av_sample_fmt_is_planar(in_fmt);out_planar=av_sample_fmt_is_planar(out_fmt);if(in_planar==out_planar){ac->func_type=CONV_FUNC_TYPE_FLAT;ac->planes=in_planar?ac->channels:1;}elseif(in_planar) ac->func_type=CONV_FUNC_TYPE_INTERLEAVE;elseac->func_type=CONV_FUNC_TYPE_DEINTERLEAVE;set_generic_function(ac);ff_audio_convert_init_arm(ac);ff_audio_convert_init_x86(ac);returnac;}intff_audio_convert(AudioConvert *ac, AudioData *out, AudioData *in){intuse_generic=1;intlen=in->nb_samples;intp;if(ac->dc){av_dlog(ac->avr,"%dsamples-audio_convert:%sto%s(dithered)\n", len, av_get_sample_fmt_name(ac->in_fmt), av_get_sample_fmt_name(ac->out_fmt));returnff_convert_dither(ac-> out
int channels
number of audio channels
VLC_TYPE(* table)[2]
code, bits
static struct twinvq_data tab
static enum AVSampleFormat sample_fmts[]
static void decode_gain_info(GetBitContext *gb, int *gaininfo)
Fill the gain array for the timedomain quantization.
static int decode(AVCodecContext *avctx, void *data, int *got_frame, AVPacket *avpkt)
static void imlt_window_float(COOKContext *q, float *inbuffer, cook_gains *gains_ptr, float *previous_buffer)
Apply transform window, overlap buffers.
static void decouple_float(COOKContext *q, COOKSubpacket *p, int subband, float f1, float f2, float *decode_buffer, float *mlt_buffer1, float *mlt_buffer2)
function decouples a pair of signals from a single signal via multiplication.
static const int vpr_tab[7]
#define AV_CH_LAYOUT_MONO
This structure stores compressed data.
void ff_free_vlc(VLC *vlc)
int nb_samples
number of audio samples (per channel) described by this frame
unsigned int channel_mask
void * av_mallocz(size_t size) av_malloc_attrib 1(1)
Allocate a block of size bytes with alignment suitable for all memory accesses (including vectors if ...
Cook AKA RealAudio G2 compatible decoderdata.
void(* saturate_output)(struct cook *q, float *out)