#include #include #include #include #ifdef _MSC_VER #include int ssse3_available() { return 1; } #else #include #include int ssse3_available() { #ifdef __APPLE__ return 1; #elif defined(__clang__) /* https://github.com/numpy/numpy/issues/8130 */ return __builtin_cpu_supports("sse4.1"); #else return __builtin_cpu_supports("ssse3"); #endif } #endif void increment_ssse3(float arr[4]) { double darr[4]; __m128d val1 = _mm_set_pd(arr[0], arr[1]); __m128d val2 = _mm_set_pd(arr[2], arr[3]); __m128d one = _mm_set_pd(1.0, 1.0); __m128d result = _mm_add_pd(val1, one); __m128i tmp1, tmp2; tmp1 = tmp2 = _mm_set1_epi16(0); _mm_store_pd(darr, result); result = _mm_add_pd(val2, one); _mm_store_pd(&darr[2], result); tmp1 = _mm_hadd_epi32(tmp1, tmp2); /* This does nothing. Only here so we use an SSSE3 instruction. */ arr[0] = (float)darr[1]; arr[1] = (float)darr[0]; arr[2] = (float)darr[3]; arr[3] = (float)darr[2]; }