SSE2饱和算法

Question

我正在编写一些音频处理软件，我需要知道如何使用SSE2双精度指令来完成饱和算术。我的值需要在-1到1之间进行规范化。有什么聪明的方法可以用SSE2内部的方法来实现这一点呢?还是我需要2组if/ each语句(每个值一个)？

Answer 1

要将双精度值剪辑到-1.0到+1.0之间，可以使用max/min操作。例如，如果您有N个double值的缓冲区，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);
}