2 * AAC encoder utilities
3 * Copyright (C) 2015 Rostislav Pehlivanov
5 * This file is part of FFmpeg.
7 * FFmpeg is free software; you can redistribute it and/or
8 * modify it under the terms of the GNU Lesser General Public
9 * License as published by the Free Software Foundation; either
10 * version 2.1 of the License, or (at your option) any later version.
12 * FFmpeg is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 * Lesser General Public License for more details.
17 * You should have received a copy of the GNU Lesser General Public
18 * License along with FFmpeg; if not, write to the Free Software
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
24 * AAC encoder utilities
25 * @author Rostislav Pehlivanov ( atomnuker gmail com )
28 #ifndef AVCODEC_AACENC_UTILS_H
29 #define AVCODEC_AACENC_UTILS_H
32 #include "aac_tablegen_decl.h"
33 #include "aacenctab.h"
35 #define ROUND_STANDARD 0.4054f
36 #define ROUND_TO_ZERO 0.1054f
37 #define C_QUANT 0.4054f
39 static inline void abs_pow34_v(float *out, const float *in, const int size)
42 for (i = 0; i < size; i++) {
43 float a = fabsf(in[i]);
44 out[i] = sqrtf(a * sqrtf(a));
49 * Quantize one coefficient.
50 * @return absolute value of the quantized coefficient
51 * @see 3GPP TS26.403 5.6.2 "Scalefactor determination"
53 static inline int quant(float coef, const float Q, const float rounding)
56 return sqrtf(a * sqrtf(a)) + rounding;
59 static inline void quantize_bands(int *out, const float *in, const float *scaled,
60 int size, float Q34, int is_signed, int maxval,
65 for (i = 0; i < size; i++) {
67 out[i] = (int)FFMIN(qc + rounding, (double)maxval);
68 if (is_signed && in[i] < 0.0f) {
74 static inline float find_max_val(int group_len, int swb_size, const float *scaled)
78 for (w2 = 0; w2 < group_len; w2++) {
79 for (i = 0; i < swb_size; i++) {
80 maxval = FFMAX(maxval, scaled[w2*128+i]);
86 static inline int find_min_book(float maxval, int sf)
88 float Q = ff_aac_pow2sf_tab[POW_SF2_ZERO - sf + SCALE_ONE_POS - SCALE_DIV_512];
89 float Q34 = sqrtf(Q * sqrtf(Q));
91 qmaxval = maxval * Q34 + C_QUANT;
92 if (qmaxval == 0) cb = 0;
93 else if (qmaxval == 1) cb = 1;
94 else if (qmaxval == 2) cb = 3;
95 else if (qmaxval <= 4) cb = 5;
96 else if (qmaxval <= 7) cb = 7;
97 else if (qmaxval <= 12) cb = 9;
102 /** Return the minimum scalefactor where the quantized coef does not clip. */
103 static inline uint8_t coef2minsf(float coef)
105 return av_clip_uint8(log2f(coef)*4 - 69 + SCALE_ONE_POS - SCALE_DIV_512);
108 /** Return the maximum scalefactor where the quantized coef is not zero. */
109 static inline uint8_t coef2maxsf(float coef)
111 return av_clip_uint8(log2f(coef)*4 + 6 + SCALE_ONE_POS - SCALE_DIV_512);
115 * Returns the closest possible index to an array of float values, given a value.
117 static inline int quant_array_idx(const float val, const float *arr, const int num)
120 float quant_min_err = INFINITY;
121 for (i = 0; i < num; i++) {
122 float error = (val - arr[i])*(val - arr[i]);
123 if (error < quant_min_err) {
124 quant_min_err = error;
132 * linear congruential pseudorandom number generator, copied from the decoder
134 static inline int lcg_random(unsigned previous_val)
136 union { unsigned u; int s; } v = { previous_val * 1664525u + 1013904223 };
140 #define ERROR_IF(cond, ...) \
142 av_log(avctx, AV_LOG_ERROR, __VA_ARGS__); \
143 return AVERROR(EINVAL); \
146 #define WARN_IF(cond, ...) \
148 av_log(avctx, AV_LOG_WARNING, __VA_ARGS__); \
152 #endif /* AVCODEC_AACENC_UTILS_H */