c++ - SSE 和 AVX 的 channel / channel 改组?

标签 c++ c vectorization sse avx

哪些 SSE/AVX 指令将 channel 从 a 打乱为 bc

float4 a = {data[0], data[1], data[2], data[3]};
float4 b = {data[1], data[2], data[3], data[0]};  // lanes shifted left
float4 c = {data[3], data[0], data[1], data[2]};  // lanes shifted right

float8 a = {data[0], data[1], data[2], data[3],
            data[4], data[5], data[6], data[7]};
float8 b = {data[1], data[2], data[3], data[4],
            data[5], data[6], data[7], data[0]};  // lanes shifted left
float8 c = {data[7], data[0], data[1], data[2],
            data[3], data[4], data[5], data[6]};  // lanes shifted right

背景:

我有一个算法需要来自它的邻居点的值;这意味着我目前正在混合对齐的负载和未对齐的负载:

(
      plate[row + 1][column]    // aligned
    + plate[row - 1][column]    // aligned
    + plate[row][column + 1]    // unaligned
    + plate[row][column - 1]    // unaligned
    + (4 * plate[row][column])  // aligned
) / 8;

这是在 SSE 中:

__m128 bottom  = _mm_load_ps(&from[row-1][column]);
__m128 left    = _mm_loadu_ps(&from[row][column-1]);
__m128 middle  = _mm_load_ps(&from[row][column]);
__m128 right   = _mm_loadu_ps(&from[row][column+1]);
__m128 top     = _mm_load_ps&from[row+1][column]);

(top + bottom + left + right + _mm_set1_ps(4.0f) * middle) * _mm_set1_ps(0.125f);

我意识到 current 和 left 或 right 中的值仅相差一个值。所以我有一个想法,而不是做 2 个未对齐的加载,我也许可以洗牌 channel ,然后插入一个不同的值。我需要检查指令的延迟/吞吐量,看看是否会更快,但我对这些类型的 SSE/AVX 指令都不熟悉。

最佳答案

在最近的 Intel CPU(Core i7 等)上,未对齐负载是一种合理的方法,但在较旧的 CPU 上,它相对昂贵。另一种方法是使用 _mm_alignr_epi8 (PALIGNR) - 通常你沿着一行迭代并保持 3 个连续的 vector - 在每次迭代之后你沿着一个 vector 打乱这些 vector 然后加载一个新 vector ,因此每次迭代只有一个负载。

__m128 va = _mm_setzero_ps();
__m128 vb = _mm_load_ps(&from[row][0]);
for (col = 0; col < N; col += 4)
{
    __m128 vc = _mm_load_ps(&from[row][col + 4]);

    __m128 centre = vb;
    __m128 left = (__m128)_mm_alignr_epi8((__m128i)va, (__m128i)vb, sizeof(float));
    __m128 right = (__m128)_mm_alignr_epi8((__m128i)vb, (__m128i)vc, 3 * sizeof(float));

    // do stuff ...

    va = vb;  // shuffle vectors along
    vb = vc;
}

由于 128 位 channel 的限制,AVX 有点棘手 - 您最好坚持使用未对齐的负载。

关于c++ - SSE 和 AVX 的 channel / channel 改组?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/19864015/

相关文章:

c - 逗号运算符在分配期间如何工作?

c++ - 是否有机会使用 SIMD 加速循环代码?

c++ - 使用混合(浮点、 double )输入 vector 执行 AVX 内积运算的最快方法

algorithm - 用于半稀疏数组中最近邻计算的 Numpy 矢量化

c++ - 删除条目 map<int, A*> 的最佳方法是什么

c++ - unsigned char 的 gdb 表示形式——是否不一致?

c++ - test.cpp :(. text._ZN1AC2Ev[_ZN1AC5Ev]+0x8)||对 `vtable for A' 的 undefined reference |

c - 扩展程序以使用三个命令实现流水线

c++ - 派生类中的显式模板静态成员实例化

c++ - 高效布局和减少虚拟2D数据(抽象)