From 96571043707be0e0e8b69a7b60c79b9afaa51ffd Mon Sep 17 00:00:00 2001 From: Yang Hau Date: Thu, 11 Jul 2024 00:57:04 +0800 Subject: [PATCH] fix: Fix type ignoring --- Makefile | 4 +- sse2neon.h | 288 ++++++++++++++++++++++++++--------------------------- 2 files changed, 146 insertions(+), 146 deletions(-) diff --git a/Makefile b/Makefile index aaee7fe4..999a3a7b 100644 --- a/Makefile +++ b/Makefile @@ -57,7 +57,7 @@ ARCH_CFLAGS := $(ARCH_CFLAGS)+$(subst $(COMMA),+,$(FEATURE)) endif endif -CXXFLAGS += -Wall -Wcast-qual -I. $(ARCH_CFLAGS) -O3 -std=gnu++14 +CXXFLAGS += -Wall -Wcast-qual -I. $(ARCH_CFLAGS) -std=gnu++14 LDFLAGS += -lm OBJS = \ tests/binding.o \ @@ -77,7 +77,7 @@ $(EXEC): $(OBJS) check: tests/main ifeq ($(processor),$(filter $(processor),aarch64 arm64 arm armv7l)) - $(CC) $(ARCH_CFLAGS) -c -O3 sse2neon.h + $(CC) $(ARCH_CFLAGS) -c sse2neon.h endif $(EXEC_WRAPPER) $^ diff --git a/sse2neon.h b/sse2neon.h index 2c179a94..cae5f09c 100644 --- a/sse2neon.h +++ b/sse2neon.h @@ -132,13 +132,13 @@ #include #include -FORCE_INLINE double recast_f64(uint64_t u64) +FORCE_INLINE double sse2neon_recast_u64_f64(uint64_t u64) { double f64; memcpy(&f64, &u64, sizeof(uint64_t)); return f64; } -FORCE_INLINE int64_t recast_i64(double f64) +FORCE_INLINE int64_t sse2neon_recast_f64_s64(double f64) { int64_t i64; memcpy(&i64, &f64, sizeof(uint64_t)); @@ -2995,10 +2995,10 @@ FORCE_INLINE __m128d _mm_add_pd(__m128d a, __m128d b) return vreinterpretq_m128d_f64( vaddq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); double c[2]; c[0] = a0 + b0; c[1] = a1 + b1; @@ -3016,9 +3016,9 @@ FORCE_INLINE __m128d _mm_add_sd(__m128d a, __m128d b) return _mm_move_sd(a, _mm_add_pd(a, b)); #else double a0, a1, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); double c[2]; c[0] = a0 + b0; c[1] = a1; @@ -3275,10 +3275,10 @@ FORCE_INLINE __m128d _mm_cmpge_pd(__m128d a, __m128d b) return vreinterpretq_m128d_u64( vcgeq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); uint64_t d[2]; d[0] = a0 >= b0 ? ~UINT64_C(0) : UINT64_C(0); d[1] = a1 >= b1 ? ~UINT64_C(0) : UINT64_C(0); @@ -3298,9 +3298,9 @@ FORCE_INLINE __m128d _mm_cmpge_sd(__m128d a, __m128d b) #else // expand "_mm_cmpge_pd()" to reduce unnecessary operations double a0, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); uint64_t a1 = vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); uint64_t d[2]; d[0] = a0 >= b0 ? ~UINT64_C(0) : UINT64_C(0); d[1] = a1; @@ -3345,10 +3345,10 @@ FORCE_INLINE __m128d _mm_cmpgt_pd(__m128d a, __m128d b) return vreinterpretq_m128d_u64( vcgtq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); uint64_t d[2]; d[0] = a0 > b0 ? ~UINT64_C(0) : UINT64_C(0); d[1] = a1 > b1 ? ~UINT64_C(0) : UINT64_C(0); @@ -3368,9 +3368,9 @@ FORCE_INLINE __m128d _mm_cmpgt_sd(__m128d a, __m128d b) #else // expand "_mm_cmpge_pd()" to reduce unnecessary operations double a0, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); uint64_t a1 = vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); uint64_t d[2]; d[0] = a0 > b0 ? ~UINT64_C(0) : UINT64_C(0); d[1] = a1; @@ -3388,10 +3388,10 @@ FORCE_INLINE __m128d _mm_cmple_pd(__m128d a, __m128d b) return vreinterpretq_m128d_u64( vcleq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); uint64_t d[2]; d[0] = a0 <= b0 ? ~UINT64_C(0) : UINT64_C(0); d[1] = a1 <= b1 ? ~UINT64_C(0) : UINT64_C(0); @@ -3411,9 +3411,9 @@ FORCE_INLINE __m128d _mm_cmple_sd(__m128d a, __m128d b) #else // expand "_mm_cmpge_pd()" to reduce unnecessary operations double a0, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); uint64_t a1 = vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); uint64_t d[2]; d[0] = a0 <= b0 ? ~UINT64_C(0) : UINT64_C(0); d[1] = a1; @@ -3461,10 +3461,10 @@ FORCE_INLINE __m128d _mm_cmplt_pd(__m128d a, __m128d b) return vreinterpretq_m128d_u64( vcltq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); uint64_t d[2]; d[0] = a0 < b0 ? ~UINT64_C(0) : UINT64_C(0); d[1] = a1 < b1 ? ~UINT64_C(0) : UINT64_C(0); @@ -3483,9 +3483,9 @@ FORCE_INLINE __m128d _mm_cmplt_sd(__m128d a, __m128d b) return _mm_move_sd(a, _mm_cmplt_pd(a, b)); #else double a0, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); uint64_t a1 = vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); uint64_t d[2]; d[0] = a0 < b0 ? ~UINT64_C(0) : UINT64_C(0); d[1] = a1; @@ -3530,10 +3530,10 @@ FORCE_INLINE __m128d _mm_cmpnge_pd(__m128d a, __m128d b) vcgeq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b)), vdupq_n_u64(UINT64_MAX))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); uint64_t d[2]; d[0] = !(a0 >= b0) ? ~UINT64_C(0) : UINT64_C(0); d[1] = !(a1 >= b1) ? ~UINT64_C(0) : UINT64_C(0); @@ -3561,10 +3561,10 @@ FORCE_INLINE __m128d _mm_cmpngt_pd(__m128d a, __m128d b) vcgtq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b)), vdupq_n_u64(UINT64_MAX))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); uint64_t d[2]; d[0] = !(a0 > b0) ? ~UINT64_C(0) : UINT64_C(0); d[1] = !(a1 > b1) ? ~UINT64_C(0) : UINT64_C(0); @@ -3592,10 +3592,10 @@ FORCE_INLINE __m128d _mm_cmpnle_pd(__m128d a, __m128d b) vcleq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b)), vdupq_n_u64(UINT64_MAX))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); uint64_t d[2]; d[0] = !(a0 <= b0) ? ~UINT64_C(0) : UINT64_C(0); d[1] = !(a1 <= b1) ? ~UINT64_C(0) : UINT64_C(0); @@ -3623,10 +3623,10 @@ FORCE_INLINE __m128d _mm_cmpnlt_pd(__m128d a, __m128d b) vcltq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b)), vdupq_n_u64(UINT64_MAX))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); uint64_t d[2]; d[0] = !(a0 < b0) ? ~UINT64_C(0) : UINT64_C(0); d[1] = !(a1 < b1) ? ~UINT64_C(0) : UINT64_C(0); @@ -3657,10 +3657,10 @@ FORCE_INLINE __m128d _mm_cmpord_pd(__m128d a, __m128d b) vceqq_f64(vreinterpretq_f64_m128d(b), vreinterpretq_f64_m128d(b)); return vreinterpretq_m128d_u64(vandq_u64(not_nan_a, not_nan_b)); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); uint64_t d[2]; d[0] = (a0 == a0 && b0 == b0) ? ~UINT64_C(0) : UINT64_C(0); d[1] = (a1 == a1 && b1 == b1) ? ~UINT64_C(0) : UINT64_C(0); @@ -3679,9 +3679,9 @@ FORCE_INLINE __m128d _mm_cmpord_sd(__m128d a, __m128d b) return _mm_move_sd(a, _mm_cmpord_pd(a, b)); #else double a0, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); uint64_t a1 = vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); uint64_t d[2]; d[0] = (a0 == a0 && b0 == b0) ? ~UINT64_C(0) : UINT64_C(0); d[1] = a1; @@ -3704,10 +3704,10 @@ FORCE_INLINE __m128d _mm_cmpunord_pd(__m128d a, __m128d b) return vreinterpretq_m128d_s32( vmvnq_s32(vreinterpretq_s32_u64(vandq_u64(not_nan_a, not_nan_b)))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); uint64_t d[2]; d[0] = (a0 == a0 && b0 == b0) ? UINT64_C(0) : ~UINT64_C(0); d[1] = (a1 == a1 && b1 == b1) ? UINT64_C(0) : ~UINT64_C(0); @@ -3726,9 +3726,9 @@ FORCE_INLINE __m128d _mm_cmpunord_sd(__m128d a, __m128d b) return _mm_move_sd(a, _mm_cmpunord_pd(a, b)); #else double a0, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); uint64_t a1 = vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); uint64_t d[2]; d[0] = (a0 == a0 && b0 == b0) ? UINT64_C(0) : ~UINT64_C(0); d[1] = a1; @@ -3746,8 +3746,8 @@ FORCE_INLINE int _mm_comige_sd(__m128d a, __m128d b) return vgetq_lane_u64(vcgeq_f64(a, b), 0) & 0x1; #else double a0, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); return a0 >= b0; #endif } @@ -3761,8 +3761,8 @@ FORCE_INLINE int _mm_comigt_sd(__m128d a, __m128d b) return vgetq_lane_u64(vcgtq_f64(a, b), 0) & 0x1; #else double a0, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); return a0 > b0; #endif @@ -3777,8 +3777,8 @@ FORCE_INLINE int _mm_comile_sd(__m128d a, __m128d b) return vgetq_lane_u64(vcleq_f64(a, b), 0) & 0x1; #else double a0, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); return a0 <= b0; #endif @@ -3793,8 +3793,8 @@ FORCE_INLINE int _mm_comilt_sd(__m128d a, __m128d b) return vgetq_lane_u64(vcltq_f64(a, b), 0) & 0x1; #else double a0, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); return a0 < b0; #endif @@ -3866,8 +3866,8 @@ FORCE_INLINE_OPTNONE __m128i _mm_cvtpd_epi32(__m128d a) #else __m128d rnd = _mm_round_pd(a, _MM_FROUND_CUR_DIRECTION); double d0, d1; - d0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(rnd), 0)); - d1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(rnd), 1)); + d0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(rnd), 0)); + d1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(rnd), 1)); return _mm_set_epi32(0, 0, (int32_t) d1, (int32_t) d0); #endif } @@ -3879,8 +3879,8 @@ FORCE_INLINE_OPTNONE __m64 _mm_cvtpd_pi32(__m128d a) { __m128d rnd = _mm_round_pd(a, _MM_FROUND_CUR_DIRECTION); double d0, d1; - d0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(rnd), 0)); - d1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(rnd), 1)); + d0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(rnd), 0)); + d1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(rnd), 1)); int32_t ALIGN_STRUCT(16) data[2] = {(int32_t) d0, (int32_t) d1}; return vreinterpret_m64_s32(vld1_s32(data)); } @@ -3896,8 +3896,8 @@ FORCE_INLINE __m128 _mm_cvtpd_ps(__m128d a) return vreinterpretq_m128_f32(vcombine_f32(tmp, vdup_n_f32(0))); #else double a0, a1; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); return _mm_set_ps(0, 0, (float) a1, (float) a0); #endif } @@ -3997,7 +3997,7 @@ FORCE_INLINE double _mm_cvtsd_f64(__m128d a) #if defined(__aarch64__) || defined(_M_ARM64) return (double) vgetq_lane_f64(vreinterpretq_f64_m128d(a), 0); #else - double _a = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double _a = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); return _a; #endif } @@ -4011,7 +4011,7 @@ FORCE_INLINE int32_t _mm_cvtsd_si32(__m128d a) return (int32_t) vgetq_lane_f64(vrndiq_f64(vreinterpretq_f64_m128d(a)), 0); #else __m128d rnd = _mm_round_pd(a, _MM_FROUND_CUR_DIRECTION); - double ret = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(rnd), 0)); + double ret = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(rnd), 0)); return (int32_t) ret; #endif } @@ -4025,7 +4025,7 @@ FORCE_INLINE int64_t _mm_cvtsd_si64(__m128d a) return (int64_t) vgetq_lane_f64(vrndiq_f64(vreinterpretq_f64_m128d(a)), 0); #else __m128d rnd = _mm_round_pd(a, _MM_FROUND_CUR_DIRECTION); - double ret = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(rnd), 0)); + double ret = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(rnd), 0)); return (int64_t) ret; #endif } @@ -4047,7 +4047,7 @@ FORCE_INLINE __m128 _mm_cvtsd_ss(__m128 a, __m128d b) vget_lane_f32(vcvt_f32_f64(vreinterpretq_f64_m128d(b)), 0), vreinterpretq_f32_m128(a), 0)); #else - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); return vreinterpretq_m128_f32( vsetq_lane_f32((float) b0, vreinterpretq_f32_m128(a), 0)); #endif @@ -4081,7 +4081,7 @@ FORCE_INLINE __m128d _mm_cvtsi32_sd(__m128d a, int32_t b) return vreinterpretq_m128d_f64( vsetq_lane_f64((double) b, vreinterpretq_f64_m128d(a), 0)); #else - int64_t _b = recast_i64((double) b); + int64_t _b = sse2neon_recast_f64_s64((double) b); return vreinterpretq_m128d_s64( vsetq_lane_s64(_b, vreinterpretq_s64_m128d(a), 0)); #endif @@ -4109,7 +4109,7 @@ FORCE_INLINE __m128d _mm_cvtsi64_sd(__m128d a, int64_t b) return vreinterpretq_m128d_f64( vsetq_lane_f64((double) b, vreinterpretq_f64_m128d(a), 0)); #else - int64_t _b = recast_i64((double) b); + int64_t _b = sse2neon_recast_f64_s64((double) b); return vreinterpretq_m128d_s64( vsetq_lane_s64(_b, vreinterpretq_s64_m128d(a), 0)); #endif @@ -4157,8 +4157,8 @@ FORCE_INLINE __m128d _mm_cvtss_sd(__m128d a, __m128 b) FORCE_INLINE __m128i _mm_cvttpd_epi32(__m128d a) { double a0, a1; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); return _mm_set_epi32(0, 0, (int32_t) a1, (int32_t) a0); } @@ -4168,8 +4168,8 @@ FORCE_INLINE __m128i _mm_cvttpd_epi32(__m128d a) FORCE_INLINE_OPTNONE __m64 _mm_cvttpd_pi32(__m128d a) { double a0, a1; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); int32_t ALIGN_STRUCT(16) data[2] = {(int32_t) a0, (int32_t) a1}; return vreinterpret_m64_s32(vld1_s32(data)); } @@ -4187,7 +4187,7 @@ FORCE_INLINE __m128i _mm_cvttps_epi32(__m128 a) // https://www.intel.com/content/www/us/en/docs/intrinsics-guide/index.html#text=_mm_cvttsd_si32 FORCE_INLINE int32_t _mm_cvttsd_si32(__m128d a) { - double _a = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double _a = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); return (int32_t) _a; } @@ -4199,7 +4199,7 @@ FORCE_INLINE int64_t _mm_cvttsd_si64(__m128d a) #if defined(__aarch64__) || defined(_M_ARM64) return vgetq_lane_s64(vcvtq_s64_f64(vreinterpretq_f64_m128d(a)), 0); #else - double _a = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double _a = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); return (int64_t) _a; #endif } @@ -4218,10 +4218,10 @@ FORCE_INLINE __m128d _mm_div_pd(__m128d a, __m128d b) return vreinterpretq_m128d_f64( vdivq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); double c[2]; c[0] = a0 / b0; c[1] = a1 / b1; @@ -4469,15 +4469,15 @@ FORCE_INLINE __m128d _mm_max_pd(__m128d a, __m128d b) vmaxq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b))); #endif #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); - uint64_t d[2]; - d[0] = a0 > b0 ? recast_i64(a0) : recast_i64(b0); - d[1] = a1 > b1 ? recast_i64(a1) : recast_i64(b1); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + int64_t d[2]; + d[0] = a0 > b0 ? sse2neon_recast_f64_s64(a0) : sse2neon_recast_f64_s64(b0); + d[1] = a1 > b1 ? sse2neon_recast_f64_s64(a1) : sse2neon_recast_f64_s64(b1); - return vreinterpretq_m128d_u64(vld1q_u64(d)); + return vreinterpretq_m128d_s64(vld1q_s64(d)); #endif } @@ -4491,9 +4491,9 @@ FORCE_INLINE __m128d _mm_max_sd(__m128d a, __m128d b) return _mm_move_sd(a, _mm_max_pd(a, b)); #else double a0, a1, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); double c[2] = {a0 > b0 ? a0 : b0, a1}; return vreinterpretq_m128d_f32(vld1q_f32((float32_t *) c)); #endif @@ -4532,14 +4532,14 @@ FORCE_INLINE __m128d _mm_min_pd(__m128d a, __m128d b) vminq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b))); #endif #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); - uint64_t d[2]; - d[0] = a0 < b0 ? recast_i64(a0) : recast_i64(b0); - d[1] = a1 < b1 ? recast_i64(a1) : recast_i64(b1); - return vreinterpretq_m128d_u64(vld1q_u64(d)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + int64_t d[2]; + d[0] = a0 < b0 ? sse2neon_recast_f64_s64(a0) : sse2neon_recast_f64_s64(b0); + d[1] = a1 < b1 ? sse2neon_recast_f64_s64(a1) : sse2neon_recast_f64_s64(b1); + return vreinterpretq_m128d_s64(vld1q_s64(d)); #endif } @@ -4553,9 +4553,9 @@ FORCE_INLINE __m128d _mm_min_sd(__m128d a, __m128d b) return _mm_move_sd(a, _mm_min_pd(a, b)); #else double a0, a1, b0; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); double c[2] = {a0 < b0 ? a0 : b0, a1}; return vreinterpretq_m128d_f32(vld1q_f32((float32_t *) c)); #endif @@ -4711,10 +4711,10 @@ FORCE_INLINE __m128d _mm_mul_pd(__m128d a, __m128d b) return vreinterpretq_m128d_f64( vmulq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); double c[2]; c[0] = a0 * b0; c[1] = a1 * b1; @@ -5007,7 +5007,7 @@ FORCE_INLINE __m128d _mm_set1_pd(double d) #if defined(__aarch64__) || defined(_M_ARM64) return vreinterpretq_m128d_f64(vdupq_n_f64(d)); #else - int64_t _d = recast_i64(d); + int64_t _d = sse2neon_recast_f64_s64(d); return vreinterpretq_m128d_s64(vdupq_n_s64(_d)); #endif } @@ -5300,8 +5300,8 @@ FORCE_INLINE __m128d _mm_sqrt_pd(__m128d a) return vreinterpretq_m128d_f64(vsqrtq_f64(vreinterpretq_f64_m128d(a))); #else double a0, a1; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); double _a0 = sqrt(a0); double _a1 = sqrt(a1); return _mm_set_pd(_a1, _a0); @@ -5318,8 +5318,8 @@ FORCE_INLINE __m128d _mm_sqrt_sd(__m128d a, __m128d b) return _mm_move_sd(a, _mm_sqrt_pd(b)); #else double _a, _b; - _a = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - _b = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + _a = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + _b = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); return _mm_set_pd(_a, sqrt(_b)); #endif } @@ -5675,10 +5675,10 @@ FORCE_INLINE __m128d _mm_sub_pd(__m128d a, __m128d b) return vreinterpretq_m128d_f64( vsubq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); double c[2]; c[0] = a0 - b0; c[1] = a1 - b1; @@ -5985,10 +5985,10 @@ FORCE_INLINE __m128d _mm_hadd_pd(__m128d a, __m128d b) return vreinterpretq_m128d_f64( vpaddq_f64(vreinterpretq_f64_m128d(a), vreinterpretq_f64_m128d(b))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); double c[] = {a0 + a1, b0 + b1}; return vreinterpretq_m128d_u64(vld1q_u64((uint64_t *) c)); #endif @@ -6023,10 +6023,10 @@ FORCE_INLINE __m128d _mm_hsub_pd(__m128d a, __m128d b) return vreinterpretq_m128d_f64( vsubq_f64(vuzp1q_f64(_a, _b), vuzp2q_f64(_a, _b))); #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); double c[] = {a0 - a1, b0 - b1}; return vreinterpretq_m128d_u64(vld1q_u64((uint64_t *) c)); #endif @@ -6824,8 +6824,8 @@ FORCE_INLINE __m128d _mm_ceil_pd(__m128d a) return vreinterpretq_m128d_f64(vrndpq_f64(vreinterpretq_f64_m128d(a))); #else double a0, a1; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); return _mm_set_pd(ceil(a1), ceil(a0)); #endif } @@ -7034,10 +7034,10 @@ FORCE_INLINE __m128d _mm_dp_pd(__m128d a, __m128d b, const int imm) vgetq_lane_f64(vreinterpretq_f64_m128d(b), 1) : 0; #else - double a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - double a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); - double b0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); - double b1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); + double a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + double a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + double b0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 0)); + double b1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(b), 1)); double d0 = (imm & 0x10) ? a0 * b0 : 0; double d1 = (imm & 0x20) ? a1 * b1 : 0; #endif @@ -7047,8 +7047,8 @@ FORCE_INLINE __m128d _mm_dp_pd(__m128d a, __m128d b, const int imm) #if defined(__aarch64__) || defined(_M_ARM64) double sum = vpaddd_f64(vreinterpretq_f64_m128d(tmp)); #else - double _tmp0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(tmp), 0)); - double _tmp1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(tmp), 1)); + double _tmp0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(tmp), 0)); + double _tmp1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(tmp), 1)); double sum = _tmp0 + _tmp1; #endif // Conditionally store the sum @@ -7140,8 +7140,8 @@ FORCE_INLINE __m128d _mm_floor_pd(__m128d a) return vreinterpretq_m128d_f64(vrndmq_f64(vreinterpretq_f64_m128d(a))); #else double a0, a1; - a0 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); - a1 = recast_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); + a0 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 0)); + a1 = sse2neon_recast_u64_f64(vgetq_lane_u64(vreinterpretq_u64_m128d(a), 1)); return _mm_set_pd(floor(a1), floor(a0)); #endif }