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