50 #define BITSTREAM_WRITER_LE
53 #define POW_TABLE_SIZE (1<<11)
54 #define POW_TABLE_OFFSET 3
55 #define OPT_SIZE ((1<<15) + 3000)
73 0, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 4, 4,
74 5, 5, 5, 6, 7, 7, 8, 8, 9, 10, 11, 11, 12, 13, 13, 14,
75 15, 15, 16, 17, 17, 18, 19, 19, 20, 21, 22, 22, 23, 24, 25, 26,
76 27, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 37, 38, 39, 40,
77 41, 41, 42, 43, 44, 45, 45, 46, 47, 48, 49, 50, 51, 52, 52, 53,
78 54, 55, 55, 56, 57, 57, 58, 59, 59, 60, 60, 60, 61, 61, 61, 62,
82 0, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 4, 4,
83 4, 5, 5, 5, 6, 6, 7, 7, 8, 8, 9, 10, 10, 11, 11, 12,
84 13, 13, 14, 15, 16, 17, 17, 18, 19, 19, 20, 21, 21, 22, 22, 23,
85 23, 24, 24, 25, 25, 25, 26, 26, 26, 26, 27, 27, 27, 27, 27, 28,
86 28, 28, 28, 28, 28, 29, 29, 29, 29, 29, 29, 29, 29, 30,
96 0, 1, 1, 2, 2, 3, 3, 4, 5, 6, 7, 8, 9, 10, 11, 11,
99 0, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 8,
100 8, 9, 10, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22,
101 22, 23, 23, 24, 24, 25, 25, 26, 26, 27, 27, 28, 28, 29, 29, 29,
104 0, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3,
105 4, 4, 4, 5, 5, 5, 6, 6, 7, 7, 7, 8, 8, 9, 9, 9,
106 10, 10, 11, 11, 11, 12, 12, 13, 13, 13, 13, 14, 14, 14, 15, 15,
107 15, 15, 16, 16, 16, 17, 17, 17, 18, 18, 18, 19, 19, 20, 20, 20,
108 21, 21, 22, 22, 23, 23, 24, 25, 26, 26, 27, 28, 29, 30, 31, 32,
109 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 42, 43, 44, 44, 45, 45,
110 46, 47, 47, 48, 48, 49, 49, 50, 50, 50, 51, 51, 51, 52, 52, 52,
111 53, 53, 53, 54, 54, 54, 55, 55, 55, 56, 56, 56, 57, 57, 57, 57,
112 58, 58, 58, 58, 59, 59, 59, 59, 60, 60, 60, 60, 60, 61, 61, 61,
116 static const float quant_lut_mul[7] = { 0.0, 0.0, 2.0, 2.0, 5.0, 12.0, 36.6 };
117 static const float quant_lut_add[7] = { 0.0, 0.0, 2.0, 7.0, 21.0, 56.0, 157.0 };
128 s->mdct_ctx.mdct_calc(&
s->mdct_ctx,
s->mdct_out,
s->in_buff);
141 if (
s->avctx->trellis) {
158 for (
int i = 1;
i < 513;
i++) {
184 av_log(avctx,
AV_LOG_ERROR,
"Nellymoser works only with 8000, 16000, 11025, 22050 and 44100 sample rate\n");
198 if (
s->avctx->trellis) {
201 if (!
s->opt || !
s->path)
210 #define find_best(val, table, LUT, LUT_add, LUT_size) \
212 LUT[av_clip ((lrintf(val) >> 8) + LUT_add, 0, LUT_size - 1)]; \
213 if (fabs(val - table[best_idx]) > fabs(val - table[best_idx + 1])) \
218 int band, best_idx, power_idx = 0;
219 float power_candidate;
223 idx_table[0] = best_idx;
227 power_candidate = cand[band] - power_idx;
229 idx_table[band] = best_idx;
234 static inline float distance(
float x,
float y,
int band)
243 int i, j, band, best_idx;
244 float power_candidate, best_val;
253 for (
i = 0;
i < 64;
i++) {
261 int idx_min, idx_max, idx;
262 power_candidate = cand[band];
263 for (q = 1000; !
c && q <
OPT_SIZE; q <<= 2) {
264 idx_min =
FFMAX(0, cand[band] - q);
267 if (
isinf(opt[band - 1][
i]) )
269 for (j = 0; j < 32; j++) {
273 if (idx >= idx_min) {
274 tmp = opt[band - 1][
i] +
distance(idx, power_candidate, band);
275 if (opt[band][idx] >
tmp) {
276 opt[band][idx] =
tmp;
291 if (best_val > opt[band][
i]) {
292 best_val = opt[band][
i];
297 idx_table[band] = path[band][best_idx];
313 int i, j, band,
block, best_idx, power_idx = 0;
314 float power_val,
coeff, coeff_sum;
327 coeff_sum +=
s->mdct_out[
i ] *
s->mdct_out[
i ]
334 if (
s->avctx->trellis) {
351 s->mdct_out[
i] *= power_val;
423 .
name =
"nellymoser",
static enum AVSampleFormat sample_fmts[]
av_cold void ff_af_queue_init(AVCodecContext *avctx, AudioFrameQueue *afq)
Initialize AudioFrameQueue.
void ff_af_queue_close(AudioFrameQueue *afq)
Close AudioFrameQueue.
void ff_af_queue_remove(AudioFrameQueue *afq, int nb_samples, int64_t *pts, int64_t *duration)
Remove frame(s) from the queue.
int ff_af_queue_add(AudioFrameQueue *afq, const AVFrame *f)
Add a frame to the queue.
#define av_assert1(cond)
assert() equivalent, that does not lie in speed critical code.
Libavcodec external API header.
#define FF_COMPLIANCE_NORMAL
static av_cold int init(AVCodecContext *avctx)
common internal and external API header
static __device__ float fabs(float a)
int ff_alloc_packet2(AVCodecContext *avctx, AVPacket *avpkt, int64_t size, int64_t min_size)
Check AVPacket size and/or allocate data.
#define AV_CODEC_FLAG_BITEXACT
Use only bitexact stuff (except (I)DCT).
#define AV_CODEC_CAP_DELAY
Encoder or decoder requires flushing with NULL input at the end in order to give the complete and cor...
#define AV_CODEC_CAP_SMALL_LAST_FRAME
Codec can be fed a final frame with a smaller size.
#define AV_LOG_ERROR
Something went wrong and cannot losslessly be recovered.
#define DECLARE_ALIGNED(n, t, v)
Declare a variable that is aligned in memory.
AVSampleFormat
Audio sample formats.
static void put_bits(Jpeg2000EncoderContext *s, int val, int n)
put n times val bit
#define FF_CODEC_CAP_INIT_THREADSAFE
The codec does not modify any global variables in the init function, allowing to call the init functi...
#define FF_CODEC_CAP_INIT_CLEANUP
The codec allows calling the close function for deallocation even if the init function returned a fai...
av_cold AVFloatDSPContext * avpriv_float_dsp_alloc(int bit_exact)
Allocate a float DSP context.
common internal API header
#define NULL_IF_CONFIG_SMALL(x)
Return NULL if CONFIG_SMALL is true, otherwise the argument without modification.
static int ff_thread_once(char *control, void(*routine)(void))
const uint16_t ff_nelly_init_table[64]
const uint8_t ff_nelly_band_sizes_table[NELLY_BANDS]
void ff_nelly_get_sample_bits(const float *buf, int *bits)
const int16_t ff_nelly_delta_table[32]
const float ff_nelly_dequantization_table[127]
The 3 alphanumeric copyright notices are md5summed they are from the original implementors.
#define NELLY_DETAIL_BITS
#define NELLY_HEADER_BITS
static float distance(float x, float y, int band)
#define find_best(val, table, LUT, LUT_add, LUT_size)
static av_cold void nellymoser_init_static(void)
static void encode_block(NellyMoserEncodeContext *s, unsigned char *output, int output_size)
Encode NELLY_SAMPLES samples.
static const uint8_t quant_lut_offset[8]
AVCodec ff_nellymoser_encoder
static av_cold int encode_init(AVCodecContext *avctx)
static float pow_table[POW_TABLE_SIZE]
pow(2, -i / 2048.0 - 3.0);
static const float quant_lut_mul[7]
static const float quant_lut_add[7]
static void get_exponent_greedy(NellyMoserEncodeContext *s, float *cand, int *idx_table)
static av_cold int encode_end(AVCodecContext *avctx)
static void get_exponent_dynamic(NellyMoserEncodeContext *s, float *cand, int *idx_table)
static void apply_mdct(NellyMoserEncodeContext *s)
static const uint8_t sf_lut[96]
static const uint8_t sf_delta_lut[78]
static const uint8_t quant_lut[230]
static int encode_frame(AVCodecContext *avctx, AVPacket *avpkt, const AVFrame *frame, int *got_packet_ptr)
static const uint16_t table[]
static void init_put_bits(PutBitContext *s, uint8_t *buffer, int buffer_size)
Initialize the PutBitContext s.
static uint8_t * put_bits_ptr(PutBitContext *s)
Return the pointer to the byte where the bitstream writer will put the next bit.
static int put_bits_count(PutBitContext *s)
static void flush_put_bits(PutBitContext *s)
Pad the end of the output stream with zeros.
void ff_init_ff_sine_windows(int index)
initialize the specified entry of ff_sine_windows
main external API structure.
int strict_std_compliance
strictly follow the standard (MPEG-4, ...).
int initial_padding
Audio only.
int sample_rate
samples per second
int flags
AV_CODEC_FLAG_*.
int channels
number of audio channels
int frame_size
Number of samples per channel in an audio frame.
const char * name
Name of the codec implementation.
This structure describes decoded (raw) audio or video data.
int nb_samples
number of audio samples (per channel) described by this frame
uint8_t * data[AV_NUM_DATA_POINTERS]
pointer to the picture/channel planes.
This structure stores compressed data.
int64_t duration
Duration of this packet in AVStream->time_base units, 0 if unknown.
int64_t pts
Presentation timestamp in AVStream->time_base units; the time at which the decompressed packet will b...
float in_buff[NELLY_SAMPLES]
float buf[3 *NELLY_BUF_LEN]
sample buffer
uint8_t(* path)[OPT_SIZE]
float mdct_out[NELLY_SAMPLES]
static const double coeff[2][5]