44#ifndef INCLUDED_volk_8ic_x2_s32f_multiply_conjugate_32fc_a_H
45#define INCLUDED_volk_8ic_x2_s32f_multiply_conjugate_32fc_a_H
55volk_8ic_x2_s32f_multiply_conjugate_32fc_a_avx2(
lv_32fc_t* cVector,
59 unsigned int num_points)
61 unsigned int number = 0;
62 const unsigned int oneEigthPoints = num_points / 8;
64 __m256i x, y, realz, imagz;
65 __m256 ret, retlo, rethi;
69 __m256i conjugateSign =
70 _mm256_set_epi16(-1, 1, -1, 1, -1, 1, -1, 1, -1, 1, -1, 1, -1, 1, -1, 1);
72 __m256 invScalar = _mm256_set1_ps(1.0 / scalar);
74 for (; number < oneEigthPoints; number++) {
80 realz = _mm256_madd_epi16(x, y);
83 y = _mm256_sign_epi16(y, conjugateSign);
86 y = _mm256_shufflehi_epi16(_mm256_shufflelo_epi16(y,
_MM_SHUFFLE(2, 3, 0, 1)),
90 imagz = _mm256_madd_epi16(x, y);
93 retlo = _mm256_cvtepi32_ps(_mm256_unpacklo_epi32(realz, imagz));
96 retlo = _mm256_mul_ps(retlo, invScalar);
99 rethi = _mm256_cvtepi32_ps(_mm256_unpackhi_epi32(realz, imagz));
102 rethi = _mm256_mul_ps(rethi, invScalar);
104 ret = _mm256_permute2f128_ps(retlo, rethi, 0b00100000);
105 _mm256_store_ps((
float*)c, ret);
108 ret = _mm256_permute2f128_ps(retlo, rethi, 0b00110001);
109 _mm256_store_ps((
float*)c, ret);
116 number = oneEigthPoints * 8;
117 float* cFloatPtr = (
float*)&cVector[number];
118 int8_t* a8Ptr = (int8_t*)&aVector[number];
119 int8_t* b8Ptr = (int8_t*)&bVector[number];
120 for (; number < num_points; number++) {
121 float aReal = (float)*a8Ptr++;
122 float aImag = (float)*a8Ptr++;
124 float bReal = (float)*b8Ptr++;
125 float bImag = (float)*b8Ptr++;
129 *cFloatPtr++ =
lv_creal(temp) / scalar;
130 *cFloatPtr++ =
lv_cimag(temp) / scalar;
137#include <smmintrin.h>
140volk_8ic_x2_s32f_multiply_conjugate_32fc_a_sse4_1(
lv_32fc_t* cVector,
144 unsigned int num_points)
146 unsigned int number = 0;
147 const unsigned int quarterPoints = num_points / 4;
158 for (; number < quarterPoints; number++) {
200 number = quarterPoints * 4;
201 float* cFloatPtr = (
float*)&cVector[number];
202 int8_t* a8Ptr = (int8_t*)&aVector[number];
203 int8_t* b8Ptr = (int8_t*)&bVector[number];
204 for (; number < num_points; number++) {
205 float aReal = (float)*a8Ptr++;
206 float aImag = (float)*a8Ptr++;
208 float bReal = (float)*b8Ptr++;
209 float bImag = (float)*b8Ptr++;
213 *cFloatPtr++ =
lv_creal(temp) / scalar;
214 *cFloatPtr++ =
lv_cimag(temp) / scalar;
220#ifdef LV_HAVE_GENERIC
227 unsigned int num_points)
229 unsigned int number = 0;
230 float* cPtr = (
float*)cVector;
231 const float invScalar = 1.0 / scalar;
232 int8_t* a8Ptr = (int8_t*)aVector;
233 int8_t* b8Ptr = (int8_t*)bVector;
234 for (number = 0; number < num_points; number++) {
235 float aReal = (float)*a8Ptr++;
236 float aImag = (float)*a8Ptr++;
238 float bReal = (float)*b8Ptr++;
239 float bImag = (float)*b8Ptr++;
243 *cPtr++ = (
lv_creal(temp) * invScalar);
244 *cPtr++ = (
lv_cimag(temp) * invScalar);
252#ifndef INCLUDED_volk_8ic_x2_s32f_multiply_conjugate_32fc_u_H
253#define INCLUDED_volk_8ic_x2_s32f_multiply_conjugate_32fc_u_H
260#include <immintrin.h>
263volk_8ic_x2_s32f_multiply_conjugate_32fc_u_avx2(
lv_32fc_t* cVector,
267 unsigned int num_points)
269 unsigned int number = 0;
270 const unsigned int oneEigthPoints = num_points / 8;
272 __m256i x, y, realz, imagz;
273 __m256 ret, retlo, rethi;
277 __m256i conjugateSign =
278 _mm256_set_epi16(-1, 1, -1, 1, -1, 1, -1, 1, -1, 1, -1, 1, -1, 1, -1, 1);
280 __m256 invScalar = _mm256_set1_ps(1.0 / scalar);
282 for (; number < oneEigthPoints; number++) {
288 realz = _mm256_madd_epi16(x, y);
291 y = _mm256_sign_epi16(y, conjugateSign);
294 y = _mm256_shufflehi_epi16(_mm256_shufflelo_epi16(y,
_MM_SHUFFLE(2, 3, 0, 1)),
298 imagz = _mm256_madd_epi16(x, y);
301 retlo = _mm256_cvtepi32_ps(_mm256_unpacklo_epi32(realz, imagz));
304 retlo = _mm256_mul_ps(retlo, invScalar);
307 rethi = _mm256_cvtepi32_ps(_mm256_unpackhi_epi32(realz, imagz));
310 rethi = _mm256_mul_ps(rethi, invScalar);
312 ret = _mm256_permute2f128_ps(retlo, rethi, 0b00100000);
313 _mm256_storeu_ps((
float*)c, ret);
316 ret = _mm256_permute2f128_ps(retlo, rethi, 0b00110001);
317 _mm256_storeu_ps((
float*)c, ret);
324 number = oneEigthPoints * 8;
325 float* cFloatPtr = (
float*)&cVector[number];
326 int8_t* a8Ptr = (int8_t*)&aVector[number];
327 int8_t* b8Ptr = (int8_t*)&bVector[number];
328 for (; number < num_points; number++) {
329 float aReal = (float)*a8Ptr++;
330 float aImag = (float)*a8Ptr++;
332 float bReal = (float)*b8Ptr++;
333 float bImag = (float)*b8Ptr++;
337 *cFloatPtr++ =
lv_creal(temp) / scalar;
338 *cFloatPtr++ =
lv_cimag(temp) / scalar;