gcc - 假人的 SSE 浮点点积

Question

我已经阅读了许多关于 SSE/SIMD 的 SO 问题（例如，SSE 入门），但我仍然对所有这些问题感到困惑。我想要的只是两个双精度浮点向量之间的点积，用 C 语言（C99 FWIW）。我正在使用 GCC。

有人可以发布一个简单而完整的示例，包括如何将双向量转换为 SSE 类型并再次转换回来？

[编辑2012-10-08]

这是我设法拼凑起来的一些 SSE2 代码，批评？

#include <emmintrin.h>

double dotprod(double *restrict a, double *restrict b, int n)
{
   __m128d aa, bb, cc, ss;
   int i, n1 = n - 1;
   double *s = calloc(2, sizeof(double));
   double s2 = 0;

   ss = _mm_set1_pd(0);

   for(i = 0 ; i < n1 ; i += 2) 
   {
     aa = _mm_load_pd(a + i);
     bb = _mm_load_pd(b + i);
     cc = _mm_mul_pd(aa, bb);
     ss = _mm_add_pd(ss, cc);
   }

   _mm_store_pd(s, ss);
   s2 = s[0] + s[1];

   if(i < n)
      s2 += a[i] * b[i];

   free(s);

   return s2;
}