使用 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);

}

我在寻找一种方法，用它我可以像 _mm_store_ps(a4, _mm_hadd_ps(a4,b4)) 一样将结果向量直接分配给 quadfloat a4 有这样的英特尔方法吗？（这是我第一次使用 SSE - 可能整个代码片段都是错误的）

Answer 1

正如彼得所建议的那样，不要使用水平求和。使用垂直总和。

例如，在伪代码中，simd width = 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

使用 intel Intrinsics 进行赋值 - 水平添加

assignment with intel Intrinsics - horizontal add

c++

sse

simd