37 uint64_t accu = 0,
round;
38 uint64_t accu0 = 0, accu1 = 0, accu2 = 0, accu3 = 0;
43 for (i = 0; i < n; i += 2) {
44 accu0 += (
int64_t)x[i + 0][0] * x[i + 0][0];
45 accu1 += (
int64_t)x[i + 0][1] * x[i + 0][1];
46 accu2 += (
int64_t)x[i + 1][0] * x[i + 1][0];
47 accu3 += (
int64_t)x[i + 1][1] * x[i + 1][1];
48 if ((accu0|accu1|accu2|accu3) > UINT64_MAX - INT32_MIN*(
int64_t)INT32_MIN || i+2>=n) {
53 while ((accu0|accu1|accu2|accu3) > (UINT64_MAX - accu) >> 2) {
61 accu += accu0 + accu1 + accu2 + accu3;
62 accu0 = accu1 = accu2 = accu3 = 0;
71 while (u < 0x80000000U) {
78 round = 1ULL << (nz-1);
79 u = ((accu +
round) >> nz);
89 for (i = 1; i < 64; i += 2)
90 x[i] = -(
unsigned)x[
i];
98 for (k = 1; k < 32; k++) {
99 z[64+2*k ] = -z[64 - k];
100 z[64+2*k+1] = z[ k + 1];
107 for (k = 0; k < 32; k++) {
116 for (i = 0; i < 32; i++) {
117 v[
i] = (
int)(0x10U + src[63 - 2*i ]) >> 5;
118 v[63 -
i] = (
int)(0x10U - src[63 - 2*i - 1]) >> 5;
126 int i = (
int)(accu >> 32);
131 while (
FFABS(i) < 0x40000000) {
138 round = 1
U << (nz-1);
139 mant = (
int)((accu + round) >> nz);
140 mant = (mant + 0x40LL)>>7;
150 int64_t accu_re = 0, accu_im = 0;
153 for (i = 1; i < 38; i++) {
154 accu_re += (uint64_t)x[i][0] * x[i+lag][0];
155 accu_re += (uint64_t)x[i][1] * x[i+lag][1];
156 accu_im += (uint64_t)x[i][0] * x[i+lag][1];
157 accu_im -= (uint64_t)x[i][1] * x[i+lag][0];
163 accu_re += (uint64_t)x[ 0][0] * x[lag][0];
164 accu_re += (uint64_t)x[ 0][1] * x[lag][1];
165 accu_im += (uint64_t)x[ 0][0] * x[lag][1];
166 accu_im -= (uint64_t)x[ 0][1] * x[lag][0];
174 accu_re += (uint64_t)x[38][0] * x[39][0];
175 accu_re += (uint64_t)x[38][1] * x[39][1];
176 accu_im += (uint64_t)x[38][0] * x[39][1];
177 accu_im -= (uint64_t)x[38][1] * x[39][0];
183 for (i = 1; i < 38; i++) {
184 accu_re += (uint64_t)x[i][0] * x[i][0];
185 accu_re += (uint64_t)x[i][1] * x[i][1];
188 accu_re += (uint64_t)x[ 0][0] * x[ 0][0];
189 accu_re += (uint64_t)x[ 0][1] * x[ 0][1];
194 accu_re += (uint64_t)x[38][0] * x[38][0];
195 accu_re += (uint64_t)x[38][1] * x[38][1];
209 const int alpha0[2],
const int alpha1[2],
210 int bw,
int start,
int end)
216 accu = (
int64_t)alpha0[0] * bw;
217 alpha[2] = (
int)((accu + 0x40000000) >> 31);
218 accu = (
int64_t)alpha0[1] * bw;
219 alpha[3] = (
int)((accu + 0x40000000) >> 31);
221 bw = (
int)((accu + 0x40000000) >> 31);
222 accu = (
int64_t)alpha1[0] * bw;
223 alpha[0] = (
int)((accu + 0x40000000) >> 31);
224 accu = (
int64_t)alpha1[1] * bw;
225 alpha[1] = (
int)((accu + 0x40000000) >> 31);
227 for (i = start; i <
end; i++) {
228 accu = (
int64_t)X_low[i][0] * 0x20000000;
229 accu += (
int64_t)X_low[i - 2][0] * alpha[0];
230 accu -= (
int64_t)X_low[i - 2][1] * alpha[1];
231 accu += (
int64_t)X_low[i - 1][0] * alpha[2];
232 accu -= (
int64_t)X_low[i - 1][1] * alpha[3];
233 X_high[
i][0] = (
int)((accu + 0x10000000) >> 29);
235 accu = (
int64_t)X_low[i][1] * 0x20000000;
236 accu += (
int64_t)X_low[i - 2][1] * alpha[0];
237 accu += (
int64_t)X_low[i - 2][0] * alpha[1];
238 accu += (
int64_t)X_low[i - 1][1] * alpha[2];
239 accu += (
int64_t)X_low[i - 1][0] * alpha[3];
240 X_high[
i][1] = (
int)((accu + 0x10000000) >> 29);
245 const SoftFloat *g_filt,
int m_max, intptr_t ixh)
250 for (m = 0; m < m_max; m++) {
251 if (22 - g_filt[m].
exp < 61) {
253 accu = (
int64_t)X_high[m][ixh][0] * ((g_filt[m].mant + 0x40)>>7);
254 Y[m][0] = (
int)((accu + r) >> (23-g_filt[m].
exp));
256 accu = (
int64_t)X_high[m][ixh][1] * ((g_filt[m].mant + 0x40)>>7);
257 Y[m][1] = (
int)((accu + r) >> (23-g_filt[m].
exp));
272 for (m = 0; m < m_max; m++) {
273 unsigned y0 =
Y[m][0];
274 unsigned y1 =
Y[m][1];
275 noise = (noise + 1) & 0x1ff;
279 shift = 22 - s_m[m].
exp;
283 }
else if (shift < 30) {
284 round = 1 << (shift-1);
285 y0 += (s_m[m].
mant * phi_sign0 +
round) >> shift;
286 y1 += (s_m[m].
mant * phi_sign1 +
round) >> shift;
292 shift = 22 - q_filt[m].
exp;
296 }
else if (shift < 30) {
297 round = 1 << (shift-1);
299 accu = (
int64_t)q_filt[m].mant * ff_sbr_noise_table_fixed[noise][0];
300 tmp = (
int)((accu + 0x40000000) >> 31);
301 y0 += (tmp +
round) >> shift;
303 accu = (
int64_t)q_filt[m].mant * ff_sbr_noise_table_fixed[noise][1];
304 tmp = (
int)((accu + 0x40000000) >> 31);
305 y1 += (tmp +
round) >> shift;
310 phi_sign1 = -phi_sign1;
static int shift(int a, int b)
static av_always_inline SoftFloat autocorr_calc(int64_t accu)
static void sbr_autocorrelate_c(const int x[40][2], SoftFloat phi[3][2][2])
Macro definitions for various function/variable attributes.
static void sbr_hf_g_filt_c(int(*Y)[2], const int(*X_high)[40][2], const SoftFloat *g_filt, int m_max, intptr_t ixh)
static av_cold int end(AVCodecContext *avctx)
#define u(width, name, range_min, range_max)
static SoftFloat sbr_sum_square_c(int(*x)[2], int n)
#define i(width, name, range_min, range_max)
#define AV_LOG_ERROR
Something went wrong and cannot losslessly be recovered.
static av_always_inline void autocorrelate(const int x[40][2], SoftFloat phi[3][2][2], int lag)
static av_always_inline av_const double round(double x)
static av_always_inline int sbr_hf_apply_noise(int(*Y)[2], const SoftFloat *s_m, const SoftFloat *q_filt, int noise, int phi_sign0, int phi_sign1, int m_max)
static void sbr_qmf_post_shuffle_c(int W[32][2], const int *z)
#define FFABS(a)
Absolute value, Note, INT_MIN / INT64_MIN result in undefined behavior as they are not representable ...
AAC definitions and structures.
static void sbr_neg_odd_64_c(int *x)
static const int16_t alpha[]
static void sbr_qmf_pre_shuffle_c(int *z)
static int noise(AVBSFContext *ctx, AVPacket *pkt)
static void sbr_qmf_deint_neg_c(int *v, const int *src)
static av_const SoftFloat av_int2sf(int v, int frac_bits)
Converts a mantisse and exponent to a SoftFloat.
static void sbr_hf_gen_c(int(*X_high)[2], const int(*X_low)[2], const int alpha0[2], const int alpha1[2], int bw, int start, int end)