Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

volk_32fc_s32f_atan2_32f: Add NaN tests for avx2 and avx2_fma code #731

Merged
merged 1 commit into from
Jan 7, 2024
Merged

volk_32fc_s32f_atan2_32f: Add NaN tests for avx2 and avx2_fma code #731

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
12 changes: 12 additions & 0 deletions kernels/volk/volk_32fc_s32f_atan2_32f.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -117,6 +117,7 @@ static inline void volk_32fc_s32f_atan2_32f_a_avx2_fma(float* outputVector,
const __m256 pi_2 = _mm256_set1_ps(0x1.921fb6p0f);
const __m256 abs_mask = _mm256_castsi256_ps(_mm256_set1_epi32(0x7FFFFFFF));
const __m256 sign_mask = _mm256_castsi256_ps(_mm256_set1_epi32(0x80000000));
const __m256 zero = _mm256_setzero_ps();

unsigned int number = 0;
unsigned int eighth_points = num_points / 8;
Expand All @@ -133,6 +134,8 @@ static inline void volk_32fc_s32f_atan2_32f_a_avx2_fma(float* outputVector,
_mm256_and_ps(y, abs_mask), _mm256_and_ps(x, abs_mask), _CMP_GT_OS);
__m256 input = _mm256_div_ps(_mm256_blendv_ps(y, x, swap_mask),
_mm256_blendv_ps(x, y, swap_mask));
__m256 nan_mask = _mm256_cmp_ps(input, input, _CMP_UNORD_Q);
input = _mm256_blendv_ps(input, zero, nan_mask);
__m256 result = _m256_arctan_poly_avx2_fma(input);

input =
Expand Down Expand Up @@ -174,6 +177,7 @@ static inline void volk_32fc_s32f_atan2_32f_a_avx2(float* outputVector,
const __m256 pi_2 = _mm256_set1_ps(0x1.921fb6p0f);
const __m256 abs_mask = _mm256_castsi256_ps(_mm256_set1_epi32(0x7FFFFFFF));
const __m256 sign_mask = _mm256_castsi256_ps(_mm256_set1_epi32(0x80000000));
const __m256 zero = _mm256_setzero_ps();

unsigned int number = 0;
unsigned int eighth_points = num_points / 8;
Expand All @@ -190,6 +194,8 @@ static inline void volk_32fc_s32f_atan2_32f_a_avx2(float* outputVector,
_mm256_and_ps(y, abs_mask), _mm256_and_ps(x, abs_mask), _CMP_GT_OS);
__m256 input = _mm256_div_ps(_mm256_blendv_ps(y, x, swap_mask),
_mm256_blendv_ps(x, y, swap_mask));
__m256 nan_mask = _mm256_cmp_ps(input, input, _CMP_UNORD_Q);
input = _mm256_blendv_ps(input, zero, nan_mask);
__m256 result = _m256_arctan_poly_avx(input);

input =
Expand Down Expand Up @@ -235,6 +241,7 @@ static inline void volk_32fc_s32f_atan2_32f_u_avx2_fma(float* outputVector,
const __m256 pi_2 = _mm256_set1_ps(0x1.921fb6p0f);
const __m256 abs_mask = _mm256_castsi256_ps(_mm256_set1_epi32(0x7FFFFFFF));
const __m256 sign_mask = _mm256_castsi256_ps(_mm256_set1_epi32(0x80000000));
const __m256 zero = _mm256_setzero_ps();

unsigned int number = 0;
unsigned int eighth_points = num_points / 8;
Expand All @@ -251,6 +258,8 @@ static inline void volk_32fc_s32f_atan2_32f_u_avx2_fma(float* outputVector,
_mm256_and_ps(y, abs_mask), _mm256_and_ps(x, abs_mask), _CMP_GT_OS);
__m256 input = _mm256_div_ps(_mm256_blendv_ps(y, x, swap_mask),
_mm256_blendv_ps(x, y, swap_mask));
__m256 nan_mask = _mm256_cmp_ps(input, input, _CMP_UNORD_Q);
input = _mm256_blendv_ps(input, zero, nan_mask);
__m256 result = _m256_arctan_poly_avx2_fma(input);

input =
Expand Down Expand Up @@ -292,6 +301,7 @@ static inline void volk_32fc_s32f_atan2_32f_u_avx2(float* outputVector,
const __m256 pi_2 = _mm256_set1_ps(0x1.921fb6p0f);
const __m256 abs_mask = _mm256_castsi256_ps(_mm256_set1_epi32(0x7FFFFFFF));
const __m256 sign_mask = _mm256_castsi256_ps(_mm256_set1_epi32(0x80000000));
const __m256 zero = _mm256_setzero_ps();

unsigned int number = 0;
unsigned int eighth_points = num_points / 8;
Expand All @@ -308,6 +318,8 @@ static inline void volk_32fc_s32f_atan2_32f_u_avx2(float* outputVector,
_mm256_and_ps(y, abs_mask), _mm256_and_ps(x, abs_mask), _CMP_GT_OS);
__m256 input = _mm256_div_ps(_mm256_blendv_ps(y, x, swap_mask),
_mm256_blendv_ps(x, y, swap_mask));
__m256 nan_mask = _mm256_cmp_ps(input, input, _CMP_UNORD_Q);
input = _mm256_blendv_ps(input, zero, nan_mask);
__m256 result = _m256_arctan_poly_avx(input);

input =
Expand Down
Loading