|
|
|
#include<simdconfig.h>
|
|
|
|
#include<simdfuncs.h>
|
|
|
|
|
|
|
|
#include<emmintrin.h>
|
|
|
|
#include<tmmintrin.h>
|
|
|
|
|
|
|
|
#ifdef _MSC_VER
|
|
|
|
#include<intrin.h>
|
|
|
|
|
|
|
|
int ssse3_available(void) {
|
|
|
|
return 1;
|
|
|
|
}
|
|
|
|
|
|
|
|
#else
|
|
|
|
|
|
|
|
#include<cpuid.h>
|
|
|
|
#include<stdint.h>
|
|
|
|
|
|
|
|
int ssse3_available(void) {
|
|
|
|
#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]) {
|
|
|
|
ALIGN_16 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];
|
|
|
|
}
|