sse - 如何将 __m256d 中的所有元素设置为另一个 __m256d 的第三个元素？

Question

有了 4 个打包的 float (__m128)，我可以使用 SSE intrinsic

__m128 X;
__m128 H = _mm_shuffle_ps(X,X,_MM_SHUFFLE(3,3,3,3));

将H的所有元素设置为X的第三个元素(这是最快的方法吗？)

现在，我想对 4 个打包的 double (__m256d) 执行相同的操作。我天真地编码

__m256d X;
__m256d H = _mm256_shuffle_pd(X,X,_MM_SHUFFLE(3,3,3,3));

但这不是正确的做法!相反，它设置 H={X[1],X[1],X[3],X[3]}。

那么，怎样做才是正确的呢？

编辑

使用 Intel(R) Xeon(R) CPU E5-2670 0 @ 2.60GHz

Answer 1

它并不总是最优的，但询问你的编译器它认为什么可能是一个很好的提示。

#include <x86intrin.h>
__m256d f(__m256d x){
  __m256i m={3,3,3,3};
  return __builtin_shuffle(x,m);
}

对于 gcc-4.8，这会生成:

vpermilpd   $15, %ymm0, %ymm0
vperm2f128  $17, %ymm0, %ymm0, %ymm0

clang 有一个不同的 builtin for shuffling，我不知道其他编译器是否有。