задание с Intel Intrinsics - горизонтальное добавление

Question

Я хочу суммировать все элементы большого вектора ary.Моя идея состояла в том, чтобы сделать это с горизонтальной суммой.

const int simd_width = 16/sizeof(float); 
float helper[simd_width];

//take the first 4 elements
const __m128 a4 = _mm_load_ps(ary);

for(int i=0; i<N-simd_width; i+=simd_width){
     const __m128 b4 = _mm_load_ps(ary+i+simd_width);
     //save temporary result in helper array
     _mm_store_ps(helper, _mm_hadd_ps(a4,b4)); //C
     const __m128 a4 = _mm_load_ps(helper);

}

Я искал метод, с помощью которого я мог бы присвоить результирующий вектор непосредственно четырехугольному плавучему элементу a4 непосредственно как _mm_store_ps(a4, _mm_hadd_ps(a4,b4))Есть ли такой метод Intel?(Я впервые работаю с SSE - возможно, весь фрагмент кода неверен)

Answer 1

Как предположил Питер, не используйте горизонтальные суммы.Используйте вертикальные суммы.

Например, в псевдокоде с шириной simd = 2

SIMD sum = {0,0}; // we use 2 accumulators
for (int i = 0; i + 1 < n; i += 2)
    sum = simd_add(sum, simd_load(x+i));
float s = horizzontal_add(sum);
if (n & 1)  // n was not a multiple of 2?
   s += x[n-1]; // deal with last element