c - SSE2 饱和算法

Question

我正在编写一些音频处理软件，我需要知道如何使用 SSE2 double 指令进行饱和运算。我的值需要在 -1 和 1 之间归一化。是否有一种巧妙的方法可以使用 SSE2 内在函数来执行此操作，或者我是否需要 2 组 if/else 语句(每个值一个)？

Answer 1

要将 double 值限制在 -1.0 到 +1.0 的范围内，您可以使用最大/最小操作。例如。如果你有一个缓冲区 buff，N 个 double 值:

const __m128d kMax = _mm_set1_pd(1.0);
const __m128d kMin = _mm_set1_pd(-1.0);

for (int i = 0; i < N; i += 2)
{
    __m128d v = _mm_loadu_pd(&buff[i]);
    v = _mm_max_pd(v, kMin);
    v = _mm_min_pd(v, kMax);
    _mm_storeu_pd(&buff[i], v);
}