31 for (k = 0; k < 64; k++) {
32 float f = z[k] + z[k + 64] + z[k + 128] + z[k + 192] + z[k + 256];
39 float sum0 = 0.0f, sum1 = 0.0f;
42 for (i = 0; i < n; i += 2)
44 sum0 += x[i + 0][0] * x[i + 0][0];
45 sum1 += x[i + 0][1] * x[i + 0][1];
46 sum0 += x[i + 1][0] * x[i + 1][0];
47 sum1 += x[i + 1][1] * x[i + 1][1];
57 for (i = 1; i < 64; i += 4) {
58 xi[i + 0].
i ^= 1
U << 31;
59 xi[i + 2].
i ^= 1
U << 31;
69 for (k = 1; k < 31; k += 2) {
70 zi[64 + 2 * k + 0].
i = zi[64 - k].
i ^ (1
U << 31);
71 zi[64 + 2 * k + 1].
i = zi[ k + 1].
i;
72 zi[64 + 2 * k + 2].
i = zi[63 - k].
i ^ (1
U << 31);
73 zi[64 + 2 * k + 3].
i = zi[ k + 2].
i;
75 zi[64 + 2 * 31 + 0].
i = zi[64 - 31].
i ^ (1
U << 31);
76 zi[64 + 2 * 31 + 1].
i = zi[31 + 1].
i;
84 for (k = 0; k < 32; k += 2) {
85 Wi[2 * k + 0].
i = zi[63 - k].
i ^ (1
U << 31);
86 Wi[2 * k + 1].
i = zi[ k + 0].
i;
87 Wi[2 * k + 2].
i = zi[62 - k].
i ^ (1
U << 31);
88 Wi[2 * k + 3].
i = zi[ k + 1].
i;
97 for (i = 0; i < 32; i++) {
98 vi[
i].
i = si[63 - 2 *
i ].
i;
99 vi[63 -
i].
i = si[63 - 2 * i - 1].
i ^ (1
U << 31);
106 for (i = 0; i < 64; i++) {
107 v[
i] = src0[
i] - src1[63 -
i];
108 v[127 -
i] = src0[
i] + src1[63 -
i];
113 float phi[3][2][2],
int lag)
116 float real_sum = 0.0f;
117 float imag_sum = 0.0f;
119 for (i = 1; i < 38; i++) {
120 real_sum += x[
i][0] * x[i+lag][0] + x[
i][1] * x[i+lag][1];
121 imag_sum += x[
i][0] * x[i+lag][1] - x[
i][1] * x[i+lag][0];
123 phi[2-lag][1][0] = real_sum + x[ 0][0] * x[lag][0] + x[ 0][1] * x[lag][1];
124 phi[2-lag][1][1] = imag_sum + x[ 0][0] * x[lag][1] - x[ 0][1] * x[lag][0];
126 phi[0][0][0] = real_sum + x[38][0] * x[39][0] + x[38][1] * x[39][1];
127 phi[0][0][1] = imag_sum + x[38][0] * x[39][1] - x[38][1] * x[39][0];
130 for (i = 1; i < 38; i++) {
131 real_sum += x[
i][0] * x[
i][0] + x[
i][1] * x[
i][1];
133 phi[2][1][0] = real_sum + x[ 0][0] * x[ 0][0] + x[ 0][1] * x[ 0][1];
134 phi[1][0][0] = real_sum + x[38][0] * x[38][0] + x[38][1] * x[38][1];
148 float real_sum2 = x[0][0] * x[2][0] + x[0][1] * x[2][1];
149 float imag_sum2 = x[0][0] * x[2][1] - x[0][1] * x[2][0];
150 float real_sum1 = 0.0f, imag_sum1 = 0.0f, real_sum0 = 0.0f;
152 for (i = 1; i < 38; i++) {
153 real_sum0 += x[
i][0] * x[
i ][0] + x[
i][1] * x[
i ][1];
154 real_sum1 += x[
i][0] * x[i + 1][0] + x[
i][1] * x[i + 1][1];
155 imag_sum1 += x[
i][0] * x[i + 1][1] - x[
i][1] * x[i + 1][0];
156 real_sum2 += x[
i][0] * x[i + 2][0] + x[
i][1] * x[i + 2][1];
157 imag_sum2 += x[
i][0] * x[i + 2][1] - x[
i][1] * x[i + 2][0];
159 phi[2 - 2][1][0] = real_sum2;
160 phi[2 - 2][1][1] = imag_sum2;
161 phi[2 ][1][0] = real_sum0 + x[ 0][0] * x[ 0][0] + x[ 0][1] * x[ 0][1];
162 phi[1 ][0][0] = real_sum0 + x[38][0] * x[38][0] + x[38][1] * x[38][1];
163 phi[2 - 1][1][0] = real_sum1 + x[ 0][0] * x[ 1][0] + x[ 0][1] * x[ 1][1];
164 phi[2 - 1][1][1] = imag_sum1 + x[ 0][0] * x[ 1][1] - x[ 0][1] * x[ 1][0];
165 phi[0 ][0][0] = real_sum1 + x[38][0] * x[39][0] + x[38][1] * x[39][1];
166 phi[0 ][0][1] = imag_sum1 + x[38][0] * x[39][1] - x[38][1] * x[39][0];
171 const float alpha0[2],
const float alpha1[2],
172 float bw,
int start,
int end)
177 alpha[0] = alpha1[0] * bw * bw;
178 alpha[1] = alpha1[1] * bw * bw;
179 alpha[2] = alpha0[0] * bw;
180 alpha[3] = alpha0[1] * bw;
182 for (i = start; i < end; i++) {
184 X_low[i - 2][0] * alpha[0] -
185 X_low[i - 2][1] * alpha[1] +
186 X_low[i - 1][0] * alpha[2] -
187 X_low[i - 1][1] * alpha[3] +
190 X_low[i - 2][1] * alpha[0] +
191 X_low[i - 2][0] * alpha[1] +
192 X_low[i - 1][1] * alpha[2] +
193 X_low[i - 1][0] * alpha[3] +
199 const float *g_filt,
int m_max, intptr_t ixh)
203 for (m = 0; m < m_max; m++) {
204 Y[m][0] = X_high[m][ixh][0] * g_filt[m];
205 Y[m][1] = X_high[m][ixh][1] * g_filt[m];
219 for (m = 0; m < m_max; m++) {
222 noise = (noise + 1) & 0x1ff;
224 y0 += s_m[m] * phi_sign0;
225 y1 += s_m[m] * phi_sign1;
228 y1 += q_filt[m] * ff_sbr_noise_table[
noise][1];
232 phi_sign1 = -phi_sign1;
237 const float *q_filt,
int noise,
244 const float *q_filt,
int noise,
247 float phi_sign = 1 - 2 * (kx & 1);
252 const float *q_filt,
int noise,
259 const float *q_filt,
int noise,
262 float phi_sign = 1 - 2 * (kx & 1);