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 = 1U << (nz-1);
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;
233 X_high[
i][0] = (
int)((accu + 0x10000000) >> 29);
235 accu = (
int64_t)X_low[
i][1] * 0x20000000;
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];
283 }
else if (
shift < 30) {
296 }
else if (
shift < 30) {
299 accu = (
int64_t)q_filt[m].mant * ff_sbr_noise_table_fixed[
noise][0];
300 tmp = (
int)((accu + 0x40000000) >> 31);
303 accu = (
int64_t)q_filt[m].mant * ff_sbr_noise_table_fixed[
noise][1];
304 tmp = (
int)((accu + 0x40000000) >> 31);
310 phi_sign1 = -phi_sign1;
Macro definitions for various function/variable attributes.
#define u(width, name, range_min, range_max)
#define FFABS(a)
Absolute value, Note, INT_MIN / INT64_MIN result in undefined behavior as they are not representable ...
#define AV_LOG_ERROR
Something went wrong and cannot losslessly be recovered.
static const int16_t alpha[]
static av_always_inline av_const double round(double x)
static int noise(AVBSFContext *ctx, AVPacket *pkt)
static void sbr_neg_odd_64_c(int *x)
static void sbr_qmf_deint_neg_c(int *v, const int *src)
static av_always_inline SoftFloat autocorr_calc(int64_t accu)
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)
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_autocorrelate_c(const int x[40][2], SoftFloat phi[3][2][2])
static void sbr_qmf_pre_shuffle_c(int *z)
static void sbr_qmf_post_shuffle_c(int W[32][2], const int *z)
static av_always_inline void autocorrelate(const int x[40][2], SoftFloat phi[3][2][2], int lag)
static SoftFloat sbr_sum_square_c(int(*x)[2], int n)
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_const SoftFloat av_int2sf(int v, int frac_bits)
Converts a mantisse and exponent to a SoftFloat.
static int shift(int a, int b)